[R] Need help to read the data file like this

2013-07-14 Thread Houhou Li
Hi, I have several really big data files in csv format like this: the first line is the header, the second to fourth lines have info about the file and are the lines I need to skip (data in 2-4th lines are not correspoding to variable names in the hearder), from the fifth line, real data begins

Re: [R] Why are coefficient estimates using ML and REML are different in lme?

2012-10-29 Thread Houhou Li
similar. Is it normal to see such difference or did I obviously miss something? Thank you!    --- On Mon, 10/29/12, S Ellison wrote: From: S Ellison Subject: RE: [R] Why are coefficient estimates using ML and REML are different in lme? To: "Houhou Li" , "r-help@r-project.org&

[R] Why are coefficient estimates using ML and REML are different in lme?

2012-10-28 Thread Houhou Li
Hi, All,   My data collection is from 4 regions (a, b, c, d). Within each region, it has 2 or 3 units. Within each unit, it has measurement from about 25 sample site. I was trying to use lme function to discribe relationship between y and a few covariates. Both y and covariates were measured at

Re: [R] Error " subscript out of bounds"

2012-03-13 Thread Houhou Li
Thank you very much David. I should realize myslef that operator precedence cause the problem. Sorry about this:-) I solved the problem. --- On Tue, 3/13/12, David Winsemius wrote: From: David Winsemius Subject: Re: [R] Error " subscript out of bounds" To: "Houhou Li

[R] Error " subscript out of bounds"

2012-03-13 Thread Houhou Li
Hello, R-users,   I have a datafile with 37313 records and each record has 5 different measurements on the same variables. The format looks like this: treeID, VIG0, VIG1, VIG2, VIG3, VIG4 I was trying to convert the one row record to 5 rows record with format like this (treeID, MEASUREMENT, VIGO

[R] extremely large negative AIC value in lme mode fitting

2010-05-17 Thread Houhou Li
Hi, Everyone, I am trying to fit a ring specific gravity dataset using the linear mixed-effects model (lme function) in R. However, the AIC value is extrememly low (such as -8000) and the loglikehihhod value is extremely big (such as 4000), which seems not correct to me. What could the result i