Hello,
I have difficulties combining boxes and lines in plot legend. I
searched previous R-posts and found this (with no solution):
http://tolstoy.newcastle.edu.au/R/help/06/07/30248.html. Is there a
way to avoid boxes behind the line legends?
x1 - rnorm(100)
x2 - rnorm(100, 2)
hist(x1, main = ,
Check out:
http://tolstoy.newcastle.edu.au/R/e2/help/07/05/16777.html
On 9/10/07, Lauri Nikkinen [EMAIL PROTECTED] wrote:
Hello,
I have difficulties combining boxes and lines in plot legend. I
searched previous R-posts and found this (with no solution):
Hello,
I would like to know what can I do if I use strplit with a string and I want
to use the middle left,I mean I have this:
strsplit(bA531F16-rep,\\-)
[[1]]
[1] bA531F16 rep
I would like to work just with bA531F16 in another variable, what could I
do?, Thank you
Why not use the more simple
xyplot(total.fat~x|variable,groups=Group,
data=tmp1,type=c(p,r))
???
See ?panel.xyplot and especially the type argument of that panel function.
Best regards
Frede Aakmann Tøgersen
Scientist
UNIVERSITY OF AARHUS
Faculty of Agricultural Sciences
Dept. of
Thanks Gabor, I got it!
For example:
x1 - rnorm(100)
x2 - rnorm(100, 2)
hist(x1, main = , col = orange,ylab = density, xlab = x, freq
= F, density = 55, xlim = c(-2, 5), ylim = c(0, 0.5))
par(new = T)
hist(x2, main = , col = green, ylab = , xlab = ,axes = F, xlim
= c(-2, 5), ylim = c(0, 0.5),
Dear Friends.
I found something very puzzling with constOptim(). When I change the
parameters for ConstrOptim, the error messages do not seem to be
consistent with each other:
constrOptim(c(0.5,0.3,0.5), f=fit.error, gr=fit.error.grr, ui=ui,ci=ci)
Error in constrOptim(c(0.5, 0.3, 0.5), f =
Dear Rosa,
please be more specific. Statistical tests for which hypothesis?
For example, some tests can be made robust using Heteroskedasticity-
*and Autocorrelation-* Consistent (HAC) covariance matrices in package
'sandwich': see
- waldtest{lmtest} for a redundant variables test much like
unlist(strsplit(bA531F16-rep,\\-))[1]
[1] bA531F16
--- Carlos Morales [EMAIL PROTECTED]
wrote:
Hello,
I would like to know what can I do if I use
strplit with a string and I want to use the middle
left,I mean I have this:
strsplit(bA531F16-rep,\\-)
[[1]]
[1] bA531F16 rep
Yogesh Tiwari wrote:
Hi,
How to plot a variance over a data point, something like error bar.
Hi Yogesh,
The generic method is to use the arrows function, and there are quite
a few variations on this theme (in alpha order):
brkdn.plot(plotrix)
dispbars(plotrix)
errbar(Hmisc and sfsmisc)
Thanks Frede
I didn't know about the r type.
Ross Darnell
-Original Message-
From: Frede Aakmann Tøgersen [mailto:[EMAIL PROTECTED]
Sent: Mon 10-Sep-07 4:45 PM
To: Ross Darnell; r-help@stat.math.ethz.ch
Subject: SV: [R] lattice panel.lmline problem
Why not use the more simple
prcomp() in stats handles matrices with n p well, IMO.
--
Bjørn-Helge Mevik
__
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
and provide
PS == Paul Smith [EMAIL PROTECTED]
on Sun, 9 Sep 2007 12:17:32 +0100 writes:
PS On 9/9/07, kevinchang [EMAIL PROTECTED] wrote:
I tried to made the matrix with this size by either matrix() or array().
However, there seems to be default limit of number for rows made. I got
sort
Hi,
I have been suggested to use Information theoretic approach instead of a
Geneal linear model for a multiple regression. Is is possible to perform
this type of anlaysis in R?
Thanks,
Adela
--
Adela González Megías
Depto. Biología Animal
Fac. Ciencias
Universidad de Granada
18071
How do I install and load the BradleyTerry add on package in R 2.5.1 in
MSWindowsXP environment?
Kalyan Roy
Indian Market Research Bureau (IMRB) International
New Delhi, India
[[alternative HTML version deleted]]
__
However you may also need to install the package brlr, since the
BradleyTerry package depends on this.
For Windows users, it's usually easiest to install packages using the
Packages menu in the RGui - any dependencies are then automatically
installed.
Heather
-Original Message-
From:
On 9/10/07, Martin Maechler [EMAIL PROTECTED] wrote:
PS On 9/9/07, kevinchang [EMAIL PROTECTED] wrote:
I tried to made the matrix with this size by either matrix() or
array().
However, there seems to be default limit of number for rows made. I
got sort
of error message
Kalyan Roy (DEL/MSG) wrote:
How do I install and load the BradleyTerry add on package in R 2.5.1 in
MSWindowsXP environment?
Hi Kalyan,
If
R CMD INSTALL
doesn't work, you can use WinZip or Zip Reader to unzip the package to:
C:\Program Files\R-2.5.1\library
or whatever your path to the
Hello,
Some weeks ago, thanks to you, I managed to install R, to connect to a local
MySQL Database and to launch some queries with a script written with Tinn-R.
My script is now ok and would like to test it with the real database.
I did the same installation of R, DBI package and RMySQL package
Hi,
I'm looking for an implementation of the corrected resampled t-test
[1] to compare to different machine learning classifiers. Any idea,
whether this is implemented in a R package?
Thanks in advance,
Christian
[1] Nadeau Bengio, 2003, Inference for the Generalization Error,
Machine
This is not quite what you want but you can try this:
legend(3, 0.45, legend = c(x1, x2, mean(x1), mean(x2)), col =
c(orange, green),pch = c(15,15,-1,-1), lty=c(-1,-1,2,2))
Although pch=22 should draw a filled square with a border - but it draws only
the border instead
Monica
I'm trying to build an extension that uses the \u characters, but
building the library fails. The Writing R Extensions guide says
that these characters should not be used, but I really can't find any
way around them. Is there any way to configure the building of a
library so that it will
Hi,
sapply(formals(readBin), mode)
con what n sizesignedendian
namename numeric logical logicalcall
returns for the mode of size logical. But in the documentation is said
that size should be integer. Does anyone know why the mode is logical?
Thanks
The following is one of the examples in the help page for histogram:
histogram( ~ height | voice.part, data = singer,
xlab = Height (inches), type = density,
panel = function(x, ...) {
panel.histogram(x, ...)
On 9/10/2007 10:26 AM, Sigbert Klinke wrote:
Hi,
sapply(formals(readBin), mode)
con what n sizesignedendian
namename numeric logical logicalcall
returns for the mode of size logical. But in the documentation is said
that size should be
You could directly use
strsplit(bA531F16-rep,\\-)[[1]][1]
and
strsplit(bA531F16-rep,\\-)[[1]][2]
Regards
Carlos Morales-2 wrote:
Hello,
I would like to know what can I do if I use strplit with a string and I
want to use the middle left,I mean I have this:
I have recently been introduced to the ggplot package by Hadley Wickham
and must say I am quite impressed so far at how easy it is to make
attractive plots, but one thing I am struggling over is how to
consolidate legends.
It's not currently possible to consolidate them (although in the
Hi,
I was trying to install the packageRWinEdt in my computer with Vista OS.
Once I finished the installation, R just cannot load this library at all,
reporting some error information. I reinstalled R and Winedit and reloaded
this package. RWinEdt seemed to register in the computer but never work
Hi, there:
just a little bit off-topic:
any algorithm in classification is good for ordered features, like all
variables used are 1, 2, 3,...; not really continuous but ordered. i
tried random forest and dlda already.
thanks.
--
Weiwei Shi, Ph.D
Research Scientist
GeneGO, Inc.
Did you
Hi,
I was trying to install the packageRWinEdt in my computer with Vista OS.
O
Use Tinn-R it works with Vista although has some minor issues (which are
probably system specific).
https://sourceforge.net/projects/tinn-r
Stefan
-=-=-
... Time is an illusion, lunchtime doubly so. (Ford
Yuchen Luo wrote:
Dear Friends.
I found something very puzzling with constOptim(). When I change the
parameters for ConstrOptim, the error messages do not seem to be
consistent with each other:
constrOptim(c(0.5,0.3,0.5), f=fit.error, gr=fit.error.grr, ui=ui,ci=ci)
Error in
On Thu, 6 Sep 2007, eugen pircalabelu wrote:
Good afternoon!
I'm trying to use the Survey package for a stratified sample which has
4 criteria on which the stratification is based. I would like to get the
corrected weights and for every element i get a weight of 1
E.g: tipping
On Thu, 6 Sep 2007, Tobias Verbeke wrote:
eugen pircalabelu wrote:
I'm trying to use the survey package to get a better point of view
for my data, but i need some piece of advice:
i have some data from a survey which has been stratified using 2
criteria: region(7 values), size of
Steve,
This example works:
x-TukeyHSD(fm1, tension, ordered = TRUE)
as.table(x$tension)
Steve Powers wrote:
So I've come across a few cases where complex outputs from functions
will not write to tables. The most recent case involves the TukeyHSD
function in the stats package. If I save
The error message about the feasible region comes from constrOptim(),
before your function is called. The error message about missing lambda1
comes from calling your function.
-thomas
On Sun, 9 Sep 2007, Yuchen Luo wrote:
Dear Friends.
I found something very puzzling with
On Sun, 9 Sep 2007, eugen pircalabelu wrote:
A short example:
stratum id weight nh Nh y sex
1 1 3 5 15 23 1
1 2 3 5 15 25 1
1 3 3 5 15 27 2
1 4 3 5 15 21 2
1 5 3 5 15 22 1
2 6
This example works fine:
test-matrix(c(1,2,'VOICIUNPETITTES',3),ncol=2,nrow=2)
write.csv(test,file='C:/xavier/test.csv')
Could you provide the same small example when it doesn't work?
kwaj wrote:
Hello,
I have a peculiar problem which I am hoping I can get help
On Mon, 10 Sep 2007, Sigbert Klinke wrote:
Hi,
sapply(formals(readBin), mode)
con what n sizesignedendian
namename numeric logical logicalcall
returns for the mode of size logical. But in the documentation is said
that size should be integer. Does
Good morning, everyone,
I am sorry for this off-topic post but think I can get great answer
from this list.
My question is what is the best OS on PC (laptop) for statistical
computing and why.
I really appreciate your insight.
Have a nice day.
__
Hi Hadley,
I just tried out your suggestion, but it does not look like the
get_legends function is working correctly. Instead of returning a grob
back to me it returns NULL.
Here is my modified code and the results of running it.
Any help would be appreciated. I believe that once I can get the
Dear R users,
Below I have written 4 functions CIT1, CIT2a and CIT2b and CIT3 which
recode a variable CLD_ISCH into 3 new variables(T1 T2 T3), I wish to use
T1, T2 and T3 based on the values of tf1 and tf2.
(NOTE:- T2a is used to create T2 in a long winded manner due to my lack
of programming
Hello,
I am new to R exploratory data analysis and plotting. Is anyone aware of a way
to overlay a set of conditional histograms with conditional PDFs? Below, I
generate a lattice plot of precipitation histograms based on different months
and stations, given a subset of the dataset:
You want whatever all the people you are working with are using
to make it as easy as possible to work together with them.
On 9/10/07, Wensui Liu [EMAIL PROTECTED] wrote:
Good morning, everyone,
I am sorry for this off-topic post but think I can get great answer
from this list.
My question is
Are there any packages in R that reproduce the package resample of S-Plus?
The sample() function in R doesn't provide equivalent flexibility of
bootstrap() and bootstrap2().
Robert A. LaBudde, PhD, PAS, Dpl. ACAFS e-mail: [EMAIL
I have 316 files storing a number of cycles of the respiratory signal
from 316 different patients. Each file record is made up of the
following data:
amplitude,phase,timestamp,validflag,ttlin,mark,ttlout
amplitude: is the position in cm relative to an arbitrary reference
(signal amplitude)
useRs,
I am looking to find the minimum positive value of some data I have.
Currently, I am able to find the minimum of data after I apply some other
functions to it:
x
[1] 1 0 1 2 3 3 4 5 5 5 6 7 8 8 9 9 10 10
sort(x)
[1] 0 1 1 2 3 3 4 5 5 5 6 7 8 8 9 9 10
Try this:
min(diff(sort(x))[diff(sort(x))0])
--
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40 S 49° 16' 22 O
On 10/09/2007, dxc13 [EMAIL PROTECTED] wrote:
useRs,
I am looking to find the minimum positive value of some data I have.
Currently, I am able to find the minimum of
On Mon, 2007-09-10 at 11:20 -0700, dxc13 wrote:
useRs,
I am looking to find the minimum positive value of some data I have.
Currently, I am able to find the minimum of data after I apply some other
functions to it:
x
[1] 1 0 1 2 3 3 4 5 5 5 6 7 8 8 9 9 10 10
Hello friends,
I loaded R 2.4.1 onto a Fedora Core 6 Linux box (taking all defaults). Then
I ran these commands from within R:
options(CRAN=http://cran.stat.ucla.edu;)
install.packages(CRAN.packages()[,1])
As a new user of R, I was shocked when I finished loading R and discovered
the following
Sorry, I should have mentioned that get_legend won't work on the plots
that you are actually plotting - you have turned their legends off!
You'll need a plot which isn't plotted, but is used to produce the
legends.
Hadley
On 9/10/07, Te, Kaom [EMAIL PROTECTED] wrote:
Hi Hadley,
I just tried
http://www.burns-stat.com/pages/Tutor/bootstrap_resampling.html
includes a synopsis of R packages that do bootstrapping.
It is brief and incomplete, but hopefully useful.
Patrick Burns
[EMAIL PROTECTED]
+44 (0)20 8525 0696
http://www.burns-stat.com
(home of S Poetry and A Guide for the Unwilling
I am having some trouble getting the persp() package to change the x
and y axis on a 3d plot. It defaults to the [0,1] interval and when I
try to change it I get errors.
Example:
This works:
D - c(1,2,3,4,5,6,7,8,9,10)
M - c(11,12,13,14,15,16,17,18,19,20)
DM - cbind(D,M)
persp(DM,
A Lenzo wrote:
Hello friends,
I loaded R 2.4.1 onto a Fedora Core 6 Linux box (taking all defaults). Then
I ran these commands from within R:
options(CRAN=http://cran.stat.ucla.edu;)
install.packages(CRAN.packages()[,1])
As a new user of R, I was shocked when I finished loading R and
On 9/10/07, Ross Darnell [EMAIL PROTECTED] wrote:
Thanks Frede
I didn't know about the r type.
For the record, this is probably what you wanted:
xyplot(...
panel = panel.superpose,
panel.groups = function(x, y, ...) {
panel.xyplot(x, y, ...)
if (length(x)
I am just trying to teach myself how to use the mle function in R because it is
much better than what is provided in MATLAB. I am following tutorial material
from the internet, however, it gives the following errors, does anybody know
what is happening to cause such errors, or does anybody know
On 11/09/2007, at 4:22 AM, Wensui Liu wrote:
Good morning, everyone,
I am sorry for this off-topic post but think I can get great answer
from this list.
My question is what is the best OS on PC (laptop) for statistical
computing and why.
I really appreciate your insight.
Have a nice day.
Linux! Mac OS is ok to me2.
On 9/10/07, Wensui Liu [EMAIL PROTECTED] wrote:
Good morning, everyone,
I am sorry for this off-topic post but think I can get great answer
from this list.
My question is what is the best OS on PC (laptop) for statistical
computing and why.
I really appreciate
Dear Professor Murdoch.
Thank you for your help!
1. I believe c(0.5,0.3,0.5) satisfies the constrain because I did the
following experiment
ui=-1*ui
ci=-1*ci
constrOptim(c(0.5,0.3,0.5), f=fit.error, gr=fit.error.grr, ui=ui,ci=ci)
The same error message pops up. Any theta ( in this case,
Terence Broderick wrote:
I am just trying to teach myself how to use the mle function in R because it
is much better than what is provided in MATLAB. I am following tutorial
material from the internet, however, it gives the following errors, does
anybody know what is happening to cause such
Hi
I was wondering if someone could help me with an rpart problem. I can see
that cross-validation is the default for tree selection in rpart -- has a
bootstrap method been implemented anywhere? I think this is a different
thing to 'bagging' or 'boosting' -- I still want 'one' tree at the end,
From ?trackObjs:
Overview of trackObjs package
Description:
The trackObjs package sets up a link between R objects in memory
and files on disk so that objects are automatically resaved to
files when they are changed. R objects in files are read in on
demand and do not
Hello All,
I have searched many help forums, message boards, etc. and I just can't
apply the comments to what I need my program to do. I am running R 2.5.1 on
an XP system, and my desire is to produce replicate datasets for a
simulation study I am running. Essentially, I have sets of
See the article Cluster in R Task Views
http://cran.at.r-project.org/src/contrib/Views/Cluster.html
It lists names of packages for clustering analysis, which you can install.
Or, go to choose another CRAN mirror closest to you from the R web site,
then click on 'Task views' in the left frame.
My sense is that R users are even split between UNIX and Windows
users so either will do in terms of the larger community.
Some R packages may not be avaliable on every platform or will
be available on one platform before another or there will be
certain platform-specific issues. So in the end
Hi Jonathan,
What exactly do you mean by replication?
Do you want to keep a1,b1,c1,... unchanged but have 30
different sets of random numbers?
Regards,
Moshe.
--- VTLT1999 [EMAIL PROTECTED] wrote:
Hello All,
I have searched many help forums, message boards,
etc. and I just can't
Here are some solutions each of which
1. has only one line,
2. x only occurs once so you can just plug in a complex expression
3. no temporary variables are left
min(sapply(x, function(z) if (z 0) z else Inf))
(function(z) min(ifelse(z 0, z, Inf))) (x)
with(list(z = x), min(z[z 0]))
local({
x-seq(-1,1,length=10)
y-seq(-1,1,length=10)
a-matrix(c(1,2,2,1),2,2)
b-matrix(c(2,1,1,2),2,2)
fv-function(x,y) {
m-x*a+y*b
t-m[1,1]+m[2,2]; d-m[1,1]*m[2,2]-m[1,2]^2
return((t-sqrt(t^2-4*d))/2)
}
gv-function(x,y) {
t-x*(a[1,1]+a[2,2])+y*(b[1,1]+b[2,2])
Its a FAQ:
http://hermes.sdu.dk/Rdoc/faq.html#Why%20does%20outer()%20behave%20strangely%20with%20my%20function%3f
On 9/10/07, Jan de Leeuw [EMAIL PROTECTED] wrote:
x-seq(-1,1,length=10)
y-seq(-1,1,length=10)
a-matrix(c(1,2,2,1),2,2)
b-matrix(c(2,1,1,2),2,2)
fv-function(x,y) {
Dear Listers,
I am a little tired of installing all packages I want every time when
I instill a new version of R.
Say, if I have a list of packages I need to use, is it possible to
tell R to install them all for me automatically rather than I install
them one by one?
Thx.
Try
?update.packages
Ross Darnell
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Wensui Liu
Sent: Tuesday, 11 September 2007 11:45 AM
To: r-help@stat.math.ethz.ch
Subject: [R] install packages automatically
Dear Listers,
I am a little tired of
Version 1.1-8 of package 'scuba' has been uploaded to CRAN.
'scuba' is a package for scuba diving calculations and decompression
models. It supports dive profiles (tables, plotting etc), analysis of dive
profiles using decompression models, gas toxicity calculations, and gas
usage calculations.
Mange tak!
FYI, this is the way it is able to run (I was going to attach
station.precip.R, but I read that attaching files is not recommended
- let me
know if you would like it)
x - dget(file=C://Documents and Settings/Bradley/My
Documents/Arizona/CourseResources/ATMO529/station.precip.R)
I am using 'xyplot' in lattice to plot some data where the x-axis is a
POSIXct date. I have data which spans a 6 month period, but when I
plot it, only the last month is printed on the right hand side of the
axis. I would have expected that at least I would have a beginning
and an ending point
Either
min(diff(sort(x))[diff(sort(x))0])
or
min(diff(sort(unique(x
--- dxc13 [EMAIL PROTECTED] wrote:
useRs,
I am looking to find the minimum positive value of
some data I have.
Currently, I am able to find the minimum of data
after I apply some other
functions to it:
x
Fiona Callaghan wrote:
I was wondering if someone could help me with an rpart problem. I can see
that cross-validation is the default for tree selection in rpart -- has a
bootstrap method been implemented anywhere? I think this is a different
thing to 'bagging' or 'boosting' -- I still want
74 matches
Mail list logo