Re: [R] Automating R for Hypothesis Testing

2012-05-16 Thread meredith
Rui- Just a quick question. I understand your comment on using ANOVA, but doesn't this only test for similarities of the mean. We are trying to see if a the same model can fit for two or three months, therefore have the similar slope and intercept. The ANOVA would only do one part of this correct

Re: [R] Automating R for Hypothesis Testing

2012-05-10 Thread Rui Barradas
Hello, I'm glad it helped. As for your second question, I don't know, but I'm not very comfortable with the way you're doing things. Why subtract the coefficients of model 1 from model 2? And why the dummy? Why set model 1 to zero? Isn't it better to use anova's F? After all, it's designed for i

Re: [R] Automating R for Hypothesis Testing

2012-05-10 Thread meredith
Rui- Thanks this definitely helps, just one quick question. How would you code the values of chi-fm and chi-fms to change based on the degrees of freedom of each model H(i)? Meredith Rui Barradas wrote > > Hello, > > Yes, it does help. Now we can see your data and what you're doing. > What f

Re: [R] Automating R for Hypothesis Testing

2012-05-09 Thread Rui Barradas
Hello, Yes, it does help. Now we can see your data and what you're doing. What follows is a suggestion on what you could do, not full solution. (You forgot to say what X1 is, but I don't think it's important to understand the suggestion.) (If I'm wrong, say something.) milwaukeephos <- read.csv(

Re: [R] Automating R for Hypothesis Testing

2012-05-09 Thread meredith
dput: http://r.789695.n4.nabble.com/file/n4620188/milwaukeephos.csv milwaukeephos.csv # Feb-march > modelH_febmarch<-lm(llfeb_march~lffeb_march) >modelHa_febmarch<-lm(llfeb_march~X1feb_mar+lffeb_march) > anova(modelHa_febmarch) > coefficients(modelH_febmarch) (Intercept) lffeb_march -2.429890

Re: [R] Automating R for Hypothesis Testing

2012-05-08 Thread Rui Barradas
Hello, I'm not at all sure if I understand your problem. Does this describe it? test first model for months 1 and 2 if test statistic less than critical value{ test second model for months 1 and 2 print results of the first and second tests? just one of them? } move on to months

[R] Automating R for Hypothesis Testing

2012-05-08 Thread meredith
R Users- I have been trying to automate a manual code that I have developed for calling in a .csv file, isolating certain rows and columns that correspond to specified months: something to the effect i=name.csv N=length(i$month) iphos1=0 iphos2=0 isphos3=0 for i=1,N if month=1 iphos1=iphos+1