R version 2.15.0 (2012-03-30)


> library(nlme)
> dataf<-read.table("C:/Documents and Settings/New user.HOME-2DCF24A8DD/Desktop/analysis/pdataf.txt",header=T)
> dataf
    tag year sex phase age   mstart     mend        mc    propmc    kmtoF1
1 PO043 2009  AF    PM  10 373.4657 475.6858 102.22008 27.370676  420.2365
2 RR435 2009  AF    PM   4 270.0599 516.3509 246.29104 91.198685 2195.4183
3 RR435 2010  AF    PM   5 331.8893 531.0843 199.19494 60.018484 1623.7955
4 WR336 2008  AF    PB  11 418.6430 507.9936  89.35066 21.342926  781.5365
5 YY189 2009  AF    PM   5 300.1781 508.0098 207.83174 69.236153 1010.0676
6 YY193 2007  AF    PB   4 346.3593 430.2214  83.86210 24.212458    0.0000
7 YY240 2007  AF    PB   4 372.1941 406.9324  34.73827  9.333375 1268.4038
  meankmtoF1 htoF1 durF1 meandurFpat totFtime tripdur propforage    farkm fx
1   394.4627   280   232    214.2222     1928    5632   34.23295 1476.228  9
2  2217.2725  1136   240    433.0000     3464    6480   53.45679 2268.101  8
3  1642.6462  1240   168    288.0000     2880    8064   35.71429 2022.378 10
4   896.8356   376   736    736.0000      736    1536   47.91667 1031.621  1
5  1011.0565   496    96    190.0000     1520    5832   26.06310 2201.082  8
6     0.0000     0     0      0.0000        0    1824    0.00000 1931.764  0
7  1324.0800   616   448    448.0000      448    2056   21.78988 1425.120  1
          mv    febpm
1 0.20312652 30.80966
2 0.25818408 48.11111
3 0.07446295 32.14286
4 1.96484928 43.12500
5 0.07896997 23.45679
6 0.00000000  0.00000
7 2.71668110 19.61089
> attach(dataf)
> fm0f<-lme(propmc~kmtoF1+htoF1+tripdur+farkm,random=~1|tag,method="ML")
Error in MEestimate(lmeSt, grps) : 
  Singularity in backsolve at level 0, block 1
> fm0f<-lme(propmc~kmtoF1+tripdur,random=~1|tag,method="ML")
> summary(fm0f)
Linear mixed-effects model fit by maximum likelihood
 Data: NULL 
       AIC      BIC    logLik
  68.36043 68.08998 -29.18022

Random effects:
 Formula: ~1 | tag
         (Intercept) Residual
StdDev: 0.0007653108 15.63737

Fixed effects: propmc ~ kmtoF1 + tripdur 
                Value Std.Error DF   t-value p-value
(Intercept) -1.619729 17.010746  5 -0.095218  0.9278
kmtoF1       0.014981  0.013920 -1  1.076211     NaN
tripdur      0.006514  0.003898 -1  1.671239     NaN
 Correlation: 
        (Intr) kmtoF1
kmtoF1  -0.269       
tripdur -0.544 -0.568

Standardized Within-Group Residuals:
        Min          Q1         Med          Q3         Max 
-1.37118825 -0.93396390  0.07986324  1.01258418  1.13408446 

Number of Observations: 7
Number of Groups: 6 
Warning message:
In pt(q, df, lower.tail, log.p) : NaNs produced
>