求助：排序以及index number 或row name

ywh19860616 发表于 2014-4-2 21:49
可以分开来抽样，你对应修改就可以

非常感谢耐心的回复！
分开抽样会有个问题，比如第一个人是30岁，收入1000元，支出300元，第二个人23岁，收入500元，支出200元。分开对于年龄，收入，支出抽样的话，子样本的年龄收入和支出就不是一一对应的关系了。

winniewang2222 发表于 2014-4-2 21:56
非常感谢耐心的回复！
分开抽样会有个问题，比如第一个人是30岁，收入1000元，支出300元，第二个人23岁， ...

x1和x2放在一起抽样，可以这样

1. x1 <- rnorm(100, mean = 0, sd = 1)
   
2. x2 <- rnorm(100, mean = 2, sd = 5)
   
3. y <- rnorm(100, mean = 1, sd = 2)
   
4. data <- cbind(x1,x2,y)
   
5. x <- data[,1:2]
   
6. y <- data[,3]
   
7. n <- dim(x)[1]
   
8. Indices <- seq(1:n)
   
9. num <- 200  #Number of bootstrap samples
   
10. for (i in 1:num) {         
    
11. Ind_sim <- sample(Indices,100,replace = TRUE)
    
12. ##Ind_sim <- c(Ind_sim)  
    
13. xsim <- x[Ind_sim,1:dim(x)[2]]         #Construct the x data
    
14. ysim <- y[Ind_sim]                      #Construct the y data
    
15. }

复制代码

网络有问题...无法评分和回复，很感谢，有些不太明白，我琢磨琢磨，我也搜到一个类似的抽样

http://stackoverflow.com/questions/20390155/writing-a-loop-for-randomly-selecting-rows-of-a-matrix-and-doing-a-linear-regres?rq=1
抽样回归的
#set.seed(23)
A <- matrix(runif(84),ncol=2)

randco <- function(A) {
  B <- A[sample(42,12),]
  lm(B[,2] ~ B[,1])$coefficients
}

t(replicate(10,randco(A)))

#      (Intercept)        B[, 1]
# [1,]   0.6018459 -0.1643174222
# [2,]   0.4411607  0.0005322798
#...
# [9,]   0.3201649  0.4848679516
#[10,]   0.5413830  0.1850853748