Minor edit > dta2 <- unstack(dta, score~Name)
is simpler. ------- David C > -----Original Message----- > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > project.org] On Behalf Of David L Carlson > Sent: Wednesday, February 13, 2013 9:40 AM > To: 'Jorge I Velez'; 'Nico Met' > Cc: 'R help' > Subject: Re: [R] Correlation with p value > > You can use package Hmisc. Assuming your data.frame is called "dta": > > > library(Hmisc) > > dta2 <- do.call(cbind, split(dta[,2], dta$Name)) > > dta2 > CTJ PKR TTK > [1,] 83.99429 88.54571 86.43714 > [2,] 86.43714 85.71714 89.70286 > [3,] 84.59600 84.90200 87.72800 > [4,] 84.74000 89.88286 89.99000 > [5,] 87.87714 91.55429 89.42000 > [6,] 83.12000 89.42000 85.69143 > [7,] 81.82400 86.81000 82.25600 > > rcorr(dta2) > CTJ PKR TTK > CTJ 1.00 0.29 0.85 > PKR 0.29 1.00 0.22 > TTK 0.85 0.22 1.00 > > n= 7 > > > P > CTJ PKR TTK > CTJ 0.5214 0.0152 > PKR 0.5214 0.6277 > TTK 0.0152 0.6277 > > ---------------------------------------------- > David L Carlson > Associate Professor of Anthropology > Texas A&M University > College Station, TX 77843-4352 > > > > -----Original Message----- > > From: r-help-boun...@r-project.org [mailto:r-help-bounces@r- > > project.org] On Behalf Of Jorge I Velez > > Sent: Wednesday, February 13, 2013 5:22 AM > > To: Nico Met > > Cc: R help > > Subject: Re: [R] Correlation with p value > > > > Nico, > > > > Check https://stat.ethz.ch/pipermail/r-help/2008-May/161725.html for > > some > > alternatives. You might have to change the structure of your data, > > though. > > > > HTH, > > Jorge.- > > > > > > On Wed, Feb 13, 2013 at 10:13 PM, Nico Met <> wrote: > > > > > Dear all, > > > > > > I have a data (bellow) and I want to make a correlation test with > p- > > value > > > > > > structure(list(Name = structure(c(3L, 3L, 3L, 3L, 3L, 3L, 3L, > > > 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = > > c("CTJ", > > > "PKR", "TTK"), class = "factor"), score = c(86.4371428571428, > > > 89.7028571428572, 87.728, 89.99, 89.42, 85.6914285714286, 82.256, > > > 83.9942857142857, 86.4371428571429, 84.596, 84.74, > 87.8771428571428, > > > 83.12, 81.824, 88.5457142857143, 85.7171428571429, 84.902, > > > 89.8828571428571, > > > 91.5542857142857, 89.42, 86.81)), .Names = c("Name", "score"), > class > > = > > > "data.frame", row.names = c(NA, > > > 21L)) > > > > > > And finally , I would like to display as matrix where I can fill up > > with > > > p-values > > > > > > > > > *TTK* > > > *CTJ* > > > *PKR* > > > *TTK* > > > ** > > > ** > > > ** > > > *CTJ* > > > ** > > > ** > > > ** > > > *PKR* > > > > > > > > > > > > > > > > > > Many thanks for your help in advance > > > > > > Nico > > > > > > [[alternative HTML version deleted]] > > > > > > ______________________________________________ > > > R-help@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. > > > > > > > [[alternative HTML version deleted]] > > > > ______________________________________________ > > R-help@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. > > ______________________________________________ > R-help@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. ______________________________________________ R-help@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.
