Dear All, I am trying to bootstrap a large data matrix 'exp' of dimension 9275x898 using package 'boot'. I am trying to calculate correlation of the obtained samples. I have done following so far:
library(boot)
load("ratio_exp.RData") #loading the data matrix
dim(ratio_exp)
[1] 9275 898
sample.cor=function(ratio_exp,d){return(cor(t(ratio_exp[d,])))} #function for obtaing correlation
boot.exp=boot(ratio_exp,sample.cor,R=50)
Error in matrix(NA, sum(R), lt0) : too many elements specified But when I reduced the no of replication that is 10 then I got following
boot.exp=boot(ratio_exp,sample.cor,R=10) str(boot.exp$t)
num [1:5, 1:86025625] 1 1 1 1 1 ... In both the cases its not correct. So I am doing something wrong in writing the function sample.cor. Please help! regards Amit