[R] Number of used observations in stats::arima with NAs

g" is not NULL, "dx" is required in other part of the code, so it is convenient to use it here as well; replacing "x" by "dx" and "xreg" by "dxreg": isna <- is.na(dx) | apply(dxreg, 1L, anyNA) n.used <- sum(!isna) Are there any reaso

Re: [R] Extending sparklyr

For versions 1.6.1 and 2.0.0 of Spark, the GaussianMixture is under the ml namespace not mllib, try this instead: envir$model <- "org.apache.spark.mllib.clustering.GaussianMixture" Best, Javier On Sun, Oct 9, 2016 at 1:47 PM, Axel Urbiz wrote: > Hi All, > > Just star

[R] metaprop

hi, i am using the funtion metaprop but when i try to save the proportion estimation by subgrup using the option "byvar", the value in TE.fixed.w are the transformed proportion, and i need the untranformed proportion. I try with the option backtranf=true, but doesn`t change anithing. I have so m

[R] R in UNIX OS

Hello, I have a silly questions: is possilbe to install R in a AIX OS Machine? Where can I found the download files? Thank you in advance Javier -- [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE

Re: [R] Heatmap

Please read R for Mac OS X FAQ (see Help) section 9 and adjust your system preferences accordingly. I looking the R for Mac document but there is not section 9.!!! Thanks for any help. Carlos. 2015-05-05 18:01 GMT-05:00 Carlos Javier Rincon Rodriguez < cjrinc...@unal.edu.co>: >

[R] Image processing using R

and I would like to complete the project, if possible, in R. Thanks in advance for your time. Javier [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo

[R] Function StructTS: covariance matrix of the initial state vector

Hi everyone: My question is about the function "StructTS" of the core package "stats", which fits by maximum likelihood the basic structural time series model. According to theory and the references given in "?StructTS", the covariance matrix of the initial state vector is a diagonal matrix. Ho

[R] Configuring Rserve for remote access on Windows Server

I am trying to configure Rserve for use with Tableau on a Windows Server 2012, however no remote telnet commands are successful on port 6311. Is there any link or guide you can recommend on configuring this package on a windows environment? Thanks! Javier Guillen [[alternative HTML

[R] Compilation fails for package 'BayesTree' on R 3.0.1 for Windows

ain this R project requiring BayesTree package. I am new to R and the error message I am getting does not provide much insight on how to solve the problem, so I have no clue what to do. Thank you in advance for any help, Javier [[alternative HTML version deleted]] _

Re: [R] fitting grid-based models

Hi Thomas, Thank you very much for your help. I will try to calculate some summary statistics and fit observed and simulated data using a time series (as your examples). With observed data, I would use two values (e.g. 0 and 10), because I do not have intermediate values. Thanks! Javier 2014-06

[R] fitting grid-based models

autocorrelation or sum of values in the landscape) and try to fit this statistics with others calculated from observed data. Thank you very much for your help! Javier [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

[R] KappaSize

Hello, good morning. I need to calculate a kappa coeficient between three raters and the scale have five categories outcome. Until now I read the command CI3Cats for the package KappaSize in order to calculate the sample size, but when i save the package into the library, R can read it and i dont

[R] Linear combination of time series for approximating another one

Dear R community, I would like to approximate a time serie as linear combination of a set of many time series, minimizing the number of the time series involved in the linear combination. Can you recommend to me any method for this? Any paper about this? I will appreciate much your suggestions.

[R] Linear combination of time series for approximating another one.

Dear R community, I would like to approximate a time serie as linear combination of a set of many time series, minimizing the number of the time series involved in the linear combination. Can you recommend to me any method for this? Any paper about this? I will appreciate much your suggestions.

[R] fitting grid-based models

e but all examples are for differential eq. models. It is mentioned in the that such functions could also be used for other kind of models. Could thus be used for grid-based ones? If that, could be used even the model does not follow the simecol syntax? Thanks for your help! J

[R] Help

ts? Thank you very much. Jonson M. Javier  [[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 a

[R] How to treat dominant alternatives & infeasible/unrealistic alternatives in the choice set

help. Jonson M. Javier [[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, mi

[R] RWeka - weighted instances

Is there any way to train Weka classifiers using weighted instances in RWeka? I have done it in Weka in Java, but I don't see any way to pass the weights to the classifier in the RWeka documentation. Thanks __ R-help@r-project.org mailing list https://s

[R] scale_y_logit not present in ggplot2 0.9?

I have updated from ggplot2 0.8.9 on Windows to ggplot 0.9.0 (and then 0.9.1) and I can't find scale_y_logit() anymore ... On Mac, I can't see scale_y_logit() in ggplot 0.9.0 either. Am I missing anything? Thanks __ R-help@r-project.org mailing list ht

[R] Spectral Angle Mapper, Spectral Information Divergence

I am looking for R packages that can calculate the Spectral Angle Mapper (SAM) or the Spectral Information Divergence (SID) between two reflectance spectrums. Thank you, -- Dr Javier Bustamante Estación Biológica de Doñana, CSIC Dept. Wetland Ecology Américo Vespucio s/n 41092 - Sevilla Spain

Re: [R] differences between survival models between STATA and R

n using two types of survival programs to analyse a data set. > > The first one is an R function called aftreg. The second one an STATA > function called streg. > > Both of them include the same analyisis with a weibull distribution. Yet, > results are very different. > > Sho

Re: [R] Estimation of the mode

arameters (bandwidth, etc.) at a minimum. > > So perhaps you need to rethink what you are trying to do and/or > perhaps get some advice from a local statistical resource. > > -- Bert > > On Fri, Jan 20, 2012 at 8:27 AM, Javier xyz wrote: > > > > Hi all, > >

[R] Estimation of the mode

Hi all, I am trying to estimate the mode of a 4-dimensional nonparametric density estimator (any) using a sample of size n=10,000. I have tried using the package 'ks' and 'np' but they are extremely slow; this is related to the estimation of the bandwidth matrix. I also checked the package

Re: [R] any DCCA function in R?

after this I still need to run the DCCA I will go to CANOCO. I can access to this software too, but as I have began to use R some months ago I would like to do everything in R, but perhaps this analysis can be an exception. Thanks. Sincerely, Javier Murillo -- View this message in context

Re: [R] any DCCA function in R?

this. Could be a good idea to post that message on the R-SIG-Ecology list too. Again thank you very much, Javier -Mensaje original- De: Gavin Simpson [mailto:gavin.simp...@ucl.ac.uk] Enviado el: martes, 20 de diciembre de 2011 9:58 Para: Juan José Santos Blanco CC: r-help@r-project.org

Re: [R] any DCCA function in R?

this. Could be a good idea to post that message on the R-SIG-Ecology list too. Again thank you very much, Javier -Mensaje original- De: Juan José Santos Blanco Enviado el: martes, 20 de diciembre de 2011 9:32 Para: r-help@r-project.org CC: Javier Murillo Pérez Asunto: any DCCA function

[R] help about fitdistr funtion

Hi, I have a variable X classified in a lot of groups and I need to run the [fitdistr] funtion for each group. I tried with the [by] or the [tapply] funtions because my data is organize in two columns (variable and the groups), but neither of these command work. If somebody have a tip to help me

Re: [R] COXPH TIME-DEPENDENT

Thanks, Javier -- View this message in context: http://r.789695.n4.nabble.com/COXPH-TIME-DEPENDENT-tp3754837p3756128.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

[R] AFT model time-dependent with weibull distribution

2- How can I solve this error? 3- Why the phreg function does not need an ID? Thanks, Javier -- View this message in context: http://r.789695.n4.nabble.com/AFT-model-time-dependent-with-weibull-distribution-tp3755079p3755079.html Sent from the R help maili

[R] Procesamiento paralelo

Hola, soy nuevo en la lista y en R, disculpen esta pregunta, pero no he encontrado información actualizada. Mi pregunta, es R capaz de realizar procesamiento paralelo? -- Atentamente Ing. Javier Alfonso Valdés Universidad de las Ciencias Informáticas

[R] Extract the names of the arguments in an "expression"

Hi everybody: I need to get the names of the arguments in an object of class "expression". I've got the following expression: > x <- expression(rho * cos(omega)) Is there any function or procedure that returns the names of the arguments (in the example: "rho" and "omega")? I tried a rough appro

[R] Problem with lm4 and lapack library

1ubuntu3 --- R CMD config BLAS_LIBS -lblas R CMD config LAPACK_LIBS -L/usr/lib64/R/lib -lRlapack R seems to be looking for LAPACK in "/usr/lib64/R/lib -lRlapack", probably the wrong place. Can somebody point me to a possible solution? -- Dr Javier Bustamant

[R] Efficient nested loops

Dear R community, I am working with huge arrays, so I spend a lot of time computing. This is my code: for (x in 1:dim(variable)[1]){ for (y in 1:dim(variable)[2]){ for (z in 1:dim(variable)[3]){ result <- max(variable[x,y,z,]) } } } Is there a more efficient p

[R] How to call to R_KalmanLike from outside StructTS

Dear all: I want to modify the 'StructTS' function from the 'stats' package. First, I am writing a working copy of the original version and got some problems. I have two versions of the function plus the original one. The first version is the same code as the 'StructTS' function: StructTS.v1 <

[R] How can I get the principal components after a varimax rotation using Varimax function?

Dear R users, I have some question about Varimax rotation of the loadings obtained from a PCA. Imagine X is a field where the rows are the observations and the columns are the variables. I obtain the loadings (L), the principal components (PC) and the percentage of the variance for each mode (PVa

Re: [R] How to plot the sample size on top of the bars of a barchart plot?

Dear Marco, I was trying using lattice barchart() but your suggestion will do the trick. Many thanks Javier Dr. Javier Perez-Barberia Population and evolutionary ecologist The Macaulay Institute Craigiebuckler Aberdeen AB15 8QH Scotland, UK Tel: +44(0)1224 395221 e-mail: j.perez-barbe

[R] How to plot the sample size on top of the bars of a barchart plot?

Dear users, Is there anyone how knows how to plot the sample size on top of the bars of a barchart plot? Many thanks in advance Javier -- View this message in context: http://r.789695.n4.nabble.com/How-to-plot-the-sample-size-on-top-of-the-bars-of-a-barchart-plot-tp2223062p2223062.html

[R] format: from list to data frame

67329" "367309" "367339" ... ..$ y: chr [1:173] "780175" "780185" "780195" "780205" ... > length(yo) [1] 44807 where x is the name of the rows and y the name of the columns (coordinates) of a matrix of 44807 values ho

[R] Error analysis for circular data

Dear R users, I would like to know if it is possible to calculate the Mean Error (ME), the Root Mean-squared error (RMSE) and absolute error (MAE) for two temporal series of directional data. Where Can I get documentation about it? Thanks in advance [[alternative HTML version deleted]

[R] how to plot single frames as a movie?

the playing speed? Sorry if this questions was already in the forum, but I couldn’t find it. Thanks in advance Javier -- View this message in context: http://n4.nabble.com/how-to-plot-single-frames-as-a-movie-tp1468628p1468628.html Sent from the R help mailing list archive at Nabble.com

[R] Doubt about cluster analysis

Dear R community, I'm a beginner with Cluster Analysis. I would like to know if there is a criterion to select the best set of clusters to do this analysis. Thanks in advance, [[alternative HTML version deleted]] __ R-help@r-project.org mailin

[R] Problem with Jri (rJava)

library that is available with the rJava distribution, but I get the same error. I work on Mac OS X with java 1.6, the java.library.path is well configured. Thank you very much! Javier This is the error: java.lang.UnsatisfiedLinkError: /../rJava/jri/libjri.jnilib: no suitable image found. Did

[R] Problem with Jri (rJava)

library that is available with the rJava distribution, but I get the same error. I work on Mac OS X with java 1.6, the java.library.path is well configured. Thank you very much! Javier This is the error: java.lang.UnsatisfiedLinkError: /../rJava/jri/libjri.jnilib: no suitable image found. Did

[R] Doubt about CCA and PCA

Dear R community, I'm working with PCA and CCA methods, and I have a theoretical question. Why is it necesary to have more temporal values than variables when the CCA O PCA are going to be used? Could you advise to me some any paper about it? Thanks in advance, [[alternative HTML versi

[R] Spatstat: xy binary data into mask type to use in owin(mask=)

x[i], y[j] (1 if the area is TRUE and 0 is the area is FALSE). How can I covert this xy file into a mask type "my_area" that I can read using owin(mask=my_area)? Many thanks in advance Javier -- View this message in context: http://www.nabble.com/Spatstat%3A-xy-binary-data-into

[R] How can I run a function to a piece of text?

ep() steps mean(script) ##function that I want to apply: doesn't work Is there any way to convert the "script" class so that the function mean() can read it as if it were text typed in the console? Thanks and have a superb day Javier -

Re: [R] Output an RWeka model via sink

Thanks. Just figured out what happened: out of quota. -Original Message- From: Coen van Hasselt Sent: Saturday, March 28, 2009 9:05 PM To: Francisco Javier Perez Caballero Cc: r-h...@stat.math.ethz.ch Subject: Re: [R] Output an RWeka model via sink I tried your code and it seems to

[R] Output an RWeka model via sink

When I sink the output of an RWeka model to a text file, the output file appears empty: library(RWeka) model = LogitBoost(Species~.,data=iris) print(model) sink("output.txt") print(model) #file output.txt is created, but it is blank sink() Am I doing anything wrong? [[alternative HTML v

[R] RODBC crashes connecting to Teradata

Hi, I'm trying to connect to a Teradata database via RODBC on a Linux 64 machine (Red Hat Enterprise Linux 5). The ODBC driver is properly configured and queries sent via unixODBC's isql tool work properly. However, this is what happens when I try to connect via RODBC: > library(RODBC) > conn = o

[R] Plotting regression lines and points from lm using lattice

del as in complex models I get confused when I have to combine the coefficients to build up each regression line. Many thanks in advance and sorry for bothering you. Javier fm<-lm(c4~SEX+AREA+c3+SEX:AREA+SEX:c3+AREA:c3,data=mydata) mygrid<-expand.grid(SEX=c("male","femal

Re: [R] Unexpected behaviour when testing for independence with multiple factors

Michael, so you're suggesting that I should do: aux <- interaction( Topology, Drift, lambda) and then fligner.test(dT~aux) Is that correct? On Thu, Sep 18, 2008 at 8:32 AM, Michael Dewey aghmed.fsnet.co.uk> wrote: > At 16:03 17/09/2008, Javier Acuña wrote: >> >> H

Re: [R] Unexpected behaviour when testing for independence, with multiple factors

>Ben Bolker ufl.edu> writes: > >I would try > >fligner.test(dT ~ Topology:Drift:lambda) > >in response to: > >Javier Acuna gmail.com> writes: > > Hi, I'm a new user of R. My background is Electrical Engineering, so > please bear with me if thi

[R] Unexpected behaviour when testing for independence with multiple factors

if I'm using the test incorrectly or something else. I would really appreciate any help. I'm using R version 2.7.2 (2008-08-25) in Windows XP. Many thanks in advance Javier ---- Javier Acuna Electrical Engineering Grad Student Univer

[R] Help with eigenvectors

Hi everybody, I have some problems with the function eigen. I have a square matrix and I want to calculate the eigenvalues and eigenvectors. I apply the function eigen and I get it, however when I solve the same problem in Statistica software, I realise that some eigenvectors are the opposite. How

[R] reshape matrices

Hello everyone, I need reshape an array. For example, if we have next array: > a <- c(1,2,3,4,5,6,7,8,9,10,11,12) > dim(a) <- c(2,2,3) > a , , 1 [,1] [,2] [1,]13 [2,]24 , , 2 [,1] [,2] [1,]57 [2,]68 , , 3 [,1] [,2] [1,]9 11 [2,] 10 12

[R] reshape matrices

Hello everyone, I need reshape an array. For example, if we have next array: > a <- c(1,2,3,4,5,6,7,8,9,10,11,12) > dim(a) <- c(2,2,3) > a , , 1 [,1] [,2] [1,]13 [2,]24 , , 2 [,1] [,2] [1,]57 [2,]68 , , 3 [,1] [,2] [1,]9 11 [2,] 10 12

[R] I need help with eofs

that I want. I would like to get the EOFS and to export it in a txt file. I would be pleased if you could help me. Thank you, and sorry for my english, Regards Francisco Javier [[alternative HTML version deleted]] __ R-help@r-project.or

[R] I need help with eofs

that I want. I have got the PCs series, but I need obtain the spatial patterns (EOFS) and to export it in a txt file. I would be pleased if you could help me. Thank you, and sorry for my english, Regards Francisco Javier [[alternative HTML version deleted]]

Re: [R] bug in file.path?

> of reading it? Or attach it as a plain text file and tell us what encoding > the file is in? > > Yes, there probably is a bug here but we do need to be able to reproduce it. > > On Thu, 24 Apr 2008, Javier Muñoz wrote: > > > > Hello, > > > > T

[R] bug in file.path?

.1252 attached base packages: [1] stats graphics grDevices datasets utils methods base other attached packages: [1] RODBC_1.2-3 ---- Is this a bug in the new file.path() or can I change something to make this work as before? Thanks. Javier. [[alternative HTML