On Mon, 13 Jun 2005, ronggui wrote:
it takes a long time to load the lme4 package.anyone else encounter this
problem?
Yes: ca 10 secs on my machine. Why do you call it a `problem', though?
system.time(library(lme4))
ÔØÈëÐèÒªµÄ³Ì¼°ü£ºMatrix
ÔØÈëÐèÒªµÄ³Ì¼°ü£ºlattice
[1] 19.90 0.30 25.56
In a plot on an X11 device, I'd like to mark a few points with
a large thin crosshairs. By mark I mean: draw the crosshair
on the plot.
The two unsatisfactory methods that I have are:
METHOD 1:
points( x, y, pch=+ ,cex=2)
METHOD 2:
line(xa, xb, y0, y0)
line(x0, x0, ya, yb)
On Sat, 11 Jun 2005, John Fox wrote:
Dear Marc,
I get the same results -- same coefficients, standard errors, and fitted
probabilities -- from multinom() and glm(). It's true that the deviances
differ, but they, I believe, are defined only up to an additive constant:
Yes. There are many
New to R, can't afford SPSS!
Why can't I assign spss variables? Here are the details.
### I've read in my collaborator's sav file using read.spss, eg
children = read.spss(filename)
### It has many variables
length(children)
[1] 347
### and I would like to assign individual variables.
Hi ,
I am attempting to format the 'identify' function labels.
I can format the colour but the 'cex' parameter appears not to work
for me.
example--
x-1:5
y-1:5
plot(x,y)
identify(x,y,cex=0.5,col=2)
[1] 3
The label is coloured red but the 'cex=0.5' does not reduce
the label size. Why
Thank you so much for all your answers.
Papers, codes, examples, methods...THANKS A LOT! :-)
P.S. Thanks to Richard R, Berton, Gabor, Roger P, Ted H et all :-)
[EMAIL PROTECTED] wrote:
On 12-Jun-05 Berton Gunter wrote:
I assume that you know the usual procedure is to 'score'
each
On Mon, 13 Jun 2005, John Wilkinson (pipex) wrote:
I am attempting to format the 'identify' function labels.
I can format the colour but the 'cex' parameter appears not to work
for me.
example--
x-1:5
y-1:5
plot(x,y)
identify(x,y,cex=0.5,col=2)
[1] 3
The label is coloured red but the
Mike Day [EMAIL PROTECTED] writes:
New to R, can't afford SPSS!
Why can't I assign spss variables? Here are the details.
### I've read in my collaborator's sav file using read.spss, eg
children = read.spss(filename)
### It has many variables
length(children)
[1] 347
### and
Le 13.06.2005 09:22, Mike Day a crit :
New to R, can't afford SPSS!
Why can't I assign spss variables? Here are the details.
### I've read in my collaborator's sav file using read.spss, eg
children = read.spss(filename)
### It has many variables
length(children)
[1] 347
### and I
Dear:
How do you do?So sorry to bother you.
Nowadays I deal with my datum by Vis5D,but the datum I have
is in text style.I find that you had post the programs to
convert the text to V5D on the WWW.Unfortunately the links to
the R2v5d.f and v5df.h are void.Can you be so friendly to send
hi,
i am developing an application on the windows platform with the win32
api. my application has to accept a file and then call R to process
that file.
i would like to know how i can do this? is there a header file that i
can use or is there any other way to do it.
i am writting my application
Mike R wrote:
In a plot on an X11 device, I'd like to mark a few points with
a large thin crosshairs. By mark I mean: draw the crosshair
on the plot.
Try using the numerical pch symbols. '3' is a crosshair, so see what
the following do:
plot(1:10,pch=3)
plot(1:10,pch=3,cex=2)
It is actually the Matrix package that is taking so long, not lme4. It is is
extremely large and is required for use with lme4. I think Doug Bates or Duncan
Murdoch can confirm Matrix() contains more than 6000 or so lines of code. But,
it is not a problem, it just is what it is.
I should have also noted there is a huge efficiency gain when using lme4 over
nlme in that the Matrix package introduces faster matrix operations using
sparse matrix algorithms.
So, while it may take a moment or two longer to load this package, it saves you
*significant* time when fitting
ronggui a crit :
it takes a long time to load the lme4 package.anyone else encounter this
problem?
system.time(library(lme4))
Matrix
lattice
[1] 19.90 0.30 25.56NANA
version
_
platform i386-pc-mingw32
arch i386
os mingw32
On Mon, 13 Jun 2005, Doran, Harold wrote:
It is actually the Matrix package that is taking so long, not lme4.
Not so: lme4 is taking longer than Matrix:
system.time(library(Matrix))
[1] 4.17 0.11 4.28 0.00 0.00
system.time(library(lme4))
Loading required package: lattice
[1] 6.55 0.08 6.64
The key to solving your problem is that read.spss per default
gives you a *list* and not a *dataframe* (can anyone explain this
choice of default?).
So most likely wou want:
children = read.spss(filename,to.data.frame=TRUE)
attach(children)
or to get things a little more handy:
children -
Dear R-helpers,
On this day 6/12/2005 10:48 AM, Muhammad Subianto wrote:
Dear All,
Many thanks to Marc Schwartz and Gabor Grothendieck who have explained
me about using expand.grid function and clearly explain how to use
JGR.
dd - expand.grid(interface = interface, screen = screen,
On 6/13/05, Muhammad Subianto [EMAIL PROTECTED] wrote:
Dear R-helpers,
On this day 6/12/2005 10:48 AM, Muhammad Subianto wrote:
Dear All,
Many thanks to Marc Schwartz and Gabor Grothendieck who have explained
me about using expand.grid function and clearly explain how to use
JGR.
Hi:
I'm using function lmer from package lme4, and I get this message:
There were 12 warnings (use warnings() to see them)
So I checked them:
Warnings 1 to 11 said:
1: optim returned message ERROR: ABNORMAL_TERMINATION_IN_LNSRCH
in: LMEoptimize-(`*tmp*`, value = structure(list(maxIter =
I need some informations about the use of gtkgraph.
I want to plot a graph with data which are in a file.
I don't want some icon like (diamond, cercle or anything else)
I have done these successives step.
gtkgraph
windows
show data
file
On 6/13/05, CENDOYA, María Gabriela [EMAIL PROTECTED] wrote:
Hi:
I'm using function lmer from package lme4, and I get this message:
There were 12 warnings (use warnings() to see them)
So I checked them:
Warnings 1 to 11 said:
1: optim returned message ERROR:
Dear all,
I am trying to create a lattice plot which consists of 1 xyplot and 2
histograms (each for x and y).
My first try was like this:
x-rnorm(1000)
y-rnorm(1000)
xy - xyplot(y~x)
hist.x - histogram(x)
hist.y - histogram(y)
print(xy, position=c(0, 0.2, 1, 1), more=TRUE)
print(hist.x,
Hi
how do I extract those components of a list that satisfy a certain
requirement? If
jj - list(list(a=1,b=4:7),list(a=5,b=3:6),list(a=10,b=4:5))
I want just the components of jj that have b[1] ==4 which in this case
would be the first and
third of jj, vizlist (jj[[1]],jj[[3]]).
How to
maybe something like this:
jj - list(list(a = 1, b = 4:7), list(a = 5, b = 3:6), list(a = 10, b
= 4:5))
###
jj[sapply(jj, function(x) x$b[1] == 4)]
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic
Check out this recent thread:
Solution using names (except you should use the 'with' implementation
from the second link here):
https://www.stat.math.ethz.ch/pipermail/r-devel/2005-June/033439.html
Solution that traverses children directly (somewhat more open
to breaking if grid changes but
On Mon, 2005-06-13 at 15:23 +0100, Robin Hankin wrote:
Hi
how do I extract those components of a list that satisfy a certain
requirement? If
jj - list(list(a=1,b=4:7),list(a=5,b=3:6),list(a=10,b=4:5))
I want just the components of jj that have b[1] ==4 which in this case
would be
On Mon, 13 Jun 2005, BXC (Bendix Carstensen) wrote:
The key to solving your problem is that read.spss per default
gives you a *list* and not a *dataframe* (can anyone explain this
choice of default?).
The reason for the default is speed.
It's possible that the default could be changed to
All-
Does anyone on the list have experience with building RODBC from source
on a Linux box for use with DB2?
I am using (all from source):
R 2.0.1
unixODBC 2.2.9
RODBC 1.1-3
For example:
[jcole]$ R CMD INSTALL RODBC_1.1-3.tar.gz 2 rodbc.log
* Installing *source* package 'RODBC' ...
checking
I have had success by downloading the zipcode (approximate) shapefiles
(the .shp files) from: http://www.census.gov/geo/www/cob/z52000.html
Then using the maptools package (rather than the maps package).
hope this helps,
Greg Snow, Ph.D.
Statistical Data Center, LDS Hospital
Intermountain
Hi,
I have some problem when I tried to install R from source:
work as root
cd /root/dls/R-patched # this is where I unzip the file:
R-patched_2005-06-08.tar.gz
make clean
# i need my R_HOME=/usr/lib/R and I need to create libR.so for
RSPython, so I did:
./configure --prefix=/usr/lib/R
This topic is sometimes called wordprinting or stylometry. The spring
2003 issue of Chance magazine had several articles on the topic.
A colleague of mine and I have been working on a perl program (along
with various graduate students) to extract many of the common statistics
used in
Weiwei Shi [EMAIL PROTECTED] writes:
Hi,
I have some problem when I tried to install R from source:
work as root
cd /root/dls/R-patched # this is where I unzip the file:
R-patched_2005-06-08.tar.gz
make clean
# i need my R_HOME=/usr/lib/R and I need to create libR.so for
Hello,
is there any implementation of Kalman filter in R?
Thanks,
Mark
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Dear all,
Again, I would like to thank Gabor Grothendieck for your help.
I can improve which you suggest with the others combination.
And thank you for your time.
Sincerely,
Muhammad Subianto
On this day 6/13/2005 2:38 PM, Gabor Grothendieck wrote:
The pattern seems to be that each row
yep! please type
?KalmanLike
or check the dse libraries
Tom
m p [EMAIL PROTECTED] wrote:
Hello,
is there any implementation of Kalman filter in R?
Thanks,
Mark
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
Can you fix your include paths? You want the sql.h and sqlext.h from
unixODBC.
On Mon, 13 Jun 2005, John B. Cole, Ph.D wrote:
All-
Does anyone on the list have experience with building RODBC from source
on a Linux box for use with DB2?
I am using (all from source):
R 2.0.1
unixODBC
Hi,
I got an error message when I run make install:
/usr/bin/install: cannot stat `index.html': No such file or directory
make[2]: *** [install] Error 1
the full description as below:
[EMAIL PROTECTED] R-patched]# make install
make[1]: Entering directory `/root/dls/R-patched/m4'
make[1]:
Weiwei Shi [EMAIL PROTECTED] writes:
Hi,
I got an error message when I run make install:
/usr/bin/install: cannot stat `index.html': No such file or directory
make[2]: *** [install] Error 1
the full description as below:
[EMAIL PROTECTED] R-patched]# make install
make[1]: Entering
never mind:
my sys admin (he is a newbie too, now i think) tells me I should run
make install then make though i thought in another way. anyway, i
found the problem: I should run make then make install. :)
On 13 Jun 2005 18:41:58 +0200, Peter Dalgaard [EMAIL PROTECTED] wrote:
Weiwei Shi
Hi R-Help,
I have just downloaded RSPerl_0.7-0.tar.gz under SuSe.
If I type
[EMAIL PROTECTED]: ~work R CMD INSTALL -c RSPerl_0.7-0.tar.gz
i get the following error message
makedir: cannot create directory '/user/lib/R/library/00LOCK' : Permission
denied
ERROR: failed to lock directory
Dear R-help folks,
I am seeing unexpected behaviour from the function mean
with option na.rm =TRUE (which is removing a whole column of a data frame
or matrix.
example:
testcase - data.frame( x = 1:3, y = rep(NA,3))
mean(testcase[,1], na.rm=TRUE)
[1] 2
mean(testcase[,2], na.rm = TRUE)
[1] NaN
Greetings all,
I'm working on a project trying to apply fourier analysis to timestamped router
logs, using R to perform the analysis. The idea is to determine if any type of
traffic (say, outgoing ICMP requests) has strong periodic features because it
may indicate a compromise somewhere on
Jim Robison-Cox wrote:
Dear R-help folks,
I am seeing unexpected behaviour from the function mean
with option na.rm =TRUE (which is removing a whole column of a data frame
or matrix.
example:
testcase - data.frame( x = 1:3, y = rep(NA,3))
mean(testcase[,1], na.rm=TRUE)
[1] 2
Jim Robison-Cox [EMAIL PROTECTED] writes:
Summary:
If I have a data frame or a matrix where one entire column is NA's,
mean(x, na.rm=T) works on that column, returning NaN, but fails using
apply, in that apply returns NA for ALL columns.
lapply works fine on the data frame.
If you
On Mon, 13 Jun 2005 11:05:46 -0600 (MDT) Jim Robison-Cox wrote:
Dear R-help folks,
I am seeing unexpected behaviour from the function mean
with option na.rm =TRUE (which is removing a whole column of a data
frame or matrix.
example:
testcase - data.frame( x = 1:3, y = rep(NA,3))
In
I see.
So apply() first changes the dataframe to a matrix, which made it
a matrix of text values, then mean made no sense for any column, so I got
all NA's.
Thanks, Peter,
Jim
On 13 Jun 2005, Peter Dalgaard wrote:
Jim Robison-Cox [EMAIL PROTECTED] writes:
Summary:
If I have a data
On Mon, 13 Jun 2005, Werner Bier wrote:
Hi R-Help,
I have just downloaded RSPerl_0.7-0.tar.gz under SuSe.
If I type
[EMAIL PROTECTED]: ~work R CMD INSTALL -c RSPerl_0.7-0.tar.gz
i get the following error message
makedir: cannot create directory '/user/lib/R/library/00LOCK' : Permission
Werner Bier wrote:
Hi R-Help,
I have just downloaded RSPerl_0.7-0.tar.gz under SuSe.
If I type
[EMAIL PROTECTED]: ~work R CMD INSTALL -c RSPerl_0.7-0.tar.gz
i get the following error message
makedir: cannot create directory '/user/lib/R/library/00LOCK' : Permission
denied
Hi All,
Two weeks ago, I received instructions from the R list about creating
pseudodatasets by bootstrapping an existing dataset 1000 times. I received
code that bootstrapped my dataset and ran an ANOVA on each pseudodataset,
producing a histogram of F-values. Here is the code I received
Probably a simple question, but I can't find the answer to it ...
In the 'ppinit' code it describes how it takes a 'file in standard format' and
creates a point process object with it. Can anyone enlighten me as to what
this 'standard format' is? The example given doesn't allow me to view the
Dear Rs,
I have tried to use conditional expressions to calculate biomass for
different life forms (trees, lianas, and palms).
Here is an example:
lifeform
dbh form
1 30 tree
2 29 tree
3 28 tree
4 27 tree
5 26 tree
6 25 tree
7 24 tree
8 23 tree
9 22 tree
10
Dear Rs,
I have tried to use conditional expressions to calculate biomass for
different life forms (trees, lianas, and palms).
Here is an example:
lifeform
dbh form
1 30 tree
2 29 tree
3 28 tree
4 27 tree
5 26 tree
6 25 tree
7 24 tree
8 23 tree
9 22 tree
10
On Mon, 13 Jun 2005 [EMAIL PROTECTED] wrote:
Probably a simple question, but I can't find the answer to it ... In the
'ppinit' code it describes how it takes a 'file in standard format' and
creates a point process object with it. Can anyone enlighten me as to
what this 'standard format'
You could use ifelse() here:
lifeform - data.frame(dbh = c(30,29,15,14,30,29),
form=factor(c(tree, tree,
liana, liana,
palm, palm)))
lifeform$biomass -
ifelse(lifeform$form==tree,
This may be of help
dat -
data.frame(dbh = c(30,
29 , 28, 27, 26, 25, 24, 23, 22,
21, 20, 15, 14, 13, 12, 11, 10,
9, 8, 7, 6, 5, 30, 29, 28, 27,
26, 25, 24, 23, 22, 21, 20),
form = c( tree, tree, tree, tree, tree, tree,
tree, tree, tree,
Prof. Ritley, Mr. Dalgaard, thank you very much for the immediate reply
and your efforts!
We investigated a littlle further, though not directly following Mr.
Dalgaard's findings
on dynamic linking. We got aware of another short notice published on
the R-devel list,
providing some hints for
R already contains a xerbla: you should not be adding one.
On Mon, 13 Jun 2005, Stefan Sobernig wrote:
Prof. Ritley, Mr. Dalgaard, thank you very much for the immediate reply
and your efforts!
We investigated a littlle further, though not directly following Mr.
Dalgaard's findings
on
I believe that FFT is not appropriate. However Lomb-Scargle periodogram
could be used.
Sincerely
Milos Zarkovic
**
Milos Zarkovic MD, Ph.D.
Associate Professor of Internal Medicine
Institute of Endocrinology
Dr Subotica 13
11000 Beograd
Suppose I fit the following model:
library(gam)
fit - gam(y~x1+x2+s(x3),family=binomial)
and then I use
fitf$coef
x1x2 s(x3)
4.1947460 2.7967200 0.0788252
are the coefficients for x1 and x2 the estimated
coefficients? what is the meaning of s(x3)? since this
is a
The first hit in RSiteSearch(html) guides you to the package R2HTML. Take
a look at HTMLInsertGraph {R2HTML}. Try RSiteSearch(pdf) for pdf devices.
Hint: savePlot(MyPlot,pdf)
Cheers
Francisco
From: Mike R [EMAIL PROTECTED]
Reply-To: r-help@stat.math.ethz.ch
To: Francisco J. Zagmutt [EMAIL
There are lots of ways to do this, and lots of sources of
information. Go to the R website and download the Introduction to R
document from the Manuals link.
Here is one simple way. Not the most succinct in terms of coding
style, perhaps, but relatively easy to understand, and relatively
good
Although I see similar, but more complex, questions addressed in the help
archive, I'm having trouble adding superscripted text to the y-axis labels of
some figures, and I can't find anything in the R documentation on this.
I have:
ylab=BA (m2/ha)
but I want the 2 to be superscripted.
Thanks in
Dear R folks:
I am trying to use SSOAP (version 0.2-2) package in R (version
2.1.0,linux) to access SOAP service on NCBI
(http://www.ncbi.nlm.nih.gov)
its WSDL file is at http://www.ncbi.nlm.nih.gov/entrez/eutils/soap/eutils.wsdl
but some errors occured:
ncbi -
See ?plotmath.
plot(1:10, 1:10, ylab=expression(BA (m^{2}/ha)))
Benjamin M. Osborne wrote:
Although I see similar, but more complex, questions addressed in the help
archive, I'm having trouble adding superscripted text to the y-axis labels of
some figures, and I can't find anything in the R
On Mon, 13 Jun 2005, Benjamin M. Osborne wrote:
Although I see similar, but more complex, questions addressed in the help
archive, I'm having trouble adding superscripted text to the y-axis labels of
some figures, and I can't find anything in the R documentation on this.
I have:
ylab=BA
Ben
Others have pointed out plotmath. However, for some superscripts (including 2)
it may be easier to use the appropriate escape sequence (at in Windows):
ylab = 'BA (m\262/ha)'
Cheers
Peter Alspach
Benjamin M. Osborne [EMAIL PROTECTED] 14/06/05 13:42:53
Although I see
Am 13 Jun 2005 um 11:40 hat Mwalili, S. M. geschrieben:
See ?layout
[...]
Thanks to Roy, Gabor and Mwalili for their helpful suggestions. As
Gabor wrote, I have to make myself familiar with the grid package and
viewports.
Bernd
__
Hi all --
I'm seeing the following commands reliably produce a crash in RGui,
version 2.0.1, for both my home and office machine:
rm(list = ls(all = TRUE));
load(dataset.R);
library(wle);
data.wle = wle.lm(abortion ~ year * lib.con + age + gender +
+ urbanism + census + income +
69 matches
Mail list logo