R中的对数函数

Question

我目前正在研究使用统计软件R的lme包的增长模型。我已经建立了具有线性和poly(Time,2)趋势的模型，而后者更适合。另外，我想测试一个对数趋势，但不知道R中的函数是什么样的或如何变换它，因为大多数指南只涵盖线性，二次，立方趋势但不包括对数趋势。

因此，如果有人可以帮助我并为我提供对数趋势的功能，我将非常感激。

为了更清楚地说明我正在做什么 - ＆gt;我目前的模特看起来像：

model.linear<-lme(DV~1+Time,random=~1|Subnum,data=dataset,na.action=na.omit,control=list(opt="optim"))

model.quadratic<-lme(DV~poly(Time,2),random=~1|Subnum,data=dataset,na.action=na.omit,control=list(opt="optim"))

Answer 1

我认为你在寻找：

model.logistic<-lme(DV~log1p(Time),random=~1|Subnum,data=dataset,na.action=na.omit,control=list(opt="optim"))

函数log1p()在记录之前为每个观察添加1，这适用于计数变量或其他下限为0和整数增量的变量。