Nothing
screen.pfc <-
function(X, fy, cutoff=0.1)
{
unipfc = function(x)
{
temp.dat<-data.frame(cbind(x, fy))
xnam<-paste("xx", 1:r, sep="")
names(temp.dat)<-c("yy", xnam)
fm.lm<- as.formula( paste("yy ~ ", paste(xnam, collapse= "+")))
fstat<-summary( lm(fm.lm, data=temp.dat))$fstatistic
return( as.vector( c(fstat[1], 1-pf(fstat[1], df1=fstat[2], df2=fstat[3]))))
}
op <- dim(X); nobs <- op[1]; npred <- op[2]; r <- dim(fy)[2]
if (is.null(colnames(X))) colnames(X) <- paste(X, 1:npred, sep="")
cnames <- c("F", "P-value", "Index")
OneXs<-data.frame(matrix(data=NA, ncol=3, nrow=npred))
colnames(OneXs) <- cnames
rownames(OneXs) <- colnames(X)
for(i in 1:npred)
{
tempX<-array(X[,i], c(nobs, 1))
OneXs[i,]<-c(round(unipfc(x=tempX), digits=5), i)
}
Xrelevant=OneXs[which(OneXs$`P-value` <= cutoff),]
o <- order(as.numeric(Xrelevant$F), decreasing = TRUE)
return(Xrelevant[o, ])
}
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.