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):
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),
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
: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Jim Lemon
Sent: 10 September 2007 12:07
To: [EMAIL PROTECTED]
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] Help in installing and loading the BradleyTerry add on
package in R
Kalyan Roy (DEL/MSG) wrote:
How do I install and load
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
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
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.
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
,
layout.pos.row = 2))
grid.rect()
p.legend - get_legends(p)
grid.draw(p.legend)
--END CODE
-Original Message-
From: hadley wickham [mailto:[EMAIL PROTECTED]
Sent: Monday, September 10, 2007 7:58 AM
To: Te, Kaom
Cc: r-help@stat.math.ethz.ch
Subject: Re
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
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
, Kaom
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] ggplot legend consolidation
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
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
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)
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
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({
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) {
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
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)
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
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 from R .To be specific,
m--matrix(ncol=3,nrow=10)
error message:[ reached
Your problem is not related to missing file permissions but to a problem
with the linker (ld). As you can see from the error message, ld cannot
find a compatible libR.so so compilation fails. You did not specify
which architecture you're on but in case it is a 64-bit platform, the
following thread
On 09/09/2007 7:01 AM, Yuchen Luo wrote:
Dear friends.
I use ConstrOptim( ) and got error message initial value not available.
My understanding of initial value not available is that one of the
following 3 cases happens:
1.The objective function is not well defined at the point of the
On 09/09/2007 7:24 AM, Werner Wernersen wrote:
Hello,
using Sweave, is there any option to preserve the
original aspect ratio of plots generated from R code?
Consider this Sweave chunk:
test,echo=F,fig=T,width=2,height=2=
x - 1:10
y - sin(x)
par(mar=c(4,4,0,4))
plot(x,y,
xlab=x
Have a look at rbind.fill in the reshape package.
Hadley
On 9/9/07, Muenchen, Robert A (Bob) [EMAIL PROTECTED] wrote:
Hi All,
If I need to stack two data frames, I can use rbind, but it requires
that all variables exist in both sets. I can make that happen, but other
stat packages would
Perfect. Thanks Hadley!
-Original Message-
From: hadley wickham [mailto:[EMAIL PROTECTED]
Sent: Sunday, September 09, 2007 10:11 AM
To: Muenchen, Robert A (Bob)
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R] stacking data frames with different variables
Have a look at rbind.fill
Wolfgang,
This looks great and certainly puts to shame my code .
Thanks,
Monica
Date: Fri, 7 Sep 2007 21:04:32 +0100 From: [EMAIL PROTECTED] To: [EMAIL
PROTECTED] CC: r-help@stat.math.ethz.ch Subject: Re: [R] confusion matrix -
better code? Dear Monica, try this: cm = table(tr
On 09/09/2007 1:36 PM, kevinchang wrote:
This is bizzare. The for loop I set seems to go infinitly because the
hourglass cursor never disappers and I never get the result. But when I
click on X of the window try to exit R, the execution finishes right after
I abort exiting when RGui asks me
I assume that you want to do the fitdistr on one of the columns of the
dataframe that you have read in. What does 'str(ONES3)' show? If the
data is in the first column, try:
fitdistr(ONES3[[1]],chi-squared)
On 9/9/07, Terence Broderick [EMAIL PROTECTED] wrote:
I am trying to fit the
Dear Professor Mordoch.
Thank you very much for your help! Your time is highly appreciated!
I do intend to optimize over 3 parameters and the way I did it is
constrOptim(c(0.5,0.3,0.5), fit.error, fit.error.grr, ui=ui,ci=ci).
Also, There is a missing value for both sigmae and ss. (The last 3rd
Many thanks for the quick reply, Duncan.
Now I see that something is going wrong in R already.
I overlooked this because the R window was so big. I
have to persuade plotViewport() not to create a square
area but one which takes the different margins into
account.
Best regards,
Werner
---
Thanks all!
On 09/09/2007, Nick Chorley [EMAIL PROTECTED] wrote:
Hi,
This may be a stupid question, but why when creating an AR model, why does
R set the first n components of the ar object's resid vector to NA (where n
is the order of the model)?
Regards,
NC
[[alternative
Hi Gabor,
Wow, this is awesome although I eventually should learn MySQL for
integrating it on web-based DB management using PHP or Perl, this is a
very helpful tool for me to start with!
Thank you very much
Gabor Grothendieck wrote:
Others have already pointed out %in% but regarding
Great. Regarding the web, note that there are actually quite a few R
web projects as well:
http://www.lmbe.seu.edu.cn/CRAN/doc/FAQ/R-FAQ.html#R-Web-Interfaces
I have used rpad (www.rpad.org) which has an integrated web server right
in the R package making setup a non-issue.
On 9/8/07,
François Pinard wrote:
[Roland Rau]
[François Pinard]
I wonder what happened, for R to hide the underlying Scheme so fully,
at least at the level of the surface language (despite there are
hints).
To further foster portability, we chose to write R in ANSI
# Create a matrix of ball locations
# You'd do this using the calls within your points function
balls - matrix(c(0,50,25,-150,-100,-50), ncol=2, byrow=F)
# Draw a line from the origin to each ball location
apply(balls, 1, function(x) lines(c(125, x[1]), c(-210, x[2]), col='red'))
A more
a - 1:3
b - 11:13
c - 21:23
names - c('a','b','c')
do.call(data.frame, list(sapply(names, function(x) get(x
runner wrote:
What I am trying to do is as follows:
- I have listed names of all wanted objects (datasets A,B,C... ) in
current workspace as a vector:
obj -
Frank E Harrell Jr wrote:
Martin Becker wrote:
Dear Cory,
I am not familiar with SAS, but is this what you are looking for?
divisionTable - matrix(c(1, New England,
2, Middle Atlantic,
3, East North Central,
4,
KateM == Katharine Mullen [EMAIL PROTECTED]
on Fri, 7 Sep 2007 20:07:41 +0200 (CEST) writes:
KateM The thread you linked to regarding Levenberg-Marquardt's supposed
lack of
KateM availability is from 2001; it has been possible to get
KateM to the MINPACK implementation of
Paul, what is the question? If the question is why you get this
warning message, the reason is that the 'lib' argument is missing
and install.packages is using '/usr/lib/R/library'.
If you want to get rid of the warning supply the 'lib' argument.
Gabor
On Sat, Sep 08, 2007 at 11:26:44AM +0100,
On 7/9/07 11:42 PM, eugen pircalabelu wrote:
I have a sample from a survey where household were interviewed. The
sample has 4 criteria on which the stratification was based: REGION,
SIZE OF HOUSEHOLD, SIZE OF LOCALITY, AGE OF HEAD OF HOUSEHOLD. Since i
don't have the whole information in
I'm wondering about experiences:
Do you know of cases where minpack.lm's nls.lm() solved a
(real) problem that nls() would have a problem with ?
In short, no. However, I looked at this question in the limited context
of fitting the parameters of a linear superposition of 2 exponentials
Luis Naver wrote:
I have a list of observations that are -1, 1 or 0. I would like to
represent them in a horizontal bar color coded based on value like a
stacked bar graph. I can achieve this in the form of a png with the
following code:
A = floor(runif(10)*3) - 1
png(width=100,
On 9/8/07, Gabor Csardi [EMAIL PROTECTED] wrote:
Paul, what is the question? If the question is why you get this
warning message, the reason is that the 'lib' argument is missing
and install.packages is using '/usr/lib/R/library'.
If you want to get rid of the warning supply the 'lib'
Hi Shubha,
genoud does not return the initial fit value. But you could easily
obtain it by passing your starting values to your function directly.
Alternatively, one can have genoud print out the entire initial
population (or the entire population as is evolves), and one can then
decide to
Dear Ken.
You could also try and use RSPerl (http://www.omegahat.org/RSPerl/). It
allows one to use R commands in Perl and vice-versa.
regards,
Paul
Pierce, Ken schreef:
Is there a way to run a simple perl script from R?
Kenneth B. Pierce Jr.
Research Ecologist
Landscape Ecology,
On 9/8/07, Peter Dalgaard [EMAIL PROTECTED] wrote:
François Pinard wrote:
[Roland Rau]
[François Pinard]
I wonder what happened, for R to hide the underlying Scheme so fully,
at least at the level of the surface language (despite there are
hints).
To further foster
will need them to display/capture what is happening.
Dirk
| -Original Message-
| From: Dirk Eddelbuettel [mailto:[EMAIL PROTECTED]
| Sent: Friday, September 07, 2007 1:11 PM
| To: Pierce, Ken
| Cc: r-help
| Subject: Re: [R] Running a PERL script from R
|
| On Fri, Sep 07, 2007 at 12:15:51PM
Michael,
Thank you very much. My code is certainly put to shame by yours. I promise to
read about factor to see how you use it and why ;-))
I really appreciate your help.
Monica
Subject: RE: [R] confusion matrix - better code? Date: Fri, 7 Sep 2007
15:36:00 -0500 From: [EMAIL PROTECTED
[Peter Dalgaard]
[François Pinard]
I meant that R might have implemented a Scheme engine [...] with
a surface language [...] which is purposely not Scheme, but could have
been. [...] one could dare dreaming that the Scheme engine in R be
completed, and Scheme offered as an alternate extension
Hi Jim,
Thanks for going to the trouble of writing the function. I'll copy and try it.
Gerard
At 06:30 PM 9/7/2007, jim holtman wrote:
This function should do it for you:
file1 - read.table(textConnection( id rx week dv1
+ 1 1 11 1
+ 2 1 12 1
+ 3 1 13 2
+ 4
On Sat, 2007-09-08 at 11:41 -0400, John Kane wrote:
Is there any convenient way to supress the x that
appears in csv export files? I would like to be able
to export a file and add a comment to it yet still be
able to read it back into R. I don't see any way to
get rid of the x that
uv wrote:
Hi. I need to use a few different clustering functions. I managed to run the
kmeans() one which is in my stats library, but I can't use any function,
such as agnes(), that is in my cluster library. Any idea how to access
other libraries?
Thanks!
Both stats and cluster are
Hi,
I am trying to install RSQLite package on my Fedora workstation. I
tried to install other packages as well, but each time I got the same
error messages saying compilation error and non zero exit status.
Do I have to specify lib=? I never specified the library path before
when I was using
Hi, try install packages whit 'sudo'.
$sudo R
--
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40 S 49° 16' 22 O
On 08/09/2007, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
Hi,
I am trying to install RSQLite package on my Fedora workstation. I
tried to install other packages as well,
Hi,
Still got the same error message. I did su R when I got the error
message for the first time. I have never seen this error message. I
will be googling for solutions as well...
Thank you.
Quoting Henrique Dallazuanna [EMAIL PROTECTED]:
Hi, try install packages whit 'sudo'.
$sudo R
Thanks Marc. It works. I had not thought of using
col.names = FALSE as I wanted to keep the colnames. I
see that I will just have to do another write.table
command to do this. Humm, actually it took a a bit of
juggling to do the names but it's looking fine now.
--- Marc Schwartz [EMAIL
Hi,
On Sep 7, 2007, at 11:00 PM, Yuchen Luo wrote:
constrOptim(c(0.5,0.3,0.5), fit.error, fit.error.grr,
ui=-1*ui,ci=-1*ci)
and I am confronted with error message initial value not feasible
I plug in the initial value of (0.5,0.3,0.5) to function fit.error and
fit.error.grr and have
On Wed, Sep 05, 2007 at 04:43:19PM -0700, Yuchen Luo wrote:
Dear friends.
I use nls() and encounter the following puzzling problem:
I have a function f(a,b,c,x), I have a data vector of x and a vectory y of
realized value of f.
Case1
I tried to estimate c with (a=0.3,
Oliver,
I am attaching an HTML document in which I have plotted -2Log(x) vs. x. If you
examine the plot you will see that -2Log(x) can be negative. Since -2Log(x) is
part of AIC and BIC, AIC and BIC can be negative.
John
John Sorkin M.D., Ph.D.
Chief, Biostatistics and Informatics
Baltimore VA
Dear Niccolo,
a) your example program
works for me
and on a Pentium 4, 3 Ghz, 1GB ram, The model fit uses
user system elapsed
146.920.76 143.32
library(geepack)
set.seed(299)
header - rep.int(seq(1:615),sample(seq(1:19),size=615,replace=T))
inr -
Others have already pointed out %in% but regarding your comment about
SQL, you can use SQL to manipulate R data frames using the sqldf package
which provides an interface to lower level RSQLite (and RMySQL in the future)
routines. The following examples use SQLite underneath:
DF -
Hi Greg,
Actually, I'm using the scatterplot function in the car namespace. My
understanding is that the scatterplot is turned into a boxplot when the x
variable is a factor, but I can't see that in the code of the function.
I've successfully modified the function to output regression stats
Your columns are factors, not character strings. Use as.is = TRUE as
an argument to read.table. Also its a bit dangerous to use T although
not wrong. Its safer to use TRUE.
On 9/7/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
Dear List,
I have a newbie question. I have read in a
A matrix is for situations where every element is of the same class
but your columns have different classes so use a data frame:
DF - data.frame(a = 11:15, b = letters[1:5], stringsAsFactors = FALSE)
subset(DF, a %in% 11:13)
subset(DF, a %in% c(0, 11:13)) # same
Suggest you review the
From: Dong-hyun Oh
Dear UseRs,
I would like to know which function is the most efficient in finding
convex hull of points in 3(or 2)-dimensional case?
Functions for finding convex hull is the following:
convex.hull (tripack), chull (grDevices), in.chull (sgeostat),
convhulln
Try this:
library(zoo)
plot(na.approx(zoo(as.matrix(data[-1]), data[,1])), plot.type = single)
See ?na.approx, ?plot.zoo, ?xyplot.zoo and vignette(zoo)
On 9/7/07, Markus Schmidberger [EMAIL PROTECTED] wrote:
Hello,
I have this kind of dataframe and have to plot it.
data - data.frame(sw=
Thanks a lot Gabor, that was very helpful. All sorted now!
Best
David
Your columns are factors, not character strings. Use as.is = TRUE as
an argument to read.table. Also its a bit dangerous to use T
although
not wrong. Its safer to use TRUE.
On 9/7/07, [EMAIL PROTECTED] [EMAIL
: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Yves Moisan
Sent: Friday, September 07, 2007 5:05 AM
To: r-help@stat.math.ethz.ch
Subject: Re: [R] text in boxplots
Hi Greg,
Actually, I'm using the scatterplot function in the car
namespace. My understanding is that the scatterplot
Markus Schmidberger wrote:
Hello,
I have this kind of dataframe and have to plot it.
data - data.frame(sw= c(1,2,3,4,5,6,7,8,9,10,11,12,15),
zehn =
c(33.44,20.67,18.20,18.19,17.89,19.65,20.05,19.87,20.55,22.53,NA,NA,NA),
zwanzig =
Phil,
If you have only a list of addresses and nothing else - i have to recognize i
am lost. But if you have more info you can use a GIS software and it will
calculate your distance automatically for all your addresses. To have your
distance in miles or Km you need to have a projection in
6. set the reset.par argument to FALSE in calling scatterplot as it
states in the documentation for scatterplot. (still look at the help on
'xpd' argment to par if you want to add text to the margins).
Hi Greg,
That seems to have done it partly. I set it to false once and then I could
see
Dear Cory,
I am not familiar with SAS, but is this what you are looking for?
divisionTable - matrix(c(1, New England,
2, Middle Atlantic,
3, East North Central,
4, West North Central,
5, South
Moisan; r-help@stat.math.ethz.ch
Subject: Re: [R] text in boxplots
The posting guide (a link at the bottom of every e-mail from
the list) suggests including a small example of what you are
trying to do.
Without an example, we have to guess at what you are trying to do.
I am still guessing
Yes but delete does not return anything so its not useful. In the devel
version of sqldf you can pass multiple command so try this using the
builtin data frame BOD noting that the record with demand = 8.3 was
removed:
library(sqldf)
Loading required package: RSQLite
Loading required package:
1 - 100 of 57901 matches
Mail list logo