[Bioc-devel] C function is wrong under Windows 7.

[Vincent Carey]

your attachment is automatically scrubbed so we cannot see it.  this
is not really a bioconductor question; you might
consider r-devel.  you would also do well to include the result of
sessionInfo().

On Mon, Oct 24, 2011 at 8:04 AM, Evarist Planet
<evarist.planet at irbbarcelona.org> wrote:
> Dear mailing list,
>
> I have a C function that gives me a wrong result when I run it under Windows
> 7.
>
> This is the code under Linux (RHEL5):
>> library(phenoTest)
>> data(epheno)
>> sign <- sample(featureNames(epheno))[1:20]
>> score <- getFc(epheno)[,1]
>> head(score)
> 1007_s_at   1053_at    117_at    121_at 1255_g_at   1294_at
> -1.183019  1.113544  1.186186 -1.034779 -1.044456 -1.023471
>> s <- which(names(score) %in% sign)
>> es.c <- .Call('getEs',score,s,PACKAGE='phenoTest')
>> head(es.c)
> [1] -0.001020408 -0.002040816 -0.003061224 -0.004081633 -0.005102041
> [6] -0.006122449
>> es.c <- .Call('getEs',score,s,PACKAGE='phenoTest')
>> head(es.c)
> [1] -0.001020408 -0.002040816 -0.003061224 -0.004081633 -0.005102041
> [6] -0.006122449
>
> As you see es.c is correct. I checked it doing the same computation with
> R. It runs without problems under Mac. I run valgrind on the same piece of
> code and got no errors.
>
> This is the same piece of code under Windows 7:
>> library(phenoTest)
>> data(epheno)
>> sign <- sample(featureNames(epheno))[1:20]
>> score <- getFc(epheno)[,1]
>> head(score)
> 1007_s_at   1053_at    117_at    121_at 1255_g_at   1294_at
> -1.183019  1.113544  1.186186 -1.034779 -1.044456 -1.023471
>> s <- which(names(score) %in% sign)
>> es.c <- .Call('getEs',score,s,PACKAGE='phenoTest')
>> head(es.c)
> [1] 1.447208e+215 1.447208e+215 1.447208e+215 1.447208e+215 1.447208e+215
> 1.447208e+215
>> es.c <- .Call('getEs',score,s,PACKAGE='phenoTest')
>> head(es.c)
> [1] 3.176615e+170 3.176615e+170 3.176615e+170 3.176615e+170 3.176615e+170
> 3.176615e+170
>
> es.c is not correct under Windows. It also gives a different result when i
> rerun the same function.
>
> I attached the C code.
>
> Could you please help me to find out what I am doing wrong?
>
> Many thanks in advance,
>
> --
> Evarist Planet
> Research officer, Bioinformatics and Biostatistics unit
> IRB Barcelona
> Tel (+34) 93 402 0553
> Fax (+34) 93 402 0257
>
> evarist.planet at irbbarcelona.org
> http://www.irbbarcelona.org/bioinformatics
>
> _______________________________________________
> Bioc-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/bioc-devel
>
>