Charles Berry ucsd.edu> writes:
>
> Scherber, Christoph gwdg.de> writes:
>
> >
> > Dear all,
> >
> > I am trying to express a multinomial GLM (using nnet) as a series of GLM
> models.
[deleted]
>
> Doing the obvious comparison:
>
> ggen.preds <-
> sapply( levels(multicats),
>
On Jul 21, 2014, at 4:53 AM, Prof Brian Ripley wrote:
> On 21/07/2014 12:24, Charles Thuo wrote:
>> I have just installed R 3.1.1 in a machine where R 3.0.1 is already
>> installed. Is it possible to use packages in the 3.0.1 on the 3.1.1.
>> version as the same are in a single workstation.
>
In R, I have a raster entitled "raster_crude" with the following details:
> raster_crude
class : RasterLayer dimensions : 320, 392, 125440 (nrow, ncol, ncell)
resolution : 0.125, 0.125 (x, y) extent : -152, -103, 30, 70 (xmin, xmax,
ymin, ymax) coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS8
Hey Shane,
Sorry you're having trouble.
The quick start is here and walks through installation: https://plot.ly/r/.
A note. If you're on Windows, you'll need Rtools to install devtools:
http://cran.rstudio.com/bin/windows/Rtools/. As Sarah noted, Plotly isn't
on CRAN.
If you're having trouble,
Dear Christoph,
If I understand correctly what you've done, the two approaches are not
equivalent and should not in general produce the same fitted probabilities.
Letting {a, b} represent logit(a vs. b) = log(Pr(a)/Pr(b)) and {ab, cd}
represent logit(a or b vs. c or d), and numbering the respon
Thanks, that was exactly it -- switching the values did the trick (and
was actually correct in terms of theory.) And of course, you are right
-- i changed the starting values to mean(x) - mean(y) for mu and
sqrt(var(x-y)) for sigma.
I also see your point about the theoretical justification for the
Here are 2 possibilities to consider.
The shiny package allows a web browser interface to R. It can run
from a server, but can also run just on a single computer. You could
set this up like you suggest, where the user would double click on an
icon on the desktop which would then run R, load the
Thanks Bert. I will post it on StatsExchange.
David
2014-07-23 2:08 GMT+08:00 Bert Gunter :
> Wrong list!
>
> Try stats.stackexchange.com for statistics questions. This list is
> about R programming related issues.
>
> Also, note that HTML does not work and should not be used here.
>
> Cheers,
Wrong list!
Try stats.stackexchange.com for statistics questions. This list is
about R programming related issues.
Also, note that HTML does not work and should not be used here.
Cheers,
Bert
Bert Gunter
Genentech Nonclinical Biostatistics
(650) 467-7374
"Data is not information. Information
Hi all r-mailling listers:
Can anyone explain the theory (or the formula) about computing Sum Sq
(color highligh below) related to regression items? The link of Wikipedia (
http://en.wikipedia.org/wiki/Partition_of_sums_of_squares) gives an
introduction on how to calculate the total, model, and r
Hi Chris,
Thanks for your reply. Certainly no need to apologize for a delayed response.
Appreciate your taking the time to answer my question.
My concern was about the value of "16". My understanding is it is based on one
expert's opinion. I wondered if this is in keeping with whatever the sta
Dear Lynn,
On Tue, 22 Jul 2014 13:59:39 +0200
Lynn Govaert wrote:
> Hi,
>
> I have a problem with doing a repeated measures ANOVA. I will first give
> you an idea of what my dataset looks like. We have 20 ponds, and for each
> pond we took some individuals (waterfleas), say 10 and tested them a
Scherber, Christoph gwdg.de> writes:
>
> Dear all,
>
> I am trying to express a multinomial GLM (using nnet) as a series of GLM
models.
>
> However, when I compare the multinom() predictions to those from GLM, I
see differences that I can´t
> explain. Can anyone help me out here?
>
> Here com
Yes, but note that such comparisons don't necessarily tell you much.
Using your l1 on my computer:
> system.time(r1 <- sapply(l1, `[`, 1))
user system elapsed
1.2 0.0 1.2
> system.time(r2<- rapply(l1,function(x)x[1]))
user system elapsed
0.810.000.81
## But
> syste
Or
rapply(l,function(x) x[1])
#[1] 1 3 7
set.seed(42)
l1 <- replicate(1e6, list(sample(1:5,sample(8),replace=T)))
system.time(r1 <- sapply(l1, `[`, 1))
# user system elapsed
# 1.324 0.000 1.326
system.time(r2 <- rapply(l1, function(x) x[1]))
# user system elapsed
# 0.736 0.004
Dear all,
I have a problem with performing statistical analyses of longitudinal data
after the imputation of missing values using mice. After the imputation of
missings in the wide data-format I convert the extracted data to the
longformat. Because of the longitudinal data participants have dup
Hi,
I have a problem with doing a repeated measures ANOVA. I will first give
you an idea of what my dataset looks like. We have 20 ponds, and for each
pond we took some individuals (waterfleas), say 10 and tested them along
two treatments (A and B). Now, for each pond we also know the fish
backgro
On Tue, 22 Jul 2014, John McKown wrote:
I like to keep the individual lines in my source files relatively
short. Mainly so that I can print them, email them, or display them on
a narrow screen without needing to shift left & right. So, for a
really long character string, such as an SQL query, I
Dear mahe,
This is, I assume, the scatterplotMatrix() function in the car package.
I don't think that you'll be able to do what you want with scatterplotMatrix(),
because it uses the pairs() function, which resets par("mfrow"). For the same
reason, using pairs() directly won't help. You should,
Dear all,
I am trying to express a multinomial GLM (using nnet) as a series of GLM models.
However, when I compare the multinom() predictions to those from GLM, I see
differences that I can´t
explain. Can anyone help me out here?
Here comes a reproducible example:
##
# set up data: (don´t care
I like to keep the individual lines in my source files relatively
short. Mainly so that I can print them, email them, or display them on
a narrow screen without needing to shift left & right. So, for a
really long character string, such as an SQL query, I do something
like:
query=paste0("select CO
This seems like a more general statistics question rather than an R
question.
You might have more luck getting a response by posting your question to
Cross Validated, http://stats.stackexchange.com/
Jean
On Mon, Jul 21, 2014 at 10:13 AM, Lenth, Russell V
wrote:
> This is a question only about
On 22 Jul 2014, at 06:04 , David Winsemius wrote:
>
> On Jul 21, 2014, at 12:10 PM, Ronald Kölpin wrote:
>
>> Dear R-Community,
>>
>> I'm trying to estimate the parameters of a probability distribution
>> function by maximum likelihood estimation (using the stats4 function
>> mle()) but can't
On 22 Jul 2014, at 08:31 , PIKAL Petr wrote:
> Hi
>
> Anyway. AFAIK your code looks OK to me, provided you want find the number of
> rows without NA in each file.
>
...however, it is not clear to any of us whether the teacher has been
instructing students to avoid expanding the vector of an
Dear Community,
Just a short and simple question, but the code does not come to my mind.
I want to plot the Resids by quarterly Data. But there are many quarters.
Just rotate them, does not make the plot clearer/ easier to read.
Is there a possibility to just show every 4th value on the x-axi
Hello Friends,
I want to put two plots of scatterplotMatrix side by side in one plot.
I have tried below command. But somehow below code doesnt join them
together side by side.
> par(mfrow=c(1,2));
> scatterplotMatrix(~ SBP + DBP + Leptin +WC+BMI+Weight |
> factor(dat2$rs3827103), data=dat2,diag
There are 2 remaining seats on the following course:
Course: Data exploration, regression, GLM & GAM with
introduction to R.
Location: Murdoch University. Murdoch. Australia
When: 28 July - 1 August, 2014
Course flyer: http://www.highstat.com/Courses/Flyer2014_08Murdoch.p
27 matches
Mail list logo
