[R-sig-ME] nlme issue

[Jeff Newmiller]

Read ?dput and use it to create an R code that will recreate your data in our R environment. 
-- 
Sent from my phone. Please excuse my brevity.

On May 10, 2017 11:09:56 AM PDT, Jochen Wirsing <jw1085 at wildcats.unh.edu> wrote:
>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
>library(nlme)
>library(broom)
>
>uncmod <- lme(fixed = year_c~1, random = ~1|Place, data = data2, method
>="ML")
>tidy(uncmod)
>summary(uncmod)
>varCorr(uncmod)
>
>
>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.
>
>Best,
>Jochen Wirsing
>
>
>
>______________
>> R.version
>               _                           
>platform       x86_64-pc-linux-gnu         
>arch           x86_64                      
>os             linux-gnu                   
>system         x86_64, linux-gnu           
>status                                     
>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
>
>locale:
> [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
> [3] LC_TIME=en_US.UTF-8        LC_COLLATE=en_US.UTF-8    
> [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
> [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
> [9] LC_ADDRESS=C               LC_TELEPHONE=C            
>[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=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
>
>_______________________________________________
>R-sig-mixed-models at r-project.org mailing list
>https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models