The fact related to Rolf's vague recollection is here:
https://stat.ethz.ch/pipermail/r-help/2015-March/426799.html
And the relevant post (Martin Maechler's answer):
https://stat.ethz.ch/pipermail/r-help/2015-March/426834.html
Should that be a FAQ?
Helios
> On 16/03/17 03:57, Bert Gunter wrot
s.
In addition, I would like to refer the package users to the development
repository in GitHub: <https://github.com/heliosdrm/phia>, where new
contributions, issues/bug reports, and other comments are welcome.
Best regards
Helios De Rosario
biomecanicamente.org
Conoce la actualidad on-line d
use alone, but you want
to distribute it. And if you want to distribute it by mainstream
repositories (e.g. CRAN), you should publish the source anyway, and
won't be allowed to include such "malicious" mechanisms.
Helios De Rosario
>>> El día 12/05/2014 a las 11:31, Ashis Deb
ot use lm(). Perhaps you are looking for ?glm (family "binomial") if it
is binomial data, or ?polr (in package MASS), ?multinom (in nnet), or
?clm and others (in ordinal) if the variable has more than two levels.
Helios De Rosario
INSTITUTO DE BIOMECÁNICA DE VALENCIA
Universidad Politécnic
versity.edu/avdrark/mokken
The ESD method for identifying outliers, also used in the paper to
handle outlier scores, is implemented (together with others) in the
package "parody":
http://www.bioconductor.org/packages/release/bioc/html/parody.html
Hope it helps
Helios De Rosario
INSTITU
--
Helios de Rosario Martínez
Researcher
>>> El día 25/07/2013 a las 11:44, "mei_yuan"
escribió:
> Hi,
>
> In the R console, I have the following:
>
>> runif(10)
> Error in runif(10) :
> '.Random.seed' is not an integer vector
kzDevice.
>
>How can I avoid to get these messages?
There are several options in ?library to control the messages that are
displayed when loading packages. However, this does not seem be able to
supress all the messages. Some messages are defined by the package
authors, because they feel ne
>>> El día 15/07/2013 a las 13:30, Duncan Murdoch
escribió:
> On 13-07-15 5:27 AM, Helios de Rosario wrote:
>> Hi,
>>
>> What should I do if I want to use "non-default" R-options for the
>> functions inside a package, but not affect the options of
disadvantages of the previous
solutions, but perhaps it is overkill, hence my question: Is there a
"clean" way of setting the R-session options in such a way that they
only affect the functions and data inside a package?
Thanks and best regards
Helios De Rosario
INSTITUTO DE BIOMECÁNICA DE
tion
testInteractions() from package "phia". For instance:
testInteractions(m1)
will give you a test of all the pairwise contrasts between color and
season. You can also test simple main effects, or other specific
contrasts by adding further arguments (see the documentation and th
A (late) update to this question:
On Fri Aug 17 07:33:29, Henrik Singmann wrote:
> Hi Diego,
>
> I am struggeling with this question also for some time and there does
> not seem to be an easy and general solution to this problem. At least
I
> haven't found one.
> However, if you have just one re
=rep(0,8))
DF2$y[DF1$x] <- DF1$y
But I'd recommend you use NA for the missing values, unless you have a
good reason to code them as zero.
> Many many thjanks in advance.
You are welcome.
Helios De Rosario
--
Helios de Rosario Martínez
Researcher
INSTITUTO DE BIOMECÁNICA
r
>
> The plural of anecdote is not data.
> ~ Roger Brinner
>
> The combination of some data and an aching desire for an answer does
not
> ensure that a reasonable answer can be extracted from a given body of
data.
> ~ John Tukey
>
>
> -Oorspronkelijk bericht-
&
cide:
> tapply(predict(mod2),shoe,FUN=mean)
S1 S2 S3
118.608 118.608 118.608
My questions are:
Is this correct? And is there an easier way of doing this?
Thanks
Helios De Rosario
--
Helios de Rosario Martínez
Researcher
INSTITUTO DE BIOMECÁNICA DE VALENCIA
Universidad P
p
page of that function, or in the web appendix to the "CAR" book:
http://socserv.mcmaster.ca/jfox/Books/Companion/appendix/Appendix-Multivariate-Linear-Models.pdf
Helios De Rosario
--
Helios de Rosario Martínez
Researcher
>>> El día 20/05/2012 a las 0:58, jacaranda tree
Hi, after re-reading I think that I misunderstood your question. You
don't provide many details, but I suppose that the "brood" effect is
nested within the fixed effects, so you don't mean a multivariate
approach for a split-plot or a repeated-measures design, but the
analysis of a multivariate mix
-Multivariate-Linear-Models.pdf
The package ez also has ezANOVA(), a wrapper to Anova() that may
simplify the task for full factorial designs.
Helios De Rosario
>>> El día 17/05/2012 a las 12:05, "David Costantini"
escribió:
> Dear All
> I would need to perform a MANOVA with both
I wrote in my previous message the following Octave code:
[Octave code:]
octave:1> x=1;
octave:2> save -ascii testdata.mat x
Forget the "-ascii". It should be "-text" or nothing ("-text" is the
default).
By the way, read.octave() does not really "fail" (it does return a
value), but the result is
RY=Spanish_Spain.1252 LC_NUMERIC=C
[5] LC_TIME=Spanish_Spain.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] foreign_0.8-49
--
Helios de Rosario Martínez
Researcher
INSTITUTO DE BIOMECÁNICA DE VALENCIA
Uni
>>> El día 21/03/2012 a las 11:27, Marco Tommasi
wrote
> To whom it may concern
>
> I made some analysis with R using the command Anova. However, I found
> some problmes with the output obtained by selecting type II o type
III
> sum of squares.
>
> Briefly, I have to do a 2x3 mixed model anova
You can tell that levels must not be dropped with scale_x_discrete():
library(ggplot2)
mtcars$cyl<-factor(mtcars$cyl)
plot1 <- ggplot(mtcars[!mtcars$cyl==4,], aes(cyl))+geom_bar()
plot1 + scale_x_discrete(drop=FALSE)
Or explicitly set the values you want in the x axis with the argument
"limits":
her report of a possible bug in facet_grid() today
that
> could be related:
>
https://groups.google.com/group/ggplot2/browse_thread/thread/5213ac35da6b36d
> 4
>
> Michael
>
> On Mon, Mar 12, 2012 at 7:16 AM, Helios de Rosario
> wrote:
>> I have encountered a p
Jianghong wrote:
> Hello,
>
> I have created two boxplots with following R code. There is one
outlier in
> B group.
> The outlier is 33. But the all other data are between 0 to 4.
>
> How can I skip y-axis around 5 to 25, and expand 0-4 for this case.
Also I
> want keep the outlier in my boxplot.
I have encountered a problem with faceted bar plots. I have tried to
create something like the example explained in the ggplot2 book (see pp.
126-128):
library(ggplot2)
mpg4 <- subset(mpg, manufacturer %in%
c("audi", "volkswagen", "jeep"))
mpg4$manufacturer <- as.character(mpg4$manufacture
--
Helios de Rosario Martínez
Researcher
>>> El día 22/02/2012 a las 11:32, Florian Weiler
escribió:
> Dear all,
>
> I have a (probably very basic) question. I am imputing data with the
mice
> package, using 10 chains. I can then write out the 10 final values of
t
That's exactly what I wanted. Thanks!
Helios
>>> El día 21/01/2012 a las 0:33, David Winsemius
escribió:
> On Jan 20, 2012, at 6:26 PM, Helios de Rosario wrote:
>
>> Hi, a quick question:
>>
>> Is there a way to retrieve the default value of a function arg
Hi, a quick question:
Is there a way to retrieve the default value of a function argument -
if it exists?
(I know I can see it if I type the function name, but I would like get
the value programaticaly.)
Thanks,
--
Helios de Rosario Martínez
Researcher
INSTITUTO DE BIOMECÁNICA DE
in.1252
[3] LC_MONETARY=Spanish_Spain.1252 LC_NUMERIC=C
[5] LC_TIME=Spanish_Spain.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
--
Helios de Rosario Martínez
Researcher
__
R-help@r-project.org mai
Dave, your situation is clearer now. You wrote (see the full context at
the end of the message):
> From this, you will see that I have 4 control sites and 7 treatment
> sites that are measured each week. All 13 locations have different
> names, and Location is a random varaible. Is Location n
Dear Dave, there are some inconsistencies in your explanation of the
problem. You said your variables are:
> CO is a continuous response variable,
>
> Week is a fixed categorical factor,
>
> Habitat is a fixed categorical factor, and
>
> Location is a random categorical factor nested within Hab
I have a doubt about the calculation of tests for fixed effects in
mixed-effects models.
I have read that, except in well-balanced designs, the F statistic that
is usually calculated for ANOVA tables may be far from being distributed
as an exact F distribution, and that's the reason why the anova
To separate the parts of a formula, use as.character
(check the examples in ?character)
Helios
22 Sep 2011 16:14:05 -0400
From: Jean-Christophe BOU?TT?
> Hello,
> You can check ?model.frame.
> I do not know however to extract only the right-hand of left-hand
part
> of a formula.
>
> JC
>
> 201
For some time I have been looking for a convenient way of performing
post-hoc analysis to Repeated Measures ANOVA, that would be acceptable
if sphericity is violated (i.e. leaving aside post-hoc to lme models).
The best solution I found was John Fox's proposal to similar requests
in R-help:
http:/
33 matches
Mail list logo
