search for: model01

...have problems with the script. I don?t know what to do, I?ve already change the languages, be sure of being working in the correct directory and doesn?t work the script is: #1. example data01<-read.table("ch01.txt",header=T) names(data01) attach(data01) FERTIL<-factor(FERTIL) model01<-(YIELD~FERTIL) summary(model01) anova(model01) all my classmates, even those with a mac operator, wrote the same, and it worked. When copying the script in the R console, it seems that there is a problem with the ~ symbol. this symbol is not in the keyboard, so I select it from spetial ch...

...ying to conduct sensitivity analysis in R using the 'sensitivity' package. Although the description of functions seem straightforward, I can?t succeed. The definition of input factors can be the problem. library(sensitivity) #A simple model with 4 input factor to test the morris function: model01=function(a1,a2,a3,a4) { Z<-numeric(10) Z[1]<-runif(1) Z[2]<-runif(1,a1,30) Z[3]<-6*runif(1,min(a1,a2),max(a1,a3)) Z[4]<-runif(1,2,5)*runif(1,min(a2,a4),max(a2,a4)) Z[5]<-0.5*runif(1,min(a3,a4),max(a3,a4)) Z[6]<-2*runif(1,min(a1,a4),max(a1,a4)) Z[7]<-runif(1) Z[8]...

...ve 4 different partition of how light is absorbed in the leaf and they all add up to 1 (part1 + part2 + part3 + part4 = 1). So I have a data frame with these colums: plant | genotype | date | temperature | part1 | part2 | part3 | part4 So the logic tells me to keep it as simple as this: *model01*<- lmer(part1,part2,part3,part4~genotype:date:temperature+(1|plant),data=data,family="binomial") However, I was just wondering how these partitions correlate. So I did a test for "Variance inflation factors" on them. Correlations of the variables Part1 p...