Question

从下面的错误信息中我猜测“步骤”功能不适用于mle，这是一个S4类？

我可以通过AIC为lm选择最佳模型的步进功能，它说它也适用于glm，但到目前为止我只看到了非检测和数据分析（NADA）包中的mle我正在使用。大约20％的结果数据（W. Nevada地下水中的砷）低于检测限。

“step”代码比我见过的用于逐个比较AIC模型的其他代码要简单得多，但我想我会看看后者是否适用于NADA包。

NADA软件包提供x2测试来比较模型，但由于某些原因，我的所有模型都有巨大的x2，p = 0，所以我不知道如何比较这种方式（例如下面的例子）。但是，我对x2测试，统计数据或R一般都不是很熟悉。 NADA作者提到在“部分对数可能性”上选择最佳模型。

谢谢！

我的代码的执行情况如此：

install.packages("NADA")
library(NADA)
As_nada_mle=with(clean2503, cenreg(Cen(As001, nd_epa)~Elev_m + Tahoe + Walker + UpperCarson + Tcan2car + TruckeeTribOrMain))
Best_nada=step(As_nada_mle)
##Error: $ operator not defined for this S4 class

我的NADA包装CENREG输出示例：

                     Value Std. Error       z        p
(Intercept)       11.03589     4.1515  2.6583 7.85e-03
Elev_m            -0.00825     0.0013 -6.3453 2.22e-10
Tahoe              0.29773     4.8417  0.0615 9.51e-01
Walker             1.36853     3.9538  0.3461 7.29e-01
UpperCarson       -2.13845     3.9506 -0.5413 5.88e-01
Tcan2car           2.56088     3.9238  0.6527 5.14e-01
TruckeeTribOrMain -2.98352     4.0807 -0.7311 4.65e-01
Log(scale)         1.36526     0.0171 79.7684 0.00e+00

Scale = 3.92 

Log Normal distribution

Loglik(model)= -12303.8   Loglik(intercept only)= -12982.3 
Loglik-r:  0.64693 

Chisq= 1357.07 on 6 degrees of freedom, p= 0 
Number of Newton-Raphson Iterations: 5 
n = 2503

Answer 1

我不具备您正在尝试的技术经验。同那个警告，这是我的想法：

似乎难以尝试运行步骤包装器你的cenreg命令。由于步骤只是一个方便的算法自动选择最好的AIC，这不是绝对必要的。从而另一种方法是制作所有可能模型的表格你想要考虑（这些应该是机械上可信的），和确定每个AIC和/或日志的可能性，然后进行比较他们并挑选最好的。这通常在＆＃34;信息中完成理论＆＃34;使用AIC的方法。这种方法有一种标准的方法这样做 - 如果你需要，我可以给你发一些例子。这也是在Burnham KP，Anderson DR中概述。 2002年。模型选择和多模推理：一种实用的信息理论方法。第2 编辑。柏林，德国：施普林格。（伯克利将此作为电子书）。一个优点是你可以看到哪些参数正在进行中似乎在各种情况下都很重要。这是唯一的情况如果你比较几十个和几十个，那真是不切实际不同种类的模特。
x2如此重要的原因可能是因为它是整体模型，您可能有大样本和几个适合的术语。我假设x2以某种方式进行比较两个对数可能性，我们看到之间的差异： Loglik（型号）= -12303.8 Loglik（仅拦截）= -12982.3
你可能意识到这一点，但现在有很多模型中的非重要术语（p> 0.1）：如果您运行另一个术语不包括它们的模型，你的logliklihood将会改善很多（即获得少负面）。
也许Dennis Helsel关于审查数据的教科书会有所帮助解释如何使用cenreg？

人们不像我那么无知，他们可以纠正或改进我的建议......

如何将模型与R中的NADA包进行比较？

1 个答案: