stefan semmelring wrote:
stefan semmelring [EMAIL PROTECTED] schrieb am 09.01.06 21:31:01:
Hallo,
I´m trying to find out how I can start WinEdt always directly when R is
started.
Does anybody know how to do so? I tried to find it out by myself but had no
fortune.
Maybe you have to
Hi
I tried to reproduce what you have told us by copy and paste
read.delim(clipboard)
but was not successful.
Even with several blank values in each column in Excel i got correct
import to R by this process. As I do not use RODBC I do not know all
possible settings and features. If
Dear R People:
In S Plus, if you have a function which calls the plot
function several times, you get several pages of graphics
output.
Is there an eqivalent in R, please?
R version 2.2.1 windows
Thanks in advance!
Sincerely,
Erin Hodgess
Associate Professor
Department of Computer and
Erin Hodgess wrote:
Dear R People:
In S Plus, if you have a function which calls the plot
function several times, you get several pages of graphics
output.
Is there an eqivalent in R, please?
Yes, for e.g. postscript() and pdf() devices the default is to plot each
new plot on a
On Tue, 10 Jan 2006, Erin Hodgess wrote:
Dear R People:
In S Plus, if you have a function which calls the plot
function several times, you get several pages of graphics
output.
Only on a graphsheet, I believe.
Is there an eqivalent in R, please?
Pretty close, on a windows() device. Look
It is nearly the same example I wrote about in
http://www.mail-archive.com/r-help@stat.math.ethz.ch/msg54238.html. I'll print
it out again:
In the first column I have numbers from 0 to 23 (hours of a day), the second
column contains the name of a weekday (Day as factor) and the third column
-- Message transféré --
From: Uwe Ligges [EMAIL PROTECTED]
To: Erin Hodgess [EMAIL PROTECTED]
Date: Sun, 08 Jan 2006 15:09:06 +0100
Subject: Re: [R] exporting methods/classes
Erin Hodgess wrote:
Dear R People:
I'm still struggling with sending methods and classes as
Dear R-Helpers,
I want to compare the results of outputs from glmmPQL and lmer analyses.
I could do this if I could extract the coefficients and standard errors
from the summaries of the lmer models. This is easy to do for the glmmPQL
summaries, using
glmm.fit - try(glmmPQL(score ~ x*type,
you can use something like:
Vcov - vcov(lmer.fit, useScale = FALSE)
betas - fixef(lmer.fit)
se - sqrt(diag(Vcov))
zval - betas / se
pval - 2 * pnorm(abs(zval), lower.tail = FALSE)
###
cbind(betas, se, zval, pval)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D.
Jean-Christophe BOUETTE wrote:
-- Message transféré --
From: Uwe Ligges [EMAIL PROTECTED]
To: Erin Hodgess [EMAIL PROTECTED]
Date: Sun, 08 Jan 2006 15:09:06 +0100
Subject: Re: [R] exporting methods/classes
Erin Hodgess wrote:
Dear R People:
I'm still struggling with sending
On 09-Jan-06 Mark Leeds wrote:
Sometimes I print out a package
and read about it and there
are sometimes nice examples
that I would like to run myself.
Is there a way to bring them
into R from the package or
are they only meant to be typed
in manually ? If manual is the
only way,
On 10-Jan-06 Ted Harding wrote:
On 09-Jan-06 Mark Leeds wrote:
Sometimes I print out a package
and read about it and there
are sometimes nice examples
that I would like to run myself.
Is there a way to bring them
into R from the package or
are they only meant to be typed
in manually
Soon on CRAN a new version of package gstat will emerge, which
has a few minor changes and possible incompatibilities w.r.t. the
previous version(s).
The new gstat (0.9-23) now:
+ depends on sp, and uses internally with Spatial* classes from sp
if data are provided in the old-fashoned way (as
Another year another project.
Does anyone have any experience / pointers in obtaining data from an
analog / digital converter hardware such as a National Instruments data
aquisition card using R?
I really need to obtain and process this data in a real time enviroment
and I'm hoping that R already
Hi all,
I think there is something I have forgotten to emphasize. I'd like to have an
xyplot (with lattice library). A trellis graphic. For every weekday I like to
have a graphic with all hours. And inside the graphic there are the two
frequencies as lines.
The problem finally is the use
On 09-Jan-06 Mark Leeds wrote:
Sometimes I print out a package
and read about it and there
are sometimes nice examples
that I would like to run myself.
Is there a way to bring them
into R from the package or
are they only meant to be typed
in manually ? If manual is the
only way,
Dear Petr,
Thank you for your help. I have tried (and succeded) to import myfile
after conversion to text and with the help of read.table (also with the
file = 'clipboard' alternative). Both methods give correct results,
albeit dateformat turns up as a factor (minor problem). Also the
We've had a small review time within R-core on this topic,
amd would like to state the following:
--
The R-core team welcomes proposals to develop an R-wiki.
- We would consider linking a very small number of Wikis (ideally
Another year another project.
Does anyone have any experience / pointers in obtaining data from an
analog / digital converter hardware such as a National Instruments data
aquisition card using R?
I really need to obtain and process this data in a real time enviroment
and I'm hoping that R already
Dear list,
I would like to perform an analysis on the following model:
aov(ampratio ~ Type * Place * agemF + Error(speakerid/Place) ,data=aspvotwork)
using the approach from http://www.psych.upenn.edu/~baron/rpsych/rpsych.html .
Now, I got the test results, wich indicate a significant
Dear R Users,
I want to test the dimensions of an incoming vector, matrix or array safely
and succinctly. Specifically I want to check if the unknown object has
exactly 2 dimensions with a specified number of rows and columns.
I thought that the following would work:
you could use: isTRUE(all.equal(dim(obj), c(3, 5)))
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax: +32/(0)16/337015
Web:
Gregory == Gregory Jefferis [EMAIL PROTECTED]
on Tue, 10 Jan 2006 14:47:43 + writes:
Gregory Dear R Users,
Gregory I want to test the dimensions of an incoming
Gregory vector, matrix or array safely
Gregory and succinctly. Specifically I want to check if
Dear ladies and gentlemen!
When I use the plot funtion how can I change the size of the title for the x and
y axes (xlab, ylab)and the size of the axes label ?
Thank you very much.
With best regards
Claudia
__
R-help@stat.math.ethz.ch mailing list
Dear r-helpers,
Happy New Year.
To a plot
xYplot(lo ~ vaR, groups = v, data = abc1.fp,aspect = xy,
xlab=expression(frac(abs( bold(v) ),abs( bold(a) ))),
ylab = grid::textGrob(expression(paste(log, frac( italic(p) ( italic
(v) ), italic(p) ( italic(a) ),
abline
Hi everyone!!
This is my first message to the list, so I hope not to disturb anyone if the
subject of my message has been already treated.
The question is that I have a tool, a GUI made with Java, connected to R
using Rserve, and I'd like to get R and Rserve in a multi-processor machine.
Now,
Hi,
Before I reinvent the wheel I wanted to kindly ask you for your opinion if
there is a simple way to do it.
I want to merge a larger number of tables from different data sources in R
and the matching criterium are country names. The tables are of different size
and sometimes the
I'm having another problem with lmer(), this time something simpler (I
think) involving the coef() function for a model with varying
coefficients. Here's the R code. It's a simple model with 2
observations per group and 10 groups:
# set up the predictors
n.groups - 10
n.reps - 2
n -
[EMAIL PROTECTED] wrote:
Dear ladies and gentlemen!
When I use the plot funtion how can I change the size of the title for the x
and
y axes (xlab, ylab)and the size of the axes label ?
Thank you very much.
With best regards
Claudia
__
At 16:06 09/01/06, Gabor Grothendieck wrote:
[snip various earlier posts]
In addition to books, the various manuals, contributed documents and
mailing list archives, all of which one should review,
the key thing to do if you want to really learn R is to read source code
and lots of it. I think
Hello,
Why not just copy and paste the examples from the help pages? No need to
type anything.
However, if you'd like to run the entire help section for a function,
check out the example function, i.e.:
example(plot)
Robert
-Original Message-
From: Mark Leeds [mailto:[EMAIL PROTECTED]
Hi people,
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 regression, you enter the dataset and have in
summary the adjusted r-square. But
Andrew Gelman [EMAIL PROTECTED] writes:
Error in rownames-(x, value) : attempt to set rownames on object with
no dimensions
Am I doing something stupid here or is it actually a (minor) bug?
The latter. There's a drop=FALSE missing in the coef() method,
specifically in the line
val -
If they were the same you could use merge. To figure out
the correspondence automatically or semiautomatically, try this:
x - c(Canada, US, Mexico)
y - c(Kanada, United States, Mehico)
result - outer(x, y, function(x,y) mapply(lcs2, x, y))
result[] - sapply(result, nchar)
# try both which.max
RobMcG == McGehee, Robert [EMAIL PROTECTED]
on Tue, 10 Jan 2006 10:55:47 -0500 writes:
RobMcG Hello,
RobMcG Why not just copy and paste the examples from the help pages? No
need to
RobMcG type anything.
RobMcG However, if you'd like to run the entire help section for a
Using R, I believe only the Haar and Daubechies wavelets are available.
wavethresh does allow for you to use different filter numbers with the
Daubechies wavelets that I think you cannot do using waveslim. You can do
cross-validation to find the best fit across different filter niumbers and
On 10 Jan 2006, [EMAIL PROTECTED] wrote:
Hi everyone!!
This is my first message to the list, so I hope not to disturb
anyone if the subject of my message has been already treated.
The question is that I have a tool, a GUI made with Java, connected
to R using Rserve, and I'd like to get R
You can also get access to the code chunks in vignettes
as shown here:
http://tolstoy.newcastle.edu.au/~rking/R/help/05/12/17822.html
On 1/9/06, Uwe Ligges [EMAIL PROTECTED] wrote:
Mark Leeds wrote:
Sometimes I print out a package
and read about it and there
are sometimes nice
Hi everyone!!
This is my first message to the list, so I hope not to disturb
anyone if the subject of my message has been already treated.
The question is that I have a tool, a GUI made with Java, connected
to R using Rserve, and I'd like to get R and Rserve in a
multi-processor machine.
R is not thread safe, so you must not use it in a
re-entrant way.
If you want to exploit multiple processors, you can write code (e.g.,
in C) called from R (e.g., through .Call or .C) that performs
parallel/threaded computations in a thread-safe way (e.g., without
calling back into R).
Another
On 1/10/06 11:41 AM, Seth Falcon [EMAIL PROTECTED] wrote:
On 10 Jan 2006, [EMAIL PROTECTED] wrote:
Hi everyone!!
This is my first message to the list, so I hope not to disturb
anyone if the subject of my message has been already treated.
The question is that I have a tool, a GUI
(2) Is there a more simpler way to build the vector of positions (versus
using append in a loop, as above), using functions such as match() and
apply()?
How about
which.na - function(x) which(is.na(x))
?
Hadley
__
R-help@stat.math.ethz.ch mailing
Hello,
I wrote a version of which.na (similar to S's) in R that returns a
vector of positions at which NA appears in the target vector v:
which.na-function(v) {
retv - c()
for (i in 1:length(v)) {
if (is.na(v[i])) {
retv-append(i, retv, after = length(retv))
}
}
Amir,
if you try the following code:
require(wavelslim)
data(doppler)
dwt(doppler, wf=la6)
you get:
Error in switch(name, haar = select.haar(), d4 = select.d4(), mb4 =
select.mb4(), : Invalid selection for wave.filter
So, it seems that instead of la8 (default) and haar you can choose also
For errors with:
download.file()
install.package()
update.package()
any other command that R uses to access the internet
Setting up a proxy in R 2.2.1 for windows.
I work at a Central Bank so our firewall and proxy are very particular. To
allow R to access the net I did the following:
1. In
Hi,
I'm having trouble with glmmPQL from the MASS package.
I'm trying to fit a model with a binary response variable, two fixed
and two random variables (nested), with a sample of about 200,000
data points.
Unfortunately, I'm getting an error message that is difficult to
understand without
Thanks for suggestions. This is a simple question in principle, but there
seem to be some wrinkles - I am always having to think quite carefully about
how to test for equality in R. I should also have said that I would like
the check to be efficient as well safe and succinct.
One suggestion
Thanks for the nice code, Gabor!
Unfortunately, it seems not to work for my purpose, confuses lots of
countries when I compare two lists of over 150 countries each.
Do you have any other suggestions?
Gabor Grothendieck [EMAIL PROTECTED] schrieb: If they were the same you
could
To r-help:
Thanks to Mat Soukup for the help:
To draw two ablines:
xYplot(lo ~ vaR, groups = v, data = abc.fp, aspect = xy, col = c
(red,blue),
xlab=expression(frac(abs( bold(v) ),abs( bold(a) ))),
ylab = grid::textGrob(expression(paste(log, frac( italic(p) ( italic
(v) ),
See
http://en.wikipedia.org/wiki/Levenshtein_distance
http://thread.gmane.org/gmane.comp.lang.r.general/31499
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Werner
Wernersen
Sent: Tuesday, January 10, 2006 2:00 PM
To: Gabor Grothendieck
Cc:
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
Stephen Upton)
Thanks to some help by Doug Bates (and the updated version of the Matrix
package), I've refined my question about fitting nested and non-nested
factors in lmer(). I can get it to work in linear regression but it
crashes in logistic regression. Here's my example:
# set up the predictors
You can improve it somewhat by first accepting all the largest
matches and removing the rows and columns for those and
repeatedly doing that with what is left.
On 1/10/06, Werner Wernersen [EMAIL PROTECTED] wrote:
Thanks for the nice code, Gabor!
Unfortunately, it seems not to work for my
If its just succint you are after then this is slightly
shorter:
identical(dim(x)+0, c(3,5))
On 1/10/06, Gregory Jefferis [EMAIL PROTECTED] wrote:
Thanks for suggestions. This is a simple question in principle, but there
seem to be some wrinkles - I am always having to think quite
Hi,
I am trying to run a batch command on Windows XP and R CMD is not
recognized. I get the error,
'R' is not recognized as an internal or external command, operable program
or batch file.
I have C:\Program Files\R\rw2010\bin in my $PATH variable and Rcmd.exe
has been installed in that
There's a gotcha in using identical() to compare dimensions -- it also
compares names, e.g.:
x - array(1:14, dim=c(rows=3,cols=5))
dim(x)
rows cols
35
identical(dim(x)+0, c(3,5))
[1] FALSE
identical(as.numeric(dim(x)+0), c(3,5))
[1] TRUE
Gabor Grothendieck wrote:
If its just
One other thing to try could be soundex. ITs normally used for
last names but it might work here too. Google to find the
soundex encoding rules. Reviewing the country names might
suggest minor modifications to the soundex algorithm to
improve it for your case.
On 1/10/06, Gabor Grothendieck
In
http://cran.r-project.org/contrib/extra/batchfiles/
there are Windows XP batch files that you can use that will
automatically locate R using the registry and run it so that you
don't have to change your path. Just place them anywhere
in your path. See the README for more info.
I use them
This is a heads up to those using the lmer function for fitting
linear mixed models or generalized linear mixed models. In the next
few days Martin and I will upload a new version of the Matrix package
that employs a new representation of linear mixed models based on the
supernodal Cholesky
I reproduce from memory my exhaustive look into this issue.
RODBC uses the Microsoft ODBC DLL's developed by Microsoft.
These DLL's perform an automatic determination of column type based on the
contents of the first N rows of cells in each column, where N [0,16]. N may be
set in the Windows
On 1/10/06, Antje Schüle [EMAIL PROTECTED] wrote:
It is nearly the same example I wrote about in
http://www.mail-archive.com/r-help@stat.math.ethz.ch/msg54238.html. I'll
print it out again:
In the first column I have numbers from 0 to 23 (hours of a day), the second
column contains the name
Let's say I have a data.frame like
A B C TS other columns
1 1 1 12345
1 1 1 56789
1 2 1 23456
1 2 2 23457
2 4 7 23458
2 4 7 34567
2 4 7 45678
and I
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
__
R-help@stat.math.ethz.ch mailing list
Leif:
Rather than trying to mimic what you might do in SAS take advantage of R's
ability to use arbitrary data structures, e.g. lists. So, one approach is:
(your.df is the data frame)
your.list-split(your.df,your.ts[,1:3],drop=TRUE)
t(sapply(your.list,function(x)x[which.max(x$TS),]))
Cheers,
See 3.3.1 in the R FAQ (** not ** the Windows FAQ).
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
The business of the statistician is to catalyze the scientific learning
process. - George E. P. Box
-Original Message-
From: [EMAIL PROTECTED]
Thanks to Leif Kirschenbaum, Brian Ripley and Petr Pikal,
who helped with this problem. Unfortunately it appears as if the problem
with columns beginning with NAs is deeply connected to Microsoft ODBC
DLL's automatic determination of column type based on the contents of
the first N rows of
Hi all,
I need some help using read.ftable to read a contingency table. My columns
are organized as follows:
order--family--species--location--number of individuals
I couldn't figure out how to change the data on my text file to be
imported into R; and after you do that, is it possible to
I am trying to use the levelplot function from a command file.
Here is the code:
library(sp)
library(gstat)
library(lattice)
gatherData - read.table(~/gather.txt, header = TRUE)
grd = makegrid(gatherData$x, gatherData$y, cell.size = 5)
k - krige(z~x+y, ~x+y, data = gatherData, newdata = grd,
Hi R Users...
I'm trying to get the Odds Ratios for Logistic Regression Coefficients and
their confidence intervals when there are interactions, or effect modifiers, in
the regression model. May anyone help me how to get them?
Thanks
André Anderson
Brasília, Brazil
I think this is an FAQ. Anyway, levelplot() is a lattice plot and therefore
must be explicitly printed when source()ed or in a function; i.e.
print(levelplot())
is required.
-- Bert
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Matt Sakals
Hi, Matt,
See FAQ 7.22.
--sundar
Matt Sakals wrote:
I am trying to use the levelplot function from a command file.
Here is the code:
library(sp)
library(gstat)
library(lattice)
gatherData - read.table(~/gather.txt, header = TRUE)
grd = makegrid(gatherData$x, gatherData$y, cell.size =
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 when the
value in data1 is not NA. Otherwise it should contain
the values in data2.
Take the following example:
a - rnorm(100)
b - trunc(3*runif(100))
g - factor(trunc(4*runif(100)),labels=c('A','B','C','D'))
y - rnorm(100) + a + (b+1) * (unclass(g)+2)
m - lm(y~a+b*g)
summary(m)
Here b is discrete but not treated as a factor. I am interested in
I've done stuff like this before but
it's been a while and I'm stuck.
Suppose I have a matrix with one
column x and another column y
and both are numeric and let the
row index of the matrix be i
Starting at index i ( i would equal on the first iteration )
when the cumulative sum of x_i+1 - x_i
I just got 2 hits for 'RSiteSearch(significant peaks in
periodogram)', the first of which was
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/61482.html;. If you
would like more help from this group, PLEASE do read the posting guide!
www.R-project.org/posting-guide.html. Anecdotal
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 regression, you enter the dataset and have in
summary the adjusted r-square. But
Hi R-help citizens,
I'm having trouble making version 0.99-6 of Matrix on FreeBSD 6.0.
The error message is:
* Installing *source* package 'Matrix' ...
** libs
gcc -I/usr/local/lib/R/include -I/usr/local/include -D__NO_MATH_INLINES -fPIC
-g -O2 -c Csparse.c -o Csparse.o
... numerous lines
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 when the
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 be supported
by
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
are organized as
André Anderson Carvalho wrote:
Hi R Users...
I'm trying to get the Odds Ratios for Logistic Regression Coefficients and
their confidence intervals when there are interactions, or effect modifiers,
in the regression model. May anyone help me how to get them?
Thanks
André Anderson
Dear R-helpers,
While plotting a hclust tree, is it possible to mark group information on it?
Suppose I have a hclust tree ClusTree, and I use cutree to cut the
tree into different groups as following
a=cutree(ClusTree,h=10)
How do I mark the grouping info stored in a while I plot ClusTree?
I use mnp to run a multinominal probit regression model, but the summary
doesn't contain the model statistics, such as the log-likelihood and degree
of freedom, for the assessment of the goodness-of-fit of the fitted model.
Is there any way that I can generate these statistics for the fitted
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 Dependent Variables is Unordered ,why not use Multinomial
Logistic Regression.see
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, but I
haven't found anything in
Just a remark:
I need some help using read.ftable to read a contingency table. My columns
are organized as follows:
order--family--species--location--number of individuals
As the first three variables are nested, I would expect this table
to contain a lot of structural zeroes. I understand
86 matches
Mail list logo