Chinaunix首页 | 论坛 | 博客
  • 博客访问: 628784
  • 博文数量: 149
  • 博客积分: 3901
  • 博客等级: 中校
  • 技术积分: 1558
  • 用 户 组: 普通用户
  • 注册时间: 2009-02-16 14:33
文章分类

全部博文(149)

文章存档

2014年(2)

2013年(10)

2012年(32)

2011年(21)

2010年(84)

分类: LINUX

2010-08-18 11:35:35


非常有用 !!

R 语言 >>>
x<-c(1,2, 3, 4, 5, 6, 7, 8, 9, 10)
y<-c(0,30,44,57,65,78,68,97,99,120)
lm.reg<-lm(y ~ 1+x)

point<-data.frame(x=11)
lm.pred<-predict(lm.reg,point,interval="prediction",level=0.95)
lm.pred[1]

输出 : 128.26666666666668




python 语言 >>>
import rpy2.robjects as robjects
r = robjects.r

robjects.globalenv["x"] = robjects.FloatVector([1,2, 3, 4, 5, 6, 7, 8, 9, 10])
robjects.globalenv["y"] = robjects.FloatVector([0,30,44,57,65,78,68,97,99,120])
lm_d9 = r.lm("y ~ 1+x")
r.predict( lm_d9,
           robjects.r(" data.frame(x=11) "),
           interval="prediction",level=0.95)[0]

输出 : 128.26666666666668



阅读(3116) | 评论(0) | 转发(0) |
给主人留下些什么吧!~~