Hi, 

I'm using rlaplace distribution to draw non normal random deviates. I'm using rmutil package for this. I know rlaplace needs location and dispersion parameters but I want to pass diag matrix as dispersion parameters and I'm not getting the correct values. Here is my code snippet:

    dispersion = c(3,.20,.10)
    dispDiag = diag(dispersion)
    location = c(0,0,0)
    output = (rlaplace(3, location, dispDiag))


Here is what I get:

 -0.8475811  0.0000000  0.0000000

I dont know why it is returning 0's at 2nd and 3rd place. If I change the values in location vector, it simply picks up the value and replace it at 2nd and 3rd place.

You help will be highly appreciated

Victor


      
	[[alternative HTML version deleted]]

