R数据分析期末作业Word文档下载推荐.docx
- 文档编号:8064481
- 上传时间:2023-05-10
- 格式:DOCX
- 页数:14
- 大小:324.18KB
R数据分析期末作业Word文档下载推荐.docx
《R数据分析期末作业Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《R数据分析期末作业Word文档下载推荐.docx(14页珍藏版)》请在冰点文库上搜索。
plot(x1,x3,xlab="
已发行股票的上期期末价值"
ylab="
当期总投资"
)
abline(lm(x3~x1))
plot(x2,x3,xlab="
上期资本存量"
abline(lm(x3~x2))
c(cor(x1,x3),cor(x2,x3))
lm1=lm(x3~x1)
lm1.summary=summary(lm1)
lm1.summary
lm2=lm(x3~x2)
lm2.summary=summary(lm2)
lm2.summary
lm=lm(x3~x1+x2)
lm.summary=summary(lm)
lm.summary
predict(lm,newdata=data.frame(x1=5593.6,x2=2226.3),interval="
prediction"
level=0.95)
运行结果:
251页1.
chen=read.table("
\\r\\251.1数据.txt"
attach(chen)
lm4=lm(x4~x5+x6,data=chen)
summary(lm4)
z<
-cbind(rep(1,length(chen[,1])),chen[,-c(1,2)])
-as.matrix(z)
eigen(t(z)%*%z)
CI<
-eigen(t(z)%*%z)$values[1]/eigen(t(z)%*%z)$values[3]
CI
install.packages("
MASS"
library(MASS)
lm5<
-lm.ridge(x4~x5+x6,data=chen)
lm5
plot(lm.ridge(x4~x5+x6,data=chen,lambda=seq(0,0.5,0.0001)))
select(lm.ridge(x4~x5+x6,data=chen,lambda=seq(0,0.5,0.0001)))
lm.ridge(x4~x5+x6,data=chen,lambda=0.0184)
3数据:
.
ting=read.table("
\\r\\252.3数据.txt"
attach(ting)
lm6=lm(x9~x7+x8,data=ting)
summary(lm6)
wj<
-resid(lm6)
ja<
-abs(wj)
summary(lm(ja~I(sqrt(x7)+sqrt(x8))))
summary(lm(ja~-1+I(sqrt(x7)+sqrt(x8))))
library(tseries)
white.test(x7+x8,x9)
4.程序:
数据:
tingan=read.table("
\\r\\253.4.txt"
attach(tingan)
lm7=lm(x12~x13)
summary(lm7)
e<
-resid(lm7)
plot(e,type="
l"
abline(h=0,col="
red"
plot(e[-1],e[-10])
orcutt"
library(orcutt)
cochrane.orcutt(lm7)
274页第3题数据
.
ganzi=read.table("
\\r\\274.3.txt"
attach(ganzi)
ganzi
lm8=lm(x15~x14)
summary(lm8)
nl.f=nls((x15)~a+b*(x14),start=list(a=11.145743,b=0.898534),trace=T)
summary(nl.f)
298页数据
\\R\\13.2.txt"
lpm=lm(ST~ARA+ASSET+ATO+ROA+GROWH+LEV+SHARE,data=ganzi)
summary(lpm)
ganz.probit=glm(ST~ARA+ASSET+ATO+ROA+GROWH+LEV+SHARE,family=binomial(link="
probit"
),data=ganzi)
summary(ganz.probit)
lmtest"
library(lmtest)
lrtest(ganz.probit)
McFa.Rsquare=function(glm.object){
deviance=glm.object$deviance
null.deviance=glm.object$null.deviance
McFa.Rsquare=1-(deviance/null.deviance)
list(McFadden.Rsquare=McFa.Rsquare)
}
McFa.Rsquare(ganz.probit)
ganz.logit=glm(ST~ARA+ASSET+ATO+ROA+GROWH+LEV+SHARE,family=binomial(link="
logit"
summary(ganz.logit)
lrtest(ganz.logit)
me=c(mean(ARA),mean(ASSET),mean(ATO),mean(ROA),mean(GROWH),mean(LEV),mean(SHARE))
pnorm(1.96,0,1)
ma=c(-16.53618,3.69994,0.66470,1.78115,-4.88107,4.10902,3.77275,0.04167)
mb=c(1,me)
mc=ma*mb
sum(mc)
pnorm(sum(mc),0,1)
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数据 分析 期末 作业