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:
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