运用R语言如何将一组数据回归1000次，并将回归系数放在一个矩阵中

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

本人刚刚接触统计分析软件最近在做一个作业
请问如何用R做一千次循环，并将系数放在一个矩阵中即本题的第四问，谢谢大家

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：回归系数 R语言 distribution Independent regression

此题已经在经管爱问中有类似的解答，可以参见：
http://ask.pinggu.org/q-45170.html

jiangbeilu 发表于 2017-10-10 11:31
此题已经在经管爱问中有类似的解答，可以参见：
http://ask.pinggu.org/q-45170.html

我的问题，就是所有的估计系数都是在summary表格中，第一步要取出来，第二步再做这个矩阵，如何将第一步取出来的系数放在下面这个矩阵中？

想要飞的小菜鸟 发表于 2017-10-10 11:55
我的问题，就是所有的估计系数都是在summary表格中，第一步要取出来，第二步再做这个矩阵，如何将第一步取 ...

回归结果是一个list,里面存有系数，你可以通过查找list数据提取来学习如何从回归结果中提取出系数

dgp <- function(theta,len, mean1, sd1, lambda, sd.epsilon){
    alpha <- theta[1]
    beta1 <- theta[2]
    beta2 <- theta[3]
    x1 <- rnorm(len, mean=mean1, sd=sd1)
    x2 <- rpois(len, lambda)
    epsilon <- rnorm(len, sd=sd.epsilon)
    y <- alpha+ beta1 * x1 + beta2 * x2 + epsilon
    data.frame(Y=y,X1=x1,X2=x2)
}

estim <- function(){
    mydata <- dgp(c(1,1,-2),100,2,1,4,1)
    tmpRes <- lm(Y~X1+X2,data=mydata)
    tmpRes$coefficient
}

tmp <- replicate(1000,estim())
results <- data.frame(t(tmp))

kantdisciple 发表于 2017-10-11 17:58
dgp

thank you for your help very much!!
I will try the codes you typed~