Hello everyone,
I have a file called "trans" which has 100 records (rows) and 60 columns.
I need to make 2 samples (with 5 records in each) out of this file with no
replacement, so I used "replicate".
The sampling part works fine and I called the outcome "result".
my problem is then I need to calculate the correlation matrix for each
sample. I used "do.call" and called the outcome "CORE". but the problem is
when I get my out put it's not separated for each sample. they are all
together and I need them to be separated. just to give U an idea I copy and
paste the out put:
616 783 805 134
217
616 1.00000000 -0.08522865 -0.15018785 -0.12462239 -0.04279605
783 -0.08522865 1.00000000 -0.02986738 -0.04731350 0.33191790
805 -0.15018785 -0.02986738 1.00000000 0.15483873 -0.08998425
134 -0.12462239 -0.04731350 0.15483873 1.00000000 0.16290948
217 -0.04279605 0.33191790 -0.08998425 0.16290948 1.00000000
607 1.00000000 -0.12662932 0.16116459 -0.03479445 -0.03479445 *(from
here is for the 2nd sample)
604 -0.12662932 1.00000000 0.01298701 -0.06168397 -0.06168397
698 0.16116459 0.01298701 1.00000000 -0.03925343 -0.03925343
639 -0.03479445 -0.06168397 -0.03925343 1.00000000 -0.01694915
989 -0.03479445 -0.06168397 -0.03925343 -0.01694915 1.00000000
does any one has any idea how can I separate the result for 2 samples?
This is my code:
result=replicate(2,trans[,sample(colnames(trans),5,replace =
FALSE)],simplify=FALSE)
CORE=do.call(rbind,lapply(result,function(x) cor(x)))
--
View this message in context:
http://www.nabble.com/Help-with-%22replicate%22-tp20722891p20722891.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
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.
