R: Pander sink stack full when printing summary lm

白昼怎懂夜的黑 提交于 2019-12-08 08:52:15

问题


I am in the middle of generating a HTML report in Rstudio via pandoc for a collaborator. However pander is hitting the sink limit in R when trying to generate the output for the following summary of a lm() object.

My R instance:

version _
platform x86_64-apple-darwin13.1.0
arch x86_64
os darwin13.1.0
system x86_64, darwin13.1.0
status
major 3
minor 1.0
year 2014
month 04
day 10
svn rev 65387
language R
version.string R version 3.1.0 (2014-04-10) nickname Spring Dance


The lm I am trying to generate the HTML from:

summary(a.lm)

Call:
lm(formula = f.t ~ ((b.c + I(a.c^2) + b.t.c + 
    c.c) * (PTGENDER + DX.bl + Al1.Al2)) + BMI.t.c + FH.t.c + 
    Hb.t.c + age.c + I(age.c^2), data = subset(d, 
    !is.na(Hb.t.c) & !is.na(BMI.t.c) & !is.na(b.t.c) & !is.na(c.c)))

Residuals:
     Min       1Q   Median       3Q      Max 
-0.72042 -0.15865 -0.01311  0.15826  0.91421 

Coefficients:
                             Estimate Std. Error t value Pr(>|t|)    
(Intercept)                 1.6247957  0.1120717  14.498  < 2e-16 ***
a.c                         0.1343279  0.1481459   0.907  0.36541    
I(a.c^2)                   -0.0312254  0.0537683  -0.581  0.56193    
b.t.c                       0.1040226  0.0782932   1.329  0.18516    
c.c                         0.2012898  0.0818699   2.459  0.01461 *  
PTGENDERFemale              0.0067323  0.0473586   0.142  0.88707    
DX.blLMCI                  -0.0121387  0.0553348  -0.219  0.82654    
DX.blAD                     0.0209534  0.0769331   0.272  0.78557    
Al1.Al23:3                  0.1658594  0.1075438   1.542  0.12426    
Al1.Al23:4                  0.3506150  0.1121648   3.126  0.00198 ** 
Al1.Al24:4                  0.1386058  0.1690946   0.820  0.41316    
BMI.t.c                    -0.0280475  0.0182582  -1.536  0.12574    
FH.t.c                      0.0031254  0.0225851   0.138  0.89005    
Hb.t.c                      0.0162085  0.0171502   0.945  0.34551    
age.c                       0.0024730  0.0028751   0.860  0.39052    
I(age.c^2)                  0.0003961  0.0002926   1.354  0.17706    
a.c:PTGENDERFemale          0.0303573  0.0474341   0.640  0.52276    
a.c:DX.blLMCI               0.1143316  0.0603606   1.894  0.05934 .  
a.c:DX.blAD                 0.1406856  0.0691761   2.034  0.04302 *  
a.c:Al1.Al23:3             -0.0369030  0.1472692  -0.251  0.80234    
a.c:Al1.Al23:4             -0.0493593  0.1469924  -0.336  0.73730    
a.c:Al1.Al24:4             -0.0688362  0.1564964  -0.440  0.66041    
I(a.c^2):PTGENDERFemale    -0.0584777  0.0324483  -1.802  0.07270 .  
I(a.c^2):DX.blLMCI          0.0298364  0.0368905   0.809  0.41940    
I(a.c^2):DX.blAD            0.0147823  0.0440503   0.336  0.73747    
I(a.c^2):Al1.Al23:3         0.0001684  0.0586422   0.003  0.99771    
I(a.c^2):Al1.Al23:4        -0.0357507  0.0621727  -0.575  0.56578    
I(a.c^2):Al1.Al24:4        -0.0328713  0.0664954  -0.494  0.62149    
b.t.c:PTGENDERFemale        0.0208930  0.0466606   0.448  0.65470    
b.t.c:DX.blLMCI            -0.0640037  0.0628767  -1.018  0.30968    
b.t.c:DX.blAD              -0.0607185  0.0705763  -0.860  0.39042    
b.t.c:Al1.Al23:3           -0.0428761  0.0727674  -0.589  0.55624    
b.t.c:Al1.Al23:4           -0.0289522  0.0729465  -0.397  0.69178    
b.t.c:Al1.Al24:4           -0.0946543  0.0981224  -0.965  0.33564    
c.c:PTGENDERFemale          0.0193275  0.0445454   0.434  0.66474    
c.c:DX.blLMCI              -0.0371738  0.0470631  -0.790  0.43034    
c.c:DX.blAD                -0.1221799  0.0722145  -1.692  0.09189 .  
c.c:Al1.Al23:3             -0.1604799  0.0798383  -2.010  0.04548 *  
c.c:Al1.Al23:4             -0.2015401  0.0871098  -2.314  0.02149 *  
c.c:Al1.Al24:4             -0.3374309  0.1520609  -2.219  0.02737 *  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Residual standard error: 0.2795 on 254 degrees of freedom
Multiple R-squared:  0.4536,    Adjusted R-squared:  0.3697 
F-statistic: 5.406 on 39 and 254 DF,  p-value: < 2.2e-16

running:

pander(summary(a.lm))

leads to the following error

Error in sink(file) : sink stack is full

Any suggestions to overcome this would be much appreciated.

Cheers


In response to Mr Flick

Yes it seems to:

Here is an example:

d <- data.frame(a = rnorm(20), b = sample(letters[1:2], 20, replace = TRUE))
pander(summary(lm(a ~ b, data = d)))
Error in sink(file) : sink stack is full
pander(lm(a ~ b, data = d))
--------------------------------------------------------------
     &nbsp;        Estimate   Std. Error   t value   Pr(>|t|) 
----------------- ---------- ------------ --------- ----------
     **bb**        -0.1637      0.3797     -0.4311    0.6715  

 **(Intercept)**   -0.1634      0.2401     -0.6805    0.5049  
--------------------------------------------------------------

Table: Fitting linear model: a ~ b

So dropping the summary call works :)


回答1:


Actually, if you're just running a standard lm regression, don't pass the summary.lm object to pander, pass the model itself (the lm object)

pander(a.lm)

That should have the information you need. The function calls summary() internally.



来源:https://stackoverflow.com/questions/24359744/r-pander-sink-stack-full-when-printing-summary-lm

标签
易学教程内所有资源均来自网络或用户发布的内容,如有违反法律规定的内容欢迎反馈
该文章没有解决你所遇到的问题?点击提问,说说你的问题,让更多的人一起探讨吧!