[Rcpp-commits] r1954 - pkg/Rcpp pkg/Rcpp/inst pkg/Rcpp/inst/include/Rcpp/stats scripts
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Mon Aug 9 11:38:31 CEST 2010
Author: romain
Date: 2010-08-09 11:38:31 +0200 (Mon, 09 Aug 2010)
New Revision: 1954
Added:
pkg/Rcpp/inst/include/Rcpp/stats/logis.h
Modified:
pkg/Rcpp/NEWS
pkg/Rcpp/TODO
pkg/Rcpp/inst/ChangeLog
pkg/Rcpp/inst/include/Rcpp/stats/stats.h
scripts/stats.R
Log:
dpq logis
Modified: pkg/Rcpp/NEWS
===================================================================
--- pkg/Rcpp/NEWS 2010-08-09 07:16:41 UTC (rev 1953)
+++ pkg/Rcpp/NEWS 2010-08-09 09:38:31 UTC (rev 1954)
@@ -1,6 +1,6 @@
0.8.6 (future)
- o new sugar functions (dpq)(unif|norm|gamma|chisq|lnorm|weibull|pois|binom|t|beta)
+ o new sugar functions (dpq)(unif|norm|gamma|chisq|lnorm|weibull|logis|pois|binom|t|beta)
o new vignette Rcpp-quickref : quick reference guide of Rcpp API
Modified: pkg/Rcpp/TODO
===================================================================
--- pkg/Rcpp/TODO 2010-08-09 07:16:41 UTC (rev 1953)
+++ pkg/Rcpp/TODO 2010-08-09 09:38:31 UTC (rev 1954)
@@ -67,7 +67,7 @@
o other statistical distribution functions : nchisq, f,
multinom, cauchy, exp, geom, hyper, nbinom, nbinom_mu,
- logis, nbeta, nf, nt, tukey, wilcox, signrank
+ nbeta, nf, nt, tukey, wilcox, signrank
o sum : deal with missing values
Modified: pkg/Rcpp/inst/ChangeLog
===================================================================
--- pkg/Rcpp/inst/ChangeLog 2010-08-09 07:16:41 UTC (rev 1953)
+++ pkg/Rcpp/inst/ChangeLog 2010-08-09 09:38:31 UTC (rev 1954)
@@ -1,7 +1,7 @@
2010-08-09 Romain Francois <romain at r-enthusiasts.com>
- * inst/include/Rcpp/stats/lnorm.h: Added d-p-q (lnorm|weibull), generated
- by the script
+ * inst/include/Rcpp/stats/lnorm.h: Added d-p-q (lnorm|weibull|logis)
+ generated by the script
2010-08-06 Douglas Bates <bates at stat.wisc.edu>
Added: pkg/Rcpp/inst/include/Rcpp/stats/logis.h
===================================================================
--- pkg/Rcpp/inst/include/Rcpp/stats/logis.h (rev 0)
+++ pkg/Rcpp/inst/include/Rcpp/stats/logis.h 2010-08-09 09:38:31 UTC (rev 1954)
@@ -0,0 +1,118 @@
+
+// -*- mode: C++; c-indent-level: 4; c-basic-offset: 4; tab-width: 4 -*-
+//
+// auto generated file (from script/stats.R)
+//
+// logis.h: Rcpp R/C++ interface class library --
+//
+// Copyright (C) 2010 Douglas Bates, Dirk Eddelbuettel and Romain Francois
+//
+// This file is part of Rcpp.
+//
+// Rcpp is free software: you can redistribute it and/or modify it
+// under the terms of the GNU General Public License as published by
+// the Free Software Foundation, either version 2 of the License, or
+// (at your option) any later version.
+//
+// Rcpp is distributed in the hope that it will be useful, but
+// WITHOUT ANY WARRANTY; without even the implied warranty of
+// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+// GNU General Public License for more details.
+//
+// You should have received a copy of the GNU General Public License
+// along with Rcpp. If not, see <http://www.gnu.org/licenses/>.
+
+#ifndef Rcpp__stats__logis_h
+#define Rcpp__stats__logis_h
+
+namespace Rcpp {
+namespace stats {
+namespace impl {
+
+ template <bool NA, typename T>
+ class DLogis : public Rcpp::VectorBase< REALSXP, NA, DLogis<NA,T> >{
+ public:
+ typedef typename Rcpp::VectorBase<REALSXP,NA,T> VEC_TYPE;
+
+ DLogis( const VEC_TYPE& vec_, double location_ = 0.0, double scale_ = 1.0 , bool log_ = false ) :
+ vec(vec_), location(location_), scale(scale_) , log(log_) {}
+
+ inline double operator[]( int i) const {
+ return ::dlogis( vec[i], location, scale , log );
+ }
+
+ inline int size() const { return vec.size(); }
+
+ private:
+ const VEC_TYPE& vec;
+ double location; double scale;
+ int log;
+
+ };
+
+ template <bool NA, typename T>
+ class PLogis : public Rcpp::VectorBase< REALSXP, NA, PLogis<NA,T> >{
+ public:
+ typedef typename Rcpp::VectorBase<REALSXP,NA,T> VEC_TYPE;
+
+ PLogis( const VEC_TYPE& vec_, double location_ = 0.0, double scale_ = 1.0 ,
+ bool lower_tail = true, bool log_ = false ) :
+ vec(vec_), location(location_), scale(scale_) , lower(lower_tail), log(log_) {}
+
+ inline double operator[]( int i) const {
+ return ::plogis( vec[i], location, scale, lower, log );
+ }
+
+ inline int size() const { return vec.size(); }
+
+ private:
+ const VEC_TYPE& vec;
+ double location; double scale;
+ int lower, log;
+
+ };
+
+ template <bool NA, typename T>
+ class QLogis : public Rcpp::VectorBase< REALSXP, NA, QLogis<NA,T> >{
+ public:
+ typedef typename Rcpp::VectorBase<REALSXP,NA,T> VEC_TYPE;
+
+ QLogis( const VEC_TYPE& vec_, double location_ = 0.0, double scale_ = 1.0 ,
+ bool lower_tail = true, bool log_ = false ) :
+ vec(vec_), location(location_), scale(scale_), lower(lower_tail), log(log_) {}
+
+ inline double operator[]( int i) const {
+ return ::qlogis( vec[i], location, scale, lower, log );
+ }
+
+ inline int size() const { return vec.size(); }
+
+ private:
+ const VEC_TYPE& vec;
+ double location; double scale;
+ int lower, log;
+
+ };
+
+} // impl
+
+template <bool NA, typename T>
+inline impl::DLogis<NA,T> dlogis( const Rcpp::VectorBase<REALSXP,NA,T>& x, double location_ = 0.0, double scale_ = 1.0, bool log = false ) {
+ return impl::DLogis<NA,T>( x, location_, scale_, log );
+}
+
+template <bool NA, typename T>
+inline impl::PLogis<NA,T> plogis( const Rcpp::VectorBase<REALSXP,NA,T>& x, double location_ = 0.0, double scale_ = 1.0, bool lower = true, bool log = false ) {
+ return impl::PLogis<NA,T>( x, location_, scale_, lower, log );
+}
+
+template <bool NA, typename T>
+inline impl::QLogis<NA,T> qlogis( const Rcpp::VectorBase<REALSXP,NA,T>& x, double location_ = 0.0, double scale_ = 1.0, bool lower = true, bool log = false ) {
+ return impl::QLogis<NA,T>( x, location_, scale_, lower, log );
+}
+
+}
+}
+
+#endif
+
Modified: pkg/Rcpp/inst/include/Rcpp/stats/stats.h
===================================================================
--- pkg/Rcpp/inst/include/Rcpp/stats/stats.h 2010-08-09 07:16:41 UTC (rev 1953)
+++ pkg/Rcpp/inst/include/Rcpp/stats/stats.h 2010-08-09 09:38:31 UTC (rev 1954)
@@ -30,6 +30,7 @@
#include <Rcpp/stats/t.h>
#include <Rcpp/stats/lnorm.h>
#include <Rcpp/stats/weibull.h>
+#include <Rcpp/stats/logis.h>
#include <Rcpp/stats/binom.h>
#include <Rcpp/stats/pois.h>
Modified: scripts/stats.R
===================================================================
--- scripts/stats.R 2010-08-09 07:16:41 UTC (rev 1953)
+++ scripts/stats.R 2010-08-09 09:38:31 UTC (rev 1954)
@@ -167,7 +167,7 @@
# "min" = list( type = "double", default = "0.0" ),
# "max" = list( type = "double", default = "1.0" )
# ) )
-#
+#
# cook( "gamma", params = list(
# "shape" = list( type = "double" ),
# "scale" = list( type = "double", default = "1.0" )
@@ -177,13 +177,18 @@
# "df" = list( type = "double" )
# ) )
#
-#cook( "lnorm", params = list(
+# cook( "lnorm", params = list(
# "meanlog" = list( type = "double", default = "0.0" ),
# "sdlog" = list( type = "double", default = "1.0" )
# ) )
+#
+# cook( "weibull", params = list(
+# "shape" = list( type = "double" ),
+# "scale" = list( type = "double", default = "1.0" )
+# ) )
-cook( "weibull", params = list(
- "shape" = list( type = "double" ),
+cook( "logis", params = list(
+ "location" = list( type = "double", default = "0.0" ),
"scale" = list( type = "double", default = "1.0" )
) )
More information about the Rcpp-commits
mailing list