楼主: 盛夏果123
6772 2

[问答] R语言用nlme包的lme函数怎么加入权重 [推广有奖]

  • 0关注
  • 1粉丝

小学生

42%

还不是VIP/贵宾

-

威望
0
论坛币
0 个
通用积分
1.0032
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
46 点
帖子
3
精华
0
在线时间
10 小时
注册时间
2016-3-16
最后登录
2016-9-13

10论坛币
最近在学习Mixed model 的理论和在R中的运用,但是碰到了一些问题,我的数据是个体*环境*区组(200*2*3)有关基因与环境互作的。如下表所示 1111.JPG ,我的分析思路是这样的,
1.首先分单个环境逐一分析,每个环境下的数据就做简单的区组、个体两项方差分析,得出两个环境各自的方差,作为权重因子进入下一步的计算;
2.接下来进行个体与环境的mixed model的分析,此时个体用的就是平均值了。
具体思路是这样的,但是代码可能有些问题,具体代码如下 11111.JPG ,我想问是不是权重的加入搞错了,因为翻阅了一些英文的教程书,也没详细写怎么加入权重,只是很粗略的写weights = varFixed(~age)这种根本不知道它在说什么玩意的东西。。 很是头大!麻烦高手们有空帮忙看看!!!谢谢谢谢谢谢!!困扰我好久好久了!!

关键词:nlme LME R语言 nlm mixed model 下一步 平均值 英文
沙发
盛夏果123 发表于 2016-3-23 17:46:20 |只看作者 |坛友微信交流群
具体的说我不知道在这个函数中如何加入权重,上面我的方法,所谓的加了一列叫env.w只是我自己YY出来的。。。
谁可以帮忙看看啊~~~~[cry][cry]

使用道具

藤椅
吴熹微 发表于 2021-9-12 16:25:28 |只看作者 |坛友微信交流群
盛夏果123 发表于 2016-3-23 17:46
具体的说我不知道在这个函数中如何加入权重,上面我的方法,所谓的加了一列叫env.w只是我自己YY出来的。。。 ...
请问这个问题你解决了吗?gls函数里的参数我也不会,如果你回了能不能分享一下呢

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-5-6 11:08