Is this a *warning* from package 'foreign'? If so it is discussed on the
relevant help page (see ?read.spss).
Please note the footer of this message, and the request for 'commented,
minimal, self-contained, reproducible code' as well as the posting guide's
request for you to do your homework a
Henrique Dallazuanna wrote:
> Try this:
>
> plot(rnorm(10), main=expression(paste("Mus., 10 ",mu," g", sep="")))
or
plot(rnorm(10), main=expression("Mus., 10 " * mu * g))
without all that paste()ing, just as one formula.
Uwe Ligges
> On 18/02/2008, Samor Gandhi <[EMAIL PROTECTED]> wrote:
>>
I am trying to calculate the distance metrics, but i keep getting this error.
I am very new to R. Am i doing something wrong.
> dim(data.sub)
[1] 1 140
> dist(data.sub, method = "euclidean", diag = FALSE, upper = FALSE, p = 2)
Error: cannot allocate vector of size 781250 Kb
In addit
On 2/16/08, hadley wickham <[EMAIL PROTECTED]> wrote:
>
> On Feb 16, 2008 9:11 AM, sigalit mangut-leiba <[EMAIL PROTECTED]> wrote:
> > Hello,
> > This is the first time i'm trying to plot in R. I want to plot estimates
> of
> > OR and their confidence limits, like a scatter plot:
> > the vertical a
You could have a look at package "epitools" (on CRAN, or at www.epitools.net).
-Heinrich.
> -Ursprüngliche Nachricht-
> Von: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] Im Auftrag von stefan lhachimi
> Gesendet: Dienstag, 19. Februar 2008 10:45
> An: [EMAIL PROTECTED]
> Betreff: [R] Co
Hello,
I am looking for a function which allows to calculate the confidence
interval for a standard mortality ratio. I do have vectors with the
number of observed and expected death. Has anybody a hint where to
look?
Best,
Stefan
__
R-help@r-project.org
-- Forwarded message --
From: sigalit mangut-leiba <[EMAIL PROTECTED]>
Date: Feb 19, 2008 12:17 PM
Subject: Re: [R] plotEst
To: r-help <[EMAIL PROTECTED]>
On 2/16/08, hadley wickham <[EMAIL PROTECTED]> wrote:
>
> On Feb 16, 2008 9:11 AM, sigalit mangut-leiba <[EMAIL PROTECTED]> w
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 the
panel.groups.
bwplot(means ~ age | scales , d
On Tuesday 19 February 2008 (12:33:21), Alfonso Pérez Rodríguez wrote:
> Hello, I'm sure that this question is too simple, but, I'm begining with R
> and I'm not able to change the NA values in a matrix by 0 values, and it's
> necessary for my work, how can I do It? Thank you.
your.matrix[is.na(yo
Hello,
I have a question about "qplot": How can I add another line to the same
plot ?
(like function "lines" in "plot").
Thank You,
Sigalit.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/l
Format your dataframe into a long format. Then use a grouping variable
to distinguish both lines.
library(ggplot2)
DF <- data.frame(X = rep(0:20, 2), Y = c(rnorm(21), runif(21)), Z =
gl(2, 21))
ggplot(data = DF, aes(x = X, y = Y, colour = Z)) + geom_line()
HTH,
Thierry
PS. I suggest you read th
Try this:
name1 <- 'sample-plot'
filename <- sprintf("%s.png", name1)
inputdir <- '/path/to/dir'
paste(inputdir, filename, sep="")
On 19/02/2008, Hyunchul Kim <[EMAIL PROTECTED]> wrote:
> Hi, all
>
> How to format and join strings ?
> For example, like following short python examples.
>
> **
Hi, all
How to format and join strings ?
For example, like following short python examples.
*
name1 = 'sample-plot'
filename = '%s.png' % name1
inputdir = '/path/to/dir'
os.path.join(inputdir, filename)
**
Best,
Hyunchul Kim
[[alternative HTML version deleted]]
__
I've figured it out by repeatedly testing. It is to use a type='term'
statement, just as used in gam.
sorry to bother.
On 2/19/08, gallon li <[EMAIL PROTECTED]> wrote:
>
> Thanks a lot, Prof Lumley.
>
> Now I can fit a model like
>
> coxfit=coxpy((time,censor)~pspline(x1)+x2+x3)
>
> but I am not
Hyunchul Kim wrote:
> How to format and join strings ?
> For example, like following short python examples.
>
> *
> name1 = 'sample-plot'
> filename = '%s.png' % name1
Two options are sprintf() and paste():
> name1 = "sample-plot"
> sprintf("%s.png", name1)
[1] "sample-plot.png"
> paste
Hello, I'm sure that this question is too simple, but, I'm begining with R and
I'm not able to change the NA values in a matrix by 0 values, and it's
necessary for my work, how can I do It? Thank you.
Alfonso Pérez Rodríguez
Instituto de Investigaciones Marinas
C/ Eduardo Cabello nº 6
C.P. 3
Ops,
paste(inputdir, filename, sep="") should be paste(inputdir, filename, sep="/")
On 19/02/2008, Henrique Dallazuanna <[EMAIL PROTECTED]> wrote:
> Try this:
>
> name1 <- 'sample-plot'
> filename <- sprintf("%s.png", name1)
> inputdir <- '/path/to/dir'
> paste(inputdir, filename, sep="")
>
>
> On
Hi all!
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?
Thank you!
Irene
[[alternative HTML version d
?RSiteSeach is useful. It gives you this:
http://tolstoy.newcastle.edu.au/R/e2/help/07/09/25536.html
and then keep clicking on "Next in thread".
Gabor
On Tue, Feb 19, 2008 at 12:33:21PM +0100, Alfonso Pérez Rodríguez wrote:
> Hello, I'm sure that this question is too simple, but, I'm begining w
Thanks a lot, Prof Lumley.
Now I can fit a model like
coxfit=coxpy((time,censor)~pspline(x1)+x2+x3)
but I am not sure how to extract the estimated function for x1 alone. I
tried to use predict function but couldn't find appropriate option to do
this.
if i only have one covariates, then the exam
Dear list members,
I have a series of plots individuated by coordinates like this:
X Y plot
0 1 001
0 2 002
..
8 8 064
Starting from a random selected plot, I'd like to obtain all the possible
combinations of t
On 18/02/2008 4:53 AM, baptiste Auguié wrote:
> Hi,
>
> I have two small issues with my R code, no big deal but curiosity
> really. Here is a sample code,
>
>
>> x <- rnorm(1:10)
>>
>> foo <- function(a = 1, b = list(x = c(1:10), y = c(1:10))){
>>
>> for (ii in seq(along=b$y)){
>>
Hallo,
I have a question to plotDensity and do not understand what stand behind. In a
density plot the x-axis is the signal intensity but for what stands than
density on the y-axis? Here I have the values 0.00-0.30 Can anyone discribe it
by his own words? I do not understand the help.
Thanks,
Dear all,
how can I perform a repeated measures ANOVA using a covariance matrix as input?
E.g., I have four repeated measures (N = 200) with mean vector tau (no
BS factor):
tau <- rbind(1.10, 2.51, 2.76, 3.52)
and covariance matrix (sigma):
sigma <- matrix(c(0.523, 0.268, 0.267, 0.211,
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
Many thanks
Bernard
-
[[alt
sun71 sun wrote:
> Dear all,
> how can I perform a repeated measures ANOVA using a covariance matrix as
> input?
> E.g., I have four repeated measures (N = 200) with mean vector tau (no
> BS factor):
>
> tau <- rbind(1.10, 2.51, 2.76, 3.52)
>
> and covariance matrix (sigma):
>
> sigma <- matrix(c(
Huilin Chen wrote:
> Dear all,
>
> I am make plots using the function windrose() from the package "circular".
> I would like to have the north on the top instead of having 0 degree on
> the right side.
> Does anybody know how to change the labels in a windrose plot?
>
>
Hi Huilin,
As it happens
Stephan Spat wrote:
> Hello!
>
> I would like to visualize the hospitalization within one year of several
> patients using a bar chart. For each patient the stay in a hospital
> should be illustrated with a dark colour a if there is a stay at home
> between 2 hopital stays, it should be illustrate
Please just use file.path(inputdir, filename).
Apart from being less error-prone, in R 2.7.0 it will be internal and so
faster (and in case you wonder how paste could be too slow, it was in some
of the package installation/checking computations).
On Tue, 19 Feb 2008, Henrique Dallazuanna wrot
On Tue, 19 Feb 2008, Keizer_71 wrote:
>
> I am trying to calculate the distance metrics, but i keep getting this error.
> I am very new to R. Am i doing something wrong.
Cross-posting: see the thread at
https://stat.ethz.ch/pipermail/r-devel/2008-February/048452.html
where your question has alr
> Format your dataframe into a long format. Then use a grouping variable
> to distinguish both lines.
>
> library(ggplot2)
> DF <- data.frame(X = rep(0:20, 2), Y = c(rnorm(21), runif(21)), Z =
> gl(2, 21))
> ggplot(data = DF, aes(x = X, y = Y, colour = Z)) + geom_line()
You can also use ggplot a
I have no idea what "numbers form the RAD7 machine" look like, but
there is a whole section on importing data in the user manual.
'read.table' is a start.
To change directories, use 'setwd'. To start a new R session on
Windows, just setup a shortcut to RGUI.exe. Most of this is covered
in the do
Hello, I need to learn to use r-software for my PhD research project
involving long timelines of radon radiation variations.
I am using windows.
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
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 how do I add the mean as a point to each boxplot in the
Hello R experts
I am just learning how to build a package so this is probably a basic
question. I have read the manual( well, to be more truthful, am reading)
the manual, and I have searched the archives to the best of my ability.
I am on Windows.
When running the command "R CMD check packagen
If someone told you to learn R, then there must a local R guru in your
neighbourhood. Therefore I would suggest that you consult this person
first. Reading an introduction into R will help you too (e.g. Peter
Dalgaard. Introductory Statistics with R. Springer, 2002. ISBN
0-387-95475-9)
HTH,
Thie
Is 'plotDensity' a specific function you are referring to? If so, in
which package? Then we can give a straight answer...
/Henrik
On Feb 19, 2008 4:10 AM, <[EMAIL PROTECTED]> wrote:
> Hallo,
>
> I have a question to plotDensity and do not understand what stand behind. In
> a density plot the
I think you need "-l":
R CMD check -l
Gabor
On Tue, Feb 19, 2008 at 02:00:38PM -, john seers (IFR) wrote:
[...]
> I think this is because I install my packages in "mylibrary" and use
> "R_LIBS=C:/PROGRA~1/R/mylibrary" in my Renviron.site file. If I move the
> package to the R main library
On Tue, 19 Feb 2008, john seers (IFR) wrote:
>
>
> Hello R experts
>
> I am just learning how to build a package so this is probably a basic
> question. I have read the manual( well, to be more truthful, am reading)
> the manual, and I have searched the archives to the best of my ability.
You cou
Hello,
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.
--
-- Yianni
[[alternative HTML version deleted]]
__
I have been asking these same questions here on this list half a year before.
You will probably find what you need by following this link and the
subsequent discussion:
https://stat.ethz.ch/pipermail/r-help/2007-July/136162.html
Donatas
On Tuesday 19 February 2008 15:14:19 jim holtman rašė:
>
I have found several resources very helpful:
http://cran.r-project.org/doc/contrib/Paradis-rdebuts_en.pdf
http://socserv.mcmaster.ca/jfox/Misc/Rcmdr/Getting-Started-with-the-Rcmdr.pdf
Donatas
On Tuesday 19 February 2008 15:14:19 jim holtman rašė:
> I have no idea what "numbers form the RAD7 ma
Many thanks, this clarifies both points perfectly.
baptiste
On 19 Feb 2008, at 11:48, Duncan Murdoch wrote:
> On 18/02/2008 4:53 AM, baptiste Auguié wrote:
>> Hi,
>> I have two small issues with my R code, no big deal but curiosity
>> really. Here is a sample code,
>>> x <- rnorm(1:10)
>>>
Peter,
thank you - that was helpful!!!
I know, there is a more elegant way of doing it, but
N <- 200
dat <- as.data.frame(X2)
m <- tau - mean(dat)
tau_sim <- cbind(rep(m[1,],N), rep(m[2,],N), rep(m[3,],N), rep(m[4,], N))
mat <- X2+tau_sim
data <- as.data.frame(mat)
mean(data)
should give me the
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 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.
___
why do you say this? Check for instance the following code:
x1 <- runif(100, -5, 5)
x2 <- rnorm(100)
x3 <- rep(0:1, each = 50)
y <- 2 + 0.5 * x1 + 2 * x2 - 3 * x3 + rnorm(100)
fit <- lm(y ~ x1 + x2 + x3)
all.equal(
as.vector(fitted(fit)),
as.vector(model.matrix(fit) %*% coef(fit))
)
I
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:length(meriter)) { print(table(meriter[[1]],
meriter[[i]])) }
Kind regards,
Kim
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 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
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
Hi,
I'm wanting to install R onto my webspace. My ISP uses CentOS release 4.6
and I'd like to know which one of the available R packages for Linux would
be appropriate.
Thanks.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo
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
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
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
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
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
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
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
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
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
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
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 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
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 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
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
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
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
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
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
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
> 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:/
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:
# =
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
--
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 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
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
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
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
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
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
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
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
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]])
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
Hi
A simple example of a linear model:
x<-1:10
y<-3*x+1
m1<-lm(y~x)
y
# [1] 4 7 10 13 16 19 22 25 28 31
fitted(m1)
# 1 2 3 4 5 6 7 8 9 10
# 4 7 10 13 16 19 22 25 28 31
The fitted and calculated values look identical to me.
Can you give an example of how your calculated values
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
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
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
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
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
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
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
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
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 +
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
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
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
1 - 100 of 149 matches
Mail list logo
