全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
14621 14
2017-08-05
数据结构如图

数据结构 第一列 代表不同的人,图中只显示了前两个人。每个人有相应的最后两列数据。我想对每个人的afp.value和time.relative进行回归,以time.relative为自变量,afp.value为因变量:
   afp.value=a+b*time.relative
然后每个人都会产生一对(a,b).
鉴于我有很多个人,要是把数据按每个人拆分,再分别回归肯定是工作量太大的。
数据文件如下:
data.zip
大小:(13.23 KB)

 马上下载


如何快速得到我想要的结果呢?
提前谢谢!
二维码

扫码加我 拉你入群

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

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

全部回复
2017-8-5 10:31:56
可惜没有论坛币下载你的数据
二维码

扫码加我 拉你入群

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

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

2017-8-5 11:34:22
###求线性回归的系数
xishu=function(a){
  attach(a)
  ###a[,1]表示数据名字的那一列,m就表示名字
  jilu1=a[which(a[,1]=="m"),]
  ###a[,3],a[,2]表示因变量与自变量
  jilu1$a1=jilu1[,3]
  jilu1$b1=jilu1[,2]
fit=lm(a1~b1,data=jilu1)
jilu2=a[which(a[,1]=="f"),]
jilu2$a1=jilu2[,3]
jilu2$b1=jilu2[,2]
fit1=lm(a1~b1,data=jilu2)
  return(c(fit$coefficients,fit1$coefficients))
  }   你试试,如果可以的话  别忘了奖励
二维码

扫码加我 拉你入群

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

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

2017-8-5 19:25:35
小R小二郎 发表于 2017-8-5 10:31
可惜没有论坛币下载你的数据
不需要论坛币啊。试了下好像有点问题诶。。In model.response(mf, "numeric") :
  using type = "numeric" with a factor response will be ignored
二维码

扫码加我 拉你入群

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

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

2017-8-5 22:17:30
D3JO)42(]B8YB(IS0Y21_3B.png     这只是一个思路  因为这个分析里只用到"aetiology        AFP.value        Time.relative
",这三列,所以我只提取了这三列的数据
###求线性回归的系数
xishu=function(a){
  
  ###a[,1]表示数据名字的那一列,m就表示名字
  jilu1=a[which(a[,1]=="HEP C"),]
  ###a[,3],a[,2]表示因变量与自变量
  fit=lm(jilu1[,3]~jilu1[,2],data=jilu1)
  #jilu2=a[which(a[,1]=="ALD + NAFLD"),]
  #jilu2$a1=jilu2[,3]
  #jilu2$b1=jilu2[,2]
  #fit1=lm(a1~b1,data=jilu2)
  return(fit$coefficients)
}     这只是一个思路,就是按人名划分,有多少人,你在往函数里面添加就行,最后修改返回值。
希望能帮到你
二维码

扫码加我 拉你入群

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

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

2017-8-5 23:24:23
小R小二郎 发表于 2017-8-5 22:17
这只是一个思路  因为这个分析里只用到"aetiology        AFP.value        Time.relative
",这三列,所以我只提取了 ...
好的谢谢!不过人实在很多,数据很多,不可能一行一行添加哈~
谢谢!
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群