Tom wrote:
On Tue, 10 Jan 2006 20:25:23 -0500, r user [EMAIL PROTECTED] wrote:
I have 2 dataframes, each with 5 columns and 20 rows.
They are called data1 and data2.I wish to create a
third dataframe called data3, also with 5 columns and
20 rows.
I want data3 to contains the values in data1
ronggui wrote:
I think it can.But if you provide more information,you will be more help.
for example,you had better give a reproducable example in you email.
2006/1/11, Naiara S. Pinto [EMAIL PROTECTED]:
Hi all,
I need some help using read.ftable to read a contingency table. My columns
Johannes Hüsing wrote:
Hello,
par(las=1) sets the orientation of the axis labels
to horizontal. That is, the tick mark labels. How
do I set the orientation of the axis label, which
annotates the variable plotted along the axis, to
horizontal?
Sorry for asking such a basic question here,
Hi there,
I hope that I am in the right forum.
I am using Exceed on a WinNT2000 machine, connected to Solaris
SunOS fluke 5.9 Generic_118558-11 sun4u sparc SUNW,Sun-Fire-480R
with Mozilla Firefox 1.0.7.
Problem: Manuals like R-exts.html, R-lang.html, R-intro.html,
R-admin.html, but *not*
--
Alexandra R. M. de Almeida wrote:
Dear list
I want to obtain the adjusted r-square given a set of coefficients
(without the intercept), and I don't know if there is a function that
does it. Exist I know that if you make a linear
You can read the code of summary.lm and adapt
The equality operator is == not =. So you need is.na(data1) == FALSE (F
is a variable, and FALSE is the non-truth value), or, clearer,
ifelse(!is.na(data1), data1, data2)
Another way is
data3 - data1
data3[is.na(data1)] - data2[is.na(data1)]
which is more efficient but less clear.
On Tue,
Gabor == Gabor Grothendieck [EMAIL PROTECTED]
on Tue, 10 Jan 2006 14:47:57 -0500 writes:
Gabor If its just succint you are after then this is slightly
Gabor shorter:
Gabor identical(dim(x)+0, c(3,5))
indeed, or, less succinct, but maybe more readable (and along the
top-level
You are mixing makes. GNU make (presumably gmake) passes on its -w
argument to sub-makes, and my guess is that make is a BSD make that does
not accept it.
The simplest way out is to have the 'make' first in your path as GNU make
whilst doing this.
BTW, this really is not the appropriate
dear all,
yes but the problem with soundex for example is that it does not work when
an error occur in the first place (Canada vs Kanada) as it keeps the fist
character. It seems that you have to look after an approximate string
matching algorithm (for example, a very good one if from
Hello Alexandra,
R2 is only defined for regressions with intercept. See a decent econometrics
textbook for its derivation.
HTH,
Bernhard
-Ursprüngliche Nachricht-
Von: Alexandra R. M. de Almeida [mailto:[EMAIL PROTECTED]
Gesendet: Mittwoch, 11. Januar 2006 03:48
An:
On Wed, Jan 11, 2006 at 08:05:22AM +, Prof Brian Ripley wrote:
You are mixing makes. GNU make (presumably gmake) passes on its -w
argument to sub-makes, and my guess is that make is a BSD make that does
not accept it.
The simplest way out is to have the 'make' first in your path as
I'm writing an publication in two column format and need to shrink some
plots. After increasing the axis labels it does not look nice at all.
The y-axis label and tick labels almost touch each other and the x-axis
tick labels expand into the plot instead of away from it. Is there a
better way
David == David Forrest [EMAIL PROTECTED]
on Mon, 9 Jan 2006 11:54:30 -0600 (CST) writes:
..
..
David Since R has such an extensive set of extensions,
David maybe we need a section in the R-intro documentation
David near
On Tue, 10 Jan 2006, McGehee, Robert wrote:
I would throw a tolower() around s1 and s2 so that 'canada' matches with
'CANADA', and perhaps consider using a Levenshtein distance rather than
the longest common subsequence.
An algorithm for Levenshtein distance can be found here (courtesy of
Andrew == Andrew Robinson [EMAIL PROTECTED]
on Wed, 11 Jan 2006 14:12:11 +1100 writes:
Andrew Hi R-help citizens,
Andrew I'm having trouble making version 0.99-6 of Matrix on FreeBSD 6.0.
Andrew The error message is:
Andrew * Installing *source* package 'Matrix' ...
Dear alltogether,
how can a delay like possible with par(ask=TRUE) be attained while using
trellis-plots within a loop or something like that?
the following draws each plot without waiting for a signal
(mouse-klick), so par() does not work for that:
library(nlme)
for(i in 1:3)
{
fitlme -
normOrder() in SuppDists
Anna Oganyan wrote:
Hello,
Could somebody point me, is there any function in R which returns
expected values of order statistics for normal distribution? I have been
looking and couldn't find it.
Thanks!
Anna
__
Hi,
I'm using gregexpr(). As a result something like this:
# starting positions of the match:
[[1]]
[1] 7 18
# length of the matched text:
attr(,match.length)
[1] 4 4
Now, I'd like to have a matrix,
74
18 4
but I don't know how to handle the attr(,match.length) ...?
The format of
Dear R-users,
We expect to create N homogenic groups of n features from an
experimentation including N*n mesures. The aim of this is to prevent from
group effects. How to do that with R functionalities. Does anyone know any
methodes enabling this ?
Best regards.
Alexandre MENICACCI
Dear Rlist,
We are working with library(Design) R 2.2.1//
When using the following fitted model:
knots - 5
lrm.1- lrm(X8~rcs(X1,5),x=T,y=T)
X8 (binary 0/1 vector)
X1, X2 explantory variables
We would like to set the probability of X8=1 to zero when the X2
variable is
Now I found a solution that seems to work OK for me:
attributes(gregexpr(expression, text)[[1]])
Petri
At 15:00 11.1.2006 +0100, Petri Palmu wrote:
Hi,
I'm using gregexpr(). As a result something like this:
# starting positions of the match:
[[1]]
[1] 7 18
# length of the matched text:
I was aware of that which is why I mentioned that it is usually used
for matching last names rather than countries and noted possible need to
modify the algorithm slightly. soundex is a relatively simple algorithm
so its not too hard. For example, one could just code the first
letter too.
On
Petri Palmu petri.palmu at geneos.fi writes:
I'm using gregexpr(). As a result something like this:
# starting positions of the match:
[[1]]
[1] 7 18
# length of the matched text:
attr(,match.length)
[1] 4 4
Now, I'd like to have a matrix,
74
18 4
something like
x1
Hi Petri,
On 11 Jan 2006, [EMAIL PROTECTED] wrote:
I'm using gregexpr(). As a result something like this:
# starting positions of the match:
[[1]]
[1] 7 18
# length of the matched text:
attr(,match.length)
[1] 4 4
Now, I'd like to have a matrix,
74
18 4
but I don't know how
Jan Verbesselt wrote:
Dear Rlist,
We are working with library(Design) R 2.2.1//
When using the following fitted model:
knots - 5
lrm.1- lrm(X8~rcs(X1,5),x=T,y=T)
X8 (binary 0/1 vector)
X1, X2 explantory variables
We would like to set the probability of X8=1 to
Alexandra,
some additional remarks taken from my past struggles with R2 :^) Without
intercept the definition is indeed problematic, as Bernhard notes.
First, to estimate a model omitting the intercept you simply have to
specify -1 in the model formula (example on an in-built dataset, for
data
... and is there also such a nice tool (like spss.get) for exporting
data frames to SPSS? write.table does not keep the data frame labels -
neither did the other exporting tools that I found.
Thanks!
Michael
[[alternative HTML version deleted]]
Uwe,
FYI:
I tried: data3 - ifelse(is.na(data1), data2, data1)
It seems to me that data3 is an array of length 100.
I do NOT end up with a dataset of 5 columns and 20 rows.
Uwe Ligges [EMAIL PROTECTED] wrote:
Tom wrote:
On Tue, 10 Jan 2006 20:25:23 -0500, r user wrote:
I
I just got 59 hits from ' RSiteSearch(space-time)'. Have you tried
this?
If you would like more help from this listserve, please read the
posting guide! www.R-project.org/posting-guide.html then submit
another question. La experiencia sugiere que las prejuntas sigiendo
t c wrote:
Uwe,
FYI:
I tried: data3 - ifelse(is.na(data1), data2, data1)
It seems to me that data3 is an array of length 100.
I do NOT end up with a dataset of 5 columns and 20 rows.
I have not read carefully enough, for a data.frame you can generalize
the approach as
On 1/11/06, Leo Gürtler [EMAIL PROTECTED] wrote:
Dear alltogether,
how can a delay like possible with par(ask=TRUE) be attained while using
trellis-plots within a loop or something like that?
the following draws each plot without waiting for a signal
(mouse-klick), so par() does not work for
I want to obtain datetime differences in mins in an other column, in front of my
datetimes. I have tried this :
T1 - c(12/31/03 23:49,1/1/04 1:14,1/1/04 0:02)
T2 - c(1/1/04 0:58,1/1/04 1:16,)
toto - data.frame(T1,T2)
toto
y - strptime(T1,%m/%d/%y %H:%M)
x - strptime(T2,%m/%d/%y %H:%M)
The following seems close to the form you were trying. It works
for matrices, not dataframes. You can use as.matrix and
as.data.frame to convert back and forth:
# test data
data1 - data2 - matrix(1:6,3)
data1[2,2] - NA
data1[] - ifelse(is.na(data1), data2, data1)
On 1/11/06, t c [EMAIL
Try difftime(x,y,unit=min) or as.numeric(difftime(x,y,unit=min))
depending on what you want.
On 1/11/06, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
I want to obtain datetime differences in mins in an other column, in front of
my
datetimes. I have tried this :
T1 - c(12/31/03 23:49,1/1/04
Is this what you want?
toto$ans - difftime(x,y)
toto
T1 T2 ans
1 12/31/03 23:49 1/1/04 0:58 69
21/1/04 1:14 1/1/04 1:16 2
31/1/04 0:02 NA
On 1/11/06, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
I want to obtain datetime differences in mins in
A much shorter (but complete) description of this is on the summary.lm
help page. It includes the definitions R (and most statistics references)
uses.
On Wed, 11 Jan 2006, Millo Giovanni wrote:
Alexandra,
some additional remarks taken from my past struggles with R2 :^) Without
intercept
On Wed, 11 Jan 2006 [EMAIL PROTECTED] wrote:
I want to obtain datetime differences in mins in an other column, in front of
my
datetimes. I have tried this :
T1 - c(12/31/03 23:49,1/1/04 1:14,1/1/04 0:02)
T2 - c(1/1/04 0:58,1/1/04 1:16,)
toto - data.frame(T1,T2)
toto
y -
Mark:
I did not see a reply to your question. Did you get one? If not,here's a
solution using a while() loop which should be fast. One could also use
recursion here in a natural way. This solution assumes that there are no
NA's anywhere -- it's a bit trickier if there are NA's in the x column.
Thanks again for your answer! I tried it out. write.foreign produces SPSS
syntax, but unfortunally this syntax tells SPSS to take the names (and not the
labels) in order to produce SPSS variable labels. The former labels get lost.
I tried a data frame produced by read.spss and one by
Dear List,
Please confirm the following:
It may be my eyes playing trick on me, but I can't seem to find functions
that do the nearest neighbor method and the variable kernel method for
kernel smoothing and density estimation corresponding to the book:
Silverman, B.W., (1986) Density
Hi,
I want to permutate the following matrix and replace permutated columns. Is
it possible to control the number of columns permutated. Let's say I only
want to permute two columns. Can i do that with the sample method or should
i any bootstrapping method ?? I'm not sure this is the best
On Wed, 11 Jan 2006, Michael Reinecke wrote:
Thanks again for your answer! I tried it out. write.foreign produces
SPSS syntax, but unfortunally this syntax tells SPSS to take the names
(and not the labels) in order to produce SPSS variable labels. The
former labels get lost.
Well, yes.
Matching regular expressions
Dear useRs!
I have the following problem. I would like to find objects in my environment
that have two strings in it. For example, I might want to find objects that
have in their names MY and TARGET. I do not care about the ordering of
these two substrings in the
Ales Ziberna [EMAIL PROTECTED] writes:
Matching regular expressions
Dear useRs!
I have the following problem. I would like to find objects in my environment
that have two strings in it. For example, I might want to find objects that
have in their names MY and TARGET. I do not care about
The version of lmer based on the supernodal Cholesky factorization,
which we will release real soon, does not crash on this example. It
does give very large estimates of the variances in that model fit, at
least for the simulation that I ran.
It is best if you use set.seed(123454321) (or
Ales Ziberna [EMAIL PROTECTED] writes:
Dear useRs!
I have the following problem. I would like to find objects in
my environment
that have two strings in it. For example, I might want to
find objects that
have in their names MY and TARGET. I do not care about
the ordering of
these two
On an improved R wiki, R-intro:
I think the issue of user-friendliness of documentation has been raised.
When I first started using R, I found the S-PLUS online documentation
very useful. It is very user-friendly and a great introduction,
organized by application. See:
S-PLUS 6 Guide to
Dear R Users,
I am trying to use its, and for that, I need to use as.POSIXct .
My dates are of the format:10 January 2006.
How do I convert this into the format acceptable to its ?
Thanks,
Tolga
==
Please
See ?strptime
Also the help desk article in RNews 4/1 contains info on dates.
On 1/11/06, Uzuner, Tolga [EMAIL PROTECTED] wrote:
Dear R Users,
I am trying to use its, and for that, I need to use as.POSIXct .
My dates are of the format:10 January 2006.
How do I convert this into the format
On Wed, 11-Jan-2006 at 10:52AM +0100, Kare Edvardsen wrote:
| I'm writing an publication in two column format and need to shrink some
| plots. After increasing the axis labels it does not look nice at all.
| The y-axis label and tick labels almost touch each other and the x-axis
| tick labels
Thank you!
This is definitely an improvement!
Best,
Ales Ziberna
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Peter
Dalgaard
Sent: Wednesday, January 11, 2006 7:24 PM
To: Ales Ziberna
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] Regular
I guess I have not been clear enough.
I want both words in the results. So if we have:
text-c(somethigMYsomthing elseTARGET another thing,MY somthing TARGET
another thing,somethig somthing elseTARGETMY another
thing,somethigMTARGETY another thing, somthingMY somthing else)
The last element
I have a problem moving from multistratum aov analysis to lmer.
My dataset has observations of ampl at 4 levels of gapf and 2 levels of bl
on 6 subjects levels VP, with 2 replicates wg each, and is balanced.
Here is the summary of this set with aov:
On 1/11/06, Dean Sonneborn [EMAIL PROTECTED] wrote:
I am using the R code listed below to create 4 smoothed lines on a
xyplot. I'm having trouble fine tuning it. First I think I may need a
black and white plot so how do I get it to plot the lines with different
characters, preferable the same
Dear R-Helpers
Given a function like
foo - function(data,var1,var2,var3) {
f - formula(paste(var1,'~',paste(var2,var3,sep='+'),sep=''))
linmod - lm(f)
return(linmod)
}
By typing
foo(mydata,'a','b','c')
I get the result of the linear model a~b+c.
How can I rewrite the
This was just discussed last week:
https://www.stat.math.ethz.ch/pipermail/r-help/2006-January/083812.html
On 1/11/06, Christian Bieli [EMAIL PROTECTED] wrote:
Dear R-Helpers
Given a function like
foo - function(data,var1,var2,var3) {
f -
Does anyone know of a simple test
in any R package that given
a series of negative ones and positive
ones ( no other values are possible in the series )
returns a test of whether the series is random or not.
( a test at each point would be good but
I can use the apply function to implement
that )
SALFORD SYSTEMS DATA MINING CONFERENCE 2006
San Diego, California, March 29-31, 2006
Focusing on the Contributions of Data Mining to Solving Real-World
Challenges
Business, Biomedical and Environmental Real-World Case Study Presentations
TOPICS INCLUDE:
Credit Risk Modeling; Targeted Marketing
Dear Martin,
That works just fine too.
Thanks for the suggestion,
Andrew
Can you try and replace 'make' by '$(MAKE)' in the following
three places, and see if it works
possibly after writing (in your shell)
export MAKE=gmake
or setenv MAKE gmake
(depending on the kind of
I would compare the Shannon entropy of your test vector with the entropy
of your expected probability distribution to see if they are close. That
is, if you're binary probability distribution is half 1 and half -1,
then if your string is long you would expect about half the numbers in
your vector
Michael Reinecke wrote:
Thanks again for your answer! I tried it out. write.foreign produces SPSS
syntax, but unfortunally this syntax tells SPSS to take the names (and not
the labels) in order to produce SPSS variable labels. The former labels get
lost.
I tried a data frame produced
From: Spencer Graves The correlation between the predictions from your
two model fits is 0.95. This suggests to me that the differences between
the two sets of answers have little practical importance, and anyone who
disagrees may be trying to read more from the results than can actually
Thank you very much for pointing me to these very useful references! Best
regards. SC
At 01:54 PM 1/11/2006, ronggui wrote:
the usage of MNP is described in MNP: R Package for Fitting the
Multinomial Probit Model
http://www.jstatsoft.org/counter.php?id=128url=v14/i03/v14i03.pdfct=1
If the
Hi,
I have a question about the p-value of mshapiro test. I simulated data
from bivariate normal 1000 times and used mshapiro test to see how many
times the test would reject the null hypothesis when the p-value is 0.05.
The answer should be around 50 since the p-value is 0.05. But I got a much
1. The function glmmPQL is in the MASS package, as can be seen by
looking at the top line in the help file for glmmPQL. To find the
maintainer, type 'help(package=MASS)'. The results say, Maintainer:
Brian Ripley [EMAIL PROTECTED].
2. It is generally NOT appropriate to
Dear mailing group,
This is my first time here. Glad to have this resource!
I am currently trying to load an Excel file into R (limma package loaded)
using the source(*name of directory*) command, but it cannot open the file.
I renamed the file as .R and .RData, to no avail. The Excel data
well, I don't know anything about the limma package and I might be
misunderstanding your apparently simple question
What I do for excel files is the following:
1. I save a copy of the file as .csv (comma separated values) in the
working directory. This format allows you work perfectly with the
no idea bout limma. but you could load excel data into R by using rodbc.
here is a sample code and hope it helpful.
library(RODBC);
###
# 1. READ DATA FROM EXCEL INTO R #
###
R won't read an Excel sheet directly. You need to export it, saving it as a
CSV or tab delimited file. You can then import using read.table. The entire
path and file have to be in double-quotes as well. Try ?read.table for more
info.
JWD
On Wednesday 11 January 2006 19:48, N. Goodacre
another options: use read.xls in gdata pcakges if you have installed
perl in you machine .
12 Jan 2006 03:48:26 +, N. Goodacre [EMAIL PROTECTED]:
Dear mailing group,
This is my first time here. Glad to have this resource!
I am currently trying to load an Excel file into R (limma
I don't know how to keep factors' levels with :
data.frame(mapply(function(x,y,z) ifelse(is.na(y), z, y),
names(D), D, D2, SIMPLIFY=FALSE))
but in that way it's ok :
data.frame(mapply(function(z,x,y) { y[is.na(y)] - x[is.na(y)] ; y },
names(D), D, D2, SIMPLIFY=F))
(?)
Uwe
Dear All,
simetimes when I load an Rdata I get this message
###
Code:
load('bladder1.RData')
Carico il pacchetto richiesto: rpart ( Bad traslastion: Load required
package-...)
Carico il pacchetto richiesto: MASS
Carico il pacchetto richiesto: mlbench
Carico il pacchetto richiesto: survival
On Wed, 11 Jan 2006, Mark Leeds wrote:
Does anyone know of a simple test
in any R package that given
a series of negative ones and positive
ones ( no other values are possible in the series )
returns a test of whether the series is random or not.
( a test at each point would be good but
I
73 matches
Mail list logo