dyn.load('irls.so')
is.loaded("irls")

Call.IRLS<-function(Xm,Yv,pv,nits)
{
  tmp=.C("irls",as.double(Xm),as.double(Yv),as.integer(pv),as.integer(nits),as.double(beta))
  return(tmp$beta)
}

library(LogRegLinux)
data(test.data)

nits<-10
pv<-1
Xm<-test.data$X1
Yv<-test.data$Y
Call.IRLS(Xm,Yv,pv,nits)