young4u 大概情况如下: 数据框dd有一系列观察对象,每个观察对象都有一个uid,其变量A沿时间T变化,其趋势大致为振荡曲线,想要拟合得出这些对象的变化模型 初步的方法是尝试 nls(A~a * sin(T) + b) 算出来系数的p值倒也都是极显著 不过我想这过程中没有用到对象编号uid这一属性 也就是忽略了一部分信息量 直接把所有变量A对T作了个拟合</p> 现在想请教的是 有没有什么办法能够把观察对象uid这一信息利用起来 这个因素对分析结果比较重要 谢谢
young4u 简单的数据形式如下 <br /> dd <- data.frame(id = rep(letters[1:3], 4), A = rnorm(12), T = rep(1:4, each = 3))<br /> 每个id对应一系列随时间变化的A值 能否让id这个因子也发挥作用呢 </p>