lme4 - Errors with simulate.merMod with new data from mixed model -


i'm trying simulate glmer model in lme4. i'm trying pass newdata getting different errors various formats. not understanding format function needs, package examples generate errors.

example package data:

# mixed model gm1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd),               data = cbpp, family = binomial)  # generate new data per example newdata <- with(cbpp, expand.grid(period=unique(period), herd=unique(herd)))  # predict works expected predict(gm1, newdata=newdata)                # works predict(gm1, newdata=newdata[1,])            # works  # simulate generates different errors simulate(gm1, newdata=newdata)               # works error, fine. # warning message: #   in wts - y : #   longer object length not multiple of shorter object length   simulate(gm1, newdata=newdata[1,])           # doesn't work # error in data.frame(sim_1 = c(" 4", "10"), check.names = false, row.names = "1") :  # 'row.names' should specify 1 of variables # in addition: warning message: #   in format.data.frame(x, digits = digits, na.encode = false) : #   corrupt data frame: columns truncated or padded nas  newdata2 = data.frame(   herd=c("1"),   period=c("1") ) simulate(gm1, newdata=newdata2)              # same error  simulate(gm1, newdata=newdata[1:2,])         # works error # warning message: #   in format.data.frame(x, digits = digits, na.encode = false) : #   corrupt data frame: columns truncated or padded nas 
r version 3.2.1 (2015-06-18) platform: x86_64-pc-linux-gnu (64-bit) running under: ubuntu 14.04.1 lts  other attached packages:  [1] ggplot2_1.0.0   zelig_3.5.4     boot_1.3-9      mass_7.3-29     hmisc_3.14-5     [6] formula_1.1-2   survival_2.37-7 lattice_0.20-29 plyr_1.8.1      rcurl_1.95-4.3  [11] bitops_1.0-6    lme4_1.1-10     matrix_1.1-2     loaded via namespace (and not attached):  [1] rcpp_0.12.1         cluster_1.14.4      munsell_0.4.2       colorspace_1.2-4     [5] minqa_1.2.4         stringr_0.6.2       tools_3.2.1         nnet_7.3-7           [9] matchit_2.4-21      gtable_0.1.2        nlme_3.1-113        latticeextra_0.6-26 [13] digest_0.6.4        reshape2_1.4        rcolorbrewer_1.0-5  nloptr_1.0.4        [17] acepack_1.3-3.3     rpart_4.1-5         labeling_0.3        scales_0.2.4        [21] foreign_0.8-59      proto_0.3-10

this fixed in development version lme4. many ben bolker.

devtools::install_github("lme4/lme4")


Comments

Post a Comment

Popular posts from this blog

html - Outlook 2010 Anchor (url/address/link) -

i'm trying link specific pdf page msoutlook'10 or msword'10, #page=... anchor particularly useful , works fine elsewhere. unfortunately when click [ctr-k] in outlook or winword, links re-formatted assuming unc path. pdf file located on shared drive. using "file:" prefix links open pdf in adobe reader (outside of browser), , "http:" open inside browser. issue seem need open in browser "file:" prefix in order anchor work. know anchors work inside browser. outlook uses windows/explorer default program (based on file extension) seems open in adobe reader neglecting anchor. is aware of solution, link specific page or bookmark of shared-pdf ms outlook and/or winword link? common references: mailermailer.com/...anchor-tags-html-emails campaignmonitor.com/...anchor-links-in-email-newsletters/ stackoverflow.com/...anchor-in-pdf-document
Read more

javascript - Why does running this loop 9 times take 100x longer than running it 8 times? -

consider code: test = function() { } t = new test(); (var = 0; < 8; i++) { result = t instanceof test; } if change number of iterations 8 9 , loop take 100 times longer complete in latest version of firefox (41.0.1). tested on 2 different pcs , magic limit 8. here jsperf test used: http://jsperf.com/instanceof-8-times-vs-9-times does have idea why might happen? seems specific instanceof . not happen if else object, example check property. note: filed bugzilla bug this. jan de mooij mozilla team has posted details in bugzilla thread . here's simplistic interpretation of highly technical answers: in i < 8 case, firefox smart enough hoist result = t instanceof test; statement out of loop (according tests, doesn't seem omit altogether). in i < 9 case, apparently doesn't optimization . why? reason not clear, related fact 9 iterations threshold above function considered "hot" enough run through jit compiler. i < 8 cas...
Read more

Getting gateway time-out Rails app with Nginx + Puma running on Digital Ocean -

following this guide setup app on digital ocean droplet, not working. site not loading , it's throwing: "gateway time-out". the puma server works, since ps aux | grep puma prints: root 15878 0.0 0.4 376832 2168 ? sl 12:50 0:00 puma 2.13.4 (unix:///var/www/myapp/shared/tmp/sockets/myapp-puma.sock) [20150930164940] root 15881 0.0 0.7 903676 3692 ? sl 12:50 0:06 puma: cluster worker 0: 15878 [20150930164940] root 24654 0.0 0.1 11748 916 pts/2 s+ 18:00 0:00 grep --color=auto puma my nginx config: user root; worker_processes 1; events { worker_connections 1024; } http { upstream puma { server unix:///var/www/myapp/shared/tmp/sockets/myapp-puma.sock; } server { listen 80; server_name myapp.com; root /var/www/myapp/current/public; access_log /var/www/myapp/current/log/access.log; error_log /...
Read more