R/apollo_combineResults.R
In byu-transpolab/apollo-byu: Tools for Choice Model Estimation and Application
Defines functions
apollo_combineResults
Documented in
apollo_combineResults
#' Write model results to file
#'
#' Writes results from various models to a single CSV file.
#'
#' @param combineResults_settings List of options. It can include the following.
#' \itemize{
#' \item modelNames: Character vector. List of names of models to combine. Use an empty vector to combine results from all models in the directory.
#' \item printClassical: Boolean. TRUE for printing classical standard errors. FALSE by default.
#' \item printPVal: Boolean. TRUE for printing p-values. FALSE by default.
#' \item printT1: Boolean. If TRUE, t-test for H0: apollo_beta=1 are printed. FALSE by default.
#' \item estimateDigits: Numeric scalar. Number of decimal places to print for estimates. Default is 4.
#' \item tDigits: Numeric scalar. Number of decimal places to print for t-ratios values. Default is 2.
#' \item pDigits: Numeric scalar. Number of decimal places to print for p-values. Default is 2.
#' \item sortByDate: Boolean. If TRUE, models are ordered by date. Default is TRUE.
#' }
#' @return Nothing, but writes a file called 'model_comparison_[date].csv' in the working directory.
#' @export
apollo_combineResults = function(combineResults_settings=NULL){
if(is.null(combineResults_settings)) combineResults_settings=list()
if(is.null(combineResults_settings[["modelNames"]])){
apollo_print("The combineResults_settings does not include an object called \"modelNames\". The apollo_combineResults function will include all models for which results have been stored in the working directory. Note that this function is not applicable for models estimated using HB.")
combineResults_settings[["modelNames"]]=list.files(pattern="*estimates.csv")
if(length(combineResults_settings[["modelNames"]])==0) stop('No model files found in the working directory!')
for(j in 1:length(combineResults_settings[["modelNames"]])){
l=nchar(combineResults_settings[["modelNames"]][j])
l=l-14
combineResults_settings[["modelNames"]][j]=substr(combineResults_settings[["modelNames"]][j],1,l)
} }
if(is.null(combineResults_settings[["printClassical"]])) combineResults_settings[["printClassical"]]=FALSE
if(is.null(combineResults_settings[["printPVal"]])) combineResults_settings[["printPVal"]]=FALSE
if(is.null(combineResults_settings[["printT1"]])) combineResults_settings[["printT1"]]=FALSE
if(is.null(combineResults_settings[["estimateDigits"]])) combineResults_settings[["estimateDigits"]]=4
if(is.null(combineResults_settings[["tDigits"]])) combineResults_settings[["tDigits"]]=2
if(is.null(combineResults_settings[["pDigits"]])) combineResults_settings[["pDigits"]]=2
modelNames=combineResults_settings[["modelNames"]]
printClassical=combineResults_settings[["printClassical"]]
printPVal=combineResults_settings[["printPVal"]]
printT1=combineResults_settings[["printT1"]]
estimateDigits=combineResults_settings[["estimateDigits"]]
tDigits=combineResults_settings[["tDigits"]]
pDigits=combineResults_settings[["pDigits"]]
if(is.null(combineResults_settings[["sortByDate"]])) combineResults_settings[["sortByDate"]]=TRUE
if(combineResults_settings[["sortByDate"]]){
details = file.info(paste(modelNames,"_estimates.csv",sep=""))
rownames(details)=modelNames
details = details[with(details, order(as.POSIXct(mtime))), ]
modelNames = rownames(details)
}
Cfile_check=paste(modelNames,"_C.csv",sep="")
outputfile_check=paste(modelNames,"_output.txt",sep="")
estimatefile_check=paste(modelNames,"_estimates.csv",sep="")
files=list.files()
if(any(Cfile_check%in%files)) stop("Your list of modelNames includes some models estimated using HB!")
if(!(all(outputfile_check%in%files))) {
cat("\nCould not find: ",outputfile_check[!outputfile_check%in%files],"\n")
stop("Model output files missing for some models!")}
if(!(all(estimatefile_check%in%files))){
cat("\nCould not find: ",estimatefile_check[!estimatefile_check%in%files],"\n")
stop("Files with estimates missing for some models!")}
estimateDigits = max(1,estimateDigits)
tDigits = max(1,tDigits)
pDigits = max(1,pDigits)
if(!is.character(modelNames)) stop("Argument 'modelNames' must be a character vector.")
estimates = list()
otheroutputs = data.frame(matrix(0,nrow=9,ncol=length(modelNames)))
rownames(otheroutputs) = c("Model name","Model description","Number of individuals","Number of modelled outcomes","Estimated parameters","LL(final)","Adj.Rho-square (0)","AIC","BIC")
values = 1 + ( 1 + printClassical ) * ( 1 + printT1) * ( 1 + printPVal )
for(j in 1:length(modelNames)){
filename=paste(modelNames[[j]],"_estimates.csv",sep="")
if(!file.exists(filename)) stop("File ",filename," not found!")
inputs = tryCatch(utils::read.csv(filename),
warning=function(w) x=FALSE,
error=function(e) x=FALSE)
if(is.logical(inputs) && inputs==FALSE) stop("Could not open file ",filename)
if(printClassical==TRUE){
if(!("t.ratio.0." %in% colnames(inputs))){
cat("\nClassical t.ratios not available in ",paste(modelNames[[j]],"_estimates.csv",sep=""))
inputs[,"t.ratio.0."]=NA
}}
if(printT1==TRUE){
if(!("Rob.t.ratio.1." %in% colnames(inputs))){
cat("\nt.ratios against 1 not available in ",paste(modelNames[[j]],"_estimates.csv",sep=""))
inputs[,"Rob.t.ratio.1."]=NA
if(printClassical==TRUE) inputs[,"t.ratio.1."]=NA
}}
if(printPVal==TRUE){
if(!("Rob.p.val.0." %in% colnames(inputs))){
cat("\np-values not available in ",paste(modelNames[[j]],"_estimates.csv",sep=""))
inputs[,"Rob.p.val.0."]=NA
if(printT1==TRUE) inputs[,"Rob.p.val.1."]=NA
if(printClassical==TRUE){
inputs[,"p.val.0."]=NA
if(printT1==TRUE) inputs[,"p.val.1."]=NA
}
}}
estimates[[j]]=round(inputs[,"Estimate"],estimateDigits)
if(printClassical==TRUE){
estimates[[j]]=cbind(estimates[[j]],round(inputs[,"t.ratio.0."],tDigits))
if(printPVal==TRUE) estimates[[j]]=cbind(estimates[[j]],round(inputs[,"p.val.0."],pDigits))
if(printT1==TRUE){
estimates[[j]]=cbind(estimates[[j]],round(inputs[,"t.ratio.1."],tDigits))
if(printPVal==TRUE) estimates[[j]]=cbind(estimates[[j]],round(inputs[,"p.val.1."],pDigits))
}
}
estimates[[j]]=cbind(estimates[[j]],round(inputs[,"Rob.t.ratio.0."],tDigits))
if(printPVal==TRUE) estimates[[j]]=cbind(estimates[[j]],round(inputs[,"Rob.p.val.0."],pDigits))
if(printT1==TRUE){
estimates[[j]]=cbind(estimates[[j]],round(inputs[,"Rob.t.ratio.1."],tDigits))
if(printPVal==TRUE) estimates[[j]]=cbind(estimates[[j]],round(inputs[,"Rob.p.val.1."],pDigits))
}
rownames(estimates[[j]])=inputs[,1]
if(j==1){
combined_names=rownames(estimates[[j]])
}else{
combined_names=c(combined_names,rownames(estimates[[j]]))
}
}
combined_names=unique(combined_names)
combined_outputs=data.frame(matrix(0,nrow=length(combined_names)+1,ncol=values*length(modelNames)))
row.names(combined_outputs)=c("",combined_names)
outputnames=c("estimate")
if(printClassical==TRUE){
outputnames=c(outputnames,"t-ratio(0)")
if(printPVal==TRUE) outputnames=c(outputnames,"p-val(0)")
if(printT1==TRUE){
outputnames=c(outputnames,"t-ratio(1)")
if(printPVal==TRUE) outputnames=c(outputnames,"p-val(1)")
}
}
outputnames=c(outputnames,"Rob.t-ratio(0)")
if(printPVal==TRUE) outputnames=c(outputnames,"Rob.p-val(0)")
if(printT1==TRUE){
outputnames=c(outputnames,"Rob.t-ratio(1)")
if(printPVal==TRUE) outputnames=c(outputnames,"Rob.p-val(1)")
}
combined_outputs[1,]=rep(outputnames,length(modelNames))
for(j in 1:length(modelNames)){
applicable_names=combined_names[combined_names%in%row.names(estimates[[j]])]
combined_outputs[applicable_names,((j-1)*values+1):(j*values)]=estimates[[j]][applicable_names,]
inapplicable_names=combined_names[!(combined_names%in%row.names(estimates[[j]]))]
combined_outputs[inapplicable_names,((j-1)*values+1):(j*values)]=NA
}
for(j in 1:length(modelNames)){
filename=paste(modelNames[[j]],"_output.txt",sep="")
if(!file.exists(filename)) stop("File ",filename," not found!")
lines = tryCatch(readLines(filename),
warning=function(w) x=FALSE,
error=function(e) x=FALSE)
if(is.logical(lines) && lines==FALSE) stop("Could not open file ",filename)
inputvar = grep("Model name", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position = gregexpr(pattern=": ",inputvar)[[1]][1]
otheroutputs[1,j]=substr(inputvar,position+1,nchar(inputvar))
} else otheroutputs[1,j] = NA
k=2
inputvar = grep("Model description", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position = gregexpr(pattern=": ",inputvar)[[1]][1]
otheroutputs[k,j]=substr(inputvar,position+1,nchar(inputvar))
} else otheroutputs[k,j] = NA
k=k+1
inputvar = grep("Number of individuals", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position=gregexpr(pattern=":",inputvar)[[1]][1]
otheroutputs[k,j]=as.double(substr(inputvar,position+1,nchar(inputvar)))
} else otheroutputs[k,j] = NA
k=k+1
inputvar = grep("Number of observations", lines)
if(length(inputvar)==0) inputvar = grep("Number of modelled outcomes", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position=gregexpr(pattern=":",inputvar)[[1]][1]
otheroutputs[k,j]=as.double(substr(inputvar,position+1,nchar(inputvar)))
} else otheroutputs[k,j] = NA
k=k+1
inputvar = grep("Estimated parameters", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position=gregexpr(pattern=":",inputvar)[[1]][1]
otheroutputs[k,j]=as.double(substr(inputvar,position+1,nchar(inputvar)))
} else otheroutputs[k,j] = NA
k=k+1
inputvar = grep("LL\\(final)", lines)
if(length(inputvar)==0) inputvar = grep("LL\\(final, whole model)", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position=gregexpr(pattern=":",inputvar)[[1]][1]
otheroutputs[k,j]=as.double(substr(inputvar,position+1,nchar(inputvar)))
} else otheroutputs[k,j] = NA
k=k+1
inputvar = grep("Adj.Rho-square \\(0\\)", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position=gregexpr(pattern=":",inputvar)[[1]][1]
tmp=substr(inputvar,position+1,nchar(inputvar))
if(tmp==" Not applicable"){
otheroutputs[k,j]=NA
} else {
otheroutputs[k,j]=as.double(tmp)
}
} else otheroutputs[k,j] = NA
k=k+1
inputvar = grep("AIC", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position=gregexpr(pattern=":",inputvar)[[1]][1]
otheroutputs[k,j]=as.double(substr(inputvar,position+1,nchar(inputvar)))
} else otheroutputs[k,j] = NA
k=k+1
inputvar = grep("BIC", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position=gregexpr(pattern=":",inputvar)[[1]][1]
otheroutputs[k,j]=as.double(substr(inputvar,position+1,nchar(inputvar)))
} else otheroutputs[k,j] = NA
}
otherouputs_new=data.frame(matrix("",nrow=9,ncol=values*length(modelNames)))
rownames(otherouputs_new)=rownames(otheroutputs)
for(j in 1:length(modelNames)){
otherouputs_new[,((j-1)*values+1)]=otheroutputs[,j]
colnames(otherouputs_new)[((j-1)*values+1)]=colnames(otheroutputs)[j]
colnames(otherouputs_new)[((j-1)*values+2)]=""
if(printClassical==TRUE){
colnames(otherouputs_new)[((j-1)*values+3)]=""
}
}
otherouputs_new=rbind(otherouputs_new,rep("",ncol(otherouputs_new)))
rownames(otherouputs_new)[10]=""
filename=paste("model_comparison_",gsub("[: -]", "" , Sys.time(), perl=TRUE),".csv",sep="")
utils::write.table(otherouputs_new, filename, sep = ",", col.names = F, append = T)
#utils::write.table("",filename, sep = ",", col.names = F, append = T)
utils::write.table(combined_outputs, filename, sep = ",", col.names = F, append = T)
cat("Outputs of apollo_combineResults saved to",filename,"\n")
return(invisible(TRUE))
}
byu-transpolab/apollo-byu documentation built on Dec. 19, 2021, 12:49 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Defines functions apollo_combineResults
Documented in apollo_combineResults
#' Write model results to file
#'
#' Writes results from various models to a single CSV file.
#'
#' @param combineResults_settings List of options. It can include the following.
#' \itemize{
#' \item modelNames: Character vector. List of names of models to combine. Use an empty vector to combine results from all models in the directory.
#' \item printClassical: Boolean. TRUE for printing classical standard errors. FALSE by default.
#' \item printPVal: Boolean. TRUE for printing p-values. FALSE by default.
#' \item printT1: Boolean. If TRUE, t-test for H0: apollo_beta=1 are printed. FALSE by default.
#' \item estimateDigits: Numeric scalar. Number of decimal places to print for estimates. Default is 4.
#' \item tDigits: Numeric scalar. Number of decimal places to print for t-ratios values. Default is 2.
#' \item pDigits: Numeric scalar. Number of decimal places to print for p-values. Default is 2.
#' \item sortByDate: Boolean. If TRUE, models are ordered by date. Default is TRUE.
#' }
#' @return Nothing, but writes a file called 'model_comparison_[date].csv' in the working directory.
#' @export
apollo_combineResults = function(combineResults_settings=NULL){
if(is.null(combineResults_settings)) combineResults_settings=list()
if(is.null(combineResults_settings[["modelNames"]])){
apollo_print("The combineResults_settings does not include an object called \"modelNames\". The apollo_combineResults function will include all models for which results have been stored in the working directory. Note that this function is not applicable for models estimated using HB.")
combineResults_settings[["modelNames"]]=list.files(pattern="*estimates.csv")
if(length(combineResults_settings[["modelNames"]])==0) stop('No model files found in the working directory!')
for(j in 1:length(combineResults_settings[["modelNames"]])){
l=nchar(combineResults_settings[["modelNames"]][j])
l=l-14
combineResults_settings[["modelNames"]][j]=substr(combineResults_settings[["modelNames"]][j],1,l)
} }
if(is.null(combineResults_settings[["printClassical"]])) combineResults_settings[["printClassical"]]=FALSE
if(is.null(combineResults_settings[["printPVal"]])) combineResults_settings[["printPVal"]]=FALSE
if(is.null(combineResults_settings[["printT1"]])) combineResults_settings[["printT1"]]=FALSE
if(is.null(combineResults_settings[["estimateDigits"]])) combineResults_settings[["estimateDigits"]]=4
if(is.null(combineResults_settings[["tDigits"]])) combineResults_settings[["tDigits"]]=2
if(is.null(combineResults_settings[["pDigits"]])) combineResults_settings[["pDigits"]]=2
modelNames=combineResults_settings[["modelNames"]]
printClassical=combineResults_settings[["printClassical"]]
printPVal=combineResults_settings[["printPVal"]]
printT1=combineResults_settings[["printT1"]]
estimateDigits=combineResults_settings[["estimateDigits"]]
tDigits=combineResults_settings[["tDigits"]]
pDigits=combineResults_settings[["pDigits"]]
if(is.null(combineResults_settings[["sortByDate"]])) combineResults_settings[["sortByDate"]]=TRUE
if(combineResults_settings[["sortByDate"]]){
details = file.info(paste(modelNames,"_estimates.csv",sep=""))
rownames(details)=modelNames
details = details[with(details, order(as.POSIXct(mtime))), ]
modelNames = rownames(details)
}
Cfile_check=paste(modelNames,"_C.csv",sep="")
outputfile_check=paste(modelNames,"_output.txt",sep="")
estimatefile_check=paste(modelNames,"_estimates.csv",sep="")
files=list.files()
if(any(Cfile_check%in%files)) stop("Your list of modelNames includes some models estimated using HB!")
if(!(all(outputfile_check%in%files))) {
cat("\nCould not find: ",outputfile_check[!outputfile_check%in%files],"\n")
stop("Model output files missing for some models!")}
if(!(all(estimatefile_check%in%files))){
cat("\nCould not find: ",estimatefile_check[!estimatefile_check%in%files],"\n")
stop("Files with estimates missing for some models!")}
estimateDigits = max(1,estimateDigits)
tDigits = max(1,tDigits)
pDigits = max(1,pDigits)
if(!is.character(modelNames)) stop("Argument 'modelNames' must be a character vector.")
estimates = list()
otheroutputs = data.frame(matrix(0,nrow=9,ncol=length(modelNames)))
rownames(otheroutputs) = c("Model name","Model description","Number of individuals","Number of modelled outcomes","Estimated parameters","LL(final)","Adj.Rho-square (0)","AIC","BIC")
values = 1 + ( 1 + printClassical ) * ( 1 + printT1) * ( 1 + printPVal )
for(j in 1:length(modelNames)){
filename=paste(modelNames[[j]],"_estimates.csv",sep="")
if(!file.exists(filename)) stop("File ",filename," not found!")
inputs = tryCatch(utils::read.csv(filename),
warning=function(w) x=FALSE,
error=function(e) x=FALSE)
if(is.logical(inputs) && inputs==FALSE) stop("Could not open file ",filename)
if(printClassical==TRUE){
if(!("t.ratio.0." %in% colnames(inputs))){
cat("\nClassical t.ratios not available in ",paste(modelNames[[j]],"_estimates.csv",sep=""))
inputs[,"t.ratio.0."]=NA
}}
if(printT1==TRUE){
if(!("Rob.t.ratio.1." %in% colnames(inputs))){
cat("\nt.ratios against 1 not available in ",paste(modelNames[[j]],"_estimates.csv",sep=""))
inputs[,"Rob.t.ratio.1."]=NA
if(printClassical==TRUE) inputs[,"t.ratio.1."]=NA
}}
if(printPVal==TRUE){
if(!("Rob.p.val.0." %in% colnames(inputs))){
cat("\np-values not available in ",paste(modelNames[[j]],"_estimates.csv",sep=""))
inputs[,"Rob.p.val.0."]=NA
if(printT1==TRUE) inputs[,"Rob.p.val.1."]=NA
if(printClassical==TRUE){
inputs[,"p.val.0."]=NA
if(printT1==TRUE) inputs[,"p.val.1."]=NA
}
}}
estimates[[j]]=round(inputs[,"Estimate"],estimateDigits)
if(printClassical==TRUE){
estimates[[j]]=cbind(estimates[[j]],round(inputs[,"t.ratio.0."],tDigits))
if(printPVal==TRUE) estimates[[j]]=cbind(estimates[[j]],round(inputs[,"p.val.0."],pDigits))
if(printT1==TRUE){
estimates[[j]]=cbind(estimates[[j]],round(inputs[,"t.ratio.1."],tDigits))
if(printPVal==TRUE) estimates[[j]]=cbind(estimates[[j]],round(inputs[,"p.val.1."],pDigits))
}
}
estimates[[j]]=cbind(estimates[[j]],round(inputs[,"Rob.t.ratio.0."],tDigits))
if(printPVal==TRUE) estimates[[j]]=cbind(estimates[[j]],round(inputs[,"Rob.p.val.0."],pDigits))
if(printT1==TRUE){
estimates[[j]]=cbind(estimates[[j]],round(inputs[,"Rob.t.ratio.1."],tDigits))
if(printPVal==TRUE) estimates[[j]]=cbind(estimates[[j]],round(inputs[,"Rob.p.val.1."],pDigits))
}
rownames(estimates[[j]])=inputs[,1]
if(j==1){
combined_names=rownames(estimates[[j]])
}else{
combined_names=c(combined_names,rownames(estimates[[j]]))
}
}
combined_names=unique(combined_names)
combined_outputs=data.frame(matrix(0,nrow=length(combined_names)+1,ncol=values*length(modelNames)))
row.names(combined_outputs)=c("",combined_names)
outputnames=c("estimate")
if(printClassical==TRUE){
outputnames=c(outputnames,"t-ratio(0)")
if(printPVal==TRUE) outputnames=c(outputnames,"p-val(0)")
if(printT1==TRUE){
outputnames=c(outputnames,"t-ratio(1)")
if(printPVal==TRUE) outputnames=c(outputnames,"p-val(1)")
}
}
outputnames=c(outputnames,"Rob.t-ratio(0)")
if(printPVal==TRUE) outputnames=c(outputnames,"Rob.p-val(0)")
if(printT1==TRUE){
outputnames=c(outputnames,"Rob.t-ratio(1)")
if(printPVal==TRUE) outputnames=c(outputnames,"Rob.p-val(1)")
}
combined_outputs[1,]=rep(outputnames,length(modelNames))
for(j in 1:length(modelNames)){
applicable_names=combined_names[combined_names%in%row.names(estimates[[j]])]
combined_outputs[applicable_names,((j-1)*values+1):(j*values)]=estimates[[j]][applicable_names,]
inapplicable_names=combined_names[!(combined_names%in%row.names(estimates[[j]]))]
combined_outputs[inapplicable_names,((j-1)*values+1):(j*values)]=NA
}
for(j in 1:length(modelNames)){
filename=paste(modelNames[[j]],"_output.txt",sep="")
if(!file.exists(filename)) stop("File ",filename," not found!")
lines = tryCatch(readLines(filename),
warning=function(w) x=FALSE,
error=function(e) x=FALSE)
if(is.logical(lines) && lines==FALSE) stop("Could not open file ",filename)
inputvar = grep("Model name", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position = gregexpr(pattern=": ",inputvar)[[1]][1]
otheroutputs[1,j]=substr(inputvar,position+1,nchar(inputvar))
} else otheroutputs[1,j] = NA
k=2
inputvar = grep("Model description", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position = gregexpr(pattern=": ",inputvar)[[1]][1]
otheroutputs[k,j]=substr(inputvar,position+1,nchar(inputvar))
} else otheroutputs[k,j] = NA
k=k+1
inputvar = grep("Number of individuals", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position=gregexpr(pattern=":",inputvar)[[1]][1]
otheroutputs[k,j]=as.double(substr(inputvar,position+1,nchar(inputvar)))
} else otheroutputs[k,j] = NA
k=k+1
inputvar = grep("Number of observations", lines)
if(length(inputvar)==0) inputvar = grep("Number of modelled outcomes", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position=gregexpr(pattern=":",inputvar)[[1]][1]
otheroutputs[k,j]=as.double(substr(inputvar,position+1,nchar(inputvar)))
} else otheroutputs[k,j] = NA
k=k+1
inputvar = grep("Estimated parameters", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position=gregexpr(pattern=":",inputvar)[[1]][1]
otheroutputs[k,j]=as.double(substr(inputvar,position+1,nchar(inputvar)))
} else otheroutputs[k,j] = NA
k=k+1
inputvar = grep("LL\\(final)", lines)
if(length(inputvar)==0) inputvar = grep("LL\\(final, whole model)", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position=gregexpr(pattern=":",inputvar)[[1]][1]
otheroutputs[k,j]=as.double(substr(inputvar,position+1,nchar(inputvar)))
} else otheroutputs[k,j] = NA
k=k+1
inputvar = grep("Adj.Rho-square \\(0\\)", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position=gregexpr(pattern=":",inputvar)[[1]][1]
tmp=substr(inputvar,position+1,nchar(inputvar))
if(tmp==" Not applicable"){
otheroutputs[k,j]=NA
} else {
otheroutputs[k,j]=as.double(tmp)
}
} else otheroutputs[k,j] = NA
k=k+1
inputvar = grep("AIC", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position=gregexpr(pattern=":",inputvar)[[1]][1]
otheroutputs[k,j]=as.double(substr(inputvar,position+1,nchar(inputvar)))
} else otheroutputs[k,j] = NA
k=k+1
inputvar = grep("BIC", lines)
if(length(inputvar)!=0){
inputvar = lines[inputvar]
position=gregexpr(pattern=":",inputvar)[[1]][1]
otheroutputs[k,j]=as.double(substr(inputvar,position+1,nchar(inputvar)))
} else otheroutputs[k,j] = NA
}
otherouputs_new=data.frame(matrix("",nrow=9,ncol=values*length(modelNames)))
rownames(otherouputs_new)=rownames(otheroutputs)
for(j in 1:length(modelNames)){
otherouputs_new[,((j-1)*values+1)]=otheroutputs[,j]
colnames(otherouputs_new)[((j-1)*values+1)]=colnames(otheroutputs)[j]
colnames(otherouputs_new)[((j-1)*values+2)]=""
if(printClassical==TRUE){
colnames(otherouputs_new)[((j-1)*values+3)]=""
}
}
otherouputs_new=rbind(otherouputs_new,rep("",ncol(otherouputs_new)))
rownames(otherouputs_new)[10]=""
filename=paste("model_comparison_",gsub("[: -]", "" , Sys.time(), perl=TRUE),".csv",sep="")
utils::write.table(otherouputs_new, filename, sep = ",", col.names = F, append = T)
#utils::write.table("",filename, sep = ",", col.names = F, append = T)
utils::write.table(combined_outputs, filename, sep = ",", col.names = F, append = T)
cat("Outputs of apollo_combineResults saved to",filename,"\n")
return(invisible(TRUE))
}
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.