Dear R-helpers,
Does anyone know if the likelihoods calculated by these two packages are
comparable in this way?
That is, is this a valid likelihood ratio test?
# Reproducable example:
library(MASS)
library(nnet)
data(housing)
polr1 = MASS::polr(Sat ~ Infl + Type + Cont, weights=Freq, data=ho
On Jul 7, 2015, at 2:45 PM, Rosa Oliveira wrote:
> Iam trying to plot 6 graphs in one single plot and I was able to, nonetheless
> I wanted that all graphs had just 1 common legend, as the legend is the same
> for all the 6 graphs and there is no sense in repeating it 6 times and even
> more,
Iam trying to plot 6 graphs in one single plot and I was able to, nonetheless I
wanted that all graphs had just 1 common legend, as the legend is the same for
all the 6 graphs and there is no sense in repeating it 6 times and even more,
the legends in each graph sometimes don’t fit the graph.
I
?merge should do it.
John Kane
Kingston ON Canada
> -Original Message-
> From: lid.z...@gmail.com
> Sent: Tue, 7 Jul 2015 09:42:56 -0500
> To: r-help@r-project.org
> Subject: [R] add a special column to a matrix
>
> Hi there,
>
> I have a two matrices which they have a common column! I
> On 07 Jul 2015, at 16:20 , Dennis Fisher wrote:
>
> R 3.2.1
> OS X
>
> Colleagues
>
> I am trying to understand ANOVA contrasts and I have encountered some
> puzzling results.
>
> I have data from five studies with four different combinations of treatments.
> I want to make the following
The fact that it works when you pass your i values "by hand" should alert you
that perhaps your loop control does not do what you think it does.
Consider:
Your version:
for(i in 1:6 - 1)
What you probably meant:
for(i in 1:(6 - 1)) print(i)
The error is created on the first iteration of your l
for (i in 1:(ncol(MD_dist) - 1)){
...
}
Even better, replace
1:(n-1)
with
seq_len(n-1)
The latter does what you want (and empty integer vector) when n is 1;
the former would give c(1,0).
Bill Dunlap
TIBCO Software
wdunlap tibco.com
On Tue, Jul 7, 2015 at 2:11 PM, David Barron wr
Try reading the Help files before posting here. That's what they're for.
?get
provides the answer in a note in the Help page.
Cheers,
Bert
Bert Gunter
"Data is not information. Information is not knowledge. And knowledge
is certainly not wisdom."
-- Clifford Stoll
On Tue, Jul 7, 2015 at
You need to put the expression on the right of the colon in your for
statement in parenthesis:
for (i in 1:(ncol(MD_dist) - 1)){
...
}
On 7 July 2015 at 19:00, Karl Schilling wrote:
> Dear All:
>
> I want to use seq() inside a for-loop and use the looping counter i as the
> "by" argument in seq(
Dear list,
Let's say we have a variable (id), whose name is dynamically constructed.
This variable represents a vector or data frame with many elements. Now I
want to specifically assign a value to one of the elements. I couldn't get
it right.
test <- 'id' # "id" is dynamically constructed throug
Dear All:
I want to use seq() inside a for-loop and use the looping counter i as
the "by" argument in seq(). Here is some exemplary data and code:
# set up some data
distances <- c(0, NA, NA, NA, NA, NA,
5, 0, NA, NA, NA, NA,
5, 2, 0, NA, NA, NA,
18
Dear Arne,
Sorry for posting my mail twice and thanks a lot for your help.
Best regards,
Maram
Sent from my iPhone
> On Jul 7, 2015, at 9:55 PM, Arne Henningsen wrote:
>
> Dear Maram
>
> Please do NOT post your message twice!
>
> The warning messages occur each time, when maxLik() tries to
Dear Maram
Please do NOT post your message twice!
The warning messages occur each time, when maxLik() tries to calculate
the logLik value for theta[1] <= 0, theta[1] + theta[2] <= 0, theta[3]
<= 0 or something similar. According to the log-likelihood function,
it seems that the parameters theta[1
Hello,
Thanks John and David. I ll look into your suggestions and I ll ensure I
post more details the next time.
Thanks,
Vyshnnavi
On Jul 7, 2015 8:19 AM, "David L Carlson" wrote:
> As John suggests, more information and some sample data would help. I
> haven't used loocv, but here is a simple e
This is a **highly technical** statistical issue, not an R-Help topic.
I strongly suggest that you post to the R-sig-mixed-models list
instead.
Cheers,
Bert
Bert Gunter
"Data is not information. Information is not knowledge. And knowledge
is certainly not wisdom."
-- Clifford Stoll
On Tue, J
Dear Max,
I think that the object is to describe clearly what CRAN wants in the
description field so that package authors don't write description fields
that are unacceptable to CRAN. Clear criteria would save both package
authors' and CRAN maintainers' time. Although a mechanical check can find
s
I am trying to fit data from 23 subjects using random effects
regression, and am comparing the results of lme and lmer. The point
estimates and the SEs are the same in both models, however the degrees
of freedom are widely different. lme reports 88 DF, lmer approximately
22. Can someone help me und
As John suggests, more information and some sample data would help. I haven't
used loocv, but here is a simple example of leave one out cross validation with
simple linear regression. This does not eliminate the loop, but uses lapply()
which makes specifying the loop simpler. The results of the
Hi All,
I have just begun working on shiny package in R and have been building a
dashboard using tutorials and documentation abvailable online. I am stuck
in the middle.
It would be of great help if you could help me out with my code. I want to
build a dashboard that would take 'date range' and
> On 7 Jul 2015, at 01:12, peter dalgaard wrote:
>
>
>> On 06 Jul 2015, at 23:19 , Duncan Murdoch wrote:
>>
>> On 06/07/2015 5:09 PM, Rolf Turner wrote:
>>> On 07/07/15 07:10, William Dunlap wrote:
>>>
>>> [Rolf Turner wrote.]
>>>
> The CRAN guidelines should be rewritten so that they s
Hi
No HTML posting please, sometimes the mail is unreadable
It is a work for merge.
mat3 <- merge(mat1, mat2, all=TRUE)
shall do the trick (untested)
Cheers
Petr
> -Original Message-
> From: R-help [mailto:r-help-boun...@r-project.org] On Behalf Of Lida
> Zeighami
> Sent: Tuesday, Ju
Thanks to all who responded. It is easy to remove NA from factor so I will go
this way. I just wanted to know that I did not miss some simple canned solution.
Best regards
Petr
> -Original Message-
> From: Jeff Newmiller [mailto:jdnew...@dcn.davis.ca.us]
> Sent: Tuesday, July 07, 2015 3:
On Tue, Jul 7, 2015 at 8:19 AM, John Fox wrote:
> Dear Peter,
>
> You're correct that these examples aren't verb phrases (though the second
> one contains a verb phrase). I don't want to make the discussion even more
> pedantic (moving it in this direction was my fault), but "Paragraph" isn't
> q
Hi there,
I have a two matrices which they have a common column! I want to add the a
column of second to matrix to the equivalent column of first matrix!
my first matrix is
head(mat1)
a b c d
fg1:23 dfgv 5
pt10:18 tgtgh 1
wq 15:123oiljk
R 3.2.1
OS X
Colleagues
I am trying to understand ANOVA contrasts and I have encountered some puzzling
results.
I have data from five studies with four different combinations of treatments.
I want to make the following comparisons:
Study 1, 3: A / B vs. C / D (i.e., mean of A and B vs. m
Sorry, now it works for me as well. And yes, must drop rows with is.na( prace )
prior to calling ggplot. Note that the levels of prace do not contain NA so
your subject line is not correct.
---
Jeff Newmiller
Dear Michael,
I am very grateful to you for the detailed explanation, now everything is
clear. In the end I will opt for the Tuckey method then.
Thanks
Best
A.
Messaggio originale
Da: meyner...@pg.com
Data: 7-lug-2015 8.58
A: "angelo.arc...@virgilio.it"
Cc: "r-help@r-project.org", "
Dear John,
thank you, I normally use Tuckey but I was told that Bonferroni is better, so I
wanted to try.
For number of tests what do you mean exactly? How do I calculate it?
In the end is it sufficient to add test=adjusted("bonferroni") as previously
suggested?
Best
A.
Messaggio or
Thank you John, so far the easiest way is to use ezanova() since I do not have
to reorder the data
and create iframe etc.
Best
Messaggio originale
Da: pda...@gmail.com
Data: 7-lug-2015 10.12
A: "John Fox"
Cc: ,
Ogg: Re: [R] Sphericity for a repeated measure ANOVA with whithin and b
I disagree with you Jeff. I can reproduce it without any problem.
I see no other option than to remove the offending missing data prior to
plotting.
Best regards,
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and
Forest
team Biometrie & Kwaliteitszor
Not reproducible.
---
Jeff NewmillerThe . . Go Live...
DCN:Basics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO#.. Playing
Resear
On Tue, Jul 7, 2015 at 2:42 AM, Lionel SPINELLI
wrote:
> Hello all,
>
>
> I am looking for information on how to install R on a network disk so that
> it can be used by several users mounting this network disk from different
> computers.
>
> We are using a cluster and shared computers in my lab.
Sent from my iPhone
Begin forwarded message:
> From: "Maram Salem"
> Date: July 6, 2015 at 2:29:56 AM GMT+2
> To: "r-help@r-project.org"
> Subject: [R] NaN produced from log() with positive input
>
> Dear All
> I'm trying to find the maximum likelihood estimator of a certain
> distribution
Dear Peter,
You're correct that these examples aren't verb phrases (though the second one
contains a verb phrase). I don't want to make the discussion even more pedantic
(moving it in this direction was my fault), but "Paragraph" isn't quite right,
unless explained, because conventionally a par
Wrong diagnosis. The actual issue is a version-compare blunder for texinfo in
the configure script. This has been mended for r-patched, but r-release
requires a workaround.
-pd
On 07 Jul 2015, at 10:33 , peter dalgaard wrote:
> Some changes were made in connection with the release in order to
In line below
On 07/07/2015 11:20, peter dalgaard wrote:
...except that there is not necessarily a verb either. What we're looking for is
something like "advertisement style" as in
UGLY MUGS 7.95.
An invaluable addition to your display cabinet. Comes in an assortment of warts
and wrinkles, c
...except that there is not necessarily a verb either. What we're looking for
is something like "advertisement style" as in
UGLY MUGS 7.95.
An invaluable addition to your display cabinet. Comes in an assortment of warts
and wrinkles, crafted by professional artist Foo Yung.
However, I'm draw
I played around with the function a bit more and could actually narrow
down the problem some more. The problem seems to lie in the scale-option
of PCAgrid. So, for the example below, this code line produces a
princomp object with which the scores of Dat can be predicted correctly:
rPCA<-PCAgri
Hello Lionel,
Lionel SPINELLI writes:
> Hello all,
>>
> I am looking for information on how to install R on a network disk so
> that it can be used by several users mounting this network disk from
> different computers.
>
> We are using a cluster and shared computers in my lab. Those machines
>
Some changes were made in connection with the release in order to have better
automation of things like this (we now centralize information on the versions
in a VERSION-INFO.dcf).
I expect that something went wrong; maybe one of the periodic CRAN housekeeping
scripts doesn't quite do what it s
Another method goes via anova.mlm(), which is pretty straight forward once you
figure out how to set up the necessary within-subject projections (contrasts).
See
@Article{Rnews:Dalgaard:2007,
author = {Peter Dalgaard},
title= {New Functions for Multivariate Analysis},
journal
Hello all,
I am looking for information on how to install R on a network disk so that it
can be used by several users mounting this network disk from different
computers.
We are using a cluster and shared computers in my lab. Those machines can mount
network disks. On those disks we have inst
Dear all.
I want to drop unused (NA) level from geom_segment plot. I tried different
combination of drop = TRUE argument in several places of following code but NA
level of prace is still present.
p<-ggplot(snimek, aes(x=cas, y=prace, xend=endtime, yend=prace, colour=typ,
size=osoba))
p+geom_s
Angelo,
the conservatism of Bonferroni aside for a moment, I don't really see what you
are after.
mcp(Emotion = "Tukey")
just defines the contrasts to test; in this case, all pairwise comparisons.
The test=adjusted("bonferroni") argument then says that these contrasts should
not be correcte
44 matches
Mail list logo
