Displaying 1 result from an estimated 1 matches for "onetail".
Did you mean:
onefail
2006 Mar 30
1
warning message in hand-made function
...o summaries are entered, they are computed
if(missing(n)) {
n <- length(controls)
mean.c <- mean(controls, na.rm=na)
sd.c <- sd(controls, na.rm=na)
}
dl <- n-1 #degrees of freedom of the test
#t.test computation
t.obs <- (patient-mean.c) / (sd.c*(((n+1)/n)^0.5))
proba.onetailed <- 1-pt(abs(t.obs), df=dl)
rar <- pt(t.obs, df=dl) #point estimate of the rarity
#confidence intervals computation on the rarity (Crawford &
Garthwaite, 2002, Neuropsychologia)
c <- (patient-mean.c)/sd.c
#finding the non central parameter of t distributions
f <- function(...