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")

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


