R中的对数函数

Question

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

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

为了更清楚地说明我目前正在做什么，我现在的模型如下所示：

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和整数增量。