Look at what 'loww' is in your function. You are passing in a vector of
length 5 and that is what is being passed into integrate for 'lower' which
is expecting an object of length 1. So determine what you want to pass in
for lower (e.g., min(loww)).
On Fri, Jan 1, 2010 at 7:44 PM, roger king <rking...@gmail.com> wrote:
> I am a novice user of "R" and I'm learning with R version 2.8.1, using
> WinEdt_1.8.1, under Widows Vista Home Version.
>
> ## The test function below, from a vector input, returns vector values:
> # and it contains an "ifelse"statement
> TEST<- function(x) {
> low<- -x^2
> up<- x^4
> ifelse(x>=0,up,low )
> }
> u<- seq(-1,1,0.5)
> TEST(u)
> ###########################################################
> # #The following function, also containing "ifelse", returns correct
> individual values
> # (confirmed independently from another program).
> # Vector input, however, generates the following message:
> # "In if (is.finite(lower)) { :the condition has length > 1 and only the
> first element
> # will be used"
> # And, in fact, correct multiple values are not returned.
> ########################################################
> ASEPF<- function(y, theta, sigma, alpha, kappa){
>
> loww<-((theta - y)/(sigma* kappa))^alpha
>
> upp<- (kappa *(y - theta)/sigma)^alpha
>
> k<- 1/(1 + kappa^2)
>
> a<- 1/alpha
>
> integrand<- function(y){y^(a-1)* exp(-y)/gamma(1/alpha)}
>
> GL<- integrate(integrand, lower = loww, upper= Inf)$value
>
> GU<- integrate(integrand, lower = upp, upper= Inf)$value
> ifelse(y>=0,1 - k*GU, (kappa^2) *k * GL)
> }
> ASEPF(u,0,1,2,0.5)
> Could someone kindly give me some advice on how to
> rewrite the function ASEPF so that it returns a vector of values?
>
> many thanks,
>
> Rik King
>
> rking...@gmail.com
>
> (ex Mathematics Department,
> University of Western Sydney, Australia)
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html<http://www.r-project.org/posting-guide.html>
> and provide commented, minimal, self-contained, reproducible code.
>
--
Jim Holtman
Cincinnati, OH
+1 513 646 9390
What is the problem that you are trying to solve?
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
