我需要一个“nls”对象的 LaTex 表示。不幸的是,stargazer 不支持这种对象类型。
网上的一些研究让我找到了 nls2 库中的 as.lm.nls 函数。它声称将 nls 对象转换为相应的 lm 对象。就我而言,它失败了。以下是示例输出:
> DNase1 <- subset(DNase, Run == 1)
> xx <- nls(density ~ SSlogis(log(conc), Asym, xmid, scal), DNase1)
> summary(xx)
Formula: density ~ SSlogis(log(conc), Asym, xmid, scal)
Parameters:
Estimate Std. Error t value Pr(>|t|)
Asym 2.34518 0.07815 30.01 2.17e-13 ***
xmid 1.48309 0.08135 18.23 1.22e-10 ***
scal 1.04146 0.03227 32.27 8.51e-14 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Residual standard error: 0.01919 on 13 degrees of freedom
Number of iterations to convergence: 0
Achieved convergence tolerance: 3.281e-06
下面的 as.lm.nls 输出与上面的实际输出不匹配:
> library(nls2)
Loading required package: proto
> xlm = as.lm.nls(xx)
> summary(xlm)
Call:
lm(formula = density ~ Asym + xmid + scal - 1, offset = fitted(xx))
Residuals:
Min 1Q Median 3Q Max
-0.033513 -0.012931 -0.001454 0.009699 0.038137
Coefficients:
Estimate Std. Error t value Pr(>|t|)
Asym -8.878e-07 7.815e-02 0 1
xmid -9.328e-07 8.135e-02 0 1
scal -3.751e-07 3.227e-02 0 1
Residual standard error: 0.01919 on 13 degrees of freedom
Multiple R-squared: 0.9996, Adjusted R-squared: 0.9995
F-statistic: 1.153e+04 on 3 and 13 DF, p-value: < 2.2e-16
由于 nls 摘要输出与 lm 摘要输出非常相似,我假设 nls 到 lm 对象的内容等价物应该可以让 stargazer 发挥其魅力。
我有相当多的 nls 模型要包含在报告中,而 as.lm.nls 的失败让我陷入了困境。我需要关于 nls 对象的乳胶输出的绝望帮助。
关于如何进行的任何指针。
问候
ķ