Dear Uwe, really sorry about this.
The code i use is the following (data attached again): setwd("C:/Users/User/Desktop/Rdatabilingualstudy2012") bilingualdataforconferences2012<-read.delim("bilingualdataforconferences2012.dat", header=TRUE) bilingualdataforconferences2012$TypeBiling<-factor(bilingualdataforconferences2012$TypeBiling, levels = c(0:1), labels = c("monolingual", "bilingual")) cypriot<-subset(bilingualdataforconferences2012, TypeBiling=="monolingual") cypriotenglish<-subset(bilingualdataforconferences2012, 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) I am trying to run an ANCOVA and a bootstrapped ANCOVA analysis on the data set attached, using the ancova and ancboot functions as above. That is, i want to see what the effect of bilingualism (variable TypeBiling, 1=bilingual, 2=monolingual) is on the dependent variable Simon when partialling out the effect of the covariate FAS (socioeconomic status). Since the data for the dependent variable Simon are not normally distributed i need to also analyse the data using a robust method of ANCOVA. Kyriakos > Date: Mon, 30 Apr 2012 15:58:34 +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 30.04.2012 15:52, Kyriakos Antoniou wrote: > > > > Sorry!! Data now attached. Note that this is a .dat file while in the R > > code i import data directly from an SPSS file (.sav). > > > Oh, come on, we do not have endless time to import your data, so > inlcuding code please. > > Moreover, you have not really described which model yopu want to fit > exactly, we can only guess from your code not knowing those functions > you used. > > Uwe Ligges > > > > > 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. > >
______________________________________________ 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.