问题描述
我的问题遵循这个简洁的介绍:
我有一个数据框,包括:每个物种的物种、性状 X 值(均值)和性状 X 的 SE;以及物种分支长度的系统发育。
我想知道是否有可能通过系统发育关系和 SE 来创建“正确”的性状 X 分布的残差??
我的基本原理是这样一个事实,即当特征 X 回归时(例如通过体型),可以创建包含 SE 和系统发育关系的残差,因此使用类似:
使用 nlme 中的 gls:
gls(trait X ~ body size,correlation=corPagel(1,phy=tree),data=data,method = "ML",weights= SE)
哪里SE <- varFixed (~1/sqrt(traiX_SE))
但是,如果我的性状 X 不能通过体型回归,我如何拯救由系统发育关系和 SE 校正的残差?有意义吗?
提前致谢!
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)