Hi,
What was I doing wrong, such that B) gave warning messages ? I want the
computation of thr to be outside the apply function.
A) Uses a simple matrix of 9 elements. No warning messages.
data2_1 <- matrix (c(1,2,3,NA,4,5,NA,NA,6), 3,3)
mean <- colMeans(data2_1, na.rm = TRUE)
sd <- sd(data2_1,
1) See ?"Memory-limits": it is almost certainly memory fragmentation.
You don't need to give the memory back to the OS (and few OSes actually do
so).
2) I've never seen this running a 64-bit version of R.
3) You can easily write a script to do this. Indeed, you could write an R
script to run m
Here is another solution. It uses only R core
functions. Note that as.Date(cut(x, "years"))
gives Jan 1 of x's year where x is of class "Date":
y <- as.Date(cut(range(my.dates), "years")) + c(0, 364)
seq(y[1], y[2], "months")
On Mon, Feb 4, 2008 at 1:00 PM, Gabor Grothendieck
<[EMAIL PROTECTED]
> I'm now using image() to show image data (in my case dumps of SOM
> weights) but would like to show RGB colour data, not just single "z"
> colour values.
You can do this fairly readily with ggplot2:
install.packages("ggplot2")
library(ggplot2)
qplot(x, y, data=mydata, fill=rgb, geom="tile") + s
On 4 February 2008 at 20:45, Doran, Harold wrote:
| I have a program which reads in a very large data set, performs some
analyses, and then repeats this process with another data set. As soon as the
first set of analyses are complete, I remove the very large object and clean up
to try and make
I have a program which reads in a very large data set, performs some analyses,
and then repeats this process with another data set. As soon as the first set
of analyses are complete, I remove the very large object and clean up to try
and make memory available in order to run the second set of an
Hello all,
I'm now using image() to show image data (in my case dumps of SOM
weights) but would like to show RGB colour data, not just single "z"
colour values.
I've currently been using seq() to skip 4 values, so I can show the R, G
or B channels separately as "z". But is there a way I can show
Hi
I am using Cox regression to identify at risk groups. How can I get the
C-index in R?
Thanks,
Bereket
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Hi
I am using Cox regression to identify at risk groups. How can I get the
C-index in R?
Thanks,
Bereket
<'[EMAIL PROTECTED]'>
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-h
I have a slight conundrum. I'm attempting to write a scrip that will
take a number of objects (lm, glm, and lmer) and return AIC scores
and weights. I've run into 3 problems, and was wondering if anyone
had any pointers.
1) is there any convenient way to extract the name of the objects?
Does this do what you want?
> your.dates <- c(as.Date('2005-05-21'), as.Date('2007-03-11'))
> x.range <- as.numeric(format(range(your.dates),"%Y"))
> x <- expand.grid(x.range[1]:x.range[2], 1:12)
> x <- sort(sprintf("%d-%02d-01", x[,1], x[,2]))
>
> x
[1] "2005-01-01" "2005-02-01" "2005-03-01" "20
There are many ways to do it. The following will place a blue point on the
boxplot at the mean, then print the mean at the bottom of the plot. In some
plots I've gone too far and included median points and values as well. You
could also put 95% CI on the same plot, but it would get perhaps too "bu
Hadley,
On Feb 4, 2008 5:03 PM, hadley wickham <[EMAIL PROTECTED]> wrote:
> [...]
> Before Christmas I started working on a solution for this -
> http://crantastic.org - a site for searching, reviewing and tagging R
> packages. Unfortunately I've run out of steam lately (and the lack of
> a 64-bi
Hi everyone,
I'm investigating the possibility of calling R directly from Java using the
Omegahat RSJava package (http://www.omegahat.org/RSJava/). RSJava is LGPL,
but R is GPL. As I understand it, RSJava talks to R via dynamic linking,
which the FSF considers to be derivative work under the terms
On Feb 4, 2008 12:17 PM, Dirk Eddelbuettel <[EMAIL PROTECTED]> wrote:
>
> On 4 February 2008 at 10:03, hadley wickham wrote:
> | Before Christmas I started working on a solution for this -
> | http://crantastic.org - a site for searching, reviewing and tagging R
> | packages. Unfortunately I've ru
> The real answer was Task Views on CRAN (most of the OQs topics *are* already
> Task Views), so crantastic is very partial. If you have a little time and want
I think crantastic and task views solve somewhat different problems
(although I agree that crantastic should mirror the task views too).
T
Is this what you want?
> x <- read.table(textConnection(" chrN start end
+ 1 chr1 11122333 11122633
+ 2 chr1 11122333 11122633
+ 3 chr3 11122333 11122633
+ 8 chr3 111273334 111273634
+ 7 chr2 12122334 12122634
+ 4 chr1 21122377 21122677
+ 5 chr2 33122355 33122655
+ 6 chr2
Gabor Csardi wrote:
>On Mon, Feb 04, 2008 at 03:21:10PM +0800, Ng Stanley wrote:
>
>
>>Hi,
>>
>>Given a test matrix, test <- matrix(c(1,2,3,NA,2,3,NA,NA,2), 3,3)
>>
>>A) How to compute the counts of each column (excluding the NA) i.e., 3, 2, 1
>>
>>
>
>apply(test, 2, function(x) sum(!is.na(x
Hi Peter
I have the following data frame with chromosome name, start and end positions:
chrN start end
1 chr1 11122333 11122633
2 chr1 11122333 11122633
3 chr3 11122333 11122633
8 chr3 111273334 111273634
7 chr2 12122334 12122634
4 chr1 21122377 21122677
5 chr2 33122355 3
On Mon, 04-Feb-2008 at 10:42AM -0700, Todd Remund wrote:
|>
|> I'm ordering a new computer to increase my ability to handle large
|> data sets. I?ve tried the dual core type and also the dual
|> processor with dual cores each, and have not been satisfied. This
|> seems to agree with all the ot
Dear All,
I have a shapefile wich contents are polygons. I sort two (here random) points
and I would like to identify if the line built with this two points cross or
not polygons and in case yes, which polygons were crossed. I tryed to figure
out the problem with a sample code got from sp/mapto
I have a time series object with two columns.
a <- ts(x, frequency=1/15)
I would like to split this into two time series of the same length
(they are the same number of observations lined up time for time
already).
--
Let's not spend our time and resources thinking about things that are
so littl
Thank you, thank you both! get() is exactly what I needed. I
apologize for not putting in reproducible code. I'm not very good at
this and was more worried that I would create problems that didn't
exist (and the files I'm working with are too big). Again, thank you.
-John
On Feb 4, 200
On 4 February 2008 at 10:03, hadley wickham wrote:
| Before Christmas I started working on a solution for this -
| http://crantastic.org - a site for searching, reviewing and tagging R
| packages. Unfortunately I've run out of steam lately (and the lack of
| a 64-bit ubuntu package for R means it
OK. Here are two zoo solutions and one using only core
functions.
In the first dd is the first of the month of Jan
and month of Dec of the first and last year respectively
so we apply seq.Dates to that.
In the second yr is Jan of the first and Jan of the
last year as "yearmon" class and we adjust
hadley wickham gmail.com> writes:
> > Also it is staggering that there are over 1200 packages for R i was
suspecting close to 1000 .
>
> And the site is a couple of weeks out of date, so there are probably even
> more.
>
The real answer was Task Views on CRAN (most of the OQs topic
Im ordering a new computer to increase my ability to handle large data sets.
Ive tried the dual core type and also the dual processor with dual cores each,
and have not been satisfied. This seems to agree with all the other postings
on the help list. I dont want to do any simulations, I j
> Your site is interesting the little example with the graphics is
> actually what i was thinking about. Now I've tagged 2 packages as "robust
> analysis" but if i press the Tags button i see only my tags - i don't see the
> graphics tags anymore. I would be extremely happy if you get enoug
The recommended citation can be found with
citation()
How exactly you cite it in the text is dependent upon
what the target journal's style manual says. I'm no
expert but the information you have given in the text
below looks quite adequate.
--- Falco tinnunculus <[EMAIL PROTECTED]> wrote:
> D
hits=-2.6 tests=BAYES_00
X-USF-Spam-Flag: NO
On Mon, 2008-02-04 at 10:48 -0500, Gabor Grothendieck wrote:
> Using zoo's yearmon class:
>
> library(zoo)
> my.dates[!duplicated(as.yearmon(my.dates))]
>
> or, although you seem to disallow this in your question,
> this would be an option:
>
> my.da
You need to create a data.frame in a different way is
all
Try this:
df1 <-
data.frame(rownames(affect),matrix(affect,nrow=3))
merge(age,df1, by.x=1, by.y=1)
--- "Boks, M.P.M." <[EMAIL PROTECTED]> wrote:
> Dear R-experts,
>
>
> I have got a dataframe:
> data
> ID disease
> V1 V2
> 1
Not precisely what you asked for but see the notch= argument to boxplot
for a graphic measure of variability. If you simply wish to print certain
statistics below the numbers already on the X axis then see:
https://stat.ethz.ch/pipermail/r-help/2008-January/152994.html
On Feb 4, 2008 10:41 AM, T
> "As for the standard deviation, are you sure you want this? Standard
> deviation only makes sense if the data are normally distributed..."
> --
>
> -- This is false, of course. What you probably meant to say is something
> like:
>
> "The sample standard deviation may not tell you
?get (maybe)
If this isn't it and you do not receive satisfactory replies, see the
Posting Guide for how to post a question with reproducible code so that
readers can understand what you mean (as opposed to me being the only one
who doesn't understand you).
-- Bert Gunter
Genentech
-Orig
"As for the standard deviation, are you sure you want this? Standard
deviation only makes sense if the data are normally distributed..."
--
-- This is false, of course. What you probably meant to say is something
like:
"The sample standard deviation may not tell you what you think
Hello all,
I've run into what I bet is a silly problem; however, I've been
trying to get around it now for a couple weeks and every time I think
I have the answer it still doesn't work. So I apologize in advance
if this is painfully obvious, but I've run out of ideas and would
really ap
Good morning,
I'm trying to write a function to assign the mean separation grouping
letters to the factors of an experiment. A client wants this so I
thought I'd write my own function for future use. But I'm having
trouble doing it because I don't understand the logic of the
assignment. As
Hadley,
Your site is interesting the little example with the graphics is actually
what i was thinking about. Now I've tagged 2 packages as "robust analysis" but
if i press the Tags button i see only my tags - i don't see the graphics tags
anymore. I would be extremely happy if you get en
>> How can I add the mean and standard deviation to each of the boxplots
using the example provided in the boxplot function?
boxplot(len ~ dose, data = ToothGrowth)
#You could add the mean as a point
points(tapply(len, dose, mean), col="red")
Alternatively, replace the median value with the me
Sorry, I wasn't very helpful. Let me try this again. I have attached a
subsample of the data which still gives me the same error as when I use the
full data file. I am trying to make a decision tree using rpart. This is my
code and output.
> data <- read.table("/Users/randygriffiths/Desktop/data",
On Feb 4, 2008 8:34 AM, Monica Pisica <[EMAIL PROTECTED]> wrote:
>
> Hi everybody,
>
>
> I know this might be very off topic and it took me quite a while to up my
> courage to post this…. But I remember a thread some time ago about how we can
> find the packages we need to do specific tasks in R
On Mon, 4 Feb 2008, Falco tinnunculus wrote:
Ft > Dear all,
Ft >
Ft > How do I make reference to R in the method section in a scientific article?
Ft > Should I state the web aderess?
Ft >
Ft > And, is this the proper way to report the lme test?
Ft >
Ft >
Ft > "The relationships were assumed to be
hits=-2.6 tests=BAYES_00
X-USF-Spam-Flag: NO
On Mon, 2008-02-04 at 10:48 -0500, Gabor Grothendieck wrote:
> Using zoo's yearmon class:
>
> library(zoo)
> my.dates[!duplicated(as.yearmon(my.dates))]
>
> or, although you seem to disallow this in your question,
> this would be an option:
>
> my.da
On Mon, 4 Feb 2008, Falco tinnunculus wrote:
Ft > Dear all,
Ft >
Ft > How do I make reference to R in the method section in a scientific article?
Ft > Should I state the web aderess?
Ft >
Ft > And, is this the proper way to report the lme test?
Ft >
Ft >
Ft > "The relationships were assumed to
Dear list,
How can I add the mean and standard deviation to each of the boxplots using
the example provided in the boxplot function?
boxplot(len ~ dose, data = ToothGrowth,
boxwex = 0.25, at = 1:3 - 0.2,
subset = supp == "VC", col = "yellow",
main = "Guinea Pigs'
Using zoo's yearmon class:
library(zoo)
my.dates[!duplicated(as.yearmon(my.dates))]
or, although you seem to disallow this in your question,
this would be an option:
my.dates[!duplicated(format(my.dates, "%Y-%m"))]
On Feb 4, 2008 10:39 AM, Gavin Simpson <[EMAIL PROTECTED]> wrote:
> hits=-2.6 t
Hi Jonathan,
I didn't know about your web page - now it is in my favorites. I am sure it
will help a lot. Google does and doesn't do the job - sometimes ;-) But R and
packages are an open source and i realize that people already put a lot of good
effort into something which is essentially fre
hits=-2.6 tests=BAYES_00
X-USF-Spam-Flag: NO
Dear List,
Say I have the following sequence of dates [*]:
start <- as.Date("2005-01-05", format = "%Y-%d-%m")
end <- as.Date("2007-10-07", format = "%Y-%d-%m")
my.dates <- seq(start, end, by = "days")
What I would like to generate is a sequence of d
Or use RSiteSearch, with restrict="functions", or use the site itself:
http://finzi.psych.upenn.edu/
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide
Perhaps:
affect <- as.data.frame(do.call('rbind', tapply(data$V2, data$V1, table)))
merge(age, affect, by.x=1, by.y=0)
On 04/02/2008, Boks, M.P.M. <[EMAIL PROTECTED]> wrote:
> Dear R-experts,
>
>
> I have got a dataframe:
> data
> ID disease
> V1 V2
> 1 p1 1
> 2 p1 3
> 3 p3 3
> 4 p3 5
> 5
Hi
I think this is a good suggestion.
And I would like to add the associated problem of deciding between
packages that do the same function which one is "better". Or similarly
packages are often superceded. I find I have to spend a lot of time
learning how to use packages to decide which one
dear list,
i'd like to do some calculation for which i need the free float
(german :streubesitzanteil) of a number of swiss firms' shares. but at
swx.com i did only find very few information. does anyone know, if and
where that information is available? best would be a website where
it's extractabl
Hello,
the lm function allows us to define some weights for the least square
optimisation. I'm looking for the same type of possibility but for ARIMA model.
Actually, I think it is not possible with the ARIMA function. Does anyone know
something about that?
Best regards,
Yannig.
Google for CRAN Task Views
On Feb 4, 2008 9:34 AM, Monica Pisica <[EMAIL PROTECTED]> wrote:
>
> Hi everybody,
>
>
> I know this might be very off topic and it took me quite a while to up my
> courage to post this…. But I remember a thread some time ago about how we can
> find the packages we nee
dt Excellent <[EMAIL PROTECTED]> wrote in
news:[EMAIL PROTECTED]:
> Does anybody know if there is an implementation of Goodman-Kruskal
> lambda measures of association in R? Also, how we can analyze
> ordered contingency tables and compute the relative measures of
> associations in R? Tha
Hi everybody,
I know this might be very off topic and it took me quite a while to up my
courage to post this…. But I remember a thread some time ago about how we can
find the packages we need to do specific tasks in R if we don’t know before
hand which ones actually do it. Now all the package
Dear R-experts,
I have got a dataframe:
data
ID disease
V1 V2
1 p1 1
2 p1 3
3 p3 3
4 p3 5
5 p5 1
>From which I extract a usefull table: affect
affect
1 3 5
p1 1 1 0
p3 0 1 1
p5 1 0 0
I want to merge this with anotherdataframe:
age
p1 23
p2 24
p3 23
p4 11
p5 45
On Feb 3, 2008 2:16 PM, <[EMAIL PROTECTED]> wrote:
> Hi there,
> I am interested in using R for machine learning (supervised classification).
> Currently, I have been investigating especially the rpart, tree, and
> randomForest package, and have achieved first results.
>
> are there any experienc
On 2/3/2008 10:09 AM, Christoph Mathys wrote:
> Dear R users,
>
> I have a linear model of the kind
>
> outcome ~ treatment + covariate
>
> where 'treatment' is a factor with three levels ("0", "1", and "2"),
> and the covariate is continuous. Treatments "1" and "2" both have
> regression coeffi
FAQ Q7.31 and ?"==" should enlighten you.
None of 0.05, 0.1 ... 0.9, 0.95 are exactly representable on a binary
computer.
On Mon, 4 Feb 2008, Eric Elguero wrote:
> Hi everybody,
>
> this is a warning more than a question.
>
> I noticed that seq produces approximate results:
>
>> seq(0,1,0.05)[1
FAQ 7.31
On 2/4/08, Eric Elguero <[EMAIL PROTECTED]> wrote:
> Hi everybody,
>
> this is a warning more than a question.
>
> I noticed that seq produces approximate results:
>
> > seq(0,1,0.05)[19]==0.9
> [1] TRUE
> > seq(0,1,0.05)[20]==0.95
> [1] FALSE
> > seq(0,1,0.05)[21]==1
> [1] TRUE
>
> > seq
Takatsugu Kobayashi indiana.edu> writes:
>
> try
>
> tmp<- slot(ex_1.7.selected, 'polygons')
> sub.tmp <- slot(tmp[[1]],'Polygons')
> [EMAIL PROTECTED]
>
> will get you there.
>
> taka
>
> Jarek Jasiewicz wrote:
> > Milton Cezar Ribeiro wrote:
> >
> >> Dear all,
> >>
> >>
> >> grd <- Gri
Hi everybody,
this is a warning more than a question.
I noticed that seq produces approximate results:
> seq(0,1,0.05)[19]==0.9
[1] TRUE
> seq(0,1,0.05)[20]==0.95
[1] FALSE
> seq(0,1,0.05)[21]==1
[1] TRUE
> seq(0,1,0.05)[20]-0.95
[1] 1.110223024625157e-16
I do not understand why 0.9 and 1 are
Thanks Bernhard for the beautiful code!!
On 2/4/08, Pfaff, Bernhard Dr. <[EMAIL PROTECTED]> wrote:
>
> Dear All,
>
>
> one can visually inspect ARCH-effects by plotting acf/pacf of the
> squared residuals from an OLS-estimation. This can be as simple as a
> demeaned series. Further one can run an
Hi,
>> I am interested in using R for machine learning (supervised
classification).
>> Currently, I have been investigating especially the rpart, tree, and
randomForest package, and have
>> achieved first results.
>>
>> are there any experiences, how the learned classificators could
>> be u
When you start R you get this information:
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
That should help answer your question.
Regards
JS
---
-Original Message-
From
Dear all,
How do I make reference to R in the method section in a scientific article?
Should I state the web aderess?
And, is this the proper way to report the lme test?
"The relationships were assumed to be linear. The response, handling time
and the explanatory variable, prey mass, were log t
gmx.de> writes:
> I am interested in using R for machine learning (supervised classification).
> Currently, I have been investigating especially the rpart, tree, and
randomForest package, and have
> achieved first results.
>
> are there any experiences, how the learned classificators could
> be
[EMAIL PROTECTED] napsal dne 02.02.2008 07:10:37:
> bgchen wrote:
> > r-help
> >
> > c(2.43, 3.22, 6.9, 3.03, 5.36, 6.9, 2.29, 6.13, 6.11, 4.25, 3.85,
> > 5.09, 7.44, 2.86, 2.82, 3.64, 3.22, 7, 2.65, 4.5, 3.73, 5.09,
> > 5.8, 7.87, 2.87, 2.9, 6.63, 6.8, 2.45, 7.68, 2.56, 2.54, 7.35,
> > 4.61,
Lassana TOURE ier.ml> writes:
> Can someone help me how to set up a directory and load additional packages
> at the startup than I do not have to do it each time I go to R session.
It depends on your system and the type of installation, but searching for
.Rprofile in the manuals and CRAN might
Does anybody know if there is an implementation of Goodman-Kruskal lambda
measures of association in R?
Also, how we can analyze ordered contingency tables and compute the relative
measures of associations in R?
Thank in advance
-
[[alterna
hits=-2.5 tests�YES_00,FORGED_RCVD_HELO
X-USF-Spam-Flag: NO
check the following options:
test <- matrix(c(1,2,3,NA,2,3,NA,NA,2), 3, 3)
# A
colSums(!is.na(test))
# B
mat <- test > rep(colMeans(test, na.rm = TRUE), each = nrow(test))
colSums(!is.na(mat) & mat)
apply(test, 2, function(x) {
m
Dear All,
one can visually inspect ARCH-effects by plotting acf/pacf of the
squared residuals from an OLS-estimation. This can be as simple as a
demeaned series. Further one can run an auxiliary regression by
regressing q lagged squared values and a constant on the squared series
itself. This tes
Hi!
Can someone help me how to set up a directory and load additional packages
at the startup than I do not have to do it each time I go to R session.
Thanks.
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://
On Mon, Feb 04, 2008 at 03:21:10PM +0800, Ng Stanley wrote:
> Hi,
>
> Given a test matrix, test <- matrix(c(1,2,3,NA,2,3,NA,NA,2), 3,3)
>
> A) How to compute the counts of each column (excluding the NA) i.e., 3, 2, 1
apply(test, 2, function(x) sum(!is.na(x)))
> B) How to compute the counts of e
Use a working directory where you have write permission.
On Sat, 2 Feb 2008, [EMAIL PROTECTED] wrote:
> Good afternoon,
>
> I've installed R on Suse 10.3 compiling the source files. I opne the program
> on shell configuration and i have a problem saving the workspace and *.R or
> *.Rdata files.
On 02-Feb-08 21:16:41, R-novice wrote:
>
> I am trying to make an array c(3,8) that contains the averages of what
> is in another array c(9,8). I want to average rows 1:3, 4:6, 7:9 and
> have a loop replace the generic 1:24 values in my array with the
> average of the three rows.
>
> The probl
What OS is this?
If Windows, see rw-FAQ Q3.2
Otherwise, see 'R Installation and Administration' Chapter 7.
On Sat, 2 Feb 2008, cris wrote:
Dear all,
i have recently update R in my computer. To my surprise the menu bar
is written in German language. This is the first time this happen to
me, u
hi,
i start working with PCA method i found VSS package it is very helpfull it
show me residual matrix with the number of components to extract now i want to
verify the normality of residual matrix and than to simulate can every one help
me please .
-
79 matches
Mail list logo
