[R] error message error

[Uwe Ligges]

On 15.02.2010 06:34, Roslina Zakaria wrote:
> Hi r-users,
>
> I hope somebody can help me to understand the error message.  Here is my code;
>
> ## Newton iteration
> newton_gam<- function(z)
> { n<- length(z)
>    r<- runif(n)
>    tol<- 1E-6
>    cdf<- vector(length=n, mode="numeric")
>    fprime<- vector(length=n, mode="numeric")
>    f<- vector(length=n, mode="numeric")
>
>    for (i in 1:1000)
>    { cdf<- integrate(fprime, lower = 0, upper = z)$value


fprime is a vector of zeros here, not a function...


>      f<- cdf - r
>      # Newton method
>      z<- z - f/fprime

... and given fprime is a vector of zeros, f/fprime this is infinite or 
undefinied.



>      if (any(f<  tol)) break
>     }
>    cbind(z,cdf)
> }
>
> alp<- 2.0165
> bt1<- 29.107 ; bt2<- 41.517
> x1<- d1d4pos[,1];x1[1:10]
> x2<- d1d4pos[,2];x2[1:10]
>
>> x1<- d1d4pos[,1];x1[1:10]


We do not have d1d4pos.

Please ......!


Uwe Ligges


>   [1] 28.4 53.6  1.3 29.5 52.1 65.9 72.6 67.6 58.7 34.5
>> x2<- d1d4pos[,2];x2[1:10]
>   [1]  43.5  56.2   0.3  16.6  71.1  86.3 172.8 111.8  89.9  70.2
>
> z<- (x1/bt1)+(x2/bt2); z
>
> newton_gam(min(z))
> newton_gam(max(z))
>
> Error in -`*tmp*` : invalid argument to unary operator
>> newton_gam(min(z))
>         z          cdf
> [1,] Inf 1.141004e-05
>> newton_gam(max(z))
> Error in integrate(fprime, lower = 0, upper = z) :
>    non-finite function value
>
>> z[100]
> [1] 7.544834
>> newton_gam(z[100])
> Error in integrate(fprime, lower = 0, upper = z) :
>    non-finite function value
>
> Thank you.
>
>
>
>
>
> 	[[alternative HTML version deleted]]
>
>
>
>
> ______________________________________________
> R-help at 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.