[R-sig-ME] nlme issue

Jochen Wirsing jw1085 at wildcats.unh.edu
Wed May 10 20:09:56 CEST 2017

Dear all,

thank you very much for you good and helpful input. I would like to
provide a subset of my data, to show the problem, but don't know whether
I could (or should) just send it as an attachment.

The code I'm using is:

### Unconditional Model

uncmod <- lme(fixed = year_c~1, random = ~1|Place, data = data2, method

And the Error I'm getting is
Error in pdFactor.pdLogChol(X[[i]], ...) : object 'logChol_pd' not found

I'd really appreciate if this issue could be fixed, so that I can keep
on using R productively in my studies.

Jochen Wirsing

> R.version
platform       x86_64-pc-linux-gnu         
arch           x86_64                      
os             linux-gnu                   
system         x86_64, linux-gnu           
major          3                           
minor          4.0                         
year           2017                        
month          04                          
day            21                          
svn rev        72570                       
language       R                           
version.string R version 3.4.0 (2017-04-21)
nickname       You Stupid Darkness      


R version 3.4.0 (2017-04-21)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 16.04.2 LTS

Matrix products: default
BLAS: /usr/lib/openblas-base/libblas.so.3
LAPACK: /usr/lib/libopenblasp-r0.2.18.so

 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods  
[7] base     

other attached packages:
 [1] broom_0.4.2     nlme_3.1-131    sjPlot_2.3.1   
 [4] sjmisc_2.4.0    dplyr_0.5.0     purrr_0.2.2    
 [7] readr_1.1.0     tidyr_0.6.2     tibble_1.3.0   
[10] ggplot2_2.2.1   tidyverse_1.1.1 MASS_7.3-47    

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.10       stringdist_0.9.4.4 lubridate_1.6.0   
 [4] mvtnorm_1.0-6      lattice_0.20-35    zoo_1.8-0         
 [7] lmtest_0.9-35      assertthat_0.2.0   digest_0.6.12     
[10] psych_1.7.5        mime_0.5           R6_2.2.1          
[13] cellranger_1.1.0   plyr_1.8.4         stats4_3.4.0      
[16] coda_0.19-1        httr_1.2.1         multcomp_1.4-6    
[19] lazyeval_0.2.0     readxl_1.0.0       minqa_1.2.4       
[22] nloptr_1.0.4       Matrix_1.2-7.1     DT_0.2            
[25] splines_3.4.0      lme4_1.1-13        stringr_1.2.0     
[28] foreign_0.8-68     htmlwidgets_0.8    munsell_0.4.3     
[31] shiny_1.0.3        compiler_3.4.0     httpuv_1.3.3      
[34] modelr_0.1.0       mnormt_1.5-5       htmltools_0.3.6   
[37] nnet_7.3-12        coin_1.1-3         codetools_0.2-15  
[40] grid_3.4.0         jsonlite_1.4       arm_1.9-3         
[43] xtable_1.8-2       gtable_0.2.0       DBI_0.6-1         
[46] magrittr_1.5       scales_0.4.1       stringi_1.1.5     
[49] reshape2_1.4.2     xml2_1.1.1         effects_3.1-2     
[52] sandwich_2.3-4     TH.data_1.0-8      blme_1.0-4        
[55] tools_3.4.0        forcats_0.2.0      sjstats_0.10.0    
[58] hms_0.3            survival_2.41-3    abind_1.4-5       
[61] parallel_3.4.0     colorspace_1.3-2   rvest_0.3.2       
[64] knitr_1.15.1       haven_1.0.0        modeltools_0.2-21 
[67] merTools_0.3.0

More information about the R-sig-mixed-models mailing list