The R Project SVN R

Rev

Rev 8431 | Rev 42302 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 8431 Rev 20449
Line 33... Line 33...
33 
    if (ISNAN(x) || ISNAN(logmean) || ISNAN(logsd))
 33 
    if (ISNAN(x) || ISNAN(logmean) || ISNAN(logsd))
34 
	return x + logmean + logsd;
 34 
	return x + logmean + logsd;
35 
#endif
 35 
#endif
36 
    if(logsd <= 0) ML_ERR_return_NAN;
 36 
    if(logsd <= 0) ML_ERR_return_NAN;
37 
 
 37 
 
38 
    if(x == 0) return R_D__0;
 38 
    if(x <= 0) return R_D__0;
39 
 
 39 
 
40 
    y = (log(x) - logmean) / logsd;
 40 
    y = (log(x) - logmean) / logsd;
41 
    return (give_log ?
 41 
    return (give_log ?
42 
	    -(M_LN_SQRT_2PI   + 0.5 * y * y + log(x * logsd)) :
 42 
	    -(M_LN_SQRT_2PI   + 0.5 * y * y + log(x * logsd)) :
43 
	    M_1_SQRT_2PI * exp(-0.5 * y * y)  /	 (x * logsd));
 43 
	    M_1_SQRT_2PI * exp(-0.5 * y * y)  /	 (x * logsd));