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)
conflm1<
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 s
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 li
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
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
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,
Wrong.
B.
On May 5, 2015, at 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?
>
> regards
>
> _
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 pa
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 wrote:
> Dear Daniel,
>
> Thank yo
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 UNSUBSC
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 decr
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 li
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 e
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 wro
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
> vari
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
On
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 after
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 im
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(lis
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 me
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.numeri
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
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 a
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
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:
>
> t<-c
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
my
On 04 May 2015, at 21:38 , T.Riedle 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
sum(psi*X[(t-1):(t-K)])
> but
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,
>
> Fi
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 acco
29 matches
Mail list logo
