[R] Problem when calling FORTRAN subroutine (dll)

[poolloopus at yahoo.com]

Hello,

I am trying to call a FORTRAN subroutine from R. The Fortran code is @:

http://lib.stat.cmu.edu/apstat/206

It performs a bivariate isotonic regression on a rectangular grid (m X n) matrix. I used the g77 compiler and successfully created a dll file and it also loads successfully from R. But somehow the programs fails to run properly. (I do get the correct result when I compile the FORTRAN code). Please see below and note that since my input matrix X is monotonically increasing in both directions, input should be equal to output since the weights are set to 1 for all elements. (also verified this in FORTRAN to make sure that the subroutine works). Output[[10]] is supposed to be the monotonically regressed output, but it does not come out correctly. I don't know if its something to do with the way I am feeding the parameters to the DLL. Thanks in advance for any help I can receive. I am not god with FORTRAN - very very infrequent FORTRAN user. I am using R 2.4.0.

 Regards

Krishna.

> rm(list=ls(all=TRUE))
> dyn.load("c:/smooth.dll")
> is.loaded("smooth")
[1] TRUE  
> NROW<-4
> NCOL<-4
> NDIM<-4
> NCYCLE<-5
> X<-array(c(0.05,0.2,0.35,0.43,0.09,0.29,0.41,0.55,0.2,0.33,0.55,0.7,0.33,0.48,0.6,0.8))
> W<-array(rep(1.0,length(X)))
> X<-matrix(X,nrow=4,ncol=4)
> X<-t(X)
> W<-matrix(W,nrow=4,ncol=4)
> A<-array(0,dim=c(NROW,NCOL,4))
> B<-array(0,dim=c(NDIM,5))
> NCYCLE<-5
> ICYCLE<-5
> G<-(X*0)+0.0
> EPS=0.001
> IFAULT<-5

# Here is how the subroutine says to input the parameters;
#SUBROUTINE SMOOTH(NROW, NCOL, NDIM, X, W, A, B, NCYCLE, ICYCLE, G, EPS, IFAULT)

# Definitions of each of these input parameters are below;

#REAL X(NROW,NCOL), W(NROW,NCOL), A(NROW,NCOL,4), B(NDIM,5), G(NROW,NCOL), EPS, ZERO, DELTA, DELC, DELR, ORD, WW, FRACT
DATA ZERO/0.0/, DELTA/0.00001/, FRACT/0.5/

# Here is my R output;

> output<-.C("smooth_",as.integer(NROW),as.integer(NCOL),as.integer(NDIM),as.real(X),as.real(W),as.real(A),as.real(B),as.integer(NCYCLE),as.integer(ICYCLE),as.real(G),as.real(EPS),as.integer(IFAULT))
> 
> output
[[1]]
[1] 4
[[2]]
[1] 4
[[3]]
[1] 4
[[4]]
[1] 0.05 0.09 0.20 0.33 0.20 0.29 0.33 0.48 0.35 0.41 0.55 0.60 0.43 0.55 0.70
[16] 0.80
[[5]]
[1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
[[6]]
[1]  0.000000e+00  0.000000e+00  0.000000e+00  -5.377600e-19  0.000000e+00
[6]  5.377600e-19  0.000000e+00  0.000000e+00  1.290268e-11  1.803836e-13
[11]  1.333175e-27  0.000000e+00  1.333175e-27  0.000000e+00  -3.434319e-15
[16]  5.092237e-315  1.000000e+00  1.000000e+00  1.000000e+00  1.000000e+00
[21]  1.000000e+00  1.000000e+00  1.000000e+00  1.000000e+00  4.999999e-02
[26]  8.999999e-02  2.000000e-01  3.300000e-01  2.000000e-01  2.900000e-01
[31]  3.300000e-01  4.799998e-01  0.000000e+00  0.000000e+00  0.000000e+00
[36]  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00
[41]  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00
[46]  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00
[51]  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00
[56]  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00
[61]  0.000000e+00  0.000000e+00  0.000000e+00  0.000000e+00
[[7]]
[1] 0.3300000 0.4799998 1.0000002 1.0000002 0.2716666 0.4799998 0.0000000
[8] 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
[15] 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
[[8]]
[1] 5
[[9]]
[1] 5
[[10]]
[1] 0.0795833 0.1145833 0.2008331 0.3100000 0.2008331 0.3100000 0.2716666
[8] 0.4799998 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000
[15] 0.0000000 0.0000000
[[11]]
[1] 0.001
[[12]]
[1] 5
> 
> 
> 



      ____________________________________________________________________________________
Be a better friend, newshound, and