When you use the 'general linear model' analysis in SPSS, the first result is
a table with all terms with F-tests and significance values for all IV's. It
uses http://joyx.joensuu.fi/~ek/anova/sstypes.txt SS Type III , which has
the advantage that the order in which the variables are added to
drop1() in R meets both your objection and that to 'Type III'.
If you really want the so-called Type III, look at function Anova() in
package car. But also consider
library(fortunes); fortune(54); fortune((55); fortune(56)
On Fri, 4 Jul 2008, jeroenooms wrote:
When you use the 'general
I do not know what 'the R Utils docs' are, but --configure-args is plural!
You are not using the correct shell syntax: you need to use quotes when
arguments contain spaces. E.g.
--configure-args='--args=val1 --arg2=val2'
(and note that quotes are used in error info you quote).
That said,
Hello,
I'm trying to superimpose a line plot onto a histogram but I'm not having
any luck. I've attached the dataset. What I did was:
hist(data,freq=F)
Now I'm trying to superimpose the following points with a line connecting
them onto the histogram:
xy
100 0.535665393824959
200
Dear friends - I have a data.frame hh with columns
Na ID fit time
135 13 134.7945 60
133 13 132.9084 120
131 13 131.0848 180
129 13 129.3372 240
127 13 127.5546 300
126 13 125.8162 360
124 13 124.1836 420
123 13 122.6077 480
When I use lattice
Dear all,
I want to use the Bland-Altman method to measure agreement with repeated
measures collected over period of time (seven periods).
How can I do this with R
Many thanks
_
o.fr
On Fri, 2008-07-04 at 15:13 +0100, Prof Brian Ripley wrote:
On Fri, 4 Jul 2008, Ulrich Leopold wrote:
Dear list,
is there a way of updating the search path when using attach() for a data
set.
I am overwriting a variable in a data frame. To update teh search path I do
the follwoing:
On Fri, 2008-07-04 at 21:14 -0700, Michael Denslow wrote:
Dear R-helpers,
I am running metaMDS in the vegan package, which uses isoMDS in MASS,
to perform Nonmetric Multidimentional Scaling (NMDS).
I have seen some authors report a p-value for the NMDS ordination
based on randomization of
Dear j.c.l.omms,
The Anova() function in the car package can compute so-called type-II and
type-III tests. When there are terms in the model that are marginal to
others (such as main effects marginal to interactions), type-III tests are
of questionable interest and have to be formulated carefully
Thanks. Here is a similar example from a book by Pinheiro and Bates (2000,
chapter 6):
library(nlme) data(Soybean)
fm1Soy.lis - nlsList( weight ~ SSlogis(Time, Asym, xmid, scal),data =
Soybean ) fm1Soy.nlme - nlme( fm1Soy.lis )
If we would like to make comparisons among the years we
How about the answer by Demitris?
Regards a lot,
miltinho
From: Dimitris Rizopoulos [EMAIL PROTECTED]
Date: Jun 16, 2008 4:05 AM
Subject: Re: [R] Superimposing Line over Histogram in Density Plot
To: Gundala Viswanath [EMAIL PROTECTED]
Cc: [EMAIL PROTECTED]
try something like this:
x -
Hello.
Somebody knows how can I generate a set of n random vectors of a normal
distribution of several variables?
For example, I want to generate n=100 random vectors of two dimensions for
a normal with mean c(0,1) and variance matrix: matrix(c(2,1,1,3),2,2).
Thanks in advance,
Arnau.
Hi, Arnau,
mvrnorm() in MASS library is what you need.
? mvrnorm to see the detail but first you need to load the MASS
library, i.e,library(MASS)
regards/
On 2008-7-6, at 上午12:21, Arnau Mir wrote:
Hello.
Somebody knows how can I generate a set of n random vectors of a
normal
On Sat, 2008-07-05 at 18:21 +0200, Arnau Mir wrote:
Hello.
Somebody knows how can I generate a set of n random vectors of a normal
distribution of several variables?
For example, I want to generate n=100 random vectors of two dimensions for
a normal with mean c(0,1) and variance matrix:
on 07/05/2008 03:17 AM kende jan wrote:
Dear all,
I want to use the Bland-Altman method to measure agreement with
repeated measures collected over period of time (seven periods).
How can I do this with R
Many thanks
The following links should be helpful from a design/methodology
There is no need to load the MASS library, since the code for
mvrnorm therein is compact and self-contained:
mvrnorm - function (n=1, mu, Sigma, tol=1e-06, empirical=FALSE)
{
p - length(mu)
if(!all(dim(Sigma) == c(p, p)))
stop(incompatible arguments)
eS - eigen(Sigma, symmetric =
HI:
After following all the instructions on how to install the SciViews package
from CRAN I still can't make the GUI show. Can someone give me a hint on how to
do this? I have tried library(svGUI), library(svDialogs) and so on with the
rest of the packages. I have also placed all the 'sv'
Hello,
We need a little more information (which platform? which R version?).
Also, the GUI is NOT provided with the SciViews package. The version on
CRAN is quite outdated and is designed to communicate with the old
SciViews R Console GUI for Windows only. We are currently working on a
new
Hi ho:
Has anybody made use of this 'editor' as a front for R v2.7.0? I'm running
xp sp2. The texmacs website mentions the possibility of it but seems a
little short on instruction or direction to explain the how of it.
--
Brian Lunergan
Nepean, Ontario
Canada
---
Hello,
I am trying to extract t and pvalues from a 1000 ttests using the by-function
but everythinhg I tried did not work. Unfortunately googling by is not very
helpful. Any help will be very appreciated.
Cheers,
Danile Stall
*creating a data set
library(MASS)
dataset - mvrnorm(160, mu,
Hi
I just discovered the answer thanks to a previous thread from Peter Daalgart
The command is:
sapply(test, [[, statistic)
Cheers, Daniel
- Original Message -
To: r-help@r-project.org
Subject: extracting values from a by function
Date: Sun, 6 Jul 2008 05:59:04 +0800
Hello,
I am
library(MASS)
attach(bacteria)
table(y)
y
n y
43 177
y-1*(y==y)
table(y,trt)
trt
y placebo drug drug+
0 12 1813
1 84 4449
library(lme4)
model1-lmer(y~trt+(week|ID),family=binomial,method=PQL)
Error in match.arg(method, c(Laplace, AGQ)) :
'arg'
hpdutra at yahoo.com writes:
library(MASS)
attach(bacteria)
table(y)
y
n y
43 177
y-1*(y==y)
table(y,trt)
trt
y placebo drug drug+
0 12 1813
1 84 4449
library(lme4)
model1-lmer(y~trt+(week|ID),family=binomial,method=PQL)
Error in
I am trying to do an exercise on Crawley's book and I get this error message. I
ran the same analysis in the past but now that I upgraded my R in my Macbook I
get this message. Is this a bug or what is going on.
It doesn't make any sense
library(MASS)
attach(bacteria)
table(y)
y
n y
43
Hi, I am trying to do a repeated measures ANOVA to determine if there is a
significant difference between two sets of timecourse data. Each individual
was given a single treatment and then measured for one variable for 10 days.
Here is made-up example of what my data would look like:
Hello.
Somebody knows how can I generate a set of n random vectors of a
normal distribution of several variables?
For example, I want to generate n=100 random vectors of two dimensions
for a normal with mean c(0,1) and variance matrix:
matrix(c(2,1,1,3),2,2).
Thanks in advance,
The verbose option gives a display like:
rf.500 -
+ randomForest(new.x,trn.y,do.trace=20,ntree=100,nodesize=500,
+importance=T)
| Out-of-bag |
Tree | MSE %Var(y) |
20 | 0.9279 100.84 |
What is the meaning of %var(y)100%? I expected that to correspond
Thanks for the quick response. But the points I'm trying to superimpose
isn't exactly the density of the data. Is there any other way to do it?
From: milton ruser [mailto:[EMAIL PROTECTED]
Sent: July 5, 2008 6:56 AM
To: Edwin Lei
Cc: r-help@r-project.org
Subject: Re: [R] trying to superimpose
Hi , Arnau
Did you ever check your mailbox? your question was answered last
night Beijing time. :)
Just read the following .
-
There is no need to load the MASS library, since the code for
mvrnorm therein is compact and self-contained:
mvrnorm - function (n=1, mu,
library(MASS)
set.seed(22134)
Y - mvrnorm(100, mu = c(0,1), Sigma = matrix(c(2,1,1,3),2,2))
colMeans(Y)
[1] 0.06461359 1.08436419
var(Y)
[,1] [,2]
[1,] 1.683963 1.300363
[2,] 1.300363 3.280612
Bill Venables
CSIRO Laboratories
PO Box 120, Cleveland, 4163
AUSTRALIA
Office Phone
What's the problem with loading the MASS library???
The MASS library is one of the recommended packages, so should be universally
available.
Pinching code like this and freezing it inside your personal scripts means that
if ever we find bug fixes or improvements you miss out on them.
Bill
31 matches
Mail list logo