Skip to content
Prev 6621 / 20628 Next

Resume terminated lmer fit if verbose=TRUE?

I think you want something like:

from_verbose = "0.519965 0.143416 0.510886 0.390266 0.0253851 0.228085
 0.00000  0.00000 0.175437 0.167589 0.000151413 0.257648  0.00000
0.00000  0.00000  0.00000  0.00000  0.00000  0.00000 0.662743 0.697672
0.0943584 0.238029 0.465242 0.264901 0.177211 -0.255314 0.297930
-0.0908047 0.234670 0.232733 0.0640540 0.0191507 -0.207110 0.370959
0.639145 0.275867 -1.61519  1.14158 -0.519320 -0.0995158 -0.00705858
-0.495069 -0.995847 -0.259104 -0.751779 -0.613690 -0.281145 -0.633951
-0.128201 -2.07103  2.11410 0.817327  1.09840 0.287138 -0.699452
0.511890 -0.227955 0.293735 0.443935 -0.107533 -0.159424 -0.380081
-0.0856487 0.271927 -0.419327 0.200206 0.0556932 -0.373988 0.219712
1.11540 -0.204753 -0.250565 0.227195 0.338191 0.0638902 0.381680
0.347678 0.656791 0.610805 0.286053 0.387693 0.0700578 0.311415
0.216640 -1.43990  8.21723 -0.0677624 -1.82260  1.98123  8.53073
-1.47762  4.62704  9.22425  16.0583  9.00298 0.610718  6.34499
0.114604 -0.268301 0.659427"
from_verbose = strsplit(from_verbose,' ')[[1]]
from_verbose = as.numeric(from_verbose[from_verbose!=''])
to_start = list(
    ST = list(
        matrix(from_verbose[1],1,1)
    )
     , fixef = from_verbose[2:length(from_verbose)]
)
dimnames(to_start$ST[[1]]) = list('(Intercept)','(Intercept)')

my_new_fit = lmer(
    ...#your data, formula, family, etc., here
    , start = to_start
)
On Tue, Aug 30, 2011 at 5:07 AM, Hans Ekbrand <hans at sociologi.cjb.net> wrote: