Displaying 1 result from an estimated 1 matches for "xmod".
Did you mean:
mod
2012 Apr 23
2
automating a script to read a file
...087277)))))))))))
if (z >= 0) return(ans) else return(-1 * ans)
}
pdf <- function(x) {
return((1 / sqrt(2 * pi)) * exp(-1 * ((x * x) / 2)))
}
cdf <- function(x) {
return( 0.5 * (1 + erf(x / sqrt(2))))
}
spdf <- function(x, skewlocation, skewscale, skewshape) {
xmod <- (x - skewlocation) / skewscale
return( 2 * pdf(xmod) * (cdf(xmod * skewshape)))
}
#### Plotting Function ################
plot.spdf <- function(xmin, xmax, skewlocation, skewscale, skewshape,
skewmax, skewtitle) {
if(missing(skewtitle)) {
plottitle <- "Skewed P...