|
1 | | -// -*- mode: C++; c-indent-level: 4; c-basic-offset: 4; tab-width: 4 -*- |
| 1 | +// -*- mode: C++; c-indent-level: 4; c-basic-offset: 4; indent-tabs-mode: nil; -*- |
2 | 2 | // |
3 | 3 | // rf.h: Rcpp R/C++ interface class library -- |
4 | 4 | // |
5 | | -// Copyright (C) 2010 - 2013 Douglas Bates, Dirk Eddelbuettel and Romain Francois |
| 5 | +// Copyright (C) 2010 - 2016 Douglas Bates, Dirk Eddelbuettel and Romain Francois |
6 | 6 | // |
7 | 7 | // This file is part of Rcpp. |
8 | 8 | // |
|
23 | 23 | #define Rcpp__stats__random_rf_h |
24 | 24 |
|
25 | 25 | namespace Rcpp { |
26 | | - namespace stats { |
| 26 | +namespace stats { |
27 | 27 |
|
| 28 | +class FGenerator_Finite_Finite : public ::Rcpp::Generator<double> { |
| 29 | +public: |
28 | 30 |
|
29 | | - class FGenerator_Finite_Finite : public ::Rcpp::Generator<double> { |
30 | | - public: |
| 31 | + FGenerator_Finite_Finite( double n1_, double n2_ ) : |
| 32 | + n1__2(n1_ / 2.0 ), n2__2(n2_ / 2.0 ), ratio(n2_/n1_) {} |
31 | 33 |
|
32 | | - FGenerator_Finite_Finite( double n1_, double n2_ ) : |
33 | | - n1__2(n1_ / 2.0 ), n2__2(n2_ / 2.0 ), ratio(n2_/n1_) {} |
| 34 | + inline double operator()() const { |
| 35 | + // here we know that both n1 and n2 are finite |
| 36 | + // return ( ::rchisq( n1 ) / n1 ) / ( ::rchisq( n2 ) / n2 ); |
| 37 | + return ratio * ::Rf_rgamma( n1__2, 2.0 ) / ::Rf_rgamma( n2__2, 2.0 ) ; |
| 38 | + } |
34 | 39 |
|
35 | | - inline double operator()() const { |
36 | | - // here we know that both n1 and n2 are finite |
37 | | - // return ( ::rchisq( n1 ) / n1 ) / ( ::rchisq( n2 ) / n2 ); |
38 | | - return ratio * ::Rf_rgamma( n1__2, 2.0 ) / ::Rf_rgamma( n2__2, 2.0 ) ; |
39 | | - } |
| 40 | +private: |
| 41 | + double n1__2, n2__2, ratio ; |
| 42 | +}; |
40 | 43 |
|
41 | | - private: |
42 | | - double n1__2, n2__2, ratio ; |
43 | | - } ; |
| 44 | +class FGenerator_NotFinite_Finite : public ::Rcpp::Generator<double> { |
| 45 | +public: |
44 | 46 |
|
| 47 | + FGenerator_NotFinite_Finite( double n2_ ) : n2( n2_), n2__2(n2_ / 2.0 ) {} |
45 | 48 |
|
46 | | - class FGenerator_NotFinite_Finite : public ::Rcpp::Generator<double> { |
47 | | - public: |
| 49 | + inline double operator()() const { |
| 50 | + // return n2 / ::rchisq( n2 ) ; |
| 51 | + return n2 / ::Rf_rgamma( n2__2, 2.0 ) ; |
| 52 | + } |
48 | 53 |
|
49 | | - FGenerator_NotFinite_Finite( double n2_ ) : n2( n2_), n2__2(n2_ / 2.0 ) {} |
| 54 | +private: |
| 55 | + double n2, n2__2 ; |
| 56 | +}; |
50 | 57 |
|
51 | | - inline double operator()() const { |
52 | | - // return n2 / ::rchisq( n2 ) ; |
53 | | - return n2 / ::Rf_rgamma( n2__2, 2.0 ) ; |
54 | | - } |
| 58 | +class FGenerator_Finite_NotFinite : public ::Rcpp::Generator<double> { |
| 59 | +public: |
55 | 60 |
|
56 | | - private: |
57 | | - double n2, n2__2 ; |
58 | | - } ; |
| 61 | + FGenerator_Finite_NotFinite( double n1_ ) : n1(n1_), n1__2(n1_ / 2.0 ) {} |
59 | 62 |
|
| 63 | + inline double operator()() const { |
| 64 | + // return ::rchisq( n1 ) / n1 ; |
| 65 | + return ::Rf_rgamma( n1__2, 2.0 ) / n1 ; |
| 66 | + } |
60 | 67 |
|
61 | | - class FGenerator_Finite_NotFinite : public ::Rcpp::Generator<double> { |
62 | | - public: |
63 | | - |
64 | | - FGenerator_Finite_NotFinite( double n1_ ) : n1(n1_), n1__2(n1_ / 2.0 ) {} |
65 | | - |
66 | | - inline double operator()() const { |
67 | | - // return ::rchisq( n1 ) / n1 ; |
68 | | - return ::Rf_rgamma( n1__2, 2.0 ) / n1 ; |
69 | | - } |
70 | | - |
71 | | - private: |
72 | | - double n1, n1__2 ; |
73 | | - } ; |
74 | | - |
75 | | - } // stats |
| 68 | +private: |
| 69 | + double n1, n1__2 ; |
| 70 | +}; |
76 | 71 |
|
| 72 | +} // stats |
77 | 73 | } // Rcpp |
78 | 74 |
|
79 | 75 | #endif |
0 commit comments