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
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
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
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
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
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)
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
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
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 =
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
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
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
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
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
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
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
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
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
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
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
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
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:
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
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
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
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
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
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
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
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
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:
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
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
33 matches
Mail list logo