Nothing
checkMIS<-function(obj,nitems,MISmodel="obj",obj.names=NULL, verbose=FALSE){
nobj<-nitems
## prepare data
ncomp<-nobj*(nobj-1)/2
if(is.character(obj)){ ## datafilename supplied
datafile <- obj
if(file.access(datafile, mode=0) == 0){
dat<-as.matrix(read.table(datafile,header=TRUE)) # datafile
} else {
stop("\ninput data file does not exist!\n")
}
} else if(is.data.frame(obj)){ ## data frame supplied
dat<-as.matrix(obj) # dataframe
dat<-apply(dat,2,as.numeric)
} else {
stop("first argument must be either datafilename or dataframe")
}
varnames<-colnames(dat)
if (ncol(dat)>ncomp) {
covnames<-varnames[(ncomp+1):ncol(dat)]
covs<-as.data.frame(dat[,(ncomp+1):ncol(dat)])
} else {
covs<-NULL
}
## vector of counts of NAs, columnwise
misv<-apply(dat,2,function(x) sum(is.na(x)) )
mism<-matrix(0,nobj,nobj)
mism[upper.tri(mism,diag=FALSE)]<-misv
mism<-mism+t(mism)
objMis<-colSums(mism)
if (is.null(obj.names)){
rownames(mism)<-colnames(mism)<-paste("o",1:nobj,sep="")
names(objMis)<-paste("o",1:nobj,sep="")
} else {
rownames(mism)<-colnames(mism)<-obj.names
names(objMis)<-obj.names
}
objMis<-colSums(mism)
if(verbose){
cat("number of missing comparisons:\n")
print(mism)
cat("number of missing comparisons for objects:\n",objMis,"\n")
}
if (MISmodel=="obj")
RET<-ifelse(objMis==0,F,T)
else if (MISmodel=="comp")
RET<-ifelse(misv==0,F,T)
else
stop('\nNo output: MISmodel not correctly specified. Use "obj" or "comp"\n')
invisible(RET)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.