The R Project SVN R

Rev

Rev 80271 | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 80271 Rev 86507
Line 62... Line 62...
62 
    if(p == R_DT_1) return ML_POSINF;
 62 
    if(p == R_DT_1) return ML_POSINF;
63 
 
 63 
 
64 
    double
 64 
    double
65 
	mu = lambda,
 65 
	mu = lambda,
66 
	sigma = sqrt(lambda),
 66 
	sigma = sqrt(lambda),
67 
	// had gamma = sigma; PR#8058 should be kurtosis which is mu^-0.5 = 1/sigma
 67 
	// had gamma = sigma; PR#8058 should be skewness which is mu^-0.5 = 1/sigma
68 
	gamma = 1.0/sigma;
 68 
	gamma = 1.0/sigma;
69 
 
 69 
 
70 
     R_DBG_printf("qpois(p=%.12g, lambda=%.15g, l.t.=%d, log=%d):"
 70 
     R_DBG_printf("qpois(p=%.12g, lambda=%.15g, l.t.=%d, log=%d):"
71 
		  " mu=%g, sigma=%g, gamma=%g;\n",
 71 
		  " mu=%g, sigma=%g, gamma=%g;\n",
72 
		  p, lambda, lower_tail, log_p, mu, sigma, gamma);
 72 
		  p, lambda, lower_tail, log_p, mu, sigma, gamma);