Ejemplo n.º 1
0
// [[Rcpp::export]]
double exp_test(NumericVector x){
  int n = x.size();
  NumericVector z = 1-exp(-x/mean(x));
  z.sort();
  double w_pom = 0.0;
  for(int i = 0; i<n; i++){
    w_pom += (2.0*i+1)*(log(z[i])+log(1-z[n-i-1]));
  }
  return -n-w_pom/n;
}
Ejemplo n.º 2
0
// [[Rcpp::export]]
NumericVector rollP3(NumericVector A, int start, int end) {
  NumericVector B = A[seq(start-1, end-1)] ;
  return B.sort() ;
}