Nothing
summary.fg <-
function(object, sign.level=0.05, ...){
if (class(object)!="fg") {stop("Object needs to be of class \"fg\"")}
p.val <- function(perm,obsv){
return(min((sum(perm>=obsv))/(sum(is.na(perm)==FALSE))*2,(sum(perm<=obsv))/(sum(is.na(perm)==FALSE))*2,1))
}
fail <- function(perm){
return(sum(is.na(perm))/length(perm)*100)
}
if(object$pairwise==FALSE){
c.table <- cbind(object$scales,object$c.list[[1]][1]-sapply(object$c.list,mean,na.rm = TRUE),object$c.list[[1]][1]-sapply(object$c.list,quantile,probs=0.025,na.rm = TRUE),object$c.list[[1]][1]-sapply(object$c.list,quantile,na.rm = TRUE,probs=0.975),sapply(object$c.list,p.val,obsv=object$c.list[[1]][1]),sapply(object$c.list,fail))
c.table <- cbind(c.table,c.table[,5]<=sign.level)
colnames(c.table) <- c("cell_size","mean","lower_limit", "upper_limit", "P_value", "% failures","Reject H0")
rownames(c.table) <- rep("",dim(c.table)[1])
}else if(object$correlate==FALSE & object$pairwise==TRUE){
m.table <- cbind(object$scales,object$m.list[[1]][1],sapply(object$m.list,mean,na.rm = TRUE),sapply(object$m.list,quantile,probs=0.025,na.rm = TRUE),sapply(object$m.list,quantile,na.rm = TRUE,probs=0.975),sapply(object$m.list,p.val,obsv=object$m.list[[1]][1]),sapply(object$m.list,fail))
m.table <- cbind(m.table,m.table[,6]<=sign.level)
colnames(m.table) <- c("cell_size","observed","mean","lower_limit", "upper_limit", "P_value", "% failures","Reject H0")
rownames(m.table) <- rep("",dim(m.table)[1])
v.table <- cbind(object$scales,object$v.list[[1]][1],sapply(object$v.list,mean,na.rm = TRUE),sapply(object$v.list,quantile,probs=0.025,na.rm = TRUE),sapply(object$v.list,quantile,na.rm = TRUE,probs=0.975),sapply(object$v.list,p.val,obsv=object$v.list[[1]][1]),sapply(object$v.list,fail))
v.table <- cbind(v.table,v.table[,6]<=sign.level)
colnames(v.table) <- c("cell_size","observed","mean","lower_limit", "upper_limit", "P_value", "% failures","Reject H0")
rownames(v.table) <- rep("",dim(v.table)[1])
}else{
c.table <- cbind(object$scales,object$c.list[[1]][1],sapply(object$c.list,mean,na.rm = TRUE),sapply(object$c.list,quantile,probs=0.025,na.rm = TRUE),sapply(object$c.list,quantile,na.rm = TRUE,probs=0.975),sapply(object$c.list,p.val,obsv=object$c.list[[1]][1]),sapply(object$c.list,fail))
c.table <- cbind(c.table,c.table[,6]<=sign.level)
colnames(c.table) <- c("cell_size","observed","mean","lower_limit", "upper_limit", "P_value", "% failures","Reject H0")
rownames(c.table) <- rep("",dim(c.table)[1])
}
cat("\n")
cat("=====================================","\n")
cat(" Floating Grid Permutation Technique","\n")
cat("=====================================","\n")
cat("\n")
cat(paste("Number of scales:",length(object$scales)),"\n")
cat(paste("Number of iterations per scale:",object$iter),"\n")
if(object$pairwise==FALSE){
cat("\n")
cat("---------------------------------------------------------","\n")
cat(" Results for univariate analyses, unexplained Moran's I:","\n")
cat("---------------------------------------------------------","\n")
cat("\n")
print(c.table, row.names=FALSE)
cat("----------------------------------------------------------------------------","\n")
cat(paste(" NOTE: Rejection of H0 based on alpha =",sign.level),"\n")
cat("\n")
}else if(object$correlate==FALSE){
cat("\n")
cat("------------------------------------------","\n")
cat(" Results for multivariate analyses, mean:","\n")
cat("------------------------------------------","\n")
cat("\n")
print(m.table, row.names=FALSE)
cat("\n")
cat(" Results for multivariate analyses, variance:","\n")
cat("----------------------------------------------","\n")
cat("\n")
print(v.table, row.names=FALSE)
cat("----------------------------------------------------------------------------","\n")
cat(paste(" NOTE: Rejection of H0 based on alpha =",sign.level),"\n")
cat("\n")
}else{
cat("\n")
cat("-----------------------------------","\n")
cat(" Results for multivariate analyses,",object$correlate,"correlation:","\n")
cat("-----------------------------------","\n")
cat("\n")
print(c.table, row.names=FALSE)
cat("----------------------------------------------------------------------------","\n")
cat(paste(" NOTE: Rejection of H0 based on alpha =",sign.level),"\n")
cat("\n")
}
}
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.