Hi,
Is there any way to use R to present t test results for three groups of
experiments, each of which involves several parallel experiment series with
groups of control vs treated. I would like to present the average fold change
of the experimental parameter (concentration of enzymes) as bars
>
> Is there an easy way to integrate the generalized gamma distribution into
> the current R survival package?
This distribution is available in my gnlr3 function in my gnlm
library. It handles left-, right-, and interval censoring, with linear
and nonlinear regression on all three parameters.
A
Hello R-user
I want to compute a multiple regression but I would to include a check for
collinearity of the variables. Therefore I would like to use a ridge
regression.
I tried lm.ridge() but I don't know yet how to get p-values (single Pr() and p
of the whole model) out of this model. Can any
You are trying to load woa.obj, and you are on Windows which should be
looking for a .dll and the error message says woa.so. You say you keep
your .exe in the rw1070 directory.
You are rather confused! dyn.load on Windows loads a DLL file from the
current working directory, *as its help says*.
I will appreciate you explain me the reason for sending me this e-mail, and
the purpose of the attached file (aplication.pi).
Best regards
Leon Villan
At 11:11 a.m. 05/06/2003 +0200, you wrote:
Please see the attached file.
__
[EMAIL PROTECTED] mailing
Dear all,
Please take my apologies if that has already been asked
- at least I couldn't find it in the archives.
When trying to specify a layout within library lattice,
i.e. using xyplot, I get an error when the prepanel
function tries to subscript the automatically generated
x.limits. This seem
Martin Biuw wrote:
Hi,
Sorry if this is an obvious one, but is there a simple way to modify the lm
function to test whether a slope coefficient is significantly different
from 1 instead of different from 0?
Thanks,
There might be an easier way, but the brute force method would be:
R> set.see
No modification needed. Fit either of
lm(y-x ~ x + z)
lm(y ~ x + z + offset(x))
and the t-test in the summary will be a test of the coefficient of x being
one.
You can also use such models to do an anova against a modle with unit
coefficient.
On Thu, 5 Jun 2003, Martin Biuw wrote:
> Sorry i
Dedi,
There is a very useful package called "stable" by Philippe Lambert and Jim
Lindsey at
http://alpha.luc.ac.be/~jlindsey/rcode.html
Look down that page to the "Probability functions and generalized regression
models for stable distributions" identifier. And if you have not read it
already, I
Hi,
Sorry if this is an obvious one, but is there a simple way to modify the lm
function to test whether a slope coefficient is significantly different
from 1 instead of different from 0?
Thanks,
Martin
Martin Biuw
SEA MAMMAL RESEARCH UNIT
Gatty Marine Laboratory
School of Environmental and E
We have published the 2003/1 issue of R News on
http://cran.R-project.org/doc/Rnews
where you can download the newsletter as PDF or Postscript file. It
will propagate to the CRAN mirrors within a day or two.
Contents of this issue:
Name Space Management for R
Converting Packages to S4
> "Dedi" == Dedi Rosadi <[EMAIL PROTECTED]>
> on Thu, 05 Jun 2003 13:53:35 +0200 writes:
Dedi> I am wondering whether anyone of you already
Dedi> implemented in R or S+ the methods of estimation the
Dedi> parameters of stable distribution. If there is, i will
Dedi> be h
On Wednesday 04 June 2003 16:44, [EMAIL PROTECTED] wrote:
> Hello all-
> I have two time series, Index1stdiff and Comps1stdiff. I regressed
> the first on the second and R returned the summary stats I expected.
> Then I looked at and plotted the residuals. I then wanted to
> assess autoco
Hi,
I am wondering whether anyone of you already implemented in R or S+ the
methods of estimation the parameters of stable distribution. If there is, i
will be happy to get the copy of the codes.
Regards
Dedi
__
[EMAIL PROTECTED] mailing list
https://
Hi
I'm doing a xyplot and I wand to reduce the number of tick marks in the
x axis. My x axis are month and I want to reduce the 12 tick marks to 4.
I used the scales argument but it doesn't seem to work, althougth it
works on y axis if I use scales=list(tick.number=4).
xyplot(land~mes|porto+arte,
Sir/Madam,
I am a graduate student in Petroleum Engineering. I have a problem in using
R and hope that you could help in solving that. I have a code in C, which
displays orthogonal arrays when I give number of rows and number of columns
to it. I have another code in R that calls this C code. I
Hello, I'm using R 1.6.2's survival2.9-6 package, and I just thought
I'd note that it took me a long time to realize that the
coxph.detail()$hazard was a centered hazard (centered at the exponentiated
average covariate times beta), not a baseline hazard. The documentation
simply says "the
>-Original Message-
>From: [EMAIL PROTECTED]
>[mailto:[EMAIL PROTECTED] On Behalf Of james taylor
>Sent: Friday, June 06, 2003 9:38 AM
>To: [EMAIL PROTECTED]
>Subject: [R] Moran's index
>
>
>Hi,
>
>I am looking for some code for Moran's I. Has anyone
>previously done this?
> I have been
> From: Uwe Ligges [mailto:[EMAIL PROTECTED]
>
> [EMAIL PROTECTED] wrote:
>
> > Hello R lovers
> > I have written a little cute function to count the number
> of missing value
> > per row in a matrix and return the percentage of missing value
> >
> > it takes a lot of time to run with a 1000 ro
Just because you get the same answer from different starting points doesn't
mean the algorithm isn't using the starting points you specified.
I tried:
> set.seed(1)
> x <- matrix(rnorm(12), 6, 2)
> kmeans(x, x[c(1,6),], 1)
$cluster
[1] 2 1 2 1 1 2
$centers
[,1] [,2]
1 0.7028106 0.6
Hi,
I am looking for some code for Moran's I. Has anyone previously done this?
I have been unable to find it in the search engines.
James
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
I'm still not completely sure what you are after, but if indeed you want to
change the color of the rectangle borders, how about this ?:
x <- rnorm(100)
y <- rnorm(100)
a <- factor(sample(1:4, 100, rep = T))
lset(list(axis.line = list(col = "grey")))
xyplot(y ~ x | a)
(Make sure you have a re
It is known that viruses now often forge who sent them, so the message
this warning is about may or may not have actually been sent by your system.
==
If you mean the rectangular boxes surrounding the panels and strips, they
are drawn by grid.rect() and not under the control of lattice.
If so, it is *not* a grid but a series of boxes, and they don't always
align perfectly, as example(xyplot) shows on my screen (and if you
had some space between
I am not using panel.grid directly. It appears to use the reference.line
component of the lattice.theme. Since I have tried changing all obvious
comonents in the lattice.theme (including reference.line) I assume I am
correct in this assertion.
I probably used the wrong word, but I wasn't sure what
many Thanks!
For another lazzy typist's - this works fine.
for (i in 1:length(tList)) {
post(tList[[i]],paste("Tree", i, "ps", sep="."), title="title",
digits=getOption("digits") - 0,use.n=TRUE)
graphics.off()
}
christian
> Yes, you need dev.off() inside the loop and to name the files diff
On Wed, 4 Jun 2003, array chip wrote:
> Hi, is it possible to use the packages from CRAN in
> SPlus? and how to do it if yes?
It's going to be possible only if the packages use the common subset of
the S language implemented in both R and S-PLUS: and the porter (you) will
have to know how to ch
[EMAIL PROTECTED] wrote:
I want upload two packages to CRAN, but I can't seem to connect to
ftp://cran.r-project.org/incoming. Do I need a username and password?
1) I can connect.
2) There's anonymous login on CRAN/incoming.
Uwe Ligges
Thanks.
Benjamin Stabler
Transportation Planning Analysis
I'm not sure of the correct term, but in essence it is the equivalent to the
frame in the basic plots or might be described as the bounding box to each
plot. That is when you look at 16 plots on a page these frames (boxes) make
up a grid.
-Original Message-
> Secondly, from a more persona
I want upload two packages to CRAN, but I can't seem to connect to
ftp://cran.r-project.org/incoming. Do I need a username and password?
Thanks.
Benjamin Stabler
Transportation Planning Analysis Unit
Oregon Department of Transportation
555 13th Street NE, Suite 2
Salem, OR 97301 Ph: 503-986-4
Hi, is it possible to use the packages from CRAN in
SPlus? and how to do it if yes?
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Perhaps ?optim or ?optimize will suffice.
George W. Gilchrist wrote:
In S-Plus, there is a function ms() [minimum sums] that estimates parameters
in a multivariate model and returns, along with the estimates, a
loglikelihood score. It's very handy, but I have not found an R version yet.
Does anyon
In S-Plus, there is a function ms() [minimum sums] that estimates parameters
in a multivariate model and returns, along with the estimates, a
loglikelihood score. It's very handy, but I have not found an R version yet.
Does anyone know if such a thing exists? Thanks.
Cheers, George
==
Thanks for your reply.
I am using logistic regression because my response variable is
categorical - and this seems to be recommended in the literature (by
Heckman, Smith and others).
The response variable is named sample because that is what it is - I am
new to R so haven't quite got into habi
Please contact your system administrator.
The scanned document was QUARANTINED.
Virus Information:
The attachment documents.pif contained the virus [EMAIL PROTECTED] and could
NOT be repaired.
-- A . S . C . A . P.--
This message, including any attachments
On Wednesday 04 June 2003 13:48, Liaw, Andy wrote:
> > From: Deepayan Sarkar [mailto:[EMAIL PROTECTED]
> >
> > On Wednesday 04 June 2003 04:26, Mulholland, Tom wrote:
> > > I am probably missing something quite obvious, but any help would be
> > > appreciated. I am continually getting people misrea
> From: Deepayan Sarkar [mailto:[EMAIL PROTECTED]
>
> On Wednesday 04 June 2003 04:26, Mulholland, Tom wrote:
> > I am probably missing something quite obvious, but any help would be
> > appreciated. I am continually getting people misreading the
> lattice plots
> > because they are expecting the
On Wed, 04 Jun 2003 10:40:29 -0700, you wrote in message
<[EMAIL PROTECTED]>:
>Dear R-Help,
>
>I've had some problems installing the R 1.7.0 binary on Windows XP
>Professional. If I launch R from the Start Menu or the Desktop, R crashes.
>However, if I launch by double-clicking on an .RData file
Sounds as if you have a corrupt .RData file in your default working
directory for R. Try renaming it, or otherwise starting with --vanilla.
On Wed, 4 Jun 2003, Simon Frost wrote:
> I've had some problems installing the R 1.7.0 binary on Windows XP
> Professional.
*Installing* it? You only d
On Wednesday 04 June 2003 04:26, Mulholland, Tom wrote:
> I am probably missing something quite obvious, but any help would be
> appreciated. I am continually getting people misreading the lattice plots
> because they are expecting the strip (with the factor names in them) to be
> below the graph.
Dear R-Help,
I've had some problems installing the R 1.7.0 binary on Windows XP
Professional. If I launch R from the Start Menu or the Desktop, R crashes.
However, if I launch by double-clicking on an .RData file, R loads up just
fine. I've never had any problems with previous versions/OS. Has
scatterplot3d in package scatterplot3d
cloud in package lattice
On Wednesday 04 June 2003 10:31, Michele Grassi wrote:
> how can i draw a 3d scatter plot? thank you.
>
> __
> [EMAIL PROTECTED] mailing list
> https://www.stat.math.ethz.ch/mailman/listinf
Michele
get the scatterplot3d package by Uwe Ligges off CRAN:
scatterplot3d: 3D Scatter Plot
Plots a three dimensional (3D) point cloud.
Version:0.3-13
Depends:R (>= 1.1.0)
Date: 2003-01-16
etc...
Gav
Michele Grassi wrote:
how can i draw a 3d scatter plot? than
Michele,
Did you use the search tools in R?
Within R:
> apropos("cooks")
[1] "cooks.distance" "cooks.distance.glm" "cooks.distance.lm"
and a similar search using the html help facility returned a single
entry for the functions influence.measures() in base, under which the
above three funct
Thanks Spencer!
That's what I was looking for. My data are not at all normal and neither are
the transformed values because The majority are at an upper limit that is not
1 (these are grouped data abstracted from a paper). There's nothing that can
be done about that, but it was good to compare
Michele Grassi <[EMAIL PROTECTED]> writes:
> where is the cook-distance in R? i can't find it!please
> help me!
> apropos("cook")
[1] "cooks.distance" "cooks.distance.glm" "cooks.distance.lm"
--
O__ Peter Dalgaard Blegdamsvej 3
c/ /'_ --- Dept. of Biostatistics
On Wed, 2003-06-04 at 16:46, Michele Grassi wrote:
> where is the cook-distance in R? i can't find it!please
> help me!
>
> __
> [EMAIL PROTECTED] mailing list
> https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Hi
Try
?cooks.distance
it's a pa
At 10:43 2003-06-04 -0300, you wrote:
David A. Paul wrote:
> I don't have the reference, but a biologist friend of mine once
> showed me a refereed journal article that purported to demonstrate
> numerical errors made by MSExcel. This would have been Excel97 or
> Excel2000... In any case, the jou
Hi: Below you find two articles about Excel.
Knusel, L. (1998) On the Accuracy of Statistical Distributions in
Microsoft Excel 97. Computational Statistics & Data Analysis. Vol. 26, No.
3, pp. 375-377
McCullough, B.D. y Wilson, D. (1999) On the Accuracy of Statistical
Procedures in Microsoft Exce
where is the cook-distance in R? i can't find it!please
help me!
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
1. If none of the numbers are 0 or 1, I might try a logit
transformation log(p/(1-p)). Then I'd make a normal probability plot of
the transformed variables to check the transformation. If that seemed
OK., then I'd do the computations on logit space and back transform the
result.
2. I
how can i draw a 3d scatter plot? thank you.
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
On Wed, 4 Jun 2003, Prof Brian Ripley wrote:
> See the FAQ, Q7.12.
>
> On Wed, 4 Jun 2003, Philipp Pagel wrote:
>
> > Every once in a while I need to convert a factor to a vector of numeric
> > values. as.numeric(myfactor) of course returns a nice numeric vector of
> > the indexes of the levels wh
Christian -
I am looking at the call to post() in the code fragment below,
and I wonder whether all of the arguments that are supplied will
be matched correctly. I am also reading help("post.rpart").
The first argument is the tree, and that should be matched okay.
The second argument in the
Hello,
I need to get a point estimate and SD for a proportion, but the subjects' data
are not binary---they are proportions (of doses received). That is, I have a
proportion for each subject. In the past I have analysed these data as a
continuous (normal) variable, but I really don't want CIs o
Hej,
just about excel -- excel 97 includ a kind of flying game (or what
ever it should be called). There was a landscape you could fly above,
and in the middle of it there were a screen with different names
(authors?). You had to put the cursor on a particular cell and click
a particular icon, A
--On woensdag 4 juni 2003 15:21 +0200 Christoph Lehmann
wrote:
> what is the easiest way to get from
> x
> x1 x2
> [1,] 2 3
> [2,] 3 2
> [3,] 1 3
> [4,] 1 4
>
> xbar1
>x1 x2
> [1,] 1.75 3
> [2,] 1.75 3
> [3,] 1.75 3
> [4,] 1.75 3
>
> with the mean of the columns of x as
Dear all,
I've now spent a couple of days trying to learn R and, in particular, the
gam() function, and I now have a few questions and reflections regarding
the latter. Maybe these things are implemented in some way that I'm not yet
aware of or have perhaps been decided by the R community to no
I need to create graph composed of different plots (bar + area). See
attachment, I am sending an example of intended design.
Is it possible to plot one set as fitted line and add another one in the
column shape, as it is shown on the example? It is intended because of
hydrological convence where to
Thee are no attached files.
- Original Message -
From: [EMAIL PROTECTED]
To: [EMAIL PROTECTED]
Sent: Monday, June 02, 2003 11:23 AM
Subject: Re: Your application
Please see the attached file.
[[alternate HTML version deleted]]
_
Hello all-
I have two time series, Index1stdiff and Comps1stdiff. I regressed the
first on the second and R returned the summary stats I expected. Then I
looked at and plotted the residuals. I then wanted to assess
autocorrelation characteristics and tried to run a Durbin-Watson using:
Paul Lemmens <[EMAIL PROTECTED]> writes:
> Hoi Christoph,
>
> --On woensdag 4 juni 2003 15:21 +0200 Christoph Lehmann
> <[EMAIL PROTECTED]> wrote:
>
> > what is the easiest way to get from
> > x
> > x1 x2
> > [1,] 2 3
> > [2,] 3 2
> > [3,] 1 3
> > [4,] 1 4
> >
> > xbar1
> >x
David A. Paul wrote:
> I don't have the reference, but a biologist friend of mine once
> showed me a refereed journal article that purported to demonstrate
> numerical errors made by MSExcel. This would have been Excel97 or
> Excel2000... In any case, the journal's scope was biological in
> natur
On Wed, 2003-06-04 at 08:09, Paul, David A wrote:
> I don't have the reference, but a biologist friend of mine
> once showed me a refereed journal article that purported
> to demonstrate numerical errors made by MSExcel. This
> would have been Excel97 or Excel2000... In any case, the
> journal's
Philipp Pagel <[EMAIL PROTECTED]> writes:
> But I just can't believe there is no R function to do this in a more
> readable way. Actually, the behaviour of as.numeric() doesn't strike me
> as very intuitive. I'm sure it has been implemented that way for a
> reason - but what is it?
One reason is
65 matches
Mail list logo