Using data frame CO2, built into R, try this:
library(lattice)
p - function(..., col) {
panel.xyplot(..., type = p, col = col[1])
panel.xyplot(..., type = smooth, col = col[2])
}
xyplot(conc ~ uptake | Treatment, CO2, panel = p, col = 1:2)
On Mon, Mar 2, 2009 at 9:13 AM, Paul
Hi,
type
is(sampA)
is(sampV)
in the R-prompt. It will show you that at least one of the vectors is not a
vector of numerical values in R's sense. Then you have to convert using the
as.numeric() function, but you will have to make sure that this really gives
you the original vectors. Since
Vivienne Ozohili wrote:
I am currently trying to implement a boostrapping approach to validation
but don't know where to start. Help please.
Take a look at Prof. Harrell's ?validate.lrm in package Design.
Regards, Mark.
Vivienne_O.Ozohili wrote:
Hi,
I was wondering whether this
Whatever it is, it is also happening on my Acer, Suse 11.1, 32 bits,
self build R.
Kees
Bo Zhou wrote:
Thanks Duncan.
I haven't got a chance to try it in windows yet. But I won't be surprised if
it's my driver's fault. I had OpenGL problems on this dell laptop before.
Date: Sun,
Hello everybody,
I'm running a panel regression with dummy variables controlling for years.
While this works fine for most of my models, in some cases I get the error
message In sqrt(diag(object$vcov)) : NaNs wurden erzeugt and missing
results for some or all of my year dummies:
Hi Juliet,
Have a look at the last section in http://had.co.nz/ggplot2/book/polishing.pdf
Hadley
On Mon, Mar 2, 2009 at 12:35 PM, Juliet Hannah juliet.han...@gmail.com wrote:
Hi, Here are three plots:
library(ggplot2)
data(diamonds)
randind - sample(nrow(diamonds),1000,replace=FALSE)
Hi,
assuming the assumption of Poisson distribution is true (or at least
approximately true), you can use MLE as mentioned by Roland using
(fitdistr of MASS, mle of stats4, MLEstimator of distrMod, probably more).
Using package distrMod you can compute minimum distance estimators (cf.
Juliet Hannah juliet.hannah at gmail.com writes:
Hi, Here are three plots:
(How to layout on one page)
library(ggplot2)
data(diamonds)
randind - sample(nrow(diamonds),1000,replace=FALSE)
dsmall - diamonds[randind,]
a= qplot(carat, data=dsmall, geom=histogram,binwidth=1)
b= qplot(carat,
Hi, I've been trying to prepare some crosstables for some survey
questions for a client. I have been using the CrossTable function in the
gmodels package. However, this command only seems to be able to create
CrossTables in text documents.
I've been trying to use odfTable in odfweave to
On 3/2/2009 2:00 PM, chaogai wrote:
Whatever it is, it is also happening on my Acer, Suse 11.1, 32 bits,
self build R.
Versions of R and rgl?
Duncan Murdoch
Kees
Bo Zhou wrote:
Thanks Duncan.
I haven't got a chance to try it in windows yet. But I won't be surprised if
it's my driver's
Oops, sorry
R version 2.8.1 (2008-12-22)
i686-pc-linux-gnu
locale:
LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=C;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8;LC_IDENTIFICATION=C
I usually have a function like this:
vplayout - function(x, y) viewport(layout.pos.row=x, layout.pos.col=y)
draw4 - function(pngname, a,b,c,d,w,h) {
png(pngname,width=w, height=h)
grid.newpage()
pushViewport(viewport(layout=grid.layout(2,2) ) )
print(a, vp=vplayout(1,1))
Hi June:
As don't have your data (yogurtnp.csv), I cannot replicate your problem.
What does
fr( c(0,0,0,1,0,-25,-0.2) )
return?
Arne
On Monday, 2. March 2009 11:55:56, June Wong wrote:
Dear Rhelpers
I have the problem with initial values, could you please tell me how to
solve it? Thank
I would appreciate any help for the following question. How to increase the
font size of labels for the x and y axes in the following code?
-
xyplot(lbxglu~lbxgh|eth,reg.dat.5, panel = p, col = 1:2,
par.strip.text = list(cex = 1.25,col=brown),
On Mon, 2009-03-02 at 11:23 -0600, maxa0...@umn.edu wrote:
I'm trying to create a bargraph of means with standard error bars using the
function bargraph.CI (in the sciplot package). Like this:
bargraph.CI(x.factor, response,data,xlab, ylab, par(family=serif),font=11)
However, an error
Try this with the built in CO2:
xyplot(uptake ~ conc | Treatment, CO2,
xlab = list(A1C, cex = 2),
ylab = list(Fasting Plasma Glucose, cex = 2))
On Mon, Mar 2, 2009 at 3:47 PM, Veerappa Chetty chett...@gmail.com wrote:
I would appreciate any help for the following question. How to increase
PJ,
Hi, I've been trying to prepare some crosstables for some survey questions
for a client. I have been using the CrossTable function in the gmodels
package. However, this command only seems to be able to create CrossTables
in text documents.
I've been trying to use odfTable in odfweave to
Hi there,
How to plot exponential axis (i.e. equally spaced .1, 1, 10, 100..)?
I want to use it with xyplot; it would be great if there is a detailed example.
Thanks a lot!
Larry
__
R-help@r-project.org mailing list
Install the Design and Hmisc package then do:
library(Design)
?validate.lrm
?calibrate
Resampling methods such as the bootstrap and cross-validation assume you
have done no model/variable selection outside the bootstrap look, i.e.,
you let the bootstrap repeat all modeling steps that used Y.
On 3/2/2009 2:48 PM, chaogai wrote:
Oops, sorry
R version 2.8.1 (2008-12-22)
i686-pc-linux-gnu
locale:
Hi,
I have a problem in which i am having trouble vectorizing a snippet of
my code to avoid a for loop. A basic example of what I'm trying to do
is below.
#simple example of problem
x = c(1,1,3,3,6,7,8,8,9)
# I create a vector that contains the unique values of x
y = unique(x)
y =
Alternatively, gives a looked in the package bpca and biplotGUI, both
available on CRAN.
HTH,
JCFaria
japal wrote:
Hello,
Something incredible (at least for me) has happen. Yesterday night I
downloaded biplot.R to edit this function and add new features I wished.
Namely I wanted to
Tinn-R was programed (intentionally) to detect any opened windows with the
word
'Console' inside of the caption as a R instance. So, for while, it is
necessary to close others
process that has the word 'Console' in the caption.
It is because the necessity of Tinn-R detects others GUIs (JCR,
Tinn-R was programed (intentionally) to detect any opened windows with the
word
'Console' inside of the caption as a R instance. So, for while, it is
necessary to close others
process that has the word 'Console' in the caption.
It is because the necessity of Tinn-R detects others GUIs (JCR,
Hello Werner,
Did you already try to use the new version of Tinn-R (2.2.0.2) and R in USB
stick?
I use it without problem!
A good idea is to have a backup of your system configuration also in the USB
stick.
Unzip it (in any place) and replace the folder Tinn-R where Tinn-R stores
all ini files
Dear Juan Carlos,
x = c(1,1,3,3,6,7,8,8,9)
table(x)
x
1 3 6 7 8 9
2 2 1 1 2 1
z-table(x)
z
x
1 3 6 7 8 9
2 2 1 1 2 1
as.vector(z)
[1] 2 2 1 1 2 1
Take a look at ?table and ?vector for more information.
HTH,
Jorge
On Mon, Mar 2, 2009 at 4:41 PM, Juancarlos Laguardia
On Mon, 2009-03-02 at 13:27 +0100, Johannes Elias wrote:
Dear R-Gurus,
I wonder why 'density' values as shown in hist or plot(density(x)) are
sometimes over 1. How can that be?
Example
hist(rnorm(1000,sd=.5),freq=FALSE)
The resulting plot shows density values below 1 on the y-axis.
Ok, thank you for your guidelines. Additionally, I have had to put the
biplotes.princomp function in the first place of the script, just after
biplotes - function(x, ...) UseMethod(biplotes). By doing this it works,
althrough I don´t really know why. I hope to acquire higher skills in R
Hi all,
Is there any way we can to power analysis for prop trend test? Many thanks!
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and
Hi list,
I want to plot state Texas using map function. I wrote the following code to
plot the grid boundary for Texas:
library(maps)
require(mapproj)
longlatLimit-c(-106.65, -93.53 , 25.93 , 36.49)
par(plt=c(0,1,0,1),cex=1,cex.main=1) #Set plotting
Hi,
My first post here and new to R so please bear with me (long time programmer
though, helping a friend with some scripts).
I've noticed a behaviour when using rownames() that I think is odd,
wondering if I'm doing something wrong.
To illustrate, say I create a very simple matrix (called
Hi list,
When I type my question in R console using ? sign (example: ?par, when I
want to query for par), the following error message popps up:
Error in
print.help_files_with_topic(C:/PROGRA~1/R/R-24~1.1/library/maps/chm/map) :
CHM file could not be displayed
I appreciate if any body
Hi all,
I need to run a program (asreml) thousands of times and each time I have
to provide a slightly different dataset. Because I have to run asreml
under Windows (DOS or scripts) environment, I have trouble to pass a
command variable (or pointer variable or %counter in the following
THe problem is that 'rownames(fred)' is NULL and therefore you will
have to define the names of all the rows:
fred-matrix(,4,2)
fred
[,1] [,2]
[1,] NA NA
[2,] NA NA
[3,] NA NA
[4,] NA NA
rownames(fred)
NULL
rownames(fred) - c(Apple,'','','')
fred
[,1] [,2]
Apple
Hi,
At the 1st step, you only assign a name for the 1st row. However, fred has 4
rows which mean you need to assign 4 rownames for these 4 rows.
At the 2nd step, you only CHANGE the 1st rowname from 1 to APPLE
If I am wrong please correct me.
Thanks
-Original Message-
From:
On Mon, 2009-03-02 at 23:15 +, Dr. Alireza Zolfaghari wrote:
Hi list,
I want to plot state Texas using map function. I wrote the following code to
plot the grid boundary for Texas:
library(maps)
require(mapproj)
longlatLimit-c(-106.65, -93.53 ,
Breslow Day has a nice three page discussion in volume 2 of their
Statistical Methods in Cancer Research. See pages 285-7. Most of the
gain in power comes from the decrease in degrees of freedom and only
if the trend is approximately linear. Alternatives that are quadratic
are not well
I am running a binomial glm with response variable the no of mites of two
species y-cbind(mitea,miteb) against two continuous variables (temperature
and predatory mites) - see below. My model shows overdispersion as the
residual deviance is 48.81 on 5 degrees of freedom. If I use quasibinomial
Hello,
first of all: I'm new to R and have only used SPSS befor this (which
can't do this at all...).
I'm trying to output some regression results to latex. The regressions
are normal OLS and I'm trying to output the results with standardized
regression coefficients and t-statistics based
Jan Schulz wrote:
Hello,
first of all: I'm new to R and have only used SPSS befor this (which
can't do this at all...).
I'm trying to output some regression results to latex. The regressions
are normal OLS and I'm trying to output the results with standardized
regression coefficients and
On Mar 2, 2009, at 9:57 AM, Brajkovic J. wrote:
spec=garchSpec(model=list(ar= 0.440270860,
omega=0.000374365,alpha=0.475446583 , mu=0, beta=0))
sim-garchSim(spec, length(dp)-1,extended=T)
I get:
spec=garchSpec(model=list(ar= 0.440270860,
omega=0.000374365,alpha=0.475446583 , mu=0,
Hi,
Using the following code I can generate the attachment Original.pdf:
latex(test, file=C:/Documents and Settings/darej/My
Documents/James/READY.tex,
title = ,
rowname = NULL,
cgroup = c(Fishstock, stocknames),
n.cgroup = c(1, rep(2,nostocks)),
colheads =
on 03/02/2009 03:11 PM Max Kuhn wrote:
PJ,
Hi, I've been trying to prepare some crosstables for some survey questions
for a client. I have been using the CrossTable function in the gmodels
package. However, this command only seems to be able to create CrossTables
in text documents.
I've
Dear R-Help,
I have written a function, which in simplified format can be represented as:
temp.fn - function(my.mean,my.sd){
Parameters - list(mean = my.mean, sd = my.sd)
curve(do.call(dnorm,c(list(x), Parameters)), from = my.mean-my.sd,
to = my.mean+my.sd)
}
This works as I want it do.
I was able to extract RGoogleDocs with 7-zip. That was when I began to
appreciate the gravity of what Prof Ripley was saying. I am on a windows xp
machine and the RGoogleDocs download does not have a windows binary version.
Instead it has only a source package. I started using the R-admin manual
Werner Wernersen pensterfuzzer at yahoo.de writes:
Hi,
I have been dreaming about a complete R environment on my USB stick for a
long time. Now I finally want to
realize it but what I am missing is a good, portable editor for R which has
tabs and syntax highlighting, can
execute code,
Download and install each of R, Rtools and MiKTeX.
Each of these have automated installers so its just a matter of
pressing Enter repeatedly.
You don't have to change your path if you place Rcmd.bat from
http://batchfiles.googlecode.com
anywhere on your path (so it can be found). To see your
On 3/03/2009, at 2:46 PM, Rob Foxall wrote:
Dear R-Help,
I have written a function, which in simplified format can be
represented as:
temp.fn - function(my.mean,my.sd){
Parameters - list(mean = my.mean, sd = my.sd)
curve(do.call(dnorm,c(list(x), Parameters)), from = my.mean-my.sd,
On Mon, Mar 2, 2009 at 7:46 PM, Rob Foxall rfoxal...@gmail.com wrote:
Dear R-Help,
I have written a function, which in simplified format can be represented as:
temp.fn - function(my.mean,my.sd){
Parameters - list(mean = my.mean, sd = my.sd)
curve(do.call(dnorm,c(list(x), Parameters)),
Try this:
temp.fn.2 - function(mean = 0, sd = 1, n = 101) {
x - seq(mean - sd, mean + sd, length = n)
plot(x, dnorm(x), type = l)
}
# test
temp.fn.2()
On Mon, Mar 2, 2009 at 8:46 PM, Rob Foxall rfoxal...@gmail.com wrote:
Dear R-Help,
I have written a function, which in simplified
Regards,
Xianming
From: Xianming Wei
Sent: Tuesday, 3 March 2009 9:22 AM
To: 'r-help@r-project.org'
Subject: how to pass a command variable in DOS to R program in R CMD
BATCH
Hi all,
I need to run a program (asreml) thousands of times and each time I have
to provide a slightly
Create a file test.R containing this single line:
print(commandArgs(TRUE))
and then try this:
Rscript test.R a b c
On Mon, Mar 2, 2009 at 6:22 PM, Xianming Wei x...@bses.org.au wrote:
Hi all,
I need to run a program (asreml) thousands of times and each time I have
to provide a
Thanks for the plug on NppToR. Yes it is portable, but a few of the
features don't work. I have it on my plan to have a launcher for the
portable apps menu. Also the website for NppToR is now
https://sourceforge.net/projects/npptor/
On my personal website I also keep, semi-up-to-date a
This does not really address my point. Yes, if the few nerds who want to do
funny stuff are the ones making the purchase, then there is a good chance (but
still not guaranteed) that they will get IML, but do all companies that buy SAS
actually think about that, or do they just see the extra
I'm using R 2.8.1 on Ubuntu 8.10. I'm writing partly to ask what's
wrong, partly to tell other users who search that there is a work
around.
The General Social Survey is a long standing series of surveys
provided by NORC (National Opinion Research Center). I have
downloaded some years of the
Here is one way:
thestring - the-actual-string
plot(0:1,0:1,type=n)
text(x=0.5, y=0.5, labels=bquote(bold(.(thestring
There is an example in the plotmath help page that shows using bquote (but that
page has so much information that I did not find it on my first scan through
even knowing
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-
project.org] On Behalf Of Prof Brian Ripley
Sent: Monday, March 02, 2009 12:38 AM
To: tedzzx
Cc: r-help@r-project.org
Subject: Re: [R] rounding problem
I nominate the following for the fortunes
for an inefficient language , it sure has dominated the predictive
analytics world for 3 plus decades.
I referred once to intellectual jealousy between newton and liebnitz.
i am going ahead and creating the R package called Anne.
It basically is meant only for SAS users who want to learn R ,
Hello all,
I'm running R2.8.1 on a Mac OS 10.4.11. While trying to install the package
gdata, I was presented with the following (error at end of report):
R install.packages(gdata)
also installing the dependency ‘gtools’
trying URL
Hi all,
I am trying to run rejection sampling for the quantity z11 in the function
below. Unfortunately I can't simplify the function further so that z11 only
appears once.
Whenever I run the algorithm, R looks as if it is running it (no error
messages or anything), but then nothing happens
can someone show me how to use a regular expression to break the string
at the bottom up into its three components :
(-0.791,-0.263]
(-38,-1.24]
(0.96,2.43]
I tried to use strplit because of my regexpitis ( it's not curable. i've
been to many doctors all over NYC. they tell me there's no cure
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of ekwaters
Sent: Monday, March 02, 2009 9:39 PM
To: r-help@r-project.org
Subject: [R] R - need more memory, or rejection sampling
algorithm doesn't work?
Hi all,
I am
Sun, Larry jingrans at amgen.com writes:
How to plot exponential axis (i.e. equally spaced .1, 1, 10, 100..)?
I want to use it with xyplot; it would be great if there is a detailed
example.
Thanks a lot!
Sound like what is usually called a logarithmic axis. See parameter log
in scales of
Thanks Kjetil. This is exactly what I wanted.
Hardi
From: Kjetil Halvorsen kjetilbrinchmannhalvor...@gmail.com
Cc: r-help r-help@r-project.org
Sent: Monday, March 2, 2009 9:45:43 PM
Subject: Re: [R] Large data set in R
install.packages(biglm, dep=TRUE)
G'day Mark,
On Tue, 03 Mar 2009 00:16:34 -0600 (CST)
markle...@verizon.net wrote:
can someone show me how to use a regular expression to break the
string at the bottom up into its three components :
(-0.791,-0.263]
(-38,-1.24]
(0.96,2.43]
I tried to use strplit because of my regexpitis
On Wed, Feb 25, 2009 at 12:51 PM, Uwe Ligges
lig...@statistik.tu-dortmund.de wrote:
R User R User wrote:
Hi guys,
I'm evaluating R for basic data exploration. I produce a bar plot of the
data, with the x axis labels aligned vertically. However, the start of
labels longer than about 10
b12=log(.4/.6)
b13=log(.85/.15)
Can you point out the syntax errors in the first one?
I have been using a couple of books for the code but I am no expert.
The choice of the runif(1,0,2) density and choice of 2 as the constant in
the envelope function was made to envelope the whole variable;
So, it seems you are using Windows. One possibilty is that your
Windows security settings are stopping you displaying CHM help files:
if you have control over these you could try to resolve that.
Another alternative is to change the default type of help used: the
settings for that are in the
Hi Johannes.
My guess is that you will get no Std. Error, if you have only 1 observation
for that year.
Could that be the case ?
On Mon, Mar 2, 2009 at 8:41 PM, Johannes Habel johannes.ha...@gmx.dewrote:
Hello everybody,
I'm running a panel regression with dummy variables controlling for
101 - 169 of 169 matches
Mail list logo