1

希望对于那些比我更有经验的人来说是一个相对容易的！

尝试使用以下代码执行 Box-Cox 转换：

fit <- lm(ABOVEGROUND_BIO ~ TREATMENT * P_LEVEL, data = MYCORRHIZAL_VARIANCE)
bc <- boxcox(fit)
lambda<-with(bc, x[which.max(y)])
MYCORRHIZAL_VARIANCE$bc <- ((x^lambda)-1/lambda)
boxplot(bc ~ TREATMENT * P_LEVEL, data = MYCORRHIZAL_VARIANCE)

但是，当我运行它时，我收到以下错误消息：

Error: object 'x' not found. (on line 4)

对于上下文，这是我的数据集的 str ：

Classes ‘spec_tbl_df’, ‘tbl_df’, ‘tbl’ and 'data.frame':    24 obs. of  14 variables:
 $ TREATMENT             : Factor w/ 2 levels "Mycorrhizal",..: 1 1 1 1 1 1 1 1 1 1 ...
 $ P_LEVEL               : Factor w/ 2 levels "Low","High": 1 1 1 1 1 1 2 2 2 2 ...
 $ REP                   : int  1 2 3 4 5 6 1 2 3 4 ...
 $ ABOVEGROUND_BIO       : num  7.5 6.8 5.3 6 6.7 7 12 12.7 12 10.2 ...
 $ BELOWGROUND_BIO       : num  3 2.4 2 4 2.7 3.6 7.9 8.8 9.5 9.2 ...
 $ ROOT_SHOOT            : num  0.4 0.35 0.38 0.67 0.4 0.51 0.66 0.69 0.79 0.9 ...
 $ ROOT_SHOOT.log        : num  -0.916 -1.05 -0.968 -0.4 -0.916 ...
 $ ABOVEGROUND_BIO.log   : num  2.01 1.92 1.67 1.79 1.9 ...
 $ ABOVEGROUND_BIO.sqrt  : num  2.74 2.61 2.3 2.45 2.59 ...
 $ ABOVEGROUND_BIO.cubert: num  1.96 1.89 1.74 1.82 1.89 ...
 $ BELOWGROUND_BIO.log   : num  1.099 0.875 0.693 1.386 0.993 ...
 $ BELOWGROUND_BIO.sqrt  : num  1.73 1.55 1.41 2 1.64 ...
 $ BELOWGROUND_BIO.cubert: num  1.44 1.34 1.26 1.59 1.39 ...
 $ TOTAL_BIO             : num  10.5 9.2 7.3 10 9.4 10.6 19.9 21.5 21.5 19.4 ...
 - attr(*, "spec")=
  .. cols(
  ..   TREATMENT = col_factor(levels = c("Mycorrhizal", "Non-mycorrhizal"), ordered = FALSE, include_na = FALSE),
  ..   P_LEVEL = col_factor(levels = c("Low", "High"), ordered = FALSE, include_na = FALSE),
  ..   REP = col_integer(),
  ..   ABOVEGROUND_BIO = col_number(),
  ..   BELOWGROUND_BIO = col_number(),
  ..   ROOT_SHOOT = col_number()
  .. )

我知道在 MYCORRHIZAL_VARIANCE 数据集中没有名为 bc 的变量，但我只是按照在执行 Box-Cox 时给我的基本说明进行操作，我想我对实际上应该将“x”表示为什么感到困惑，因为我以为'x'是在第3行定义的？有关如何解决此错误的任何建议？

提前致谢！

4

1 回答 1

0

我以为'x'是在第3行定义的？

第 3 行是lambda<-with(bc, x[which.max(y)]). 它不定义x，它定义lambda。它确实使用，它 在环境中x寻找它。如果你从包中使用，确实应该包含和组件，所以不应该给你同样的错误信息。我预计替换长度会出错。因为...bcboxcox()MASSbcxybc$x

bc$x是尝试的潜在 lambda 值boxcox- 您使用的是默认值seq(-2, 2, 1/10)，如果您的数据有 41 行的倍数，在将 41 个值分配给新列时不会出错，这不太可能是巧合。

第 3 行挑选出使可能性最大化的 lambda 值，因此您不应该再需要其余的值bc。我希望您使用该lambda值来转换您的响应变量，因为这就是 Box Cox 转换的用途。((x^lambda)-1/lambda)没有任何统计或程序意义。改用这个：

MYCORRHIZAL_VARIANCE$bc <- (MYCORRHIZAL_VARIANCE$ABOVEGROUND_BIO ^ lambda - 1) / lambda

（请注意，我还更正了括号。你想要(y ^ lambda - 1) / lambda，而不是(y ^ lambda) - 1 / lambda。）

于 2020-05-13T02:24:21.080 回答