R : 루프를 사용하여 1에서 n을 리샘플링 (y)

Question

나는 각각 100 개의 행을 갖는 n 개의 변수가 있습니다. 1에서 nrows로 다시 샘플링하려면 다음 코드를 사용하여 예상 결과를 얻지 만, 지루하고 비실용적입니다.

y<-rnorm(n=5, mean=10, sd=2) 
R=1000 #number of resamplings 
boot.means = numeric(R) 
for (i in 1:R) { boot.sample = sample(y, 1, replace=T) 
boot.means[i] = mean(boot.sample) } 
m1<-mean(boot.means) 
d1<-sd(boot.means) 
cv1 =(d1*100)/m1 

R=1000 #number of resamplings 
boot.means = numeric(R) 
for (i in 1:R) { boot.sample = sample(y, 2, replace=T) 
boot.means[i] = mean(boot.sample) } 
m2<-mean(boot.means) 
d2<-sd(boot.means) 
cv2 =(d2*100)/m2 

R=1000 #number of resamplings 
boot.means = numeric(R) 
for (i in 1:R) { boot.sample = sample(y, 3, replace=T) 
boot.means[i] = mean(boot.sample) } 
m3<-mean(boot.means) 
d3<-sd(boot.means) 
cv3 =(d3*100)/m3 


R=1000 #number of resamplings 
boot.means = numeric(R) 
for (i in 1:R) { boot.sample = sample(y, 4, replace=T) 
boot.means[i] = mean(boot.sample) } 
m4<-mean(boot.means) 
d4<-sd(boot.means) 
cv4 =(d4*100)/m4 


R=1000 #number of resamplings 
boot.means = numeric(R) 
for (i in 1:R) { boot.sample = sample(y, 5, replace=T) 
boot.means[i] = mean(boot.sample) } 
m5<-mean(boot.means) 
d5<-sd(boot.means) 
cv5 =(d5*100)/m5 

CV.OK<-(c(cv1,cv2,cv3,cv4,cv5)) 
plot(CV.OK) 

나는 다음과 같은 코드를 사용하고자하는

하지만주는 예상치 못한 결과 : 상황을 재현하려면 y는 5 개 행이 suposse을 할 수 있습니다. 제발, 누군가 도울 수있어. 감사.

R = 1000 #number of resamplings 
boot.sample=seq(1,5, by=1) 
boot.means = numeric(R) 
boot.sd = numeric(R) 
m = 5 
d = 5 
for (i in 1:5) { 
    for (j in 1:R) { 
    boot.sample[i] = sample(y, i, replace=T) 
    boot.means[j] = mean(boot.sample[i]) 
    boot.sd[j] = sd(boot.sample[i]) 
    m[i]=mean(boot.means[j]) 
    d[i]=mean(boot.sd[j]) 
    } 
} 
CV.Fail<-(d*100)/m 

Answer 1

난 당신이 원하는 생각 :

y<-rnorm(n=5, mean=10, sd=2) 
R = 1000 #number of resamplings 
CVs <- numeric(5) 
for (i in 1:5) { 
    boot.means = numeric(R) 
    for (j in 1:R) { 
    boot.sample = sample(y, i, replace=T) 
    boot.means[j] = mean(boot.sample) 
    } 
    m=mean(boot.means) 
    d=sd(boot.means) 
    CVs[i] = (d*100)/m 
} 
plot(CVs) 

Answer 2

R에, 당신은 그들이 매우 느린 때문에 루프를 피하려고한다. 문제를 올바르게 이해하고 다른 점에서 시작해야하는 약간의 기능을 작성했기를 바랍니다.

library(plyr) 
library(dplyr) 

# dummy data set 
data_set = data.frame(value = runif(200), group = rep(c("a", "b"), each=100)) 

# create a function that takes the sample size as an argument 
iterative_sample = function(sample_size, data){ 
# group the data (your 'n' equals the number of groups- 
# here thats 'a' and 'b' 
    sample_temp = dplyr::group_by(data, group) %>% 
    # take x (sample size) samples from each group 
    sample_n(sample_size, replace=T) %>% 
    # compute summary stats for each group 
    summarize(mean = mean(value), sd = sd(value)) %>% 
    # attach the sample size to keep track 
    mutate(sample_size = sample_size) 
    # we must return a dataframe to uses ldply later on 
    return(sample_temp) 
} 

# thats the vector we are going to iterate over using ldply 
sample_vect = c(1:2) 

# ldplyr (plyr package) takes a list or vector and returns a dataframe and our custom 
# function -checkout the manpage 
# ?ldply 

# ... 
# 
# 
# .data: list to be processed 
# 
#  .fun: function to apply to each piece 
# 
#  ...: other arguments passed on to ‘.fun’ 
# 
# ... 
# 

ldply(.data = sample_vect, .fun = iterative_sample, data_set)