Dear Uwe. Many thanks for your reply. I attach the data set that i use for the analyses. Hopefully this will help you identify the problem. I'm really new in using r so please bear with me.
All the best, Kyriakos > Date: Mon, 30 Apr 2012 15:20:04 +0200 > From: lig...@statistik.tu-dortmund.de > To: antonio...@hotmail.com > CC: r-help@r-project.org > Subject: Re: [R] Error in if (nuhat < 2) stop("The degrees of freedom must be > greater than or equal to 2") : missing value where TRUE/FALSE needed > > > > On 29.04.2012 16:28, Kyriakos Antoniou wrote: > > > > Hi, > > > > i am trying to run an ANCOVA and a bootstrapped ANCOVA analysis on a > > specific data set. I am using the ancova and ancboot functions as in the > > following code: > > > > setwd("C:/Users/User/Desktop/Rdatabilingualstudy2012") > > > > bilingualismdata<-read.spss("bilingualdataforconferences2012.sav", > > use.value.labels = TRUE, to.data.frame = TRUE) > > > > bilingualismdata$TypeBiling<-factor(bilingualismdata$TypeBiling, levels = > > c(0:1), labels = c("monolingual", "bilingual")) > > > > cypriot<-subset(bilingualismdata, TypeBiling=="monolingual") > > cypriotenglish<-subset(bilingualismdata, TypeBiling=="bilingual") > > > > covGrp1<-cypriot$FAS > > dvGrp1<-cypriot$Simon > > covGrp2<-cypriotenglish$FAS > > dvGrp2<-cypriotenglish$Simon > > > > source("http://www-rcf.usc.edu/~rwilcox/Rallfun-v16") > > > > ancova(covGrp1, dvGrp1, covGrp2, dvGrp2) > > ancboot(covGrp1, dvGrp1, covGrp2, dvGrp2, nboot = 2000) > > > > > > However, i always get the following errors (note that i get these errors > > when i am running the analysis on a specific data set. The same code works > > fine with other data sets): > > > >> ancova(covGrp1, dvGrp1, covGrp2, dvGrp2) > > [1] "NOTE: Confidence intervals are adjusted to control the probability" > > [1] "of at least one Type I error." > > [1] "But p-values are not" > > Error in if (nuhat< 2) stop("The degrees of freedom must be greater than > > or equal to 2") : > > missing value where TRUE/FALSE needed > > Aleady the way the Note is printed shows the code is poorly written, the > non-working error message supports that. > Since we do not have the data, we cannot help, probably lm() is capable > of the ancova you have in mind. > > Uwe Ligges > > > > > > In addition: Warning messages: > > 1: In min(sub[vecn>= 12]) : > > no non-missing arguments to min; returning Inf > > 2: In max(sub[vecn>= 12]) : > > no non-missing arguments to max; returning -Inf > > 3: In near(x1, x1[isub[i]], fr1) : NAs introduced by coercion > > 4: In near(x2, x1[isub[i]], fr2) : NAs introduced by coercion > > 5: NAs introduced by coercion > > > >> ancboot(covGrp1, dvGrp1, covGrp2, dvGrp2, nboot = 2000) > > [1] "Note: confidence intervals are adjusted to control FWE" > > [1] "But p-values are not adjusted to control FWE" > > [1] "Taking bootstrap samples. Please wait." > > Error in if (nuhat< 2) stop("The degrees of freedom must be greater than > > or equal to 2") : > > missing value where TRUE/FALSE needed > > In addition: Warning messages: > > 1: In min(sub[vecn>= 12]) : > > no non-missing arguments to min; returning Inf > > 2: In max(sub[vecn>= 12]) : > > no non-missing arguments to max; returning -Inf > > 3: In near(x1, x1[isub[i]], fr1) : NAs introduced by coercion > > 4: In near(x2, x1[isub[i]], fr2) : NAs introduced by coercion > > 5: In near(x1, x1[isub[i]], fr1) : NAs introduced by coercion > > 6: In near(x2, x1[isub[i]], fr2) : NAs introduced by coercion > > > > Any help with how i van solve this would be appreciated as i am really > > stack and don't really know what to do. > > > > Best, > > Kyriakos > > [irrelevant part deleted] > > > ______________________________________________ > > 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. [[alternative HTML version deleted]] ______________________________________________ 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.