search for: betamean

Displaying 1 result from an estimated 1 matches for "betamean".

Did you mean: beaman
2007 Jan 04
1
Parameter changes and segfault when calling C code through .Call
...= &sigmaView.matrix; 1: double logans = gsl_MB_mvnorm_logpdf(x, mu, sigma, k); // <-call logpdf here PROTECT(res=allocVector(REALSXP,1)); REAL(res)[0] = logans; UNPROTECT(1); return(res); } The logpdf function is here double gsl_MB_mvnorm_logpdf(gsl_vector * beta, gsl_vector * betaMean, gsl_matrix * sigma, int k) { // computes density of multivariate normal vector at vector beta, with mean betaMean and cov sigma double logdetSigma = 0; double res; double * kern; int i, err; // pointer to Cholesky decomp of sigma gsl_matrix * sigmaChol = gsl_matrix_alloc(k, k); //...