Hi R-users,
I have a data set like this (first ten rows):
id patient_id date code class eala ID1564262 1562 6.4.2006 12:00 1
NA ID1564262 1562 6.4.2006 12:00 1 NA ID1564264 1365 14.2.2006 14:35
1 50 ID1564265 1342 7.4.2006 14:30 2 50 ID1564266 1648 7.4.200614:30
2 50
Dear Listers,
I would be interested in representing a trend surface including an
interaction term z = f(x,y,x.y) - eg the type of chart obtained with
persp() or wireframe(), then adding datapoints as a cloud, ideally with
dots which are under the surface in a color, and those who are above in
jim holtman wrote:
matches - sapply(result, function(x){
+ .comma - strsplit(x, ',')[[1]] # get the fields
+ any(my.list %in% .comma)
+ })
Thanks for that!
Joh
__
R-help@stat.math.ethz.ch mailing list
Hi All-
This seems like such a pathetic problem to be posting about, but I have no
idea why this testcase does not work. I have tried this using R 2.4.1,
2.4.0, 2.3.0, and 2.0.0 on several different computers (Mac OS 10.4.8,
Windows XP, Linux). Below the signature, you will find my test case R
[EMAIL PROTECTED] wrote:
Hello all,
I am low down on the learning curve of R but so far I have had little
trouble using most of the packages. However, recently I have run into
a wall when it comes to bootstrapping a Levene's test (from the car
package) and thought you might be able to
FAQ Q7.31
On Tue, 20 Feb 2007, Todd A. Johnson wrote:
Hi All-
This seems like such a pathetic problem to be posting about, but I have no
idea why this testcase does not work. I have tried this using R 2.4.1,
2.4.0, 2.3.0, and 2.0.0 on several different computers (Mac OS 10.4.8,
Windows
This is a function of your data and the tuning parameters you chose to
use. See ?tree.control.
On Tue, 20 Feb 2007, stephenc wrote:
Hi
I am trying to use tree() to classify movements in a futures contract. My
data is like this:
diff dip dim adx
I want to calculate the probability that a group will include a particular
point using the squared Mahalanobis distance to the centroid. I understand
that the squared Mahalanobis distance is distributed as chi-squared but that
for a small number of random samples from a multivariate normal
--- amna khan [EMAIL PROTECTED] wrote:
Sir I am very new user of R. I am not understanding
the how to write the
plot description in box outside the plot.
Have a look at
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/68585.html
You should read up on par
?par
I hope this helps
First I have to say I am sorry because I have not been so clear in my
previous e-mails. I will try to explain clearer what it is my problem.
I have the following model:
lnP=Sc+Ag+Ag2+Var+R+D
In this model the variable Sc is endogenous and the rest are all objective
exogenous variables. I
Haven't quite learned to 'cast' yet, but I have always used the 'apply'
functions for this type of processing:
x - id patient_id date code class eala
+ ID1564262 1562 6.4.200612:00 1 NA
+ ID1564262 1562 6.4.200612:00 1 NA
+ ID1564264 1365 14.2.200614:35 1 50
+
Dear R-list,
I have to design the validation of a score (ordinal values between 0 and 6)
reputed to separate 4 groups of patients with known frequencies in the
population. I think the more accurate is to calculate sample size under median
test. Is there a function for that in R (not in the pwr
I don't see why making copies of the columns you need inside the loop is
better memory management. If the data are in a matrix, accessing
elements is quite fast. If you're worrying about speed of that, do what
Charles suggest: work with the transpose so that you are accessing
elements in the
Haven't quite learned to 'cast' yet, but I have always used the 'apply'
functions for this type of processing:
x - id patient_id date code class eala
+ ID1564262 1562 6.4.200612:00 1 NA
+ ID1564262 1562 6.4.200612:00 1 NA
+ ID1564264 1365 14.2.200614:35 1 50
+
At 15:21 19/02/2007, Paolo Accadia wrote:
Hi all,
I have a problem to estimate Std. Error and
t-value by âpolrâ in library Mass.
They result from the summary of a polr object.
I can obtain them working in the R environment with the following statements:
temp - polr(formula =
Liaw, Andy wrote:
I don't see why making copies of the columns you need inside the loop is
better memory management. If the data are in a matrix, accessing
elements is quite fast. If you're worrying about speed of that, do what
Charles suggest: work with the transpose so that you are
SAULEAU Erik-André wrote:
Dear R-list,
I have to design the validation of a score (ordinal values between 0 and 6)
reputed to separate 4 groups of patients with known frequencies in the
population. I think the more accurate is to calculate sample size under
median test. Is there a
Hi everybody!
Is there any way to read a postscrit file into R?
All the best to you
Ralf Finne
SYH University of Applied Sciences
Vasa Finland
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Thanks guys! Jim got it right what I was looking for. I think that the
reason why I don't get this right with cast-function is that I don't know
how to write a proper stats-function. SQL is still much easier for me.
Lauri
2007/2/20, hadley wickham [EMAIL PROTECTED]:
Haven't quite learned to
Hello!
I written several implementations in R of Rémy's algorithm for
generating random ordered strictly binary trees on 2n+1 vertices.
One implementation involves manufacturing character strings like:
X - list(0,list(0,0))
for the case n=2. If I perform the following two steps:
cmd - X -
DearList,
I m trying to draw ONE table that summarize SEVERAL categorical variables
according to one classification variable, say sex. The result would look
like several contingency tables appended one to the other. All the variables
belong to a data frame.
The summary.formula in Hmisc
On Mon, 2007-02-19 at 15:11 -0600, Ranjan Maitra wrote:
The error is different now. It now cannot find Xext library. Do a yum search
on this and install that.
yum provides libXext
which will give you the package Xext which needs to be installed.
You may also need to install the
David LEVY wrote:
DearList,
I ‘m trying to draw ONE table that summarize SEVERAL categorical variables
according to one classification variable, say “sex”. The result would look
like several contingency tables appended one to the other. All the variables
belong to a data frame.
The
On 2/20/2007 9:44 AM, Rick Bilonick wrote:
On Mon, 2007-02-19 at 15:11 -0600, Ranjan Maitra wrote:
The error is different now. It now cannot find Xext library. Do a yum search
on this and install that.
yum provides libXext
which will give you the package Xext which needs to be installed.
On Tue, 20 Feb 2007, Ralf Finne wrote:
Hi everybody!
Is there any way to read a postscrit file into R?
See http://www.r-project.org/useR-2006/Slides/Murrell.pdf
page 4, the grImport package, now on CRAN, with further notes on Paul
Murrell's home page:
http://www.stat.auckland.ac.nz/~paul/
On 20-Feb-07 Ralf Finne wrote:
Hi everybody!
Is there any way to read a postscrit file into R?
All the best to you
Ralf Finne
SYH University of Applied Sciences
Vasa Finland
Well, yes ... since a PostScript file is ASCII text you could
use readline() ... but what you'd do with it after
At 14:41 20/02/2007, you wrote:
Please do not just reply to me,
1 - I might not know
2 - it breaks the threading
Hi
here there is an example extracted from polr help in MASS:
The function could be:
temp - function(form, dat) {
house.plr - polr(formula =
form,
Hi!
I would like to use an existing R-code in an .net project. Now I realize
that there exists a DCOM interface which enables me to use the Statconnector
object in my c# project. That way I can get the results of terms using the
.evaluate() method.
My question is: Is it somehow possible to
Hi all,
The problem is that when you try to use the function summary of a polr object
in a function, it does not work.
The problem is not related to the formula or the structure of data involved.
It is probably related to the use of the function vcov in the code of summary
for polr, and the
That's not it (the function is 'coef' not 'coeff', and R can tell
functions and lists apart).
If you read the help page for polr you will see you could have used
Hess=TRUE. It works then. THAT is why we needed an example, to see how
you used the function.
On Tue, 20 Feb 2007, Michael
I want to calculate the probability that a group will include a particular
point using the squared Mahalanobis distance to the centroid. I understand
that the squared Mahalanobis distance is distributed as chi-squared but that
for a small number of random samples from a multivariate normal
Dear r-helpers,
In basic graphics, I have a figure with x and y axes suppresed. I
would like to move the xlab and the ylab closer to the axes. Do I
have to use mtext()?
_
Professor Michael Kubovy
University of Virginia
Department of Psychology
USPS: P.O.Box
Hello
Instead of testing against 0 i would like to test regression slopes against -1.
Any idea if there's an R script (package?) available.
Thanks for any hint.
Cheers
Lukas
°°°
Lukas Indermaur, PhD student
eawag / Swiss Federal Institute of Aquatic Science and Technology
ECO -
Thank you very much
Paolo
Prof Brian Ripley [EMAIL PROTECTED] 20/02/07 4:00 PM
That's not it (the function is 'coef' not 'coeff', and R can tell
functions and lists apart).
If you read the help page for polr you will see you could have used
Hess=TRUE. It works then. THAT is why we needed
Hello,
I'm using the sem package to do a confirmatory factor analysis on data
collected with a questionnaire. In the model, there is a unique factor G
and 23 items. I would like to calculate the standardized residual
variance of the observed variables. Sem only gives the residual
variance
On Tue, 2007-02-20 at 11:31 -0500, Michael Kubovy wrote:
Dear r-helpers,
In basic graphics, I have a figure with x and y axes suppresed. I
would like to move the xlab and the ylab closer to the axes. Do I
have to use mtext()?
Michael,
You could set the first value in par(mgp) in the
Summarizing:
I'm running R 2.4.1 on a current FC6 32-bit system. In order to have the
rgl R package install, I needed to install both mesa-libGLU-devel (FC6
version is 6.5.1-9) and libXext-devel (FC6) rpm packages. Thanks to
everyone who commented.
Rick B.
two options are to use an offset term or the linear.hypothesis()
function from package car, e.g.,
y - rnorm(100, 2 - 1 * (x - runif(100, -3, 3)), 3)
fit0 - lm(y ~ 1 + offset(-x))
fit1 - lm(y ~ x)
anova(fit0, fit1)
library(car)
linear.hypothesis(fit1, c(x = -1))
I hope it helps.
On Tue, 20 Feb 2007, Federico Calboli wrote:
Liaw, Andy wrote:
I don't see why making copies of the columns you need inside the loop is
better memory management. If the data are in a matrix, accessing
elements is quite fast. If you're worrying about speed of that, do what
Charles
Charles C. Berry wrote:
This is a bit different than your original post (where it appeared that
you were manipulating one row of a matrix at a time), but the issue is
the same.
As suggested in my earlier email this looks like a caching issue, and
this is not peculiar to R.
Viz.
Hi Bernd,
It seems to me that the sharpe function uses the diff to calculate the first
differnce of the input which
is a cumalative return series. This is consistant. see sharpe
Basically you need the rate of returns ((priceFinal - Price
Initial)/priceInitial) which you used. In your case r = 0,
Sorry to answer to myself. A solution was trivial with lattice... (as
often !)
library(lattice)
modbusetmp-lm(IKA_buse ~ Ct *Cc,data=dtbuse)
G1-cloud(IKA_buse~Ct*Cc,type=h,data=dtbuse)
G2-cloud(IKA_buse~Ct*Cc,data=dtbuse)
seqCc-seq(min(dtbuse$Cc),max(dtbuse$Cc),l=20)
Hi Duncan,
I don't know if this will list all the dependencies for your documentation,
since Rick's error messages did not involve libraries and header files already
installed by him for something else, perhaps.
Just a thought.
Ranjan
On Tue, 20 Feb 2007 11:59:24 -0500 Rick Bilonick [EMAIL
Hello all,
I use R on both windows and a mainframe linux installation (RedHat
enterprise 3.0, which they tell me is soon to be upgraded to 4.0). On
windows I installed the package gplots without trouble, and it works fine.
When I attempted to do the same on the unix computer, the following
On 2/20/2007 1:11 PM, Ranjan Maitra wrote:
Hi Duncan,
I don't know if this will list all the dependencies for your documentation,
since Rick's error messages did not involve libraries and header files
already installed by him for something else, perhaps.
No, but it's a start: I'm
install gtools package firstly
Aimin
At 12:17 PM 2/20/2007, Randy Zelick wrote:
gtools'
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
well, it's complaining because you don't have gtools installed.
how about:
install.packages(gplots, dep=T)
?
b
On Feb 20, 2007, at 1:17 PM, Randy Zelick wrote:
Hello all,
I use R on both windows and a mainframe linux installation (RedHat
enterprise 3.0, which they tell me is soon to be
I am pretty pleased with baselines I fit to chromatograms using the
runquantile() function in caTools(v1.6) when its probs parameter is
set to 0.2 and its k parameter to ~1/20th of n (e.g., k ~ 225 for n ~
4500, where n is time series length). This ignores occasional low-
side outliers, and,
On Tue, 2007-02-20 at 10:17 -0800, Randy Zelick wrote:
Hello all,
I use R on both windows and a mainframe linux installation (RedHat
enterprise 3.0, which they tell me is soon to be upgraded to 4.0). On
windows I installed the package gplots without trouble, and it works fine.
When I
I have made some progress with the user defined splitting function and I got
a lot of the things I needed to work. However, I am still stuck on accessing
the node data. It would probably be enough if somebody could tell me, how I
can access the original data frame of the call to rpart.
So if the
On Tue, 20 Feb 2007, Marc Schwartz wrote:
On Tue, 2007-02-20 at 10:17 -0800, Randy Zelick wrote:
Hello all,
I use R on both windows and a mainframe linux installation (RedHat
enterprise 3.0, which they tell me is soon to be upgraded to 4.0). On
windows I installed the package gplots without
On Tue, 20 Feb 2007, Ranjan Maitra wrote:
Hi Duncan,
I don't know if this will list all the dependencies for your
documentation, since Rick's error messages did not involve libraries and
header files already installed by him for something else, perhaps.
It will not, and most of this is
I am not surprised at slowness of runquantile, since it is trying to
perform n=4500 partial sorts of k=225 elements. Here are some thoughts
at speeding it up:
1) playing with different endrule settings can save some time, but
usually results with undesirable effects at first and last 112 values.
Thanks for al the tips, it was the readline() function who did the trick for
me, so thanks for all of your input
Kind Regards
Bart
Original message:
snip
But unfortunately R will do first the calculations and then
afterwards return the strings.
Is there a way around?
On Tue, 20 Feb 2007, Steven Finch wrote:
Hello!
I written several implementations in R of Rémy's algorithm for
generating random ordered strictly binary trees on 2n+1 vertices.
One implementation involves manufacturing character strings like:
X - list(0,list(0,0))
for the case n=2. If I
On Tue, 2007-02-20 at 18:59 +, Prof Brian Ripley wrote:
On Tue, 20 Feb 2007, Marc Schwartz wrote:
On Tue, 2007-02-20 at 10:17 -0800, Randy Zelick wrote:
Hello all,
I use R on both windows and a mainframe linux installation (RedHat
enterprise 3.0, which they tell me is soon to be
Hi Mark,
thanks for your email.
I used your formula for cumul. returns and plugged them into sharpe:
mysharpe - function(x){
+ return(sharpe(cret(x), r=0, scale=1))
+ }
whereby cret is my cumul. returns function as defined by:
cret
function(x){
cumprod(diff(log(x))+1)-1
}
For the index
Duncan Murdoch schrieb:
On 2/20/2007 1:11 PM, Ranjan Maitra wrote:
Hi Duncan,
I don't know if this will list all the dependencies for your documentation,
since Rick's error messages did not involve libraries and header files
already installed by him for something else, perhaps.
No, but
The problem is that rgl is apparently written for GNU make, and has (as
shipped)
ifdef MAKINGAGL
[EMAIL PROTECTED]@ -Iext
[EMAIL PROTECTED]@
else
PKG_CPPFLAGS= -If:/R/R-2.4.1/src/extra/zlib -DHAVE_PNG_H
-If:/R/R-2.4.1/src/gnuwin32/bitmap/libpng -Iext
PKG_LIBS=-lgdi32 -lopengl32 -lglu32
Hello,
I'm looking for a package in R that performs, analysis of correlation
matrices: cross-classified by 2 factors.
The orginal reference to this method is by CJ Brien Biometrica (1998)
75(3):469-76.
THank you,
Agnes
-
E. Agnes Richards, Ph.D.
Hello Tom,
the problem is because R has assumed that pop and rep are integers,
not factor levels. Try:
test - read.table(test.txt,header=T)
sapply(test, class)
test$pop - factor(test$pop)
test$rep - factor(test$rep)
then try fitting the models. Also, there has been substantial
discussion
On 2/20/2007 4:45 PM, Prof Brian Ripley wrote:
The problem is that rgl is apparently written for GNU make, and has (as
shipped)
ifdef MAKINGAGL
[EMAIL PROTECTED]@ -Iext
[EMAIL PROTECTED]@
else
PKG_CPPFLAGS= -If:/R/R-2.4.1/src/extra/zlib -DHAVE_PNG_H
In the example R script below, horizontal gray gridlines are drawn at y
coordinates where the points are drawn with the code:
panel.abline(h=y, v=xScale, col.line=gray)
How do I change this so that the horizontal gray gridlines are drawn at y
coordinates where the y labels are drawn? The
Am 20 Feb 2007 um 20:52 hat Rene Braeckman geschrieben:
From: Rene Braeckman [EMAIL PROTECTED]
To: r-help@stat.math.ethz.ch
Date sent: Tue, 20 Feb 2007 20:52:29 -0800
Subject:[R] Different gridlines per panel in xyplot
In the
Thanks for the reply. Trunc would only work when the truncated y values
result in the desired y coordinates for the grid lines as in this simple
data set. Since my real dataset contains many more points, it does not give
a general solution.
Rene
-Original Message-
From: Bernd Weiss
I solved my own problem. Suddenly remembered the list of very useful
functions under Accessing Auxiliary Information During Plotting in the
help pages.
Here is the line that did the trick:
panel.abline(h=yScale[[panel.number()]], v=xScale, col.line=gray)
Rene
-Original Message-
From:
Hi,
Kindly let me know if I posted a wrong question in the forum.
I want to draw a splom plot with different symbols in plot. My command is as
follows:
splom(~ log10(splomData[2:3]), groups = programs, data = splomData,
panel = panel.superpose,
key = list(title = paste(splomLoop,Programs of Hog
67 matches
Mail list logo