Displaying 1 result from an estimated 1 matches for "unnorm1".
Did you mean:
unnorm
2003 Apr 16
2
trying to plot function using curve
...xp(theta*x)-exp(-theta*x))/x)
}
return(unnorm(x)/
integrate(unnorm,lower = -theta, upper = theta,
subdivisions=1000)$value)
}
}
if(len==3)
{
if((pos==1) || (pos==3))
{
unnorm1 <- function(y) ##here y is 2-dim
{
ifelse(y[2] == 0,2*theta,
(1/y[2])*(exp(y[1]*y[2] + y[2]*theta)
- exp(y[1]*y[2] - y[2]*theta)))
}
unnorm2 <- function(y)
{...