Displaying 2 results from an estimated 2 matches for "notna".
Did you mean:
nota
2007 Jan 25
1
poly(x) workaround when x has missing values
...lowed in 'poly'
>
> lm( y ~ poly(x,2) , subset=!is.na(x)) # This does not help?!?
Error in poly(x, 2) : missing values are not allowed in 'poly'
The following function seems to be an okay workaround.
Poly<- function(x, degree = 1, coefs = NULL, raw = FALSE, ...) {
notNA<-!is.na(x)
answer<-poly(x[notNA], degree=degree, coefs=coefs, raw=raw, ...)
THEMATRIX<-matrix(NA, nrow=length(x), ncol=degree)
THEMATRIX[notNA,]<-answer
attributes(THEMATRIX)[c('degree', 'coefs', 'class')]<- attributes(answer)[c...
2005 Jan 16
1
p.adjust(<NA>s), was "Re: [BioC] limma and p-values"
...methods) {
method <- match.arg(tolower(method),p.adjust.methods)
if(method=="fdr") method <- "bh"
if(is.data.frame(p)) {
if(length(p)) for(i in 1:length(p)) p[[i]] <-
Recall(p[[i]],method=method)
return(p)
}
porig <- p
notna <- !is.na(p)
p <- as.vector(p[notna])
n <- length(p)
if (n == 1) return(porig)
if (n == 2 && method=="hommel") method <- "hochberg"
porig[notna] <- switch(method,
holm = {
i <- 1:n
o &...