Hi guys,
I want to transfer the following code from R into .Call compatible form. How
can i do that?
Thanks!!!
INT sim;
for(i in 1:sim){
if(i>2) genemat <- genemat[,sample(1:ncol(genemat))]
ranklist[,1] <- apply(genemat, 1, function(x){
(mean(x[cols]) -
mean(x[-cols]))/sd(x)})
ranklist <- ranklist[order(ranklist[,1]),]
if(ranklist[1,2]==1) score[i] <- 1/Ns
if(ranklist[1,2]==0) score[i] <- -(1/(N-Ns))
for(j in 2:nrow(ranklist)){
phit <- sum(rep(1/Ns, sum(ranklist[1:j,2]==1)))
pmiss <- sum(rep(1/(N-Ns), sum(ranklist[1:j,2]==0)))
if((phit-pmiss)>score[i]) score[i] <- phit - pmiss
}
}
I tried a little bit, but not enough knowledge in C.
#include <stdio.h>
#include <R.h>
#include <Rdefines.h>
#include <math.h>
SEXP ESscore(SEXP Rgeneset, SEXP Rgenemat, SEXP Rranklist, SEXP sim)
{
int nc = ncols(Rgenemat);
double *geneset = NUMERIC_DATA(Rgeneset);
double *genemat = NUMERIC_DATA(Rgenemat);
SEXP Rscore;
PROTECT(Rscore=NEW_NUMERIC(sim));
double *score = NUMERIC_DATA(Rscore);
for(i=1; i<=sim; i++){
if(i>2) {genemat <- genemat[,sample(1:nc)];}
for(k=1; k<nc )
}
}
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
