r语言 - 如何使用Box cox传输数据?

Question

我有一个数据如下。

数据：

temperature stage Replicate week Nematode_Number
T20   All         1    1            60.0
T20   All         2    1            72.7 
T20   All         3    1            69.3
T20   All         4    1            45.3
T20   All         5    1            40.7

我喜欢找到线性模型，然后根据 Box cox 对这些数据进行归一化数据，但有一个错误。

Script: 
> mode (all)
[1] "list"
> is.data.frame (all)
[1] TRUE
> getClass(class(all))
Class "data.frame" [package "methods"]
Slots:

Name:                .Data               names           row.names            .S3Class
Class:                list           character data.frameRowLabels           character
Extends: 
Class "list", from data part
Class "oldClass", directly
Class "vector", by class "list", distance 2
> a <- boxcox(lm(Nematode_Number ~ week, data = all))
Error in terms.formula(formula, data = data) : 'data' argument is of the wrong type

Answer 1

data

df <- data.frame(temperature=rep("T20",5),
stage=rep("All",5),
Replicate=1:5, week=rep(1,5),
Nematode_Number=c(60.0, 72.7, 69.3, 45.3, 40.7),
stringsAsFactors=F)

这是你想要的吗？

这行得通

lm(Nematode_Number ~ week, data=df)

这行得通

library(MASS)
boxcox(Nematode_Number ~ week, data=df)

这行得通

temp <- lm(Nematode_Number ~ week, data=df)
boxcox(temp, data=df)

Answer 2

boxcox应该可以工作。好吧，boxcox也接受公式对象。因此，如果可以的话，请尝试清除您的环境，或者rm(all)并再次重读它。然后尝试下面的代码，如果此代码在上述之前根本不起作用。

a=boxcox(Nematode_Number~week,data=all)

希望这有帮助