The R Project SVN R

Rev

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

Rev 68947 Rev 69855
Line 48... Line 48...
48 
	if (a == 0 || a/b == 0) // point mass 1 at 0 ==> P(X <= x) = 1, all x > 0
 48 
	if (a == 0 || a/b == 0) // point mass 1 at 0 ==> P(X <= x) = 1, all x > 0
49 
	    return R_DT_1;
 49 
	    return R_DT_1;
50 
	if (b == 0 || b/a == 0) // point mass 1 at 1 ==> P(X <= x) = 0, all x < 1
 50 
	if (b == 0 || b/a == 0) // point mass 1 at 1 ==> P(X <= x) = 0, all x < 1
51 
	    return R_DT_0;
 51 
	    return R_DT_0;
52 
	// else, remaining case:  a = b = Inf : point mass 1 at 1/2
 52 
	// else, remaining case:  a = b = Inf : point mass 1 at 1/2
53 
	if (x < 0.5)
 - 
 
54 
	    return R_DT_0;
 - 
 
55 
	// else,  x >= 0.5 :
 - 
 
56 
	    return R_DT_1;
 53 
	if (x < 0.5) return R_DT_0; else return R_DT_1;
57 
    }
 54 
    }
58 
    // Now:  0 < a < Inf;  0 < b < Inf
 55 
    // Now:  0 < a < Inf;  0 < b < Inf
59 
 
 56 
 
60 
    double x1 = 0.5 - x + 0.5, w, wc;
 57 
    double x1 = 0.5 - x + 0.5, w, wc;
61 
    int ierr;
 58 
    int ierr;