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.
