Message-ID: <AANLkTinVWBfcpSjncPCUA24xHAOyoaqym1vofDsdgAcY@mail.gmail.com>
Date: 2010-08-26T22:26:39Z
From: Douglas Bates
Subject: R console not responding after defining variable as factor
In-Reply-To: <3154E0B9605B3140AFAD0A3A97DA29CF832426@SRV603.tudelft.net>
On Thu, Aug 26, 2010 at 3:29 PM, Sabira el Messlaki
<S.elMesslaki at student.tudelft.nl> wrote:
> Dear list members,
> I have a problem running the lme model. The problem is caused by the random factor TransJaarCod. TransJaarCod is a factor with 12 levels. If I run the underneath model without defining TransJaarCod (R recognize it by default as interval) ?as a factor the model works fine, but when I define TransJaarCod as a factor the R console turns not responding. I thought it was because of a lack of memory on my notebook but I tried the same on my desktop and the same happens there.
> model3.7.fit <- lme(y ~ x1+LnBankt500+Perceelkwaliteit+LnNS_Afstand+PGebouwd+Bouwjaar_cat4+LnBanIndDistr500+LeegstandJaar_L2, random = ~TransJaarCod|GebID,data, method = "REML")
It's not that the console has frozen - it's that the R process is
going to be very, very busy for a long time (and possibly run out of
memory space). When you include a random effect specification with a
continuous covariate you need to estimate 3 variance components - 2
variances and a covariance. When you have a factor with 12 levels you
need to estimate 12 variances and 66 covariances. That will take a
very long time, if it could be done at all.
You should reconsider the form of your model.
> Can someone please help me out.
>
>
>> str(data)
> 'data.frame': ? 695 obs. of ?35 variables:
> ?$ nr ? ? ? ? ? ? ? ? ? ? ?: int ?59 72 74 89 129 130 145 146 157 168 ...
> ?$ GebID ? ? ? ? ? ? ? ? ? : Factor w/ 165 levels "AFC-01","ALA-01",..: 16 20 21 23 32 32 36 36 38 41 ...
> ?$ TransJaarCod ? ? ? ? ? ?: int ?1 1 1 1 1 1 1 1 1 1 ...
> ?$ L_3L750 ? ? ? ? ? ? ? ? : int ?2100 500 0 2000 42100 42100 41250 41250 10910 42100 ...
> ?$ LnRHuurM2 ? ? ? ? ? ? ? : num ?4.72 5.27 5.17 5.21 5.07 ...
> ?$ LnReis_opaf_m ? ? ? ? ? : num ?0.945 1.446 0.236 1.187 0.896 ...
> ?$ LnNS_Afstand ? ? ? ? ? ?: num ?6.46 7.28 6.89 7.24 6.6 ...
> ?$ LnIC_afstand ? ? ? ? ? ?: num ?7.98 7.28 6.89 7.46 6.6 ...
> ?$ LnBankt500 ? ? ? ? ? ? ?: num ?5.54 6.93 7.68 6.71 8.46 ...
> ?$ LnBanIndDistr500 ? ? ? ?: int ?152 218 566 114 1776 1776 1175 1175 327 3607 ...
> ?$ LnVoorzDagel500 ? ? ? ? : num ?2.77 3 0 1.95 1.1 ...
> ?$ PGebouwd ? ? ? ? ? ? ? ?: int ?0 0 0 0 0 0 0 0 0 0 ...
> ?$ Bouwjaar_cat4 ? ? ? ? ? : int ?1 1 3 2 3 3 3 3 3 3 ...
> ?$ LogoBedrijf ? ? ? ? ? ? : int ?0 0 0 1 1 1 0 0 0 1 ...
> ?$ Perceelkwaliteit ? ? ? ?: int ?0 1 1 1 1 1 0 0 0 1 ...
> ?$ KwaliteitStraatmeubilair: int ?0 0 0 0 0 0 0 0 0 0 ...
> ?$ X_Coord ? ? ? ? ? ? ? ? : int ?117591 119543 123932 120624 125236 125236 124636 124636 122515 124799 ...
> ?$ Y_Coord ? ? ? ? ? ? ? ? : int ?486259 482125 482698 482028 479864 479864 480391 480391 486650 479899 ...
> ?$ LeegstandJaar_L2 ? ? ? ?: num ?9.78 9.78 9.78 9.78 9.78 ...
> ?[list output truncated]
>
> ? ? ? ?[[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>