On 7/9/07, Alex Baugh [EMAIL PROTECTED] wrote:
I am executing a Repeated Measures Analysis of Variance with 1 DV
(LOCOMOTOR
RESPONSE), 2 Within-Subjects Factors (AGE, ACOUSTIC CONDITION), and 1
Between-Subjects Factor (SEX).
Does anyone know whether the between-subjects factor (SEX)
On Mon, 9 Jul 2007, Thomas Lumley wrote:
On Mon, 9 Jul 2007, Jim Lemon wrote:
Hi folks,
I thought I recalled a request for turning a character string into an
object name as in:
Yes. It's a FAQ.
There is an FAQ about turning character strings into objects (Q7.21), but
this seems a bit
Hello.
I'm playing with Kalman filter.
library(sspir)
m1 - SS(y=t(t(sin(1:20 # dummy data
plot(m1$y, type=l)
m1.f - kfilter(m1)
m1.s - smoother(m1.f)
lines(m1.f$m, lty=dotted,col=blue)
lines(m1.s$m, lty=dotted,col=red)
I was wondering how it's possible to forecast using sspir library.
I
Hi,
I have some problems with my repeated measures analysis. When I compute it
with SPSS I get different results than with R. Probably I am doing something
wrong in R.
I have two groups (1,2) both having to solve a task under two conditions
(1,2). That is one between subject factor (group) and
In the absence of a data set, it may help to read the help file carefully:
?hist
Note, in particular, that the argument freq defaults to TRUE if and only if
breaks are equidistant (and probability is not specified).
Regards,
Mark.
Sarah Goslee wrote:
Well, how about an example of what you
confirm 03122ac6f26ed0e761c70bde1bb41c3063b82211
[EMAIL PROTECTED] wrote: The results of your email command are provided below.
Attached is your
original message.
- Results:
Ignoring non-text/plain MIME parts
- Unprocessed:
Regards
Dimple
[EMAIL PROTECTED] wrote: Mailing list
We don't know what you are after (or what you did in SPSS), but
dat - read.table(mb2.dat, header=TRUE,
colClasses=c(rep(factor,3), double))
summary(aov(Score ~ Group * Task + Error(Id), dat))
would seem to be the sort of thing your description indicates.
If you tell us
mb2 wrote:
Hi,
I have some problems with my repeated measures analysis. When I compute it
with SPSS I get different results than with R. Probably I am doing something
wrong in R.
I have two groups (1,2) both having to solve a task under two conditions
(1,2). That is one between subject
dimple thyagarajan wrote:
...
Building a website is a piece of cake.
But sending an email is another matter altogether.
Jim
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Hello Yihsu,
have a look at ?cajorls. With this function a VECM is estimated, whence
the cointegration rank has been determined (ca.jo). For further
analysis, you might want to consider the function vec2var in package
vars and methods irf, fevd and predict, as well as the diagnostic tests
that
Hi Sir
In iid.test package the function iid.test() has been developed for N sites
having same number of observations. Can we use this function for a single
site? Because I have sites with different numer of observations.
Regards
--
AMINA SHAHZADI
Department of Statistics
GC University Lahore,
Hello
The excellent adapt package integrates over multi-dimensional
hypercubes.
I want to integrate over a multidimensional simplex. Has anyone
implemented such a thing in R?
I can transform an n-simplex to a hyperrectangle
but the Jacobian is a rapidly-varying (and very lopsided)
function and
Hello,
is it possible to obtain type III sums of squares for a nested model as
in the following:
lmod - lm(resp ~ A * B + (C %in% A), mydata))
I have tried
library(car)
Anova(lmod, type=III)
but this gives me an error (and I also understand from the documentation
of Anova as well as from a
Dear all
Thanks in advance for replies
I am trying to make a barplot out of this data which I read from a file
tb - read.table(tmp.dat
, na.string=c(-))
tmp.dat:
(the file is much longer and includes NAs as -)
#A1A2A3B1B2B3C1C2C3D1
Carsten Jaeger wrote:
Hello,
is it possible to obtain type III sums of squares for a nested model as
in the following:
lmod - lm(resp ~ A * B + (C %in% A), mydata))
I have tried
library(car)
Anova(lmod, type=III)
but this gives me an error (and I also understand from the documentation
Dear all,
I have been trying to build R-2.5.1 on the Interix-6.0 Unix subsystem
that ships with Vista and everything looks fine during the configure
except towards the end when sed throws an error:
sed: 1: s/\*/\\\*/g: invalid command code
A few lines later I then get:
./configure: : bad
The message from this cute little data set is very clear. Consider
fm - aov(resp ~ A*B + A/C, mydata)
drop1(fm, test = F)
Single term deletions
Model:
resp ~ A * B + A/C
Df Sum of Sq RSS AIC F value Pr(F)
none 65.540 47.261
A:B 216.132
Hello,
I think that to apply the Neuman-Keuls test under R it's possible with the
TukeyHSD function and I know that to do a Neuman-Keuls test I have to sort the
means but I can't use the TukeyHSD function I don't understand how to do to
apply this function to my data.
In the R help they use the
Oakeley, Edward wrote:
Dear all,
I have been trying to build R-2.5.1 on the Interix-6.0 Unix subsystem
that ships with Vista and everything looks fine during the configure
except towards the end when sed throws an error:
sed: 1: s/\*/\\\*/g: invalid command code
A few lines later I then
Hi,
Can anyone help me with repeated meausres MANOVA in R ? For repeated measures
ANOVA I used function aov. Is there something like this exists for MANOVA?
Thanks,
Deepa
-
[[alternative HTML version deleted]]
See
?summary.manova
--
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40 S 49° 16' 22 O
On 10/07/07, deepa gupta [EMAIL PROTECTED] wrote:
Hi,
Can anyone help me with repeated meausres MANOVA in R ? For repeated
measures ANOVA I used function aov. Is there something like this exists
Hey all,
I'm doing some work with machine learning on R (I'm a fairly new user of R),
and I have a question about generating new tables from existing tables. I'm
currently using a table of measurements I read in from a CSV file to
generate training and validation data set tables for future use
barchart(Titanic, stack=F) produces a very nice horizontal barchart.
Each panel has four groups of two bars.
barchart(Titanic, stack=F, horizontal=F) doesn't produce the results I
would have expected, as it produces this warning message:
Warning message:
y should be numeric in:
You will get more useful answers if you specify exactly how you want
to overlay the boxplots (overlay them on what?). You can certainly do
this with the ggplot2 package, or lattice or base graphics.
Hadley
On 7/10/07, Hao Liu [EMAIL PROTECTED] wrote:
hi, All:
I need to overlay two boxplot, I
Greg,
You are certainly welcome to use my Perl script as the basis for a
write.xls() function for gdata. You can even change the name of the
script to csv2xls.pl if you wish, for consistency with the existing
function.
If you might want to use the script largely 'as is', I won't have time
for a
Hi all,
I would like to plot part of the emperical CDF. Suppose the variable is x, I
just need the part when x1,therefore, I am using the following codes.
tail - x1
plot(ecdf(x[tail]), do.points=FALSE, verticals=TRUE)
The x value starts from 1, but the yaxs still begins from 0, not the
Michael Hoffman said the following on 7/10/2007 7:06 AM:
barchart(Titanic, stack=F) produces a very nice horizontal barchart.
Each panel has four groups of two bars.
barchart(Titanic, stack=F, horizontal=F) doesn't produce the results I
would have expected, as it produces this warning
Hi Robin,
A Monte-Carlo approach could be attempted, if one could generate samples that
are either uniformly distributed over the simplex. There is a small section in
Luc Devroye's book (Generation of Non-uniform random deviates) on random
uniform sampling from a simplex, if I remeber
Thanks... I just realized using add=TRUE will work...
Best
Hao
hadley wickham wrote:
You will get more useful answers if you specify exactly how you want
to overlay the boxplots (overlay them on what?). You can certainly do
this with the ggplot2 package, or lattice or base graphics.
Hadley
Date Fo Co6/27/2007 57.1 13.96/28/2007 57.7 14.3
6/29/2007 57.8 14.36/30/2007 57 13.97/1/2007 57.1 13.9
7/2/2007 57.2 14.07/3/2007 57.3 14.17/4/2007 57.6 14.2
7/5/2007 58 14.47/6/2007 58.1 14.57/7/2007 58.2 14.6
hi, All:
I need to overlay two boxplot, I played around with points() but found
it does not seem to work with boxplot, it works fine with other. Is
there a way to overlay two boxplot (using different color) in R?
There was a thread talking about using ggplot package, however, I don't
think
On 7/10/2007 10:36 AM, livia wrote:
Hi all,
I would like to plot part of the emperical CDF. Suppose the variable is x, I
just need the part when x1,therefore, I am using the following codes.
tail - x1
plot(ecdf(x[tail]), do.points=FALSE, verticals=TRUE)
The x value starts from 1, but
Felipe Carrillo said the following on 7/10/2007 7:58 AM:
Date Fo Co6/27/2007 57.1 13.96/28/2007 57.7 14.3
6/29/2007 57.8 14.36/30/2007 57 13.97/1/2007 57.1 13.9
7/2/2007 57.2 14.07/3/2007 57.3 14.17/4/2007 57.6 14.2
7/5/2007 58
Sundar Dorai-Raj wrote:
It seems that barchart.table doesn't allow the horizontal = FALSE
argument. With a slight modification to barchart.table this can be
accomplished.
Thanks for supplying that.
Also, I don't get a warning with your original code using
R-2.5.1 and lattice 0.16-1.
Dear users,
please help to define the following data structure:
I would like to have a matrix, where every element is a container of
different size , containing real numbers. The containers (bins) are
addressed by an index pair [i,j] (i is number of corresponding row of
the matrix, j
Try this:
m - matrix(list(1, 1:2, 1:3, 1:4), 2)
m[[1,1]]
[1] 1
m[[2,1]]
[1] 1 2
m
[,1] [,2]
[1,] 1 Integer,3
[2,] Integer,2 Integer,4
On 7/10/07, Balazs Torma [EMAIL PROTECTED] wrote:
Dear users,
please help to define the following data structure:
I would like to
Hi,
I am writing an R package with two functions in C++. So far
everything works.
Now, i would like to write a third function which would use a pointer
(it is a pointer to a class object) created by first function.
I tried placing this pointer outside of the function
I'm trying to figure out why when I use as.character() on one row of a
data.frame, I get factor numbers instead of a character vector. Any
suggestions?
See the following code:
a-c(Abraham,Jonah,Moses)
b-c(Sarah,Hannah,Mary)
c-c(Billy,Joe,Bob)
df-data.frame(a=a,b=b,c=c)
#Suppose I'm interested
I nominate the following 2 pieces from Bill's reply for fortunes
(probably 2 separate fortunes):
All this becomes even more glaring if you take the unusal
step of plotting the data.
and
What sort of editor would overlook this clear and
demonstrable message leaping out from the data in
Thank you very much.
Duncan Murdoch-2 wrote:
On 7/10/2007 10:36 AM, livia wrote:
Hi all,
I would like to plot part of the emperical CDF. Suppose the variable is
x, I
just need the part when x1,therefore, I am using the following codes.
tail - x1
plot(ecdf(x[tail]),
Hi,
Some progress: I am using
SEXP retty;
book=Calloc(1,int);
*book=10;
PROTECT(retty=R_MakeExternalPtr(book,R_NilValue,R_NilValue));
then UNPROTECTING and returning retty.
In a another function,
foo(SEXP s){
int*
Hello all,
I would like to plot the emperical CDF, normal CDF and pareto CDF in the
same graph and I amusing the following codes. z is a vector and I just
need the part when z between 1.6 and 3.
plot(ecdf(z), do.points=FALSE, verticals=TRUE,
xlim=c(1.6,3),ylim=c(1-sum(z1.6)/length(z), 1))
x -
Maybe this is what you want:
you are right about the re-allocating the tables, but you can subset your
table into a new one:
selection - which(device_Prob_Vector 0.5)
# or via sample: selection - sample(num_Devices)
training_Set - measurements[selection]
validation_Set -
Dear All,
I have about 50 native ArcGis GRID maps and I need read them on R.
After that I need to cross these maps with a set of XY coordinates which are
stored on a table. When XY coordinates of my table match with the pixels of my
maps, I would like to store the values of the pixels as a
On 7/10/2007 6:57 AM, Robin Hankin wrote:
Hello
The excellent adapt package integrates over multi-dimensional
hypercubes.
I want to integrate over a multidimensional simplex. Has anyone
implemented such a thing in R?
I can transform an n-simplex to a hyperrectangle
but the Jacobian
Original Message
Subject: [R] ECDF, distribution of Pareto, distribution of Normal
From: livia [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Date: Tue Jul 10 2007 18:35:04 GMT+0200
Hello all,
I would like to plot the emperical CDF, normal CDF and pareto CDF in the
same
Hi. I'm trying to export a dataframe from R into Stata to use a statistical
function I have there. I attached library write.foreign and renamed my
variables to get them to match Stata's required format, and now have the
following error: file /tmp/Rtmps7rmrM/file1c06dac8.raw not found Other
Andrew:
As you haven't received a reply yet ...
?factor,?UseMethod, and An Introduction to R may help. But it's a bit
subtle.
Factors are objects that are integer vectors (codes) with a levels attribute
that associates the codes with levels as character names. So
df[df$a==Abraham,] is a
Hello all --
I would like to thicken the borders between panels -- or more generally,
all borders -- in a plot generated using lattice (specifically,
levelplot). Something similar perhaps to box() function in graphics.
I haven't been successful in reviewing available documentation.
The
I am not sure what you are doing there but what you need is
library(foreign)
and
write.dta()
see
?write.dta once you have loaded the foreign package
Stefan
Original Message
Subject: [R] Help with write.foreign (exporting data to Stata)
From: kdestler [EMAIL PROTECTED]
To:
On Jul 10, 2007, at 12:33 PM, Saptarshi Guha wrote:
Hi,
Some progress: I am using
SEXP retty;
book=Calloc(1,int);
*book=10;
PROTECT(retty=R_MakeExternalPtr(book,R_NilValue,R_NilValue));
then UNPROTECTING and returning retty.
In a another
On 7/10/07, Sundar Dorai-Raj [EMAIL PROTECTED] wrote:
Michael Hoffman said the following on 7/10/2007 7:06 AM:
barchart(Titanic, stack=F) produces a very nice horizontal barchart.
Each panel has four groups of two bars.
barchart(Titanic, stack=F, horizontal=F) doesn't produce the
On 7/10/07, Hil Lyons [EMAIL PROTECTED] wrote:
Hello all --
I would like to thicken the borders between panels -- or more generally,
all borders -- in a plot generated using lattice (specifically,
levelplot). Something similar perhaps to box() function in graphics.
I haven't been successful
On Tue, 10 Jul 2007, Stefan Grosse wrote:
I am not sure what you are doing there but what you need is
library(foreign)
and
write.dta()
write.foreign should also work, though.
My guess is that Kate used tempfile() to specify the filenames, and that the
data file would then have been
I have a data frame, lets call it dat,
with 3 columns ( mc, yr, ret) which represent market
cap, year, and return. mc is a factor, mc, and ret are
real numbers.
I want to add a column to the data calculated as
follows.
For each year, I want to split the data by mc decile,
then calculate the mean
This should work (with x containing the dataframe):
x$Id=factor(x$Id)
x$Group=factor(x$Group)
x$Task=factor(x$Task)
str(x)
'data.frame': 48 obs. of 4 variables:
$ Id : Factor w/ 24 levels 1,2,3,4,..: 1 2 3 4 5 6 7 8 9
10 ...
$ Group: Factor w/ 2 levels 1,2: 1 1 1 1 1 1 1 1 1 1 ...
I'm running windows xp, R 2.3.1 with maptools 0.6-6, I guess.
When plotting from a large SpatialLinesDataFrame and using xlim ylim to
reduce the area, the plot axes automatically have the same scale size, even if
xlim and ylim ranges differ.
E.g.:
tmp - readShapeLines(filepath)
Jose,
The Gehan-Breslow test provides a generalization of the Kruskal-Wallis test for
censored data. As an alternative, try using survdiff with rho=1. This method
uses weights w(ti) = S(ti) (where S is the Kaplan-Meier estimate of survival)
which yields Fleming and Harrington's version of
We assume the Fo and Co represent the same data except in different
units (this seems to be approximately the case) so there is really
only one variable being measured here. If that's not the case let
me know. Below we read the data, define enough padding around plot
to do what we want, call
Dear useRs,
Version 0.2 of package relations appeared on CRAN and is currently
propagating to the mirrors. In addition to some bug fixes, the new
release includes:
o an introductory vignette showing the main features;
o new SD fitters for the C (complete) and A (antisymmetric)
I second the nomination!
Simon.
On Tue, 2007-07-10 at 10:02 -0600, Greg Snow wrote:
I nominate the following 2 pieces from Bill's reply for fortunes
(probably 2 separate fortunes):
All this becomes even more glaring if you take the unusal
step of plotting the data.
and
What
Hello all,
I would like to use the lp function (lpSolve package) on a linux system. Using
sample code from the lp function help, I can solve the linear program with no
errors or problems. However, when I copy the exact same code to R in linux
(after loading the lpSolve package), I get the
Hi list,
has anyone implemented the Gap statistic for clusters based on
categorical/mixed data?
Slmisc (and SAGx) only work on numerical data.
Any suggestions welcome
Thanks
Herry
__
R-help@stat.math.ethz.ch mailing list
63 matches
Mail list logo