Check out this thread:
https://stat.ethz.ch/pipermail/r-help/2008-February/153638.html
On Feb 19, 2008 10:47 PM, Rolf Turner <[EMAIL PROTECTED]> wrote:
>
> For reasons best known only to myself ( :-) ) I wish to create a data
> frame with 0 rows and 9 columns.
>
> The best I've been able to come
On 2/14/08, Deepayan Sarkar <[EMAIL PROTECTED]> wrote:
> On 2/14/08, Felix Andrews <[EMAIL PROTECTED]> wrote:
> > You can tell Lattice to stop when an error occurs, like this:
> >
> > lattice.options(panel.error="stop")
> > xyplot(1:10 ~ 1:10, panel=function(...) stop("foo"))
> > # -> Error in p
G'day Sean,
On Fri, 15 Feb 2008 09:12:22 +0800
"Hung-Hsuan Chen (Sean)" <[EMAIL PROTECTED]> wrote:
> Assume I have 3 distributions, x1, x2, and x3.
> x1 ~ normal(mu1, sd1)
> x2 ~ normal(mu2, sd2)
> x3 ~ normal(mu3, sd3)
> y1 = x1 + x2
> y2 = x1 + x3
>
> Now that the data I can observed is only y
check out the 'approx' function.
On Feb 19, 2008 12:44 PM, Dani Valverde <[EMAIL PROTECTED]> wrote:
> Hello,
> I have two vectors, one with 13112 points and the other one with 10909.
> I wonder if there is a way to interpolate the data so the shorter
> vectors has the same number of points as the
G'day Juliet,
On Tue, 19 Feb 2008 21:35:11 -0500
"Juliet Hannah" <[EMAIL PROTECTED]> wrote:
> I implemented a try() statement that looks like:
>
> <- function(index)
> {
>
>reduced_model <- try(glm.fit(X4,n,family=poisson(link="log")))
>full_model <- try(glm.fit(X5,n,family=poisson(lin
A one-day workshop on the R programming environment
has been organised by the Centre for Mathematics and its
Applications (Australian National University) and
Geoscience Australia in Canberra on Monday 3-March 2008.
Registration fee is AUD$150 (students: $100).
The workshop will be conducted by
For reasons best known only to myself ( :-) ) I wish to create a data
frame with 0 rows and 9 columns.
The best I've been able to come up with is:
junk <- as.data.frame(matrix(0,nrow=0,ncol=9))
Is there a sexier way?
cheers,
Rolf
##
Hi
I want to do use AIC for model selection on mixed model. However, before
going deeper in the model selection, I want to assess is there is
overdispersion with the full model in order to decide if I should use QAIC
instead of AIC.
Is there a way to compute a chisquare goodness of fit test for
Dear R Users,
I implemented a try() statement that looks like:
<- function(index)
{
reduced_model <- try(glm.fit(X4,n,family=poisson(link="log")))
full_model <- try(glm.fit(X5,n,family=poisson(link="log")))
if (inherits(reduced_model,"try-error") ||
inherits(full_model,"try-error")) r
I'm using the monoreg function (with weights) from the fdrtool package.
How can I calculate the R square for this type of regression?
Thanks for your help,
Thierry
--
View this message in context:
http://www.nabble.com/R-square-for-Monotone-regression-tp15580803p15580803.html
Sent from the
Assuming the time values are in ascending order within date:
ts$prices <- ave(ts$prices, dates, FUN = function(x) tail(x, 1))
On Feb 19, 2008 8:41 PM, Edwin Hoyle <[EMAIL PROTECTED]> wrote:
> My code below makes a data frame with columns for
> date, time, and price. Time on each date runs from
My code below makes a data frame with columns for
date, time, and price. Time on each date runs from 1
to 4.
I'd like to add a new column "ts$closingprice", which
would have the closing price for that date. To find
the closing price, I'd like to take the price in the
row having the greatest time
Kenneth
See %in%. In your example:
Data2 <- subset(Data1, !Identity%in%c(1,2,3))
HTH ..
Peter Alspach
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Kenneth Takagi
> Sent: Wednesday, 20 February 2008 2:28 p.m.
> To: r-help@r-project.org
Hello,
I have a dataset which consists of 9 columns (variables) and 35 rows
(observations). I am doing a simple linear regression of one variable
on the other. There are some observations that are outliers and I would
like to remove them based on another variable (it's a unique, numeric
varia
The output from xy.coords is not being sent to Axis since if it were then
Gavin's code earlier in this thread, which does work, would not.
On Feb 19, 2008 4:08 PM, Mark Difford <[EMAIL PROTECTED]> wrote:
>
> Hi Stiffler,
>
> >> I was wondering why the plot() command ignores the datatype when
> >>
Distance matrices are not usually and end in themselves but a means to
some other end. Rather than ask what is the best way to calculate such
a huge distance matrix, maybe the question you should ask yourself is
what are you going to do with it if ever you did manage to calculate it.
Maybe you ca
On Tue, 2008-02-19 at 14:40 -0800, Stiffler wrote:
>
>
> Mark Difford wrote:
> >
> >>> I was wondering why the plot() command ignores the datatype when
> >>> displaying axis labels...
> >
> > plot() doesn't ignore the datatype:
> > [...]
> > plot(x,y) calls xy.coords(), which recasts x as: x =
There were two queries recently regarding removing
rows or columns that have all NAs.
Three respondents suggested combinations of apply() with
any() or all().
I cringe when I see apply() used unnecessarily.
Using rowSums() or colSums() is much faster, and gives more readable
code. (Two responden
?colSums
On Feb 19, 2008, at 6:39 PM, Saurav Pathak wrote:
Hi,
I cannot seem to figure out how to sum over an index of a array.
For example, let A be a 3 dimensional array. I want to, say, find
the sum over the first dimension. That is
S_jk = Sum_i A_ijk
where now S is a 2-dim matr
Hi,
I cannot seem to figure out how to sum over an index of a array.
For example, let A be a 3 dimensional array. I want to, say, find
the sum over the first dimension. That is
S_jk = Sum_i A_ijk
where now S is a 2-dim matrix. I dont want to use a loop.
Thanks,
--
saurav
__
***reading in data**
data<-read.table("microarray.txt",header=T, sep="\t")
head(data)
dim(data)
attach(data)
***creating matrix and calculating variance across probesets
x<-1:2
y<-2:141
data.matrix<-data.matrix(data[,y])
variableprobe<-apply(data.matri
First, a clarification. The subject line suggests that the
Welch procedure is not an ordinary student t-test.
That is incorrect.
There are two common two-sample t-tests:
non-pooled variance (Welch version)
pooled variance
I would refer to the non-pooled version just as a "two-sam
Kathy Gerber wrote:
>
> Earlier today I sent a question to Frank Harrell as an R developer with
> whom I am most familiar. He suggested also that I put my questions to
> the list for additional responses. Next month I'll be giving a talk on
> R as an example of high quality open source soft
Dear R gurus,
To start, let me confess to not being an experienced programmer, although I
have used R fairly
extensively in my work as a
graduate student in statistics.
I wish to find the root of a function of two variables that is defined by an
integral which must be
evaluated numerically.
On 19/02/2008 5:40 PM, Stiffler wrote:
>
>
> Mark Difford wrote:
I was wondering why the plot() command ignores the datatype when
displaying axis labels...
>> plot() doesn't ignore the datatype:
>> [...]
>> plot(x,y) calls xy.coords(), which recasts x as: x = as.double(x), which
>> is f
Mark Difford wrote:
>
>>> I was wondering why the plot() command ignores the datatype when
>>> displaying axis labels...
>
> plot() doesn't ignore the datatype:
> [...]
> plot(x,y) calls xy.coords(), which recasts x as: x = as.double(x), which
> is fine, since x is (also/primarily) numeric.
>
Gavin Simpson wrote:
>
>> PS what's the right way to get integer labels?
>
> Do them by hand, if they are (numeric) integers
>
>> plot(x,y, axes = FALSE)
>> axis(2)
>> axis(1, at = x)
>> box()
>
> You could try writing your own Axis.integer function if doing the extra
> steps is a pain - som
Hi Conny,
It still isn't clear what your question is, but a density plot "simply"
shows you the distribution of your data, say a set of measurements of
something. Think of it as a modern replacement for the histogram.
See
http://en.wikipedia.org/wiki/Density_estimation
for greater insight.
H
This is your original problem:
f <- function(x, y) as.vector(outer(x, y, "/"))
mapply(f, as.data.frame(t(X)), as.data.frame(t(Y)))
so just replace "/" with whatever function of two
variables you like. See ?outer and be sure your
function replacing "/" is vectorizable as noted there.
On Feb 19,
Sorry to disturb,
I managed to plot 2 together with 'layer' like this:
qplot(se, or, min=lcl1, max=ucl1, data=df1, geom="pointrange")+layer(data =
df2, mapping = *aes*(x = se, y = OR2,min=lcl2,max=ucl2), geom =
"pointrange")+geom_line()
My only problem is that It doesn't plot the line (geom_line
Greg,
thank you for the nice overview, which is very helpful.
Today I tested the HTML-Tool (Hmisc-library). First I had to install
the HeVea translater (it´s a little bit tricky). I could produce a little
html-table (it made me happy).
I will be out of office until monday - it would be nice to con
Hi Yianni,
This just proves that you should be using R as your calculator, and not the
other one!
Regards, Mark.
gatemaze wrote:
>
> Hello,
>
> on a simple linear model the values produced from the fitted(model)
> function
> are difference from manually calculating on calc. Will anyone have
Aah, didn't notice this one.I was using the sunspot data and couldnt
think of a way to show it in this example.
Nevertheless, axs="i" worked perfectly.
Thank you
Saptarshi
On Feb 19, 2008, at 4:17 PM, Deepayan Sarkar wrote:
>>
>>In this case there is some padding between the endpoints an
Hi Stiffler,
>> I was wondering why the plot() command ignores the datatype when
>> displaying axis labels...
plot() doesn't ignore the datatype:
> x <- as.integer(c(1,2,3))
> y <-x
> typeof(x)
[1] "integer"
> mode(x)
[1] "numeric"
plot(x,y) calls xy.coords(), which recasts x as: x = as.double
On 2/19/08, Saptarshi Guha <[EMAIL PROTECTED]> wrote:
> Hello,
> I need to plot a line graph using lattice so that the endpoints are
> flush with the left and right vertical axes.
> So the scale end points and horizontal axes length should coincide.
> This works
>
>
I'm trying to implement a recursive function using integrate, and I
suspect I need a Vectorize somewhere,
but I can't suss it out. Any help would be appreciated. I've tried
traceback() and various debugging ideas to no avail (most likely due to
my inexperience with these tools.)
Here's what I have
Hello Keith,
Based on what you have done, please try this:
M<-matrix(nrow=num.rows, ncol=num.x.col*num.y.col)
for( i in 1:num.rows){M[i,]<-rep(Y[i,], each=num.y.col)}
Z<-Z+M
not elegant, but works.
HTH
YH Deng
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Beha
Sorry,
I want to plot 2 lines (or pointrange like I used before) together, with
confidence limits of the estimates.
I want to plot another curve as I did for the first one (with qplot like you
suggested).
I know how to do it with ggplot *without* confidence limits:
SE<-c(0.6,0.7,0.8,0.9)
OR<-c(2.
"Kathy Gerber" <[EMAIL PROTECTED]> wrote in message
news:[EMAIL PROTECTED]
> Spencer,
>
> I believe this is the first mention of pricing that I've seen.
> Several additional points have been made about the comparison of R to
> Octave, some off list.
> -- Matlab did not alienate developers all th
Hi André
>> I'm trying to use bart() ... but I can't predict my test data.
There is no predict() method, as such, but read ?pdbart very carefully, and
work through the examples. The fitted/predicted values from a fit are in
bartobject$yhat.train.mean
It would be a good idea to read the two PD
On 2/19/2008 2:22 PM, Benjamin Zuckerberg wrote:
>
> I am encountering an error when I attempt to reference a glm model
> within a function. The function uses the segmented.glm command
> (package = segmented). Within the segmented.glm command one specifies
> an object, in this case a logis
I am encountering an error when I attempt to reference a glm model
within a function. The function uses the segmented.glm command
(package = segmented). Within the segmented.glm command one specifies
an object, in this case a logistic regression model, and specifies a
starting threshold
On Tuesday 19 February 2008 (19:51:15), TLowe wrote:
> Hey Folks,
>
> Could somebody show me how to loop through a list of dataframes? I want to
> be able to generically access their elements and do something with them.
>
> For instance, instead of this:
>
> df1<- data.frame(x=(1:5),y=(1:5));
> df
I have a problem I cannot run rJava properly, I have a winXP, below is the
error log I hope someone can help me because I need rJava to run other
programs like RLadayBug that I need for my research. Mr Hohle was kind
enough to assist me but I think the problem is with my PC not being able to
run r
Use lchoose and use logarithms throughout.
> x=666
> y=1287
> lchoose(x+y,x)-(x+y)*log(2)
[1] -104.4265
> Pxy = exp(lchoose(x+y,x)-(x+y)*log(2))
[1] 4.447787e-46
Joe
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Hyojin Lee
Sent: Monday, February 18,
Hello,
I need to plot a line graph using lattice so that the endpoints are
flush with the left and right vertical axes.
So the scale end points and horizontal axes length should coincide.
This works
u=newdata[1:100,] #u[,1] is from 1700-1795
xyplot
On 2/19/2008 1:51 PM, TLowe wrote:
> Hey Folks,
>
> Could somebody show me how to loop through a list of dataframes? I want to
> be able to generically access their elements and do something with them.
>
> For instance, instead of this:
>
> df1<- data.frame(x=(1:5),y=(1:5));
> df2<- data.frame(
Thank you ahead of time for help with this.
I have two matrices
X
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 660 693.00 726.0 759.00 792.0 825.00
[2,] 548 575.40 602.8 630.20 657.6 685.00
[3,] 676 709.80 743.6 777.40 811.2 845.00
[4,] 763 801.15 839.3 877.45 915.6 953.75
[5,] 768 80
Hey Folks,
Could somebody show me how to loop through a list of dataframes? I want to
be able to generically access their elements and do something with them.
For instance, instead of this:
df1<- data.frame(x=(1:5),y=(1:5));
df2<- data.frame(x=(1:5),y=(1:5));
df3<- data.frame(x=(1:5),y=(1:5));
Check out:
https://stat.ethz.ch/pipermail/r-help/2006-March/101812.html
On Feb 19, 2008 1:41 PM, Nitin Jain <[EMAIL PROTECTED]> wrote:
> Dear R-help members,
>
> I am using logistic regression on a high throughput data and would like
> to capture warning messages, if generated for a particular pr
On 19-Feb-08 18:09:18, Monica Pisica wrote:
> take a look at:
>
> Du, 2002, Master Thesis, http://www.math.mcmaster.ca/peter/mix/Rmix.pdf
>
> Macdonald, P., 2003, RMIX routine for R,
> http://www.math.mcmaster.ca/peter/mix/mix.html
>
> I don't think this package was actually posted on CRAN (the
Matt, I know you are probably busy with work, but I cannot help buy asking
you these R questions. If it is bothersome, please let me know and I will
stick with the R help... but
I have two matrices
X
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 660 693.00 726.0 759.00 792.0 825.00
[2,] 548
Dear R-help members,
I am using logistic regression on a high throughput data and would like
to capture warning messages, if generated for a particular probe. The
way I am approaching it currently is:
myResult <- data.frame(matrix(NA, nrow = 1000, ncol = 4))
colnames(myResult) <- c("intercep
Hi all,
sorry for my english, but I don't speak yours language.
I'm trying to use bart() and rbf(). The package I'm using now is
"BayesTree" and "neural", respectively. I could create the models, but I
can't predict my test data.
Does anyone have such an experience? Any advice is appreciated
perfect. many thanks.
-Original Message-
From: Gabor Grothendieck [mailto:[EMAIL PROTECTED]
Sent: Tuesday, February 19, 2008 1:26 PM
To: Sung, Iyue
Cc: r-help@r-project.org
Subject: Re: [R] Extracting original variable list from lm object
Try:
all.vars(formula(my.model))
On Feb 19, 2
Try:
all.vars(formula(my.model))
On Feb 19, 2008 1:10 PM, Sung, Iyue <[EMAIL PROTECTED]> wrote:
> Fellow R users,
>
> I have an lm object, from which I would like to extract the list of
> original variables.
> The problem I have is the formula includes functions of the covariates.
>
> I tried us
Fellow R users,
I have an lm object, from which I would like to extract the list of
original variables.
The problem I have is the formula includes functions of the covariates.
I tried using "attr", but the result stores the transformed variable
name. For example:
> my.model<-lm(y ~ a + log(b +
On Tue, 2008-02-19 at 09:17 -0800, Stiffler wrote:
> Hello,
>
> I was wondering why the plot() command ignores the datatype when displaying
> axis labels. More specifically, if the data points are integers then the
> axis labels should intuitively also be integers, right?
>
> > x <- as.integer(c
take a look at:
Du, 2002, Master Thesis, http://www.math.mcmaster.ca/peter/mix/Rmix.pdf
Macdonald, P., 2003, RMIX routine for R,
http://www.math.mcmaster.ca/peter/mix/mix.html
I don't think this package was actually posted on CRAN (the mix package on CRAN
is a different one as far as i reme
I've used a similar three-stage approach. Depending on the format of your
map (mine came from an ArcView shape file) you may be able to specify
col="transparent" for the foreground (ie land) in plot as I did (tip, use
border= to get colours for the country borders). This works where alpha
channels
I don't do "why" answers. Only how. Occasionally.
?plot.default ##with the axes=FALSE argument. Then
?axis ## note the labels and at arguments.
Bert Gunter
Genentech Nonclinical Statistics
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of Stiffler
Sen
Hi Rthoughts,
Yes, I see now that they truly are (just) Rthoughts;) but take courage, for
we are getting closer (to the start). You still need to read the basic
documentation, and you will get used to the command line.
What I think you need is a package called Rcmdr. So, start R using your
des
Hello,
I have two vectors, one with 13112 points and the other one with 10909.
I wonder if there is a way to interpolate the data so the shorter
vectors has the same number of points as the longer one.
Best,
Dani
--
Daniel Valverde Saubí
Grup de Biologia Molecular de Llevats
Facultat de Veteri
The mixreg package may help. Just fit an intercept only model to your
data (may need to have x be a column of 1's and tell it not to include
an intercept) and see if it finds your mixtures.
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL
Hi Kathy,
maybe this article could be also of use for you?
Ihaka, R., and Gentleman, R. (1996)," R: A Language for Data Analysis
and Graphics," The Journal of Computational and Graphical Statistics, 5,
299-314
Best,
Roland
Kathy Gerber wrote:
> Earlier today I sent a question to Frank Harrell
Hello,
I was wondering why the plot() command ignores the datatype when displaying
axis labels. More specifically, if the data points are integers then the
axis labels should intuitively also be integers, right?
> x <- as.integer(c(1,2,3))
> y <-x
> typeof(x)
[1] "integer"
> plot(x,y)
>
The ax
The approach that I usually use for things like this is:
Plot the map to get the aspect ratio and limits correct.
Add the image (obscuring the original map).
Add the map again on top of the image using a light grey color.
This seems to work fine for me. Another aproach that you may try (untested
Hi Rthoughts,
Don't be discouraged I'm learning R on/off for classes and when I was
under pressure to get work done, learning R likes pulling my hair (and
teeth) off of frustration. But, this forum is great; I got so much help
from this forum.
I use R on windows as well. After install R, i
There are several options available to you depending on your knowledge and
workflow.
Others have mentioned using Excel to format the table and to copy that into
word, one thing along those lines that has not been mentioned yet is that if
you have your data in a matrix or data frame then (on win
You might want to look at the datasets package
as an example:
https://svn.r-project.org/R/trunk/src/library/datasets/
The other way is to save() them in an .rda file and then
just load() the file.
On Feb 19, 2008 10:20 AM, Birgit Lemcke <[EMAIL PROTECTED]> wrote:
> I am using the recent R versio
Thanks for your suggestion, Henrique
Here I attach a link to some notes by Robert Bivand, with a SAR model, see
page 30
www.bias-project.org.uk/ASDARcourse/unit6_slides.pdf
The difference with CAR is in the covariance structure, I run and example:
# =
Hi Mark,
Thank you for the reply.
I meant the command prompts to start an R file.
To be followed on by importint data I can then use to practise the software
with. The installation did put an icon on teh desktop. I am a very skilled
user of computers but command lines for many programs is somet
Hi Rthoughts,
It isn't clear what you mean. When you install R, the installation program
usually puts an icon on your desktop that you can click on to run the
program. So, if you don't have that, but have installed R, and what you
mean is, "How do I start the R program?" or "How do I run R?" th
Irene Mantzouni difres.dk> writes:
> I am trying to run a regression where the predictor values are not real
> data but each is estimated from a different model. So, for each value I
> have a mean and variance.
>
> Which package/function should I use in this case?
See the Dobson example in th
On 2/19/08, Tom Cohen <[EMAIL PROTECTED]> wrote:
> Dear list,
>
> With the below codes, I got 8 bwplots but I would like to put 2 bwplots in
> one single graph so that instead of 8 separate bwplots I would have 4 graphs,
> each contains 2 bwplots. How can I do that?
>
> Another question is h
Thank you very much for your help - that fixed it. Sorry I did not see
it in the manual.
Regards
JS
---
-Original Message-
From: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
Sent: 19 February 2008 14:14
To: john seers (IFR)
Cc: r-help@r-project.org
Subject: Re: [R] Building a pa
On 2/19/08, Tom Cohen <[EMAIL PROTECTED]> wrote:
> Dear list,
>
> I have following plot, where I have set the color (red and green) and lines
> (lty=2:3) in the panel.groups but can't not figure out how change the lines
> and color of the legend in the "key" to the same lines and color as in th
Will Holcomb wrote:
> I have been attempting some basic power calculations using R and I am not
> getting the results I expect. I had a homework assignment in SAS, but I want
> to learn R as well, so I was attempting to reproduce my result. (No one else
> in the class is doing R, so there's no need
Hi Mark,
Thank you for your reply. There is one link I haven't come across, the last
one. I have seen them but I couldn't find where 'how to start R' is
explained for Windows platforms.
I will look further into them.
As for everyone else who sent e-mails, thank you. I have printed them out
and
Will,
Your SAS input indicates that within standard deviation is 9, not the
variance. If you use within.var=81 in your R statement you will get the
answer matching SAS.
Cheers,
Andy
__
Andy Jaworski
518-1-01
Process Laboratory
3M Corporate Research Laboratory
--
There is no plotDensity() function in 'limma', but plotDensities().
FYI, there are 1000+ CRAN packages, several hundred Bioconductor
packages, and probably another 1000 packages elsewhere, so please be
precise in order to avoid confusion (and time waste). Take a minute
or two to read the followin
Maybe one option should be:
ecdf2 <- function (x)
{
x <- sort(x, decreasing = TRUE) # Changed Line
n <- length(x)
if (n < 1)
stop("'x' must have 1 or more non-missing values")
vals <- unique(x)
rval <- approxfun(vals, cumsum(tabulate(match(x, vals)))/n,
method
Have you looked at Sweave? It allows you to put all your commands into
a template file along with details on where you want graphs and any
other markup. You can also include other text (notes, explanations,
etc) that will show up in the output, but not be processed by R. The
original version of
I have been attempting some basic power calculations using R and I am not
getting the results I expect. I had a homework assignment in SAS, but I want
to learn R as well, so I was attempting to reproduce my result. (No one else
in the class is doing R, so there's no need to obsfucate the answer, th
Hi, if I may further add it would be very helpful if you get a book on R
from your library... the introductory chapters should be very much helpful
as all of them start on how to simply start R, load your data files, etc
etc The deepness and power of R is well beyond any other software, but
yo
Hi, all
ecdf function (Empirical Cumulative Distribution Function) in "stats"
package counts from smaller values to larger values.
However, I want to draw it by counting from larger value to smaller values
and I couldn't find options for this purpose.
How can I draw ecdf or ecdf like graph by co
Hi,
well if there is no binary it seems you have to recompile from source...
>From the R FAQ:
2.5.1 How can R be installed (Unix)
If R is already installed, it can be started by typing R at the shell prompt
(of course, provided that the executable is in your path).
If binaries are available for
> Another question: Can I add another line to the same plot with qplot?
> (like function "lines" in "plot").
Yes.
(if you want more details, tell us what you're trying to do!)
Hadley
--
http://had.co.nz/
__
R-help@r-project.org mailing list
https:/
Thank you all for your help. Apologies for not giving an example.
model.matrix was useful as comparing that table with the one from the
spreadsheet showed the "mispell".
On 19/02/2008, Douglas Bates <[EMAIL PROTECTED]> wrote:
>
> On Feb 19, 2008 8:41 AM, <[EMAIL PROTECTED]> wrote
>
> > on a simpl
Hi Rthoughts,
>> I am currently discouraged by the use of r. I cannot figure out how to
>> use it despite
>> extensive searches. Can anyone help me with getting started? How can
>> import
>> a txt file with series...
There are piles of documents that you could (and should) read. I am
surprised
Hi,
Marc Bernard wrote:
> Dear all,
>
> I am looking for a good reference on "Survival analysis". I am looking for
> a booking containing both applications and Maths. Explaining different
> methods in survival analysis
>
since I just answered in another thread with a book recommen
Hello all,
I am working at the FORTH institute in Crete and it's been a long now that I
am trying to reproduce the results of the paper :
"Gene expression profiling predits clinical outcome of breast cancer", by
Van't Veer et al. It has been published in NATURE, vol 415,
31 January 2002.
http://ww
Hallo,
I just loaded the limma package.
Conny
Original-Nachricht
> Datum: Tue, 19 Feb 2008 06:41:18 -0800
> Von: "Henrik Bengtsson" <[EMAIL PROTECTED]>
> An: [EMAIL PROTECTED]
> CC: [EMAIL PROTECTED]
> Betreff: Re: [R] plotDensity
> Is 'plotDensity' a specific function you are
I am using the recent R version on a G4 PowerBook with Mac OS X 10.4.11.
I have a bunch of datasets that I would like to put all together in a
R library, so that I only have to type:
>library("name")
and all datasets are loaded. Is this like building a package or is
there an easier way
Howdee,
I am able to fit a 4-parameter logistic growth curve to a dataset which
comprise many individuals (using R v. 2.3.1). Yet, if I want to obtain the
parameters for each individual (i.e., for each 'id') using nlsList, then I
obtain an Error message which I have trouble interpreting. Any advic
Hi,
Xing Yuan wrote:
> Dear List,
>
> Does anybody no how to compare mean survival times for two (more) groups in
> R? What test statistics should I use?
my answer is less of an R answer than a literature answer:
John P. Klein and Melvin L. Moeschberger devote section 4.5 in their
book "Surviv
On Tue, Feb 19, 2008 at 10:04:13AM -0500, Duncan Murdoch wrote:
> On 2/19/2008 9:24 AM, Hans Ekbrand wrote:
[...]
> > I tried the following small code snippet which I copied from the
> > "Introduction to R":
> >
> >> for (i in 2:length(meriter)) { table(meriter[[1]], meriter[[i]]) }
>
> Where d
On Tue, Feb 19, 2008 at 04:52:19PM +0200, K. Elo wrote:
> Hi,
>
> Hans Ekbrand wrote (19.2.2008):
> > I tried the following small code snippet which I copied from the
> >
> > "Introduction to R":
> > > for (i in 2:length(meriter)) { table(meriter[[1]], meriter[[i]]) }
>
> Try:
> for (i in 2:lengt
On 2/19/2008 9:24 AM, Hans Ekbrand wrote:
> Hi list
>
> I have a data frame I would like to loop over. To begin with I would
> like crosstabulations using the first variabel in the data frame,
> which is called "meriter".
>
>> table(meriter[[1]], meriter[[3]])
>
On Feb 19, 2008 8:41 AM, <[EMAIL PROTECTED]> wrote
> on a simple linear model the values produced from the fitted(model) function
> are difference from manually calculating on calc. Will anyone have a clue...
> or any insights on how fitted function calculates the values? Thank you.
___
1 - 100 of 149 matches
Mail list logo
