On Mon, 13 Feb 2006, Deepayan Sarkar wrote:
On 2/10/06, simon chou [EMAIL PROTECTED] wrote:
Hi,
I was trying to stack a topographic map readed in from esri shapefile with a
contour map and a vector map. However, the plotMap(maptools) and
contourplot(lattice) do not seem to work well on
On Mon, 13 Feb 2006, Denis Chabot wrote:
Hi,
This is probably documented, but I cannot find the right words or
expression for a search. My attempts failed.
I have a data frame of 3 vectors (x, y and z) and would like to
transform this so that I could use persp. Presently I have
Hi R users
I have a simple question to ask
x - c(a,b,c)
x
[1] a b c
I like to have abc instead of having a b c
I tried to use paste and other functions related to character.
Is there any function (command) that enable me to combine elements of a
character vector into one character ?
Thanks
At 09:35 AM 2/10/2006, Gregor Gorjanc wrote:
Hello!
I would like to get MLE for parameter lambda of Poisson distribution. I
can use fitdistr() for this. After looking a bit into the code of this
function I can see that value for lambda and its standard error is
estimated via
estimate - mean(x)
Hi R-help members,
I have read a lot in the Archive about the Type I vs Type III sum
of square. I think I have read confusing post so
I want to have a clear idea of the problem.
Here is an example.
I have 3 groups of subjects of unequal sample size (G1 (n=7), G2
(n=7), G3 (n=4)).
for Each
probably you missed the 'collapse' argument of paste(); try this:
paste(x, collapse = )
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:
Taka Matzmoto wrote:
Hi R users
I have a simple question to ask
x - c(a,b,c)
x
[1] a b c
Try
paste(x, collapse = )
DW
I like to have abc instead of having a b c
I tried to use paste and other functions related to character.
Is there any function (command) that enable me to combine
Bernardo Rangel tura wrote:
At 09:35 AM 2/10/2006, Gregor Gorjanc wrote:
Hello!
I would like to get MLE for parameter lambda of Poisson distribution. I
can use fitdistr() for this. After looking a bit into the code of this
function I can see that value for lambda and its standard error is
On Mon, 13 Feb 2006, Roger Levy wrote:
Prof Brian Ripley wrote:
On Mon, 13 Feb 2006, Roger Levy wrote:
Hi,
Is it possible to index via a character matrix? For example, I would
like to do the following:
cont.table - table(df1$A,df1$B) # df1 a data frame with factors A and B
Hello sir:
A question about the usage of list:
The data of a list:
[[1]]
genename comp diff lwr uprp.adj
g2-g1 NewU1-R g2-g1 -242.6667 -337.9602 -147.3732 1.754308e-04
g3-g1 NewU1-R g3-g1 -266.6667 -361.9602 -171.3732 8.869357e-05
g4-g1 NewU1-R g4-g1 440.6667
It is too late now, but this would have been better
asked on the R-sig-finance list.
On 2005-09-22 there was an answer to a question
on R-help MGARCH estimation that gave the address
of a package for BEKK estimation.
There is a working paper on the Burns Statistics website
that explains how to
XinMeng wrote:
Hello sir:
A question about the usage of list:
The data of a list:
[[1]]
genename comp diff lwr uprp.adj
g2-g1 NewU1-R g2-g1 -242.6667 -337.9602 -147.3732 1.754308e-04
g3-g1 NewU1-R g3-g1 -266.6667 -361.9602 -171.3732 8.869357e-05
g4-g1
Court Arning wrote:
This is my first time making a request so please feel free to school
me on how to improve on how to make an inquiry. I am running R-1.9.1
in UNIX via a PHP script with a system command to generate a *.jpg
image from a list created by a perl script.
Terry W. Schulz wrote:
Hello,
I am new to R, so forgive my ignorance on what is probably simple.
I find package simpleboot quite useful for LOESS bootstrapping and
generation of plots.
I want to calculate the standard error for x=60 of the 100 grid points
and 50 bootstraps.
The code
Dear Colleagues,
I've been looking for a function that can perform a tree analysis, similar to
CHAID or QUEST. The key point is that in this case the variables are not binary
but nominal (10 different values), so tree (from the tree package) won't
work. Does anybody know help?
Hi Scott,
It is difficult to debug your issue without more information. Would
it be possible to email me code of a simple example?
Cheers,
Jas.
===
Jasjeet S. Sekhon
Associate Professor
On Tue, 14 Feb 2006, Markus Preisetanz wrote:
I've been looking for a function that can perform a tree analysis,
similar to CHAID or QUEST. The key point is that in this case the
variables are not binary but nominal (10 different values), so tree
(from the tree package) won't work. Does
More to the point, you are confusing multistratum AOV with single-stratuam
AOV. For a good tutorial, see MASS4 (bibliographic information in the R
FAQ). For unbalanced data we suggest you use lme() instead.
On Tue, 14 Feb 2006, WPhantom wrote:
Hi R-help members,
I have read a lot in the
Dear all,
I tried to plot cumulative hazard H(t) against t on a log-log sclale by
putting the optional function fun ='cloglog' in the plot function.
However, I came across the following error message
survcs20fit2-survfit(Surv(start,end,status)~sx,data=ALLDPcs20)
plot(survcs20fit2,fun=cloglog)
There's a periodic cubic spline smoother built into package mgcv's gam
function (although it happens not to use a B-spline basis). This can be
used in tensor product smooth surfaces. e.g.
gam(y~te(x,z,bs=c(cc,cc)))
best,
Simon
Hi..,
is there any funcrion in R to fit a periodic B-spline
On Tue, 14 Feb 2006, Markus Preisetanz wrote:
Dear Colleagues,
I've been looking for a function that can perform a tree analysis, similar to
CHAID or QUEST.
`ctree' in package `party' implements unbiased recursive partitioning in a
way inspired by CHAID. The methodological details are
Hello Erin,
Have you tried changing the font to a large, bold face font in the GUI
preferences? This may take care of the resolution issues without needing to
use power point, and give you the flexibility of a live R session.
Best,
Randy
On 2/14/06 1:46 AM, Erin Hodgess [EMAIL PROTECTED] wrote:
Diethelm Wuertz wrote:
Thanks a lot that works fine!
Next problem, if I would have my own package, and the file test.csv
would be located in the data directory
How to use the function data to get
data(test)
Also put in the data subdirectory the file test.R
with the commands to read
Thank you Berton, I will look into that method.
Roger
Berton Gunter [EMAIL PROTECTED] writes:
?unlist
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
__
R-help@stat.math.ethz.ch mailing list
Arun Kumar Saha wrote:
Dear aDVISOR,
Hope I am not disturbing you. Can you tell me how I can perform Multivariate
Garch analysis in R. Also please, it is my humble request let me know some
resource materials on Multivariate Garch analysis itself.
You could tryhelp.search(garch) or
Dear Mr Graves,
this seems work on my data. Thank you very much for your help.
Sincerely
Michael Coeurdassier
Spencer Graves a écrit :
The following appears to be an answer to your question, though
I'd be pleased to receive critiques from others. Since your example
is NOT self
I defined three functions:
fun0 - function( x=1:5, y=1:5, ... ) xyplot( y ~ x, ... )
fun1 - function( x=1:5, y=1:5, ... ) fun2( y ~ x, ... )
fun2 - function( ... ) xyplot( ... )
The call of fun0() works as expected.
The call of fun1() causes the following error:
'Error in eval(expr,
On Feb 14, 2006, at 1:46 AM, Erin Hodgess wrote:
I'm using R in a time series class. ... I have decided to put
together Power Point slides for the teaching. ... I am currently
saving the R screen as WMF files and inserting them into
PowerPoint. While this works, it seems that there
On 2/14/2006 8:56 AM, Wolfram Fischer wrote:
I defined three functions:
fun0 - function( x=1:5, y=1:5, ... ) xyplot( y ~ x, ... )
fun1 - function( x=1:5, y=1:5, ... ) fun2( y ~ x, ... )
fun2 - function( ... ) xyplot( ... )
The call of fun0() works as expected.
The call of fun1()
Dear all
I imported a Stata .dta file with the read.dta-function from the
foreign-package. The dataframe's dimensions are
dim(d.apc)
[1] 15806 1300
Importing needs up to 15 min and calculations with these data are rather
slow (although I subset the data before starting analyses).
My
On 2/14/2006 9:08 AM, Michael Kubovy wrote:
On Feb 14, 2006, at 1:46 AM, Erin Hodgess wrote:
I'm using R in a time series class. ... I have decided to put
together Power Point slides for the teaching. ... I am currently
saving the R screen as WMF files and inserting them into
On 2/14/06, Duncan Murdoch [EMAIL PROTECTED] wrote:
On 2/14/2006 8:56 AM, Wolfram Fischer wrote:
I defined three functions:
fun0 - function( x=1:5, y=1:5, ... ) xyplot( y ~ x, ... )
fun1 - function( x=1:5, y=1:5, ... ) fun2( y ~ x, ... )
fun2 - function( ... ) xyplot( ... )
The
Dear All,
I found Andy Liaw's suggestion about using a NUMA (instead of SMP) kernel when
running R on amd64 with 1 CPU
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/35109.html
A couple of questions:
1. Is this still the case with the newer dual-core opterons (e.g., the 275 et
al.,
On 2/14/2006 9:38 AM, Gabor Grothendieck wrote:
On 2/14/06, Duncan Murdoch [EMAIL PROTECTED] wrote:
On 2/14/2006 8:56 AM, Wolfram Fischer wrote:
I defined three functions:
fun0 - function( x=1:5, y=1:5, ... ) xyplot( y ~ x, ... )
fun1 - function( x=1:5, y=1:5, ... ) fun2( y ~ x, ... )
I think it is well worth the effort to start using a database system like e.g.
MySql for such purposes.
If you look at http://gbi.agrsci.dk/~sorenh/misc/R-SAS-MySql/R-SAS-MySql.html
then you'll find a short - and rudimentary - description of how to use MySql in
connection with R and SAS (on
do you want to use the dataset as it is? if not, it might be better save the
data into a database, such as sqlite, and then odbc to it? Then you can use
sql to only fetch the data you want for analysis.
On 2/14/06, Christian Bieli [EMAIL PROTECTED] wrote:
Dear all
I imported a Stata .dta file
Soren's suggestion is the right way to go provided you don't need all the
data all the time. Another thing to try is once the data is imported,
convert the numeric part of the data frame to a matrix, as calculations on a
matrix are much faster than calculations on a data frame.
I'm too lazy to
Thank you for the perfect solution.
Uwe Ligges wrote:
Terry W. Schulz wrote:
Hello,
I am new to R, so forgive my ignorance on what is probably simple.
I find package simpleboot quite useful for LOESS bootstrapping and
generation of plots.
I want to calculate the standard error for x=60
I want to numerically maximize a function with optim (maximization over several
arguments). optim() needs a function which returns a scalar only. However, it
could be nice to be able to take other things out from the function as well.
I'tried to create an attribute to the scalar with what I
Thanks Brian for the reference.
I just discover that it is available in our
library so I going to take it read it soon.
Actually, I don't even know the difference
between a multistratum vs a single-stratum AOV. A
quick search on google returned me the R materials so that I imagine
that these
Hello,
I have some Inf values in a matrix, but I don't want to replace them with
some value but rather remove the rows that contain the Inf values. Also I
would like to record the rows which were removed. Is there an easy way to
do this instead of writing loops over the matrix? Thanks.
Ita
Erin,
From an Rgui graphics window (windows() device), Ctrl-W will save the
current graph to the clipboard as a metafile; Ctrl-C will save as a
bitmap. In PPt, Ctrl-V will paste either into a blank spot on a slide.
The metafile is a Windows vector spec that will be sharper and smaller.
Søren Højsgaard wrote:
I want to numerically maximize a function with optim (maximization
over several arguments). optim() needs a function which returns a
scalar only. However, it could be nice to be able to take other
things out from the function as well. I'tried to create an attribute
to
Hi Duncan,
On Feb 14, 2006, at 9:26 AM, Duncan Murdoch wrote:
On 2/14/2006 9:08 AM, Michael Kubovy wrote:
On Feb 14, 2006, at 1:46 AM, Erin Hodgess wrote:
I'm using R in a time series class. ... I have decided to put
together Power Point slides for the teaching. ... I am currently
On 2/14/06, Duncan Murdoch [EMAIL PROTECTED] wrote:
On 2/14/2006 9:08 AM, Michael Kubovy wrote:
On Feb 14, 2006, at 1:46 AM, Erin Hodgess wrote:
I'm using R in a time series class. ... I have decided to put
together Power Point slides for the teaching. ... I am currently
saving the R
Erin Hodgess [EMAIL PROTECTED] wrote in message
news:[EMAIL PROTECTED]
I am currently saving the R screen as WMF files and inserting them
into PowerPoint. While this works, it seems that there
might be a simpler method.
Does anyone have any suggestions for the Power Point, please?
Instead
On 2/14/2006 11:12 AM, Henrik Bengtsson wrote:
On 2/14/06, Duncan Murdoch [EMAIL PROTECTED] wrote:
On 2/14/2006 9:08 AM, Michael Kubovy wrote:
On Feb 14, 2006, at 1:46 AM, Erin Hodgess wrote:
I'm using R in a time series class. ... I have decided to put
together Power Point slides for the
Erin,
Another option. I'm saving R screens as postscript files, opening them
in Ghostview and copying and pasting to blank slides in Powerpoint.
I can't detect any deterioration in quality. Convenient when you have a
lot of slides to make in a hurry. Ghostview and
Ghostscript (required to
try the following:
mat - matrix(rnorm(40), 10, 4)
mat[sample(20, 5)] - Inf
mat
##
index - rowSums(!is.finite(mat)) = 1
mat. - mat[!index, ]
which(index)
mat.
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic
WPhantom [EMAIL PROTECTED] writes:
Thanks Brian for the reference.
I just discover that it is available in our
library so I going to take it read it soon.
Actually, I don't even know the difference
between a multistratum vs a single-stratum AOV. A
quick search on google returned me
Hi
see
?is.infinite
?which
something like
y - x[-which(is.infinite(x), arr.ind=T)[,1],]
y.removed -x[which(is.infinite(x), arr.ind=T)[,1],]
shall make the trick.
On 14 Feb 2006 at 17:01, [EMAIL PROTECTED] wrote:
To: r-help@stat.math.ethz.ch
From: [EMAIL
[EMAIL PROTECTED] wrote:
Hello,
I have some Inf values in a matrix, but I don't want to replace them with
some value but rather remove the rows that contain the Inf values. Also I
would like to record the rows which were removed. Is there an easy way to
do this instead of writing
On 2/14/06, Duncan Murdoch [EMAIL PROTECTED] wrote:
On 2/14/2006 9:08 AM, Michael Kubovy wrote:
On Feb 14, 2006, at 1:46 AM, Erin Hodgess wrote:
I'm using R in a time series class. ... I have decided to put
together Power Point slides for the teaching. ... I am currently
saving the R
Hi Ross --
Not a direct answer to your question, but here's an idea.
I guess you've got some script in a file batch.sh
#! /bin/bash
R --no-save --no-restore --gui=none `hostname` 21 BYE
# some mpi commands
# other mpi commands
BYE
and you do something like
mpiexec -np 10 batch.sh
Instead,
I had the problem to do this.
__
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
Hi
it appears to me that read.table is very slow for reading large data files
(mine are 200,000 rows). Is there a better way?
Thanks!
Max
--
Maximilian O. Kauer, Ph.D.
Department of Genetics, White lab
333 Cedar St, NSB 386
PO Box 208005
New Haven, CT 06510
?scan is much faser. Also, read.table has a colClasses optional argument
which can be used to speed up the reading of large files significantly.
read.table has a pretty good help section well worth reading.
read.table(file, header = FALSE, sep = , quote = \', dec = .,
Dear all,
I would be pleased if anyone could help me.
The Rhelp description for the figs parameter is a two-element vector
describing
the number of rows and colunns in a screen matrix.
So, why does my code (below) produce a 2x1 screen matrix instead of
a 1x2 one?
Thanks in advance,
rodrigo.
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to
announce 2-day R/S-plus Fundamentals and Programming
Techniques in San Francisco: www.xlsolutions-corp.com/Rfund.htm
Raleigh,April 3-4
Miami, April 10-11
Houston,
List,
Is it possible to add a color legend to a heatmap , similar to the one in
levelplot and filled.contour plot. The legend will represent the colors used in
the heatmap along with values for each color range. Can this be done? Please
help.
Thanks,
Svakki.
On Tue, 14 Feb 2006, Rodrigo Tsai wrote:
Dear all,
I would be pleased if anyone could help me.
The Rhelp description for the figs parameter is a two-element vector
describing
the number of rows and colunns in a screen matrix.
So, why does my code (below) produce a 2x1 screen matrix
Dear Prof Brian Ripley,
Many thanks for your help. Before trying to split.screen (1x2) I had set the
mar parameter in par( ) in order to move each graph of another (2x1)
splitted figure. This setting might have caused the problem. Now it works.
Best regards,
Rodrigo.
On 2/14/06, Prof Brian
You may find that ?layout is an easier way to handle multiple plots on the
screen (I do), if you are not already aware of it.
-- 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.
I need to install several (too many) packages from local *.zip
files.
is there any form to do it without clicking?
I'm asking for a R code that allow me perform these task.
---
Mario Alfonso Morales Rivera.
Profesor Asistente.
Departamento de Matemáticas y
On Tue, 14 Feb 2006, Mario Alfonso Morales Rivera wrote:
I need to install several (too many) packages from local *.zip
files.
is there any form to do it without clicking?
I'm asking for a R code that allow me perform these task.
See ?install.packages.
You seem to be using Windows (you
I need to install several (too many) packages from local *.zip
files.
Are you aware of being able to select multiple .zip files at one time
using SHIFT-click or CTRL-click?
is there any form to do it without clicking?
I'm asking for a R code that allow me perform these task.
You could
Is it assumption Tobit regression like OLS ? :
1. Normal
2. No autocorrelation
3. Homogeneity of variance
4. No Multicolinearity
please make sure me !
I need answer because i want to Judisium/ kompre/thesis test
__
R-help@stat.math.ethz.ch mailing list
On Wed, 15 Feb 2006 05:15:28 +0700 bambang pramono wrote:
Is it assumption Tobit regression like OLS ? :
1. Normal
2. No autocorrelation
3. Homogeneity of variance
4. No Multicolinearity
please make sure me !
I need answer because i want to Judisium/ kompre/thesis test
This is the third
how do i perform the stuff i just recorded?!
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Hello,
I'm trying to create a simple network diagram like that on
http://www.pauck.de/marco/photo/infrared/comparison_of_films/diagram3.gif.
Unfortunately I've not yet found any function that could do this.
I think this should not be a difficult task since I only need two
vectors and perhaps
Dear R-helpers,
There's a method named as in this message's subject that has been proposed
by:
Chen Y, Jackson DA, Harvey HH (1992) A comparison of von Bertalanffy and
polynomial functions in modelling fish growth data. Canadian Journal of
Fisheries and Aquatic Sciences 49:1228-1235.
which is
Hi R users:
I am not sure how do I fix it.
But I erase an old R library installation (2.2.0), and also I erase all
file directories on the
R.2.2.1 library directory, and then I try library(Rcmdr) again and it
works!.
Any way, how could I trace or debug, so I get a hint whats going on
Check out ?stars
On 2/14/06, Rainer Hahnekamp [EMAIL PROTECTED] wrote:
Hello,
I'm trying to create a simple network diagram like that on
http://www.pauck.de/marco/photo/infrared/comparison_of_films/diagram3.gif.
Unfortunately I've not yet found any function that could do this.
I think
Hi. I was wondering if anyone might have some suggestions about how I can
overcome a problem of iteration limit reached without convergence when
fitting a mixed effects model.
In this study:
Outcome is a measure of heart action
Age is continuous (in weeks)
Gender is Male or Female (0 or
Dear Kenneth,
It's very difficult for me to tell from the information that you've given
what the problem might have been. I'm glad, however, that things are working
now.
One thought: Do you have the R_LIBS environment variable set? If so, it's
possible that you were picking up an old version of
75 matches
Mail list logo