[R-sig-ME] glmmADMB troubles
Ben Bolker
bbolker at gmail.com
Thu Feb 28 05:05:16 CET 2013
Ross Boylan <ross at ...> writes:
>
> It seems it's one thing after another. Ben suggested glmmADMB for
> poisson--actually for truncated poisson, but for now I'm just trying to
> get anything going. Some of my count variables have 1, not 0 as a
> minimum, and some have a maximum too, e.g., all values >10 are reported
> as 10.
If you have enough data to make it worthwhile, you might try the
'ordinal' package ...
>
> I was able to install it from r-forge, with a warning that it was build
> for 2.15.3. But ...
>
> > r <- glmmadmb(sexActs~(1|id), sexpartner)
> Error in II[, ii] = II[, ii] + REmat$codes[[i]] :
> number of items to replace is not a multiple of replacement length
> In addition: Warning messages:
> 1: In glmmadmb(sexActs ~ (1 | id), sexpartner) :
> NAs removed in constructing fixed-effect model frame: you should probably
remove them manually, e.g.
> with na.omit()
> 2: In II[, ii] + REmat$codes[[i]] :
> longer object length is not a multiple of shorter object length
>
> Figuring that its missing value handling might be off, I tried removing the
missing values:
> > r <- glmmadmb(sexActs~(1|id), sexpartner[!is.na(sexpartner$sexActs),])
> Error in system(cmd, intern = intern, wait = wait | intern,
show.output.on.console = wait, :
> 'C:/Program' not found
> > traceback()
> 4: system(cmd, intern = intern, wait = wait | intern, show.output.on.console =
wait,
> ...)
> 3: shell(cmd, invisible = TRUE, intern = !verbose)
> 2: run_bin(platform, bin_loc, file_name, cmdoptions, run, rm_binary =
!use_tmp_dir,
> debug = debug, verbose = verbose)
> 1: glmmadmb(sexActs ~ (1 | id), sexpartner[!is.na(sexpartner$sexActs),
> ])
Hmmm. Can you give
* the results of sessionInfo()
* the results of trying your command with debug=TRUE ?
* the results of glmmADMB:::get_bin_loc() ?
Ben Bolker
Ben
More information about the R-sig-mixed-models
mailing list