Re: [R-es] Muestreo de bases de datos

2015-05-05 Thread Fernando Fernández
Hola, Para mí el principal error de la anécdota es intentar hacer una traducción literal de lo que se tiene programado en SAS porque no se tiene ni idea de cómo funciona R. Muchas veces en SAS se hacen determinadas cosas simplemente porque se puede sin pensar en si realmente son necesarias o son

Re: [R] A problem with string handling to make a time duration

2015-05-05 Thread gavinr
Thanks guys. The first solution with the gsub / lapply works perfectly. The solution using substrings would work if the times were in a consistent format, but without the leading zeros and with some parts of the string absent completely it would need some extra logic to apply. I need something to

Re: [R] computing AIC when distribution is fitted by lmom or lmomco

2015-05-05 Thread Vito M. R. Muggeo
dear Milan, I think you should consult a (local) statistician for your analyses. R (and packages) are only software, and you need theoretical background. best, vito Il 05/05/2015 14.25, Milan Cisty ha scritto: Dear list members, Is it possible to compute in R AIC, when model was fitted by

[R] computing AIC when distribution is fitted by lmom or lmomco

2015-05-05 Thread Milan Cisty
Dear list members, Is it possible to compute in R AIC, when model was fitted by methods of L-moments (packages lmom or lmomco)? Or what could be used to compare different candidate distribution models fitting to data? Sorry I have no reproducible example, because I do not know how this (AIC) could

Re: [R] help - hoslem.test

2015-05-05 Thread Luciane Maria Pilotto
Hi Thanks for the help! I used the following commands to run the test hoslem: rl-id3[,c(q06, q08, q10, q11, q12, q13, q14, sexo, edcat, rendacat, cpod, getario, bwgr_et)]#select only variables of interest dim(rl)#[1] 2284312 attach(rl) q131 - as.numeric(cut(rl$q13,c(0,1.5,5)))-1

Re: [R] Drawing the regression line and the 95% confidence intervals

2015-05-05 Thread Jim Lemon
Hi Sacha, The line you have requested is off the plot. The following will produce what I think you are asking, but I cannot speak for whether it means anything sensible. plot(GDP.per.head, fitted(LinearModel.1),ylim=c(200,800)) devlm1-lm(fitted(LinearModel.1)~GDP.per.head) abline(devlm1)

[R] [R-pkgs] New release of MAVIS (Meta Analysis via Shiny) v1.1 now on CRAN

2015-05-05 Thread Kyle Hamilton
Dear R users, I'm pleased to announce that the first major update to the MAVIS package (Meta Analysis via Shiny) has been released to CRAN. MAVIS: Meta Analysis via Shiny v1.1 Smiling Fox has expanded greatly over the last couple of months to include the following. - User interface has been

Re: [R] Plot by FIPS Code using Shapefiles

2015-05-05 Thread Shouro Dasgupta
Hello, Thank you for your reply. My original data has 3109 FIPS codes. Is there a way to merge only this data into the shapefiles? I hope I am clear. Thank you also for the link, I am trying to do something like this: https://gist.github.com/reubano/1281134. Thanks again! Sincerely, Shouro

[R] Plot by FIPS Code using Shapefiles

2015-05-05 Thread Shouro Dasgupta
I am trying to plot data by FIPS code using county shapes files. library(data.table) library(rgdal) library(colourschemes) library(RColorBrewer) library(maptools) library(maps) library(ggmap) I have data by FIPS code which looks like this: dput(head(max_change)) structure(list(FIPS =

Re: [R] Plot by FIPS Code using Shapefiles

2015-05-05 Thread Anthony Damico
hi, after running each individual line of code above, check that the object still has the expected number of records and unique county fips codes. it looks like length( shapes$GEOID ) == 3233 but nrow( merged_data ) == 3109. the way for you to debug this is for you to go through line by line

Re: [R] Transformation of the circular variable into linear variable

2015-05-05 Thread David Winsemius
On May 5, 2015, at 3:59 AM, Nasr Al-Dhurafi wrote: Dear All I am working with a circular data ( wind direction ) and i have problem in transforming the circular variable into linear variable. I have found an equation that can be used to convert the circular variable into linear variable

[R] Drawing the regression line and the 95% confidence intervals

2015-05-05 Thread varin sacha
Hi, Dear R-helpers, Here below you will find a reproducible fictitious example working except the abline function. First thing : I try to draw the regression line (multiple linear regression). I try the abline function but it does not work. I don't get any error message but the straight line

Re: [R] order by which the eigenvalues are presented

2015-05-05 Thread William Dunlap
You can sort the eigenvalues in the order you want with o - order(Re(e$values), decreasing = TRUE) or o - order(abs(Re(e$values)), decreasing = TRUE) followed by e$values[o] where 'e' is the object that eigen returns. The main argument to order() is what you want to sort by and the

[R] predict function in regression analysis

2015-05-05 Thread li li
Hi all, I have the following data in which there is one factor lot with six levels and one continuous convariate time. I want to fit an Ancova model with common slope and different intercept. So the six lots will have seperate paralell regression lines.I wanted to find the upper 95% confidence

Re: [R] Plot by FIPS Code using Shapefiles

2015-05-05 Thread Anthony Damico
so check the unique number of fips codes in the objects before and after merged_data - merge(shapes_fips,max_change,by=FIPS,all.X=T, all.y=T) also note that all.X should be all.x and you might want to use FALSE for one or both of those On Tue, May 5, 2015 at 11:40 AM, Shouro Dasgupta

Re: [R-es] Muestreo de bases de datos

2015-05-05 Thread Patricio Fuenmayor Viteri
Bueno ... voy a compartirles mi experiencia en esto.Freddy, yo tuve la misma experiencia que su amigo, es mas, aun la sigo teniendo... y en resumen y sin alargar el cuento ... me quedo con R.Porque mi decisi�n... por varios temas:1. Rapidez. Uno debe saber que es lo que se quiere hacer y que

[R] update.packages() ask for mirror several times

2015-05-05 Thread Nordlund, Dan (DSHS/RDA)
I recently installed an R-3.2.0 binary for Win7 64-bit from a CRAN table. Now when using the menu to update packages, it asks me to select a mirror several times before proceeding to successfully download and install the updated packages. It is only a minor annoyance, was wondering if anyone

Re: [R] update.packages() ask for mirror several times

2015-05-05 Thread Duncan Murdoch
On 05/05/2015 12:13 PM, Nordlund, Dan (DSHS/RDA) wrote: I recently installed an R-3.2.0 binary for Win7 64-bit from a CRAN table. Now when using the menu to update packages, it asks me to select a mirror several times before proceeding to successfully download and install the updated

Re: [R] does segfault mean (always) a bug?

2015-05-05 Thread Boris Steipe
Wrong. B. On May 5, 2015, at 2:54 PM, lejeczek pelj...@yahoo.co.uk wrote: hi eveybody I'm trying something simple (Biocunductor packages), so simple I believe it's example from docs but I get segfault. I don't suppose incorrect scripting can cause segfault, right? regards

Re: [R] update.packages() ask for mirror several times

2015-05-05 Thread Nordlund, Dan (DSHS/RDA)
Thanks, I will do that. Dan Daniel J. Nordlund, PhD Research and Data Analysis Division Services Enterprise Support Administration Washington State Department of Social and Health Services -Original Message- From: Duncan Murdoch [mailto:murdoch.dun...@gmail.com] Sent: Tuesday, May

[R] does segfault mean (always) a bug?

2015-05-05 Thread lejeczek
hi eveybody I'm trying something simple (Biocunductor packages), so simple I believe it's example from docs but I get segfault. I don't suppose incorrect scripting can cause segfault, right? regards __ R-help@r-project.org mailing list -- To

Re: [R] Need Help!

2015-05-05 Thread Boris Steipe
For crying out loud Dawood - can you please stop posting screenshots? Do you not know how to copy your output and paste it into the body of a (not HTML formatted!) email? Let us know if you need help with that. B. On May 5, 2015, at 11:05 AM, dawood ahmad dawoodma...@hotmail.com wrote:

Re: [R] does segfault mean (always) a bug?

2015-05-05 Thread Duncan Murdoch
On 05/05/2015 2:54 PM, lejeczek wrote: hi eveybody I'm trying something simple (Biocunductor packages), so simple I believe it's example from docs but I get segfault. I don't suppose incorrect scripting can cause segfault, right? In R, a segfault always indicates a bug. What's not so

Re: [R] Problem in r help me

2015-05-05 Thread PIKAL Petr
From: Ghada Almousa [mailto:ghada.f...@gmail.com] Sent: Monday, May 04, 2015 9:20 PM To: PIKAL Petr Subject: Re: Problem in r help me I used internet search to get info about methods but I can't find any info about the comparison between those methods The comparison in the number of iterations

[R] Sort data and symbol change in Jitter plot (ggplot2)

2015-05-05 Thread Luis Fernando García
Dear R experts, First than all I want to thank your expertise and for sharing your knowledge with the people who are starting. Recently I have been working on a new plot style (Jitterplot) but have still some issues when making two basic functions. First than all I want to change the symbol

Re: [R] Sort data and symbol change in Jitter plot (ggplot2)

2015-05-05 Thread PIKAL Petr
Hi answeres in line -Original Message- From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Luis Fernando García Sent: Tuesday, May 05, 2015 9:34 AM To: r-help@r-project.org Subject: [R] Sort data and symbol change in Jitter plot (ggplot2) Dear R experts, First than

Re: [R] predict function in regression analysis

2015-05-05 Thread Bert Gunter
I think you want CI's for intercepts, not means (what is a mean for a line??). If so, the ?confint function will give this to you for the lot effect estimates when applied to a model fitted without an intercept: myfit -lm(y~ lot-1+time) confint(myfit) Further discussion should be directed to a

[R] order by which the eigenvalues are presented

2015-05-05 Thread Luis Borda de Agua
I’m using R 3.1.3 in OS X 10.10.3 (Yosemite) I use the function “eigen” to calculate the eigenvalues of matrix where each element is sampled from a given distribution (normal, beta, etc.). According to the information provided: values a vector containing the p eigenvalues of x, sorted in

Re: [R] using filter() to sum up

2015-05-05 Thread peter dalgaard
On 04 May 2015, at 21:38 , T.Riedle tr...@kent.ac.uk wrote: Hi everybody, I am trying to create a code for the formula in the attachment. I first tried following code: ltau - m + theta*sum(psi*X[t-k]) That's not going to work. It might work with something like

[R] append to all warnings

2015-05-05 Thread Rory Wilson
Hello all,This may be a newbie type question but I am rather stuck. I would like to alter my R options such that anytime a warning is generated, the standard warning message is generated, as well as additional details that I would specify (ie, time/date, etc.) It seems to me something like 

Re: [R] A problem with string handling to make a time duration

2015-05-05 Thread Franklin Bretschneider
Hello gavinr, I have a character string that represents a time duration. It has an hours minutes seconds structure(ish) but with letters denoting units (H,M or S) no leading zeros and no placeholder at all where one or other of the units are not required. It looks like this:

Re: [R] Transformation of the circular variable into linear variable

2015-05-05 Thread Nasr Al-Dhurafi
Dear All I am working with a circular data ( wind direction ) and i have problem in transforming the circular variable into linear variable. I have found an equation that can be used to convert the circular variable into linear variable which is included in this paper * ARMA based approaches for

Re: [R-es] Muestreo de bases de datos

2015-05-05 Thread Gilsanz, Jose Luis
Hola Desviándome un poco de la cuestión que plantea Freddy, no estoy muy de acuerdo con lo que comenta acerca de que ... R no está diseñado (operaciones en la base de datos)... Precisamente la mayor parte de mi trabajo en R se realiza vinculadolo a un servidor Microsoft SQL Server y lo