[R] bootstrapping problem
Uwe Ligges
ligges at statistik.tu-dortmund.de
Mon May 2 10:29:25 CEST 2011
On 29.04.2011 09:01, mimato wrote:
> I want to classify bipolar neurons in human cochleas and have data of the
> following structure:
>
> Vol_Nuc Vol_Soma
> 1 186.23 731.96
> 2 204.58 4370.96
> 3 539.98 7344.86
> 4 477.71 6939.28
> 5 421.22 5588.53
> 6 276.61 1017.05
> 7 392.28 6392.32
> 8 424.43 6190.13
> 9 256.41 3850.51
> 10 249.17 3118.14
> 11 276.97 3037.29
> 12 295.30 3703.76
> 13 314.43 5265.97
> 14 301.15 5781.73
>
> I already worked with Matlab (I´m not a programmer) and created nice
> colourcoded dendrograms and also made some verifications of them. I started
> now working with R and bootstrapped data with a library named pvclust. It
> worked and R computed ...
>
> here is the code:
>
> library(pvclust)
>
> data =
> data.frame(Vol_Nuk=c(186.23,204.58,539.98,477.71,421.22,276.61,392.28,424.43,256.41,249.17,276.97,295.3,314.43,301.15),
> Vol_Soma=c(731.96,4370.96,7344.86,6939.28,5588.53,1017.05,6392.32,6190.13,3850.51,3118.14,3037.29,3703.76,5265.97,5781.73))
>
> plot(data)
> result<-pvclust(data,nboot=100)
> plot(result)
>
> It is also not working using following commands:
>
> cluster.bootstrap<- pvclust(Raw, nboot=1000, method.dist="abscor")
> plot(cluster.bootstrap)
> pvrect(cluster.bootstrap)
>
> I always get the following problem:
>
> mistake in plot.hclust(x$hclust, main = main, sub = sub, xlab = xlab, col =
> col, :
> invalid input for Dendrogram
>
> Does anyone has an idea whats wrong...
Yes: You have only 2 variables and when clustering these 2 variables it
makes no sense to plot this. My I guess you want to cluster observations
rather than variables? Then transpose your data before applying pvclust.
Best,
Uwe Ligges
> Thanks a lot!!
>
> --
> View this message in context: http://r.789695.n4.nabble.com/bootstrapping-problem-tp3483068p3483068.html
> Sent from the R help mailing list archive at Nabble.com.
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
More information about the R-help
mailing list