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 .
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 problem I
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,
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. I
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 each
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
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
-
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) {
mus
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
Lassana TOURE lassana.toure at 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
[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, 2.58, 3.27,
sich at 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
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
Takatsugu Kobayashi tkobayas at 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 - GridTopology(c(1,1),
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 coefficients
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
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? Thank in
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 experiences, how
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 p5 1
From
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 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
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
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(0,1,0.05)[20]-0.95
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:
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 p1 1
2 p1
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 if we
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'
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
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 other
hadley wickham h.wickham at 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
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 enough
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,
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
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
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
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
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
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:
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 linear.
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 run out
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).
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
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?
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-bit
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
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
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 all
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
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'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) +
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] wrote:
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
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,
53 matches
Mail list logo