Skip to content

lme4Eigen's bootMer & installing latest svn

5 messages · Michael Lawrence, Ken Knoblauch, Ben Bolker

Original
Michael Lawrence
# 
Hi folks,

I'm playing with lme4Eigen (version 0.9996875-8, running on Mac OS
10.7.3 using R 2.14.2) and am quite excited by the new bootMer()
function. However, when I try to run it, regardless of what fit I
provide for the argument "x" or function I provide for the argument
"FUN" (including running the examples), I get the error:

Error in envRefInferField(x, what, getClass(class(x)), selfEnv) :
  "resp" is not a valid field or method name for reference class ?lmerResp?

I presume that this is why the bootMer documentation example section
says "## Not run: %%--- fails for now --- FIXME"? I just thought I'd
double-check.

Also, I thought I'd make sure the devs know that the latest svn
version doesn't build on mac; when I try to do so, I get the error:

glmFamily.cpp: In member function ?virtual const Eigen::ArrayXd
glm::negativeBinomialDist::variance(const Eigen::ArrayXd&) const?:
glmFamily.cpp:228: error: ?Rcout? is not a member of ?Rcpp?

Cheers,

Mike

--
Mike Lawrence
Graduate Student
Department of Psychology
Dalhousie University

Looking to arrange a meeting? Check my public calendar: http://goo.gl/BYH99

~ Certainty is folly... I think. ~
Ben Bolker
# 
Mike Lawrence <Mike.Lawrence at ...> writes:
More recently (version 12) this should work ...
Thanks for the heads-up.  I will check into it and try to see about
getting new binary versions of RcppEigen and lme4Eigen up on the
repository -- although possibly not before Monday.

  cheers
    Ben
Michael Lawrence
#
On Fri, Mar 9, 2012 at 6:44 PM, Ben Bolker <bbolker at gmail.com> wrote:
While I'm not super familiar with C++, I took a look at the
"glmFamily.cpp" file indicated by the error message, and it seems that
all that has happened is that a few lines (used for debugging?)
starting "Rcpp::Rcout" need to be be commented out to "//Rcpp::Rcout"
in "glmFamily.cpp" and also "optimizer.cpp". Once this is done, the
current svn builds fine on mac.

Mike
2 days later
Ken Knoblauch
# 
Ben Bolker <bbolker at ...> writes:
I was able to get lme4Eigen_0.9996875-13  to compile this (Monday)
morning on my Mac from source after first compiling version 0.2.0
of RcppEigen, which is available on CRAN but I don't see it on Rforge.
(R version 2.14.2 Patched (2012-02-29 r58552)
Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)
Lion 10.7.3).

On an issue from a while back (and going off subject), 
I'm pleased to say that the links from
the psyphy package that allow non-zero lower asymptotes for binomial
families seem to  work with lme4Eigen, at least for binomial aggregated
data.  So far, I get errors for binary responses.  For example to use the
mafc.probit link for a 4-alternative forced-choice experiment, where
one might want to limit the lower asymptote of the link function to 0.25,
I first do the following:

Bi4 <- glmFamily$new(family = binomial(mafc.probit( 4 )))

and then use the argument 

family = Bi4$family

in the arguments to glmer.

It seems to produce promising results in simulated data when I aggregate the
binary responses but when I try it with a binary response variable, I get:

Error in FUN(1:3[[1L]], ...) : Downdated VtV is not positive definite

I would be happy to share the simulation script, if anyone is interested.

Thanks.
1 day later
Ben Bolker
# 
Ken knoblauch <ken.knoblauch at ...> writes:
[snip]
[snip]
[snip]

  For what it's worth, binaries of lme4Eigen...13  and Rcpp 0.2.0 should
be available now from http://lme4.r-forge.r-project.org/repos ...

  Ben