[R-sig-ME] Resume terminated lmer fit if verbose=TRUE?

Mike Lawrence Mike.Lawrence at dal.ca
Tue Aug 30 10:59:52 CEST 2011


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:
> On Wed, Jul 27, 2011 at 03:39:59PM -0300, Mike Lawrence wrote:
>> Thanks to Harold for pointing out the start argument and indicating
>> that the verbose output can be supplied as its value.
>
> Hi Mike,
>
> I am just about to restart an interation that did not converge on 300
> iterations. I have saved both the fit-object and the output of
> glmer(verbose=TRUE). I don't really get how to format the input to the
> start value parameter. Does it suffice to provide the last iteration
> as a character string? Here is what I got for the last iteration:
>
> 300:     83537.084: 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 0.0675373 0.223978!
>  -0.314364 0.242165 -0.0398588 0.00348154 0.390582 0.205796 0.314690 -0.0210720 0.359508 0.0102215 0.0133513 -0.0249899 0.00303007 0.00104005 -0.00199693 -0.0162271 -0.00910787 -0.000721493 -0.00571186 -0.0221484 0.00132736 -0.00938054 -0.0559509 0.000625754 -0.0151419 -0.0342182 -0.0301707 -0.00874410 -0.0166841 -0.0118311 -0.0150344 -0.0455952 0.334382  1.62474 0.592823 -0.0838088 -0.0622300 -0.0477605 0.568942 -0.379506 -0.957889 -0.671352 -0.158491 -0.330145 0.269097 -0.316426 -1.04228 0.196337 -0.639169 -0.646029 -0.246746  1.11005  1.86385 0.428263  3.28910 0.565427 -0.285010 0.206211 -0.353546 0.913265 -0.534973 0.445224 -1.18672 -0.372053 -0.804367 -0.418456 -0.681271 0.0606201 0.0683290 0.685260 0.201087 -0.736270 -0.444255 -0.509914 -0.286142 -0.484370 -0.395121 -0.482115 0.300839 -0.444407 -0.444848 -0.0528026 -0.751176 0.670575 -0.0569850 -0.0616133 0.343210 -3.71098  1.26284  1.37018  1.13882  1.86587  1.14523  1.35873  1.81718  2.69451  2.54993  1.81382  1.974!
>  16  2.14052  1.94535 0.989847  2.04325  3.65850  1.65551  3.04955
>
> Or do you have to format it in any way? How did you do it?
>
> Kind regards,
>
> Hans Ekbrand
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>




More information about the R-sig-mixed-models mailing list