search for: propto

Hi, I have to draw samples from an asymmetric-Laplace-Normal distribution: f(u|y, x, beta, phi, sigma, tau) \propto exp( - sum( ( abs(lo) + (2*tau-1)*lo )/(2*sigma) ) - 0.5/phi*u^2), where lo = (y - x*beta) and y=(y_1, ..., y_n), x=(x_1, ..., x_n) -- sorry for this huge formula -- A WinBUGS Gibbs sampler and the HI package arms sampler were used with the same initial data for all parameters. I compared the mean...

...ing a bit into the code of this function I can see that value for lambda and its standard error is estimated via estimate <- mean(x) sds <- sqrt(estimate/n) Is this MLE? With my poor math/stat knowledge I thought that MLE for Poisson parameter is (in mixture of LaTeX code) l(\lambda|x) \propto \sum^n_{i=1}(-\lambda + x_iln(\lambda)). Is this really equal to (\sum^n_{i=1} x_i) / n -- Lep pozdrav / With regards, Gregor Gorjanc ---------------------------------------------------------------------- University of Ljubljana PhD student Biotechnical Faculty Zootechnical Department...

Hi, I have count data that have been meddled with enough to make them non integers. Using glm(poisson) returns a "non integer" error but glm(quasipoisson) does not. Just wondering if anyone knows if I am violating the assumptions of a quasipoisson error structure by using these non-integer response data? Thanks! I'd welcome your thoughts and/or references... Mark

...t;spadesuit", 170, + "biarrow", 171, + "leftarrow", 172, + "uparrow", 173, + "rightarrow", 174, + "downarrow", 175, + "degree", 176, + "plusminus", 177, + "gteq", 179, + "times", 180, + "propto", 181, + "partial", 182, + "bullet", 183, + "divide", 184, + "notequal", 185, + "equiv", 186, + "approx", 187, + "ellipsis", 188, + "aleph", 192, + "Imaginary", 193, + "Real", 194, +...

I want to fit some p-values to a beta distribution. But the problem is some of the values have 0s and 1's. I am getting an error if I use the MASS function to do this. Is there anyway to get around this? -- Thanks, Jim. [[alternative HTML version deleted]]

...set of log(Effort) as in your second model is the more > standard way to approach this problem; it corresponds to assuming that > catch is strictly proportional to effort. Adding log(Effort) as a > covariate (as illustrated below) tests whether a power-law model (catch > propto (Effort)^(b+1), b!=0) is a better description of the data.? (In > this case it is not, although the confidence intervals on b are very > wide, indicating that we have very little information -- this is not > surprising since the proportional range of effort is very small >...

Hi I am having problems with the rendering of scientific symbols (mu and degree) in my plots. Whenever I use these symbols they are rendered changed (mu is changed to the proportionality symbol and degree is changed to something resembling a gamma) in the X-device. If I make a pdf of the plot and open the file in Evince or Okular symbols are also rendered wrong, however if I open the file

Hi I am having problems with the rendering of scientific symbols (mu and degree) in my plots. Whenever I use these symbols they are rendered changed (mu is changed to the proportionality symbol and degree is changed to something resembling a gamma) in the X-device. If I make a pdf of the plot and open the file in Evince or Okular symbols are also rendered wrong, however if I open the file

Hi I am having problems with the rendering of scientific symbols (mu and degree) in my plots. Whenever I use these symbols they are rendered changed (mu is changed to the proportionality symbol and degree is changed to something resembling a gamma) in the X-device. If I make a pdf of the plot and open the file in Evince or Okular symbols are also rendered wrong, however if I open the file