全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
1269 1
2014-11-29
有三列数据(如下):月份、SVL和性腺体积,我想通过线性模型分析不同月份性腺体积的变化情况(排除svl对性腺体积的影响,将svl作为协变量),请问如何做?
(最初的分析方法:通过对svl与volume的回归分析,发现svl对volume有影响,为了剔除这种影响,通过ln(volume)/ln(svl)的方法转化得到一个值,然后再通过方差分析每个月这个值的差异,但是老板认为分析方法太复杂,让我用最简单的线性模型去分析)
Month Body length(mm)(SVL) testis volume 1
4 57 48.131
4 57.75 43.854
4 59.05 41.005
4 56.75 50.697
4 57.44 41.729
4 56.41 43.593
4 60.35 29.609
4 56.28 59.218
4 57.68 42.439
4 58.82 36.616
5 55.67 40.301
5 61.81 25.661
5 57.01 38.136
5 57.99 37.472
5 58.29 35.794
5 59.72 32.076
5 58.62 34.215
5 59.07 34.899
5 57.04 39.478
5 58.11 36.294
5 59.94 31.691
5 53.74 49.052
5 54.65 42.729
5 60.34 28.628
5 60.07 31.649
5 55.85 39.478
6 57.9 15.545
6 57.3 17.321
6 63.32 14.002
6 58.5 15.377
6 55.1 17.502
6 56.15 18.423
6 60.06 14.804
6 61.28 14.377
7 59.04 25.296
7 57.98 27.003
7 55.79 27.316
7 52.71 29.675
7 60.96 22.703
7 59.49 23.687
7 53.77 28.306
7 61.44 22.029
二维码

扫码加我 拉你入群

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

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

全部回复
2014-11-30 11:25:53
Month SVL volume
4 57 48.131
4 57.75 43.854
4 59.05 41.005
4 56.75 50.697
4 57.44 41.729
4 56.41 43.593
4 60.35 29.609
4 56.28 59.218
4 57.68 42.439
4 58.82 36.616
5 55.67 40.301
5 61.81 25.661
5 57.01 38.136
5 57.99 37.472
5 58.29 35.794
5 59.72 32.076
5 58.62 34.215
5 59.07 34.899
5 57.04 39.478
5 58.11 36.294
5 59.94 31.691
5 53.74 49.052
5 54.65 42.729
5 60.34 28.628
5 60.07 31.649
5 55.85 39.478
6 57.9 15.545
6 57.3 17.321
6 63.32 14.002
6 58.5 15.377
6 55.1 17.502
6 56.15 18.423
6 60.06 14.804
6 61.28 14.377
7 59.04 25.296
7 57.98 27.003
7 55.79 27.316
7 52.71 29.675
7 60.96 22.703
7 59.49 23.687
7 53.77 28.306
7 61.44 22.029
以上的数据存在一个名为data2.txt下。
data2<-read.table(file="D:/R workplace/data2.txt",sep=" ",header=TRUE)

plot(data2[,2],data2[,3])
abline(lm(formula=volume~SVL,data=data2))

sum(data2[,1]==4)
sum(data2[,1]==5)
sum(data2[,1]==6)
sum(data2[,1]==7)

plot(data2[1:10,2],data2[1:10,3])
abline(lm(formula=volume~SVL,data=data2[1:10,]))

plot(data2[11:26,2],data2[11:26,3])
abline(lm(formula=volume~SVL,data=data2[11:26,]))

plot(data2[27:34,2],data2[27:34,3])
abline(lm(formula=volume~SVL,data=data2[27:34,]))

plot(data2[35:42,2],data2[35:42,3])
abline(lm(formula=volume~SVL,data=data2[35:42,]))

用以上的代码可以看出,按月做线性回归的话,还是比整体好。

你可以分月去做。 具体的信息,你用
summary(lm(formula=volume~SVL,data=data2[35:42,]))去看。
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

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