search for: smonth

Hi, i would like check a repeated measurement dataset, whether cases (which (id) could be more than one time included and not all same times included.) have more than one times the same SMONTH!? #This pseudo code didn't work with a for loop, because the [i+1] isn't known. How i could refer to [i+1] ? if(ID[i] == ID[i+1] & SMONTH[i] !=[i+1]) {res[i] <- 0 } if(ID[i] == ID[i+1] & SMONTH[i] == [i+1]) {res[i] <- 1 } many thanks, christian

....ch Message-ID: <41920F0D.9030900 at web.de> Content-Type: text/plain; charset=ISO-8859-1; format=flowed Hi, i would like check a repeated measurement dataset, whether cases (which (id) could be more than one time included and not all same times included.) have more than one times the same SMONTH!? #This pseudo code didn't work with a for loop, because the [i+1] isn't known. How i could refer to [i+1] ? if(ID[i] == ID[i+1] & SMONTH[i] !=[i+1]) {res[i] <- 0 } if(ID[i] == ID[i+1] & SMONTH[i] == [i+1]) {res[i] <- 1 } many thanks, christian -----------------------...

...have made a mistake somewhere. It is also difficult to alter models if I want to change a term. For example, below are the set of models I would like to run. Is there a way to specify the full model and have R generate the rest? I.e. specify m1234567<-glm.convert(glm.nb(mantas~site*year+cosmonth+sinmonth+coslunar+sinlunar+plankton, data=mydata)) and have R run all the other models. library(MASS) #Intercept only m0<-glm.convert(glm.nb(mantas~1,data=mydata)) #One term - 7 models #Manta abundance is greater at one of the two sites m1<-glm.convert(glm.nb(mantas~site,data=myda...