如何获得“nls”对象的漂亮乳胶输出

我需要“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 对象的乳胶输出的绝望帮助。

有关我如何继续的任何指示。

regards

K

对于 LaTeX 中的自定义输出表，我的经验是xtable套餐是最方便、最灵活的解决方案。

让我们考虑一下您的情况nlsobject 和 stargazer 不支持它（也许有一天会支持？）。

### Estimate model
model <- nls(Y ~ <your function>(theta1, theta2, theta3))

### generate various parts for output
sum_model = summary(model)
mat_model = sum_model$coefficients

### generate coefficients, se, t-stat and p values 
df_model = as.data.frame(mat_model)
colnames(df_model) <- c("Coef.", "Std. error", "t-stat.", "p")

mat <- data.frame(t(df_model))
tbl <- xtable(mat)
print(tbl, only.contents=TRUE, include.rownames=T, 
      include.colnames=F, floating=F,
      hline.after=NULL,
      file="summary_nls.tex")

现在请注意，您可以使用 R 包根据您的需求定制输出xtable。例如，我省略了整个 LaTeXtabular一直无法摆脱的环境stargazer我喜欢用threeparttable与 LaTeX 包中的命令一起打包booktabs。对于给定的玩具模型，简洁的输出如下所示：

\begin{table}[t]
\centering
\begin{threeparttable}
\captionabove{Regression results for nonlinear model.}
\begin{tabular}{lccc}
      \toprule
       & Theta1 & Theta2 & Theta3\\
       & $\theta 1$ & $\beta$ & $\sigma$ \\
      \midrule
      \input{summary_nls}
      \bottomrule
   \end{tabular}
\label{tab:nls_summary}
\end{threeparttable}
\end{table}

我将上面的代码保留为模板，只为更大、更复杂的模型添加行和列。搭配knitr这对我来说就像一个魅力。