You do

  for(i in 1:n){
    a[i] = x[i,1]

where a has length n and x has iter rows, hence you get NA values for all i > iter...

Best,
Uwe Ligges



On 01.03.2016 07:39, smart hendsome via R-help wrote:
Hi R-users,
I have problem regarding my function. My function as below:
gen.m <- function(n,itr){

     set.seed(1234)

     m <- matrix(nrow = n, ncol=4)
        a <- matrix(nrow = n, ncol = 1)
        b <- matrix(nrow = n, ncol = 1)
        C <- matrix(nrow = n, ncol = 1)
        d <- matrix(nrow = n, ncol = 1)
   lambda <- 0.342
     day0 <- 0.1134

       x = matrix(runif(n*itr, 0, 1),nrow = n, ncol = itr)

   y = 0
   for(i in 1:n){
     a[i] = x[i,1]
     b[i] = y - (log (a[i])) / lambda
        y =  b[i]
     }

   if (a[1] < day0){
     C[1] = 0
     }else{
     C[1] = 1}

     for (i in 1:n){
        w = as.integer(b[i])
        if (w <= n){
          C[w] = 1
          }else{
          C[w] = 0}
      }


     for (i in 2:n)
        rain2 <- cbind(a,b,C,d)
        y <- cbind(a,b,C)
        return(y)
    }


   gen.m(31,10)

When i run my output, it gave the result NA. I want the NA values is zero. Hope 
anyone can help me. Thanks so much.

        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.


______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.

Reply via email to