Re: [R] analyzing results from Tuesday's US elections

This thread strikes me as pretty far off-topic for a forum dedicated to software support on R. https://www.r-project.org/mail.html#instructions "The ‘main’ R mailing list, for discussion about problems and solutions using R, announcements (not covered by ‘R-announce’ or ‘R-packages’, see above), a

Re: [R] Mac/PC differences in lmer results

To check whether the data are being read in appropriately, what happens when you plot the distribution of each of the independent variables on the respective systems? -A On Wed, 5 Jun 2019 12:32:28 +0200 Olivier Crouzet wrote: > Hi, > > 32bit vs. 64bit systems? > > Another thing I would look

[R] Help reinstalling rgdal (Ubuntu 16.04)

Hi all, Following some updates to R that I received via Synaptic Package Manager on Ubuntu 16.04 (looks like I now have R 3.4.4-1xenial0), I have been unable to reinstall rgdal, and I need help. Initially I was getting error messages about dependencies on GDAL 1.11.4, but after following instruct

[R] working with ordinal predictor variables?

I'm trying to develop a linear model for crop productivity based on variables published as part of the SSURGO database released by the USDA. My default is to just run lm() with continuous predictor variables as numeric, and discrete predictor variables as factors, but some of the discrete variable

[R] Dynamic forecast

Good afternoon. We have two typesof forecast using time series: statistical forecast and dynamic forecast. I would like to make dynamic forecastsimilar to what uses Eviews software. Does anyone know any package in R? Please. [[alternative HTML version deleted]] __

[R] Graphing a subset of data

alid 'envir' argument of type 'logical' Does anyone know why this is, or if there is something wrong with my code? Thanks! -- *Alexandra Hua * Yale University | MPH Candidate Class of 2016 Chronic Disease Epidemiology [[alternative HTML version deleted]] _

Re: [R] Finding values in a dataframe at a specified hour

igher than the 5*sigma value, I would think it would just pass to the next loop, yet it doesn't. I will keep trying! Thanks, Alexandra On Fri, Apr 10, 2015 at 3:43 PM, Jim Lemon wrote: > Hi Alexandra, > The error probably comes from the first iteration of i in 0:23. As indexing > in R

Re: [R] Finding values in a dataframe at a specified hour

the day 0:23. Thanks, Alexandra On Fri, Apr 10, 2015 at 1:07 PM, Alexandra Catena wrote: > Hello, > > I have a large dataframe (windHW) of wind speeds (ws) at each hour > from many days over a set of years. Some of these values are > obviously wrong (600 m/s) and I want to get ri

[R] Finding values in a dataframe at a specified hour

t for specific hours: windHW$ws[windHW$ws>=spring$sigma5] <- NA I imported the data using readLines and into the dataframe windHW. I also have R version 3.1.1 Any help would be appreciated! Thanks, Alexandra __ R-help@r-project.org mailing list -- To

[R] Plotting using tapply function output

d') I also tried to make a dataframe of the tapply output but it ends up using the hours as the column names instead of putting it into the dataframe. Please help!! I have R version 3.1.1 Thanks a lot, Alexandra __ R-help@r-project.org mailing l

[R] Error with using windRose function from the open air package

ingful for factors In addition: Warning messages: 1: In Ops.factor(mydata[[wd]], 10) : %% not meaningful for factors 2: In Ops.factor(mydata[[wd]], angle) : / not meaningful for factors Can anyone tell me what this means/what I'm doing wrong? Also, I have R

[R] Replacing 9999 and 999 values with NA

seem to do anything. > ltdl.fix.df(windMV, zero2na = FALSE, coded = 999) n = 9432 by p = 4 matrix checked, 0 NA(s) present 0 factor variable(s) present 5675 value(s) coded 999 set to NA 0 -ve value(s) set to +ve half the negative value I have R version 3.1.1 Than

[R] Help with looping

from the last year. Also, I tried Data[i] = readLines(sprintf('file/%4i,i)) but it says: "number of items to replace is not a multiple of replacement length" How do I get it to not replace each year of data? I have R version 2.15.1 Thanks, Alexandra [[alternative

[R] How to unzip a .gz file

Hello, Can someone help me with unzipping a .gz file. I used: readLines(gzfile('/home/file.gz')) I also found that I could use gunzip, but after trying to install it, it says: "package ‘gunzip’ is not available (for R version 2.15.1)" Thanks, Alexandra [[alte

Re: [R] How to download and unzip data in a loop

2005:2015? Thanks, Alexandra On Thu, Feb 5, 2015 at 3:11 AM, Jon Skoien wrote: > In addition to following Jim's suggestion, you should probably also use > full.names = TRUE, otherwise you will try to open a connection to files in > your current directory, not in tmpdir. > Another

[R] How to download and unzip data in a loop

quot;/tmp/RtmpKvB4Wz/724927-23285-2015.gz" next to 2015. and files returns 724927-23285-2015.gz. However, when I try to unzip the gz file using the last line, it says it cannot open the connection and the probable reason is that there is no such file or directory. T

[R] Gaussian mixture model based clustering

Hey everybody, I am looking for a R implementation of uni- and multivariate GMM clustering which is both fast and can provide me with estimates of the respective components' means and covariance matrices, not just the observations' labels for the components which I am not interested in. As it is me

[R] Particular case of a factorial experiment

Hi, Is there any package related to Analysis of a factorial experiment with more than one control treatment? if there is no package to do this, then, could you give me references of books about it? For help to create a program for it. Thanks. Fiorella Flores [[alternative HTML version d

Re: [R] bug in dummy.coef?

Thank you for the reply! I need to use it in the package for constructing the contrast matrices. Renaming a variable like y2 <- y^2 will not be a solution for me... Will try to look whether the function dummy.coef can be modified to work this issue. With best regards, Alexandra -Origi

Re: [R] bug in dummy.coef?

ts_2.11.3 grid_3.0.1 gtools_3.0.0 minqa_1.2.1nlme_3.1-109 [8] rpart_4.1-1tools_3.0.1 With best wishes, Alexandra From: r-help-boun...@r-project.org [r-help-boun...@r-project.org] on behalf of Alexandra Kuznetsova [a...@dtu.dk] Sent: Thurs

[R] bug in dummy.coef?

gdata_2.13.2 gplots_2.11.3 grid_3.0.1 gtools_3.0.0 minqa_1.2.1nlme_3.1-109 [8] rpart_4.1-1tools_3.0.1 With best wishes, Alexandra [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.

[R] Drawing asymmetric error bars

Hello, I have data which I have arcsin transformed to analyse. I want to plot my data with error bars however as my data is back-transformed my standard errors are uneven. Is there a simple way to draw these asymmetric error bars in R? Thanks for your help. [[alternative HTML version del

[R] fitted values with locfit

s this example? thanks a lot in advance I'm very desperate. Alexandra [[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-proje

[R] multivariate locfit regression

eights as this example? thanks a lot in advance. Alexandra [[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/post

[R] Why does the order of terms in a formula translate into different models/ model matrices?

umber of columns: ncol(X1) ## 9 > ncol(X2) ## 11 ## rank of design matrices: qr(X1)$rank ## 9 qr(X2)$rank ## 10 Will be very grateful if someone could help me here! Thanks, Alexandra -- View this message in context: http://r.789695.n4.nabble.com/Why-does-the-order-of-terms-in-a-f

[R] inverse of the weight function

Dear members, I'm trying to do a local polynomial estimation where the weight function is the inverse of the kernel (using a gaussian kernel). Does anybody know how I can do it? When I work with a weigth function usual I use: locfif(y~lp(x,deg=1,h=0.75),kern="gauss") Does anybody know

[R] varying coefficients model

Dear members, I'm trying to estimate a varying coefficients model using the local polynomial estimation method in two case (univariate and bivariate) and with two smooth functions. In the univariate case I use: m1<-smooth.lf(x=lp(z1,by=x1,deg=1,h=0.7,ev=z1)+lp(z3,by=expl2,deg=1,h=0.7,ev=z1

[R] kernel weight

his is not what I have to do. Please, can somebody help me? Thanks a lot. Alexandra [[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 ht

[R] multivariate bandwidth: regband

any other way to get a bivariate bandwidth in a local regression polynomial? I'll appreciate any query. I'm very lost and very desperate. I cannot find a solution. Thank you a lot in advance. Alexandra [[alternative HTML versi

[R] KernSmooth: dpill

t if I state "x" is a matrix with 2 columns the program also works. Please, somebody can help me? Thanks a lot Alexandra [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/list

[R] automathic bandwidth for multivariate local regression

to find a way to do it in the bivariate case. Please, help me. Thanks you in advance Alexandra [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the pos

[R] Bandwith selectors for multivariate local regression

Dear all, I'm trying to do a multivariate local regression whit the locfit instruction, such as:Y=m(Z,X,W)+u However, I have a problem at the moment of calculte the bandwith of the regression. If I had a univarite local regression model I could use the instruction "regband". However, t

[R] Generate 2 variables with normal distribution correlated betweenthem and with serial correlation

Hello, I want to create 2 variables with normal distribution that the correlation between them is equal to 0.8 and also each one has a serial correlation equal to 0.6. To generate variables with correlation between them I use: t1<-10 N<-45 sigma2<-matrix(c(1,0.8*sqrt(1),0.8*sqrt(1),1),

[R] Create two uniformly random variables correlated

w how to do it with variables uniformly distributed. Thanks a lot. Alexandra [[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://

[R] How to store the output of a loop into a matrix??

450x1). Thanks a lot for your help. Alexandra [[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

Re: [R] Ranking submodels by AIC (more general question)

rfect for making comparisons! Thanks again to everybody who's tried to help me out on this! Alexandra On Thu, 2011-06-23 at 21:29 +0200, Jan van der Laan wrote: > Alexandra, > > Have a look at add1 and drop1. > > Regards, > Jan > > > On 06/23/2011 07:32 PM, Ale

[R] Ranking submodels by AIC (more general question)

been able to find is stepAIC(), which as far as I can tell only gives back the best submodel, not a ranking of all submodels. Thanks in advance, Alexandra __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read t

Re: [R] [cleaned code] RE: AIC() vs. mle.aic() vs. step()?

On Thu, 2011-06-23 at 09:29 -0400, Alexandra Thorn wrote: > Ok, here's some example code showing how I get different output for AIC > vs. mle.aic(). Now that I've taken another look at the independent > variables, I'm wondering whether missing values in one of the vari

[R] [example code] RE: AIC() vs. mle.aic() vs. step()?

avior changes when I remove those... Alexandra #Code with outputs R> require(wle) R> xA # 1st independent variable (categorical) [1] Diffuse Diffuse Diffuse Diffuse Diffuse Diffuse Diffuse Diffuse Diffuse [10] Diffuse Diffuse RingDiffuse Diffuse RingDiffuse Diffuse Diffuse [19] Diffu

Re: [R] AIC() vs. mle.aic() vs. step()?

The packages is wle. I'll put together some code that shows the behavior I'm talking about, and send it to the list. Alexandra On Thu, 2011-06-23 at 13:51 +0200, Rubén Roa wrote: > I don't find the mle.aic function. Thus it does not ship with R and it's in > some c

[R] AIC() vs. mle.aic() vs. step()?

best models and also compare to far inferior models? Failing that, could someone point me to an appropriate resource that might help me understand? Thanks in advance, Alexandra __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r

[R] power of the comparison of two ROC curves

Hello, I would like to find the power of the comparison of two ROC curves (paired data). Is there an R function that does this? Thank you, Alexandra [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https

Re: [R] Robust regression error: Too many singular resamples

"You could try rlm in the MASS package; it doesn't use he resampling step." That seems to do the trick. Thank you! -- View this message in context: http://r.789695.n4.nabble.com/Robust-regression-error-Too-many-singular-resamples-tp2286585p2287468.html Sent from the R help mailing list archiv

[R] Robust regression error: Too many singular resamples

Hello. I've got a dataset that may have outliers in both x and y. While I am not at all familiar with robust regression, it looked like the function lmrob in package robustbase should handle this situation. When I try to use it, I get: Too many singular resamples Aborting fast_s_w_mem() Looki

[R] Modified R/S statistic

Hi R users, Does anyone know if Lo's modified R/S statistic is implemented in R? Thank you very much, Alexandra Almeida -- Alexandra R M de Almeida [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

[R] plotting fitted lme values as a smooth line

I am trying to plot fitted lme values as a smooth line of a graph showing the exponential relationship between temperature and soil respiration. In the plot, the x-axis has temperature, and the y-axis has soil respiration. When I try to add a line showing temperature versus the fitted va

[R] Help

Hi everybody! I´m with a problem that probably is easy for you but I really don´t know how to solve. On the following script: for(j in 1:length(limiares)) { excessos<-limiares[j]-estacao[estacao__ R-help@r-project.org mailing list https://stat.ethz.c

[R] Help

Hi everybody! I´m with a problem that probably is easy for you, but I really don´t know how to solve. On the following script: for(j in 1:length(limiares)) { excessos<-limiares[j]-estacao[estacao__ R-help@r-project.org mailing list https://stat.ethz.

[R] IDE

People, I'm a ubunto user and I used to write my scipts in "Java Gui for R", but it is a very slow tool to run my scripts... Do you know some efficient IDE for R? Thankssss!!! Alexandra Almeida -- Alexandra R M de Almeida [[alternative HTML

[R] Goodness of fit tests

. Can anyone help me? Alexandra [[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 comm

[R] test for the variance of a normal population

Hi, Does anyone know an instruction to perform a test (or to construct a confidence interval) for the variance of a normal population (distribution). The test statistic is (n-1)*S´^2/sigma^2 and has chi square (n-1) distribution. Thank you, Alexandra [[alternative HTML

[R] analysing mixed effects/poisson/correlated data

I am attempting to model data with the following variables: timepoint - n=48, monthly over 4 years hospital - n=3 opsn1 - no of outcomes total.patients skillmixpc - skill mix percentage nurse.hours.per.day Aims To determine if skillmix affects rate (i.e. no.of.outcomes/total.patien