标签：20 基础 gpa ols act 操作 hat col

OLS一元回归（GPA on act）伍德里奇案例2.3

#data input
gpa=c(2.8,3.4,3.0,3.5,3.6,3.0,2.7,3.7)#读入数据：GPA因变量
act=c(21,24,26,27,29,25,25,30)#读入数据：自变量
student=1:8#生成1到8编号
data=cbind(student,gpa,act)#按列合并，生成3列
#data=rbind(student,gpa,act)#按行合并，生成3行

#simple linear regresion
ols=lm(gpa~act)#lm是linear model，做gpa对act的回归
summary(ols)

u_hat=ols$residuals#$是抽取这一数据到新变量里，默认是列向量
gpa_hat=ols$fitted.values
beta_hat=ols$coefficients

plot(act,gpa,col=1)#plot是新建一个表格，points是在已有表格画点
                      #col指颜色
points(gpa_hat~act,col=4)#~将gpa_hat和act联系起来，否则x轴是1-8学生编号

plot(u_hat)#无~，输出x轴是1-8学生编号
sum(u_hat)#E(u|x)很小约等于0，满足OLS条件

gpa_20=t(beta_hat) %*% c(1,20)#%*%是矩阵乘法，c(1,20)指beta0*1,beta1*20,是在act=20下gpa的表现
#默认列向量，t(beta_hat)转置成行向量

两个好玩的包：
praise 和 rainbow，一个夸夸小可爱，一个彩虹小可爱

标签：20,基础,gpa,ols,act,操作,hat,col
来源： https://blog.csdn.net/g19ly/article/details/117471420

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。