Re: [R] Writing R package that call Fortran codes

2016-06-16 Thread Kodalore Vijayan, Vineetha W
Thanks for the comment. It resolved now. I have edited my namespace file. -V On Thu, Jun 16, 2016 at 3:20 PM, Uwe Ligges wrote: > > > On 16.06.2016 20:17, Kodalore Vijayan, Vineetha W wrote: > >> Hi, >> >> Its running in R separately using dyn.load(). Thanks!

Re: [R] maptools readShapeSpatial not working in R 3.3.0

2016-06-16 Thread MacQueen, Don
You might want to ask about this on R-sig-geo. I use readOGR() from the rgdal package to read shapefiles, and it is currently working for me in R 3.3.0. -Don -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062 On 6/16/16, 12:34

[R] maptools readShapeSpatial not working in R 3.3.0

2016-06-16 Thread Bravo, Mercedes
Hello, I have some R code that uses readShapeSpatial in the maptools package to read a shapefile from ArcGIS. This code worked as recently as May 2016, when I was using R version 3.2.3. Then, R was updated to version 3.3.0. Now, when I run the same code that worked three weeks ago, I get an

Re: [R] dplyr's arrange function - 3 solutions received - 1 New Question

2016-06-16 Thread Muhuri, Pradip (AHRQ/CFACT)
Hello David, Your revisions to the earlier code have given me desired results. library("gtools") mydata[ mixedorder(mydata$prevalence_c, decreasing=TRUE), c("indicator", "prevalence_c") ] Thanks, Pradip Pradip K. Muhuri, AHRQ/CFACT 5600 Fishers Lane # 7N142A, Rockville, MD 20857 Tel:

Re: [R] Writing R package that call Fortran codes

2016-06-16 Thread Kodalore Vijayan, Vineetha W
Thank you Don. I did send an email to R-sig-mac. -Vineetha On Thu, Jun 16, 2016 at 10:36 AM, MacQueen, Don wrote: > You might want to take this question to R-sig-mac. > > -- > Don MacQueen > > Lawrence Livermore National Laboratory > 7000 East Ave., L-627 > Livermore, CA

Re: [R] Writing R package that call Fortran codes

2016-06-16 Thread Kodalore Vijayan, Vineetha W
Thanks for the suggestion Berend. I just tried with numeric(n) and its working for x and y. But my issue is different and its not resolved yet. -Vineetha On Thu, Jun 16, 2016 at 1:19 PM, Berend Hasselman wrote: > > > On 16 Jun 2016, at 21:04, Eduardo M. A. M.Mendes

Re: [R] Writing R package that call Fortran codes

2016-06-16 Thread Berend Hasselman
> On 16 Jun 2016, at 21:04, Eduardo M. A. M.Mendes wrote: > > >> On Jun 16, 2016, at 3:50 PM, Berend Hasselman wrote: >> >> >>> On 16 Jun 2016, at 20:02, Eduardo M. A. M.Mendes >>> wrote: >>> >>> Thanks Bill for pointing this

Re: [R] Writing R package that call Fortran codes

2016-06-16 Thread Eduardo M. A. M.Mendes
> On Jun 16, 2016, at 3:50 PM, Berend Hasselman wrote: > > >> On 16 Jun 2016, at 20:02, Eduardo M. A. M.Mendes >> wrote: >> >> Thanks Bill for pointing this out. I haven’t noticed it. >> >> Vineetha, try as.double(rep(0,n)) or as.matrix(rep(0,n)), >>

Re: [R] Writing R package that call Fortran codes

2016-06-16 Thread Berend Hasselman
> On 16 Jun 2016, at 20:02, Eduardo M. A. M.Mendes wrote: > > Thanks Bill for pointing this out. I haven’t noticed it. > > Vineetha, try as.double(rep(0,n)) or as.matrix(rep(0,n)), > Why not simply numeric(n)? Berend Hasselman

Re: [R] Writing R package that call Fortran codes

2016-06-16 Thread Kodalore Vijayan, Vineetha W
Hi, Its running in R separately using dyn.load(). Thanks! But my original issue is still there. I still can not build my package. My source code file name in SRC folder is "randomepi.f95" , the package name is "NEpidemic". and the wrapper function in the R folder is "randomepi.r". Not sure if

Re: [R] Writing R package that call Fortran codes

2016-06-16 Thread Kodalore Vijayan, Vineetha W
Hi, @Eduardo: I do have statements like, double preciosion :: y(n) in my fortran code. Will change the code and run again. @ William : Yes I was initializing X ,Y,and tau, each have a vector of length "n". Thanks, Vineetha On Thu, Jun 16, 2016 at 11:54 AM, William Dunlap

Re: [R] Writing R package that call Fortran codes

2016-06-16 Thread Eduardo M. A. M.Mendes
Thanks Bill for pointing this out. I haven’t noticed it. Vineetha, try as.double(rep(0,n)) or as.matrix(rep(0,n)), Ed > On Jun 16, 2016, at 2:54 PM, William Dunlap wrote: > > > 1: .Fortran("dataxy", n = as.integer(n), tmax = as.integer(tmax), > > alpha =

Re: [R] Writing R package that call Fortran codes

2016-06-16 Thread William Dunlap via R-help
> 1: .Fortran("dataxy", n = as.integer(n), tmax = as.integer(tmax), alpha = as.double(alpha), beta = as.double(beta), x = as.double(0, length = n), y = as.double(0, length = n), tau = as.integer(0, length = n)) Are you expecting that as.integer(0, length=n) and as.double(0, length =

Re: [R] extracting coefficients from ar() output

2016-06-16 Thread peter dalgaard
> On 16 Jun 2016, at 17:07 , William Dunlap via R-help > wrote: > > help(ar) should tell you how to get the coefficients. If, like me, you > don't > read help files, you can use str() to look at the structure of ar's output. Also notice that the output of rollapply is

Re: [R] inverse table

2016-06-16 Thread Michael Friendly
That function, expand.dft is in the vcdExtra package On 6/15/2016 12:42 PM, Marc Schwartz wrote: There is a function called expand.dft(), which I posted some years ago, which is a modification of a prior version, posted a few years before that. The updated version is here:

Re: [R] Writing R package that call Fortran codes

2016-06-16 Thread Eduardo M. A. M.Mendes
Hi Vineetha It seems that there are some problems in the Fortran code as far as allocating memory is concerned. If you have statements such as double preciosion :: y(n) you will need to allocate memory for the vector within the Fortran code or in the R wrapper-function. It is hard to

Re: [R] Writing R package that call Fortran codes

2016-06-16 Thread Kodalore Vijayan, Vineetha W
Hi Eduardo, Thanks for your comments. I haven't tried the way you told me. Now when I tried, got the following error: *** caught segfault *** address 0x0, cause 'memory not mapped' Traceback: 1: .Fortran("dataxy", n = as.integer(n), tmax = as.integer(tmax), alpha = as.double(alpha), beta =

Re: [R] dplyr's arrange function - 3 solutions received - 1 New Question

2016-06-16 Thread David Winsemius
> On Jun 16, 2016, at 6:12 AM, Muhuri, Pradip (AHRQ/CFACT) > wrote: > > Hello, > > I got 3 solutions to my earlier code. Thanks to the contributors. May I > bring your attention to a new question below (with respect to David's > solution)? > > 1) Thanks to

Re: [R] Writing R package that call Fortran codes

2016-06-16 Thread MacQueen, Don
You might want to take this question to R-sig-mac. -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062 On 6/15/16, 1:20 PM, "R-help on behalf of Kodalore Vijayan, Vineetha W"

[R] Mplus in R

2016-06-16 Thread Birhanu Worku
Hi everyone why this code is not running ? here is it; library(MplusAutomation) res <- mplusModeler( mplusObject( MODEL = "mpg ON wt hp; wt WITH hp;", rdata = mtcars), "mtcars.dat", run = 1L) summary(res) note that mtcars is the data in r The out put; "No R variables to use specified.

Re: [R] dplyr's arrange function - 3 solutions received - 1 New Question

2016-06-16 Thread Muhuri, Pradip (AHRQ/CFACT)
Hello, I got 3 solutions to my earlier code. Thanks to the contributors. May I bring your attention to a new question below (with respect to David's solution)? 1) Thanks to Daniel Nordlund for the tips - replacing leading space with a 0 in the data. 2) Thanks to David Winsemius for his

Re: [R] extracting coefficients from ar() output

2016-06-16 Thread William Dunlap via R-help
help(ar) should tell you how to get the coefficients. If, like me, you don't read help files, you can use str() to look at the structure of ar's output. > str(a <- ar(sin(1:30), aic=TRUE)) List of 14 $ order : int 2 $ ar : num [1:2] 1.011 -0.918 $ var.pred: num 0.0654 $

Re: [R] Package rms: c-statistic from lrm function with weights

2016-06-16 Thread David Winsemius
> On Jun 15, 2016, at 6:44 PM, Marie-Pierre Sylvestre > wrote: > > Dear list, > > I am using the lrm function from the rms package to estimate a logistic > model with weights. The c-statistic (or area under the curve) is part of > the lrm output. > > To understand how

Re: [R] Y in Kohonen xyf function

2016-06-16 Thread ch.elahe via R-help
Hi Kimmo, Thanks for your reply, Here is a part of my df: 'data.frame': 562 obs. of 128 variables $ TE :int 37 37 35 34 37 37 35 33 32 ... $ TR :int 11 11 8 13 11 8 15 12 8 . $ BW :int 150 191 128 145 200 191 $speed :int 4 4

Re: [R] Building a binary vector out of dichotomous variables

2016-06-16 Thread Tom Wright
Does this do what you want? as.numeric(with(ds_example,paste(1,year2013,year2014,year2015,sep=''))) On Thu, Jun 16, 2016 at 8:57 AM, wrote: > Hi All, > > I need to build a binary vector made of a set of dichotomous variables. > > What I have so far is: > > -- cut -- > >

Re: [R] help for fine mappting

2016-06-16 Thread PIKAL Petr
Hi Did you test my suggestions? If not, why not? If yes, in what respect they did not work? sel <- map$POS >= ref$POS_start[1] & map$POS < ref$POS_end[1] result1 <- sum( map$post_prob[sel] ) result2 <- length( map$post_prob[sel] ) result3 <- min( map$p[sel] ) should give you desired values. It

Re: [R] Installing Caret

2016-06-16 Thread Max Kuhn
The problem is not with `caret. Your output says: > installation of package ‘minqa’ had non-zero exit status `caret` has a dependency that has a dependency on `minqa`. The same is true for `RcppEigen` and the others. What code did you use to do the install? What OS and version or R etc? On

Re: [R-es] Levantar archivos de Qlikview (.qvd) en R

2016-06-16 Thread Mauricio Monsalvo
Muchas gracias, Javier. Mi problema es que yo tampoco usé nunca qlikview y me quiero "salir" de él para poder trabajar los datos de la BD en R.- Como si me conectara con R a un Servidor SQL y pudiera hacer consultas y crear mis propias tablas. Pero por lo que pude ver también en la red, la mejor

[R] Building a binary vector out of dichotomous variables

2016-06-16 Thread G . Maubach
Hi All, I need to build a binary vector made of a set of dichotomous variables. What I have so far is: -- cut -- ds_example <- structure( list( year2013 = c(0, 0, 0, 1, 1, 1, 1, 0), year2014 = c(0, 0, 1, 1, 0, 0, 1, 1), year2015 = c(0, 1, 1, 1, 0, 1,

Re: [R] Y in Kohonen xyf function

2016-06-16 Thread K. Elo
Hi! Some sample data could help us to help you... But have you read '?xyf' in order to ensure that your 'Y' is what 'xyf' expects it to be? What kind of error messages do you get? Regards, Kimmo 16.06.2016, 15:13, ch.elahe via R-help wrote: Is there any answer? Hi all, I have a df and I

[R] Y in Kohonen xyf function

2016-06-16 Thread ch.elahe via R-help
Is there any answer? Hi all, I have a df and I want to use supervised Self Organizing Map to do classification. I should use Kohonen library and xyf function from it. As you know the xyf function looks like this and I have problem defining my Y:

[R] extracting coefficients from ar() output

2016-06-16 Thread T.Riedle
Hi everybody, I am trying to run an AR1 model using the ar() function as shown below. > rollingarma<-rollapply(data,width=36,function(data) ar(data,aic=TRUE)) > head(rollingarma,50) order arvar.pred x.mean aicn.used order.max partialacf resid methodseries

[R] Mixed modelling course in Sydney

2016-06-16 Thread Highland Statistics Ltd
We would like to announce the following statistics course: Course: Introduction to Linear Mixed Effects Models and GLMM with R. Frequentist and Bayesian approaches. Where: UNSW, Sydney, Australia When: 18-22 July 2016 Course website: http://www.highstat.com/statscourse.htm Course flyer:

Re: [R] inverse table

2016-06-16 Thread Patrizio Frederic
Thank you all, David's solution is the one that matches my taste the most. Patrizio On Wed, Jun 15, 2016 at 8:04 PM, David L Carlson wrote: > After converting the table to a data frame, replicate each row by the > number of observations: > > > ddd.df <- as.data.frame(ddd) #

[R] Installing Caret

2016-06-16 Thread TJUN KIAT TEO
I am trying to install the package but am I keep getting this error messages installation of package ��minqa�� had non-zero exit status 2: In install.packages("caret", repos = "http://cran.stat.ucla.edu/;) : installation of package ��RcppEigen�� had non-zero exit status 3: In

Re: [R] help for fine mappting

2016-06-16 Thread PIKAL Petr
Hi From posted ref and map you cannot obtain final file need, they have nothing in common. answers see in line > -Original Message- > From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of greg holly > Sent: Wednesday, June 15, 2016 5:21 PM > To: r-help@r-project.org >