Nothing
# R package 'abremout'
# output methods for the abrem object
# May 2013, Jurgen Symynck
# Copyright 2013, Jurgen Symynck
#
# For the latest version of this file, check the Subversion repository at
# http://r-forge.r-project.org/projects/abernethy/
#
# Disclaimer:
# The author is not affiliated with Dr. Abernethy or Wes Fulton - CEO of
# Fulton Findings(TM) and author of the software package SuperSMITH
#-------------------------------------------------------------------------------
# This program is free software: you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation, either version 3 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program. If not, see <http://www.gnu.org/licenses/>.
#
# For more info on this software and its predecesser, the "weibulltoolkit",
# consult following documents:
#
# - "Weibull analysis using R, in a nutshell",
# (Jurgen Symynck, Filip De Bal, 2010)
# - "Monte Carlo pivotal confidence bounds for Weibull analysis
# with implementations in R",
# (Jurgen Symynck, Filip De Bal, 2011)
#
# +-----------------------------------+
# | execute this program with R: |
# | http://www.r-project.org/ |
# +-----------------------------------+
#
legendconf <- function(fit,conftype,opa,opp){
if(!is.null(fit$options.abremplot)){
oppfit <- modifyList(opp,fit$options.abremplot)
}else{oppfit <- opp}
#opafit <- modifyList(opa,fit$options.abremplot)
if(identical(tolower(conftype),"blives")){
if(!is.null(fit$conf$blives)){
fun0 <- function(blicon){
le <- list()
le[[1]] <- paste0(conftype,", type = ",
ifelse(is.null(blicon$type),"NA",
paste0("\"",blicon$type,"\"")))
le[[2]] <- paste0(" CL = ",
ifelse(is.null(blicon$cl),"NA",
paste0(signif(blicon$cl*100,4)," [%]")),
ifelse(is.null(blicon$S),"",
paste0(", S = ",blicon$S)))
params <- unlist(list(beta=fit$beta,eta=fit$eta,lambda=fit$lambda,
meanlog=fit$meanlog,sdlog=fit$sdlog,rate=fit$rate))
le[[3]] <- lapply(opa$blives,Blifestring,blicon,opp$signif,params)
na.omit(unlist(le))
}
lapply(fit$conf$blives,fun0)
}else{NA}
}
# NA
}
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.