Subject: Creating a barplot--advice needed
Advice needed: I am preparing a computer program to do a barchart. Advice
needed: Should I learn the lattice package, or try to? As a
non-statistician, much of the terminology is unfamiliar to me. grouping
variable, object of class trellis,
Hi all,
I recently updated my to R 2.6.0 and tseries package tseries version: 0.10-11.
When i was using R Version 2.3.1 (2006-06-01) with tseries 'tseries' version:
0.10-7, the code
garch(dflnRCLC1)
* ESTIMATION WITH ANALYTICAL GRADIENT *
Call:
garch(x = dflnRCLC1)
Coefficient(s):
I have got R GUI working on SuSE 10.3
The workaround is to set an environment variable that prevents the java
related error.
See http://en.opensuse.org/Xlib.lock
I tried to install 3 packages that I need through the R GUI.
When I installed cluster I noticed on R console the error messages I'm
Hi there,
I've used R before, but I just got a new computer with Leopard and I'm
having trouble installing. I download the installer, but when I reach the
end of the install process, the install button remains gray and
unclickable. I've tried both standard and custom, as well as installing on
I answered a similar question this morning.
Refer to R Mac SIG list, the mpkg installer in 2.6 binary is broken in Leopard.
Until an updated version is release, you can install R by install each
package in Packages folder one by one.
Hope it helps.
apply(matrix, 2, cumsum)
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature
and Forest
Cel biometrie, methodologie en kwaliteitszorg / Section biometrics,
methodology and
On Thu, 25-Oct-2007 at 04:25PM +0100, Omar Baqueiro wrote:
| Hello,
|
| I am trying to plot a chart similar to the typical barley yield dotplot:
| dotplot(variety ~ yield | site, data = barley,
| groups = year, auto.key = list(space = right),
| xlab = Barley Yield (bushels/acre)
This does not answer your question but here is a shorter
version of your code:
xx - xtabs(Count ~ Outcome + Product, x)
yy - sum(xx) * xx / outer(rowSums(xx), colSums(xx))
subset(as.data.frame(yy), Freq != 0)
On Nov 5, 2007 3:18 AM, francogrex [EMAIL PROTECTED] wrote:
Hi, my question is
Hi
If your barplot is simple, why not to use plain barplot from basic
graphics. Try to look at
?barplot
Petr
[EMAIL PROTECTED]
[EMAIL PROTECTED] napsal dne 04.11.2007 01:15:20:
Subject: Creating a barplot--advice needed
Advice needed: I am preparing a computer program to do a
Hello,
I have a vector containing strings of simple operation methods.
methods - c(mean,sd,median)e.g.
Now, I'd like to apply these methods to my data and I thought, I could make a
do.call. But these methods don't need any further arguments... Do I have to
create an empty argument list?
Antje wrote:
Hello,
I have a vector containing strings of simple operation methods.
methods - c(mean,sd,median)e.g.
Now, I'd like to apply these methods to my data and I thought, I could make a
do.call. But these methods don't need any further arguments... Do I have to
create an
Thank you very much! This is exactly what I need... I did not know that there
is such a simple way :)
Antje
Peter Dalgaard schrieb:
Antje wrote:
Hello,
I have a vector containing strings of simple operation methods.
methods - c(mean,sd,median)e.g.
Now, I'd like to apply these
ANUSHA AIYALU KANNAN wrote:
Respected Sir,
I am working on something in R. I am getting the following error.
Warning messages:1: the standard deviation is zero in: cor(x, y, na.method,
method == kendall)
I like to handle this error.
I like to know how to do error handling in
Hello,
I have a problem with the checking and building of R packages that contain C or
FORTRAN code.
I have implemented R 2.3.1 (2006-06-01) and the corresponding toolset under
Windows XP and Windows 2000.
The MinGW compiler components are (exactly as recommended for R 2.3.1)
?try
?tryCatch
worked for me
David Gouache
Arvalis - Institut du Végétal
Station de La Minière
78280 Guyancourt
Tel: 01.30.12.96.22 / Port: 06.86.08.94.32
-Message d'origine-
De : ANUSHA AIYALU KANNAN [mailto:[EMAIL PROTECTED]
Envoyé : dimanche 4 novembre 2007 04:50
À : [EMAIL
Please note the R posting guide: none of the helpers will still be on R
2.3.1. This looks like a bug in the shell you have 'make' set up to use:
it is not an OS-specific issue as that should be sh.exe from the Rtools
distribution (and that recommended for R 2.3.1 is no longer available).
On
On 11/04/07 16:34, Gilbert G wrote:
Dear R people:
I wish to switch from SPSS to R, but there is one particular type of
ANOVA design that cannot be done in R. Or more likely, it can be
done, but it is nowhere documented.
The problem is typical for psychologists:
You have a repeated
MR. ANDREW L. J. BARENDS
Standard Bank Plc.
Johannesburg, South Africa1
ID Number: 0092-042.
UNTERSTÜTZUNG.
Dringend GESCHAEFTSANGEBOT
Mein Name ist Andrew L. J. Barends und ich bin der
Leiter des Corporate Affairs Committee in der Standard
Bank of South-Africa PLC in Südafrika. Ich
Hello,
I try to make several zoom in chosen region.
So i use identify to indicate in which area to zoom.
i want o make it several time.
But once I zoom one time, the identify fonction gives bad results.
There is a short example, when i try to identify in Fnew it sems that there
is an offset ,
francogrex wrote:
I see that some people ask general statistics questions here sometimes. As I
understand this list is only for technical help with R-related issues. But
does anyone know of a list/forum where we can ask help/questions regarding
general statistics problems (and where solutions
Hello again,
I was wondering whether there is a command with which I can print the structure
of a list, containing list of list of ...
If you have any hint let me know...
Antje
__
R-help@r-project.org mailing list
Hello,
I have an optimization problem which I would like to solve in R but
am not sure how to approach this.
I have a cost function C which is quadratic in rx and ry (if it
helps, the cost function is separable in rx and ry) and is always
positive for all rx,ry.
I
Hi all,
I have a series of measurements made over different types of road surfaces
and capacity that I want to summarize in tables via LaTeX. Ideally I'd like
a table similar to this (
http://www.andy-roberts.net/misc/latex/tutorial4/multirow.png ) with mean/SD
of the measurement broken down
I would like to fit a power model of type:
y=b*x(^z)
to my data (having y as the response and x as the predictor variables).
I could not use linear models like y~x since R should estimate the parameter
z and not only the coefficient b.
Thus I need a package and a function! Any suggestions?
Thanks
Hi Antje,
I know you have a solution, but
l - list(list(L=1,list(L=2)))
f - function() { i=1; function(...) {cat(level, i, :, ...); i-i+1 }}
r=rapply(l, f(), ...=\n)
level 1 : 1
level 2 : 2
lets you use R scoping rules and a relatively new addition to the
'apply' family of functions. It
Dear list,
order(x,y,z) returns a permutation to order x, ties broken by y, remaining
ties broken by z. (And so on.)
What I'd like to do is
order(X), where X is a matrix (or a list or data frame if necessary) of
unspecified size, which orders X[,1], ties broken by X[,2], remaining ties
broken
System:
Linux Ubuntu 7.10 Gibbon
kernel 2.6.22-14-generic
Emacs ver. 22.1.1
ESS ver. 5.3.0
R version 2.6.0 (2007-10-03)
Colleagues
This is a follow-up note to my earlier post under this header, giving my
solution.
First, with some trepidation I
Hi All,
I'm preparing to co-teach a class on applied MCMC using R as the
computational engine in the spring. As a consequence, I'm starting
to look at MCMC textbooks. Please send along your recommendations!
I will collate the responses and post a summary, so please reply to
me
Dear R Community,
I am trying to create a heat map where, instead of normal text labels, I
will use small PNG images.
(I know there is a function heatmap in the stats package (and also a
heatmap.plus package), but it is unclear how to modify it to accomplish
the task.)
Would anyone have any
On Mon, 5 Nov 2007, Christian Hennig wrote:
Dear list,
order(x,y,z) returns a permutation to order x, ties broken by y, remaining
ties broken by z. (And so on.)
What I'd like to do is
order(X), where X is a matrix (or a list or data frame if necessary) of
unspecified size, which orders
See
?ftable
and do
RSiteSearch(latex ftable)
to find some posts relevant to what you need.
Chuck
On Mon, 5 Nov 2007, A Friedman wrote:
Hi all,
I have a series of measurements made over different types of road surfaces
and capacity that I want to summarize in tables
That solved it!
Thank you very much!
On Mon, 5 Nov 2007, Prof Brian Ripley wrote:
On Mon, 5 Nov 2007, Christian Hennig wrote:
Dear list,
order(x,y,z) returns a permutation to order x, ties broken by y, remaining
ties broken by z. (And so on.)
What I'd like to do is
order(X), where X
Here are a couple of options that may get you started:
library(TeachingDemos)
plot(1:10, 1:10, type='n')
abline(0,1)
rhist - function(x){
tmp - hist(x, plot=FALSE)
barplot(tmp$counts, horiz=TRUE, axes=FALSE, space=0)
}
subplot( rhist(rnorm(100)), 2,2, hadj=0 )
subplot(
On Mon, 2007-11-05 at 17:50 +0100, Duccio - wrote:
I would like to fit a power model of type:
y=b*x(^z)
to my data (having y as the response and x as the predictor variables).
I could not use linear models like y~x since R should estimate the parameter
z and not only the coefficient b.
Thus
Hi,
I have been trying to use the function cochran.test from the Outliers
package to test for homogeneity of variance. This works well except when
I use transformed data. Would anyone have an idea why it doesn't work
and how I could do the cochran test on transformed data?
Thanks,
Stephanie
Stephanie Bernard wrote:
Hi,
I have been trying to use the function cochran.test from the Outliers
package to test for homogeneity of variance. This works well except when
I use transformed data. Would anyone have an idea why it doesn't work
and how I could do the cochran test on
On 4/11/2007, at 1:15 PM, Thomas L Jones, PhD wrote:
Subject: Creating a barplot--advice needed
Advice needed: I am preparing a computer program to do a barchart.
Advice
needed: Should I learn the lattice package, or try to? As a
non-statistician, much of the terminology is unfamiliar
If you look at the latex output associated with the page you
posted you need to:
- eliminated all repeated elements in column 1
- replace each non-repeated element, say X, which spans 3 rows, say, to
\hline\multirow{3}{*}{X}
- add an \hline after the last row
Using the first 10 rows of the
## I would use latex() for this example.
library(Hmisc)
team.dat - '
GK Paul Robinson
LB Lucus Radebe
DC Michael Duberry
DC Dominic Matteo
RB Didier Domi
MC David Natty
MC Eirik Bakke
MC Jody Morris
FW Jamie Mcmaster
ST Alan Smith
ST Mark Viduka
'
team - read.table(textConnection(team.dat),
You need to supply both x and y to the identify function. The suprising
thing is that the first call worked.
So try:
plot(F$x,F$y)
I=identify(F$x,F$y)
Fnew=F[I[1]:I[2],]
plot(Fnew$x,Fnew$y)
I=identify(Fnew$x,Fnew$y)
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data
Benson,
Apparently it's a break with the installer for Leopard. See this:
Refer to R Mac SIG list, the mpkg installer in 2.6 binary is broken in Leopard.
Until an updated version is release, you can install R by install each
package in Packages folder one by one.
Hope it helps.
I have looked through all the ROracle documentation, R-help, R-dev, Googled
the world and still cannot figure out this issue. Everything I read it is
about connecting to an Oracle database and pulling data, processing it, and
then pushing it back to Oracle. I am running in pl/sql code so I
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
You might want to take a look at PL/R or REmbeddedPostgres
(the former being a reimplementation of the latter).
AFAIK, nobody has done this for Oracle, but we did do it
as an experiment for MySQL and Postgres and so the ideas are
there but there will
jah wrote on 11/05/2007 02:41 PM:
I have looked through all the ROracle documentation, R-help, R-dev, Googled
the world and still cannot figure out this issue. Everything I read it is
about connecting to an Oracle database and pulling data, processing it, and
then pushing it back to Oracle. I
Dear Irene,
May be I misunderstand your problem, but are you not describing here a
classical hierarchical model ?
Do you know the book recently published, entitled 'Data Analysis Using
Regression and Multilevel/Hierarchical Models', written by Gelman and Hill ?
I think the section 12 can provide
I have looked through all the ROracle documentation, R-help, R-dev,
Googled the world and still cannot figure out this issue. Everything I
read it is about connecting to an Oracle database and pulling data,
processing it, and then pushing it back to Oracle. I am running in
pl/sql code
)]
An object of class its
1
20071105 1
20071105 2
Slot dates:
[1] 2007-11-05 16:03:09 Eastern Standard Time
[2] 2007-11-05 16:03:10 Eastern Standard Time
# zoo
library(zoo)
z - zoo(rep(1:2, 5), Sys.time() + 1:10)
z[!duplicated(z)]
2007-11-05 16:03:10 2007-11-05 16:03:11
1
Dear all
I have a small routine in S-Plus that I would like to use in R. There are some
small changes I've been able to adapt, but I don't know how to translate the
line starting with g1-sort.col and the one Size- c(Size...) (I also
add some lines before)
(...)
Size-NULL
Josh Kalish wrote:
Hi,
Has anyone used R to send out an email via an SMTP server? I run R from a
script and would like to be able to send out an email from there also. I
could shell out to an email program, but was hoping maybe someone had a
package for R or a smooth way to do
Hi all,
R newbie, but have been reading the boards and learning lots. I have
read all documents list-serve responses I could about R and memory
but still no answer.
QUESTION:
When the sizes of your objects exceeds your available RAM, R switches
to Virtual Memory right? If so, why does it so
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Lux Zhang wrote:
On 06/11/2007, Adrian Dusa [EMAIL PROTECTED] wrote:
On Monday 05 November 2007, you wrote:
On 04/11/2007, Andrew Perrin [EMAIL PROTECTED] wrote:
http://www.omegahat.org/RSPerl/
Hello,
What I am trying to do is:
Generate a density plot of a population of data. This data has a bimodal
distribution so I've isolated a couple of possible sub-populations and I
want to overlay these two density plots over the first to see whether they
are contributing to the bimodal
It would help if we knew what that command was supposed to do. My guess
is that you want the Tr$size with the smallest value of Tr$dev, in which
case
Tr$size[which.min(Tr$dev)]
would seem to do it. However, there are details like what is to happen if
there are ties for the smallest value.
On 06/11/2007, Adrian Dusa [EMAIL PROTECTED] wrote:
On Monday 05 November 2007, you wrote:
On 04/11/2007, Andrew Perrin [EMAIL PROTECTED] wrote:
http://www.omegahat.org/RSPerl/
--
Andrew J Perrin - andrew_perrin
Dear all,
Want to update to 2.6 on a RedHat EL4 64 bit system. When the RPM tries to
install, I get a message that R2.5.1-dev is necessary. I don't have
2.5.1-dev installed so the 2.6 install quits. I couldn't find 2.5.1-dev at
CRAN so I'm stuck. Any help would be appreciated.
Thanks in
Hoping someone can offer me some assistance. I'm trying to execute a script
and I keep getting this error message about Error: element 12 is empty.
I'm wondering if my syntax is incorrect within legend.list. If anyone has
any suggestions to sees something obvious that I am missing, I would
Even when it does start using virtual memory (which you never want it
to do from a performance standpoint), there is still the physical
limitation (e.g., a 2-3GB on Windows depending on how you have it
configured) that you can not exceed. Also the maximum memory is set
to:
--max-mem-size=N
Not having run your script, it looks to me like you have an extra
comma after the final element of legend.list...
On 11/5/07, Patrick Richardson [EMAIL PROTECTED] wrote:
Hoping someone can offer me some assistance. I'm trying to execute a script
and I keep getting this error message about
Aha. Thanks Jim you rock. So my system is probably hitting that upper
limit of 2-3GB (a 32 bit issue I assume). Looks like I need to go 64
bit.
BTW, --max-mem-size looks like a command line command. Can I execute
something similar from within .Rprofile so that it does this each time
I start R?
Graham,
I'm not sure what the issue is, but I would suggest using OpenOffice
directly with the odfWeave package.
Max
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Graham Smith
Sent: Monday, November 05, 2007 3:02 AM
To: [EMAIL PROTECTED]
Subject: [R]
Do
?Memory-limits
and then look at
?memory.limit
I have an icon that starts R and have that as the command line
parameter. You should be able to set that up. You should be able to
use
memory.limit(2000) # set limit to 2GB
in .Rprofile.
On 11/5/07, Johan Jackson [EMAIL PROTECTED] wrote:
RUsers,
I am trying (with little success) to determine how to combine unique and
min to do some data frame manipulations and I could use a little advice.
I have tried various combinations of 'unique', 'min' and 'apply' with
no luck.
Given this simple example data (x, y)
x-c(10, 10, 12, 13,
Here is one way of doing it:
z - as.data.frame(cbind(x, y))
a - by(z, z$x, function(b) b[which.min(b$y),])
do.call('rbind', a)
x y
10 10 112
12 12 111
13 13 117
15 15 114
16 16 113
17 17 119
On 11/5/07, Bret Collier [EMAIL PROTECTED] wrote:
RUsers,
I am trying (with little success)
Greetings,
I would like to use the boot function to generate a bootstrap confidence
interval for the slope in a SLR that has a zero intercept. My attempt to do
this is shown below. Is this the correct implementation of the boot function
to solve this problem? In particular, should I be doing
I recently bought a new Intel Mac with OSX Leopard. I reinstalled the
latest version of R and the Xcode tools. I'm trying to compile a C
source code with the R compiler in X11, and I receive the following
error message:
[Macintosh:School Work 07-08/36-711/HW 3] dajmcdon% R CMD SHLIB McDonaldHW3.c
Hi,
I finally figured out where the problem is. It was because the
user account has a 10GB quota for the systerm that I used.
Once the space limit is reached, R will have this error whenever you
want to save or load.
Gosh! It took me so long to realize this.
Hongxiao
Dear all,
I have a puzzle regarding the estimation of Neg. Binomial event count
model in R. I would greatly appreciate if anyone could shed some light
on my puzzle.
Using the glm.nb command, or the zelig command developed by Gary King
et. al., I obtain the same point estimates in R as well as
Dear UseRs,
I would like to know the created time and date of specific object.
Is there any function for it?
Thank you in advance.
Sincerely,
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the
Max,
Thanks, a few people have suggested this offlist and I am having a look at
this, but I just trying to quickly get results into a word processor in a
more easily read format than a straight copy and paste.
Open Office Writer was ideal in the past, because a straight paste produced
a
Hello everybody,
I'm using SSOAP and RCurl to call webservice. I got success in calling
webservice with string as
an input. Now i want to pass array as an input(can be single or
multidimensional). When i do that i get
error.
unable to find an inherited method for function toSOAP, for
On Tue, 6 Nov 2007, Dong-hyun Oh wrote:
Dear UseRs,
I would like to know the created time and date of specific object.
Is there any function for it?
There isn't even the concept. Most objects in R are a collection of
SEXPRECs created at different times. Suppose you create a data frame out
On Mon, Nov 05, 2007 at 10:51:08AM -0700, Greg Snow wrote:
Does the following do what you want (or at least start you in the correct
direction)?
mydata - data.frame( job=c(Ambassadör,Läkare,Domare,
Professor,Advokat,Pilot,Verkställande direktör,Forskare,
Civilingenjör,Statsråd),
72 matches
Mail list logo