On Nov 11, 2009, at 7:14 AM, Bogaso wrote:


Hi all,

I wanted to have a seasonality study like whether a particular month has
significant effect as compared to others. Here is my data :

0.10499 0       0       1       0       0       0       0       0       0       
0       0
0.00259 0       0       0       1       0       0       0       0       0       
0       0
-0.06015        0       0       0       0       1       0       0       0       
0       0       0
0.10721 0       0       0       0       0       1       0       0       0       
0       0
0.03597 0       0       0       0       0       0       1       0       0       
0       0
0.10584 0       0       0       0       0       0       0       1       0       
0       0
0.02063 0       0       0       0       0       0       0       0       1       
0       0
-0.03509        0       0       0       0       0       0       0       0       
0       1       0
-0.03485        0       0       0       0       0       0       0       0       
0       0       1
0.01632 0       0       0       0       0       0       0       0       0       
0       0
0.06844 1       0       0       0       0       0       0       0       0       
0       0
-0.01766        0       1       0       0       0       0       0       0       
0       0       0
0.00989 0       0       1       0       0       0       0       0       0       
0       0
0.11673 0       0       0       1       0       0       0       0       0       
0       0
0.01789 0       0       0       0       1       0       0       0       0       
0       0
-0.00323        0       0       0       0       0       1       0       0       
0       0       0
0.06811 0       0       0       0       0       0       1       0       0       
0       0
-0.01292        0       0       0       0       0       0       0       1       
0       0       0
-0.12244        0       0       0       0       0       0       0       0       
1       0       0
-0.06645        0       0       0       0       0       0       0       0       
0       1       0
-0.03355        0       0       0       0       0       0       0       0       
0       0       1
0.02308 0       0       0       0       0       0       0       0       0       
0       0
-0.11711        1       0       0       0       0       0       0       0       
0       0       0
0.06116 0       1       0       0       0       0       0       0       0       
0       0
0.02832 0       0       1       0       0       0       0       0       0       
0       0
0.01441 0       0       0       1       0       0       0       0       0       
0       0
-0.04412        0       0       0       0       1       0       0       0       
0       0       0
0.05558 0       0       0       0       0       1       0       0       0       
0       0
0.08363 0       0       0       0       0       0       1       0       0       
0       0
-0.01063        0       0       0       0       0       0       0       1       
0       0       0

i.e. all explanatory variables here are dichotomous. However once I run lm()
i got following error :

lm(
Error in model.frame.default(formula = dat[, 1] ~ dat[, -1],
drop.unused.levels = TRUE) :
 invalid type (list) for variable 'dat[, -1]'

Assuming that dat really is a data.frame, then the error message is saying you have attempted to pass most of it (since it is a list) to the rhs of the formula and the lm procedure was expecting a vector or vectors. What do you get with:

lm( dat[,1]~. , data = dat)

(It would be preferable for readability and possibly even for code correctness to refer to the column name for data[,1].)

Can anyone please tell me what to do?

Better to post the output of str(dat) , or even better, the output of dput(dat).


Best,

--

David Winsemius, MD
Heritage Laboratories
West Hartford, CT

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to