I haven't finished with the notebook
http://nbviewer.ipython.org/github/dmbates/JuliaWorkshop/blob/master/ParallelGLM.ipynb
but, unfortunately, must turn my attention to something else. I'll get
back to it tomorrow. In the meantime, if someone has a hint as to where
the memory allocation is
Did you try julia --track-allocation=user (or all, if necessary)? See the
docs:
http://docs.julialang.org/en/latest/stdlib/base/#Base.clear_malloc_data
--Tim
On Thursday, October 30, 2014 02:16:37 PM Douglas Bates wrote:
I haven't finished with the notebook
I tried calling XtWXXtWZ! with a Distribution object e.g. Bernoulli()
instead of a type, likewise for Link, as in
XtWXXtWz!(XtWX,XtWz,Xt,β,y,wt,Bernoulli(),LogitLink())
and changing the appropriate lines to
μ = linkinv(*typeof*(L),η)
μη = dμdη(*typeof*(L),η)
W