Unfortunately Duncan's suggestion to restart is probably the only way to
go here. I've done similar thing myself too. What I've learned was
that it is clever to include a so called hot-patch mechanism in your
code, which will load R source code found in a certain directory, say,
hot/, once
Hi
adiag() from the magic package does what you want:
library(magic)
Loading required package: abind
a - matrix(1,2,2)
b - matrix(6,2,3)
adiag(a,b)
[,1] [,2] [,3] [,4] [,5]
[1,]11000
[2,]11000
[3,]00666
[4,]006
Hello,
when I work in R, I write code in a text file that I run with the
source(filename) command. In R2.1.1 the file is read-only while the source
command is executed. This was not the case in R2.0.1. Is this a bug-fix or is
it possible not to have the file read-only when executed?
Best
Hi,
I am using robust regression, i.e. model.robust-ltsreg(MXD~ORR,data=DATA).
My question:- is there any way to determine the Robust Multiple R-Squared
(as returned in the summary output in splus)? I found an equivalent model in
the rrcov package which included R-square, residuals etc in it's
Thanks for everyone who replied to my question.
I tried newmat myself, seems to be working well. What I am interested in
is something like
* Fast element-wise operations. You know, it may be slow in R.
* (Some) control over memory allocation. I would like to specify when
the matrix should
Aeacides
[[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! http://www.R-project.org/posting-guide.html
Le 01.09.2005 10:32, Ott Toomet a écrit :
Thanks for everyone who replied to my question.
I tried newmat myself, seems to be working well. What I am interested in
is something like
* Fast element-wise operations. You know, it may be slow in R.
* (Some) control over memory allocation. I
Is there any function to plot the standard deviation with the barplots
_ _
| Deviation
|---
| -|
| |
| | Barplot
| |
| |
| |
with regards
Knut Krueger
__
On 9/1/05 6:32 AM, Knut Krueger [EMAIL PROTECTED] wrote:
Is there any function to plot the standard deviation with the barplots
_ _
| Deviation
|---
| -|
| |
| | Barplot
| |
| |
| |
I think barplot2 in
Hi Trevor,
please excuse my late reply; your e-mail was sent when I started
travelling (to Seattle and the DSC and Bioconductor workshops there).
Also, I hope you don't mind if I follow this up on R-help, since
there, the thread started and this does related to it.
Trevor == Trevor Hastie
Thank you all for the answers. I will look into both methods for doing so.
Jean
On Thu, 1 Sep 2005, Henrik Bengtsson wrote:
Unfortunately Duncan's suggestion to restart is probably the only way to
go here. I've done similar thing myself too. What I've learned was
that it is clever to
Romain == Romain Francois [EMAIL PROTECTED]
on Thu, 01 Sep 2005 11:19:59 +0200 writes:
Romain Le 01.09.2005 10:32, Ott Toomet a écrit :
Thanks for everyone who replied to my question.
I tried newmat myself, seems to be working well. What I
am interested in is
Hello,
Iâm currently working with Oracle under windows and Iâd like to use the
ROracle package which exists only for linux/unix. Is there any possibility to
port this package under windows?
L. Tessier
[[alternative HTML version deleted]]
Hello,
I am trying to draw a convex-polygon (envelope) connecting the outer points
(of a particular class) in a scatter plot
[which I want to do to illustrate the range of a particular type of samples
in an ordination plot].
Is there any function to plot this kind of polygon?
Kind regards
Hi,
ever thought about using the built-in help of R?
try
example(chull)
cu
detlef
On Thu, 1 Sep 2005 13:51:42 +0200
Zlatko Petrin [EMAIL PROTECTED] wrote:
Hello,
I am trying to draw a convex-polygon (envelope) connecting the outer points
(of a particular class) in a scatter plot
Hi, if I do
Z-rnorm(50)
Followed by
Dim(Z)
I get NULL. Is this correct? Shouldn't I get 50 instead?
TIA
Chris
DISCLAIMER:\ L'utilizzo non autorizzato del presente messagg...{{dropped}}
__
R-help@stat.math.ethz.ch mailing list
yes it is correct! look at ?dim() for more info. In this case you need
length(), i.e.,
Z - rnorm(50)
length(Z)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35,
Check out the dim vs length for vectors thread:
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/50720.html
This thread goes through the bug-or-feature discussion which is always
entertaining from a socio-R perspective.
Also, note Dim with a capital D doesn't exist.
HTH, Andy
-Original
Laurent TESSIER wrote:
Hello,
I’m currently working with Oracle under windows and I’d like to use the
ROracle package which exists only for linux/unix. Is there any possibility to
port this package under windows?
It is ported, but perhaps you need some additional files:
See
Christian Prinoth wrote:
Hi, if I do
Z-rnorm(50)
Followed by
Dim(Z)
I get NULL. Is this correct? Shouldn't I get 50 instead?
No, because Z has no dim attribute. ?dim tells you:
For an array (and hence in particular, for a matrix) dim retrieves the
dim attribute of the object. It is
Hjellvik Vidar wrote:
Hello,
when I work in R, I write code in a text file that I run with the
source(filename) command. In R2.1.1 the file is read-only while the source
command is executed. This was not the case in R2.0.1. Is this a bug-fix or is
it possible not to have the file
Should the Value: section of the help page read:
For an array (and hence in particular, for a matrix) dim retrieves the
dim attribute of the object. It is NULL for a vector of mode integer.
NOTE: for not or in the second sentence.
Uwe Ligges wrote:
Christian Prinoth wrote:
Hi, if I do
[EMAIL PROTECTED] wrote:
Dear listers,
I would like to know whether a maximum size is set for data sets in R.
Well, in theory, there are some limits re. size of some integers, but
you are far away, as far as I understand.
In practice, you will somewhere hit the limit of available
Hi,
I have built R 2.1.1 on Solaris and OSF1, and used the operating
system's native compiler set. This worked fine. However, when I try
using install.packages() it insists on using g77 etc., which doesn't
exist on those systems.
How can I specify which compilers to use? What about other
Chuck Cleland wrote:
Should the Value: section of the help page read:
For an array (and hence in particular, for a matrix) dim retrieves the
dim attribute of the object. It is NULL for a vector of mode integer.
NOTE: for not or in the second sentence.
NO!!! The help page is correct!
The
If you are indeed using lme and not lmer then the needed function is
VarCorr(). However, 2 recommendations. First, this is a busy list and
better emails subject headers get better attention. Second, I would
recommend using lmer as it is much faster. However, VarCorr seems to be
incompatible with
Nam-Ky Nguyen wrote:
Dear Rexpert,
I would like to thank those who spend time answering my email on the
burning of a CD with all R binary files for Windows and Linux. I have
tried a couple of suggestions but have not been successful. I will pass
these suggestions to our system
On Thu, 01 Sep 2005 08:25:14 -0400 Chuck Cleland wrote:
Should the Value: section of the help page read:
For an array (and hence in particular, for a matrix) dim retrieves the
dim attribute of the object. It is NULL for a vector of mode integer.
This statement is true, but less general than
On Thu, 1 Sep 2005, Trond Hasle Amundsen wrote:
Hi,
I have built R 2.1.1 on Solaris and OSF1, and used the operating
system's native compiler set. This worked fine. However, when I try
using install.packages() it insists on using g77 etc., which doesn't
exist on those systems.
How can
Achim Zeileis wrote:
On Thu, 01 Sep 2005 08:25:14 -0400 Chuck Cleland wrote:
Should the Value: section of the help page read:
For an array (and hence in particular, for a matrix) dim retrieves the
dim attribute of the object. It is NULL for a vector of mode integer.
This statement is
The OS is Windows XP.
I use to work on the file while it's executed and save changes continually. It
doesn't seem to have any effect on the current execution. I just find it
annoying not to be able to do it. It's not crucial, but if there is some easy
way around it I would like to know
Roger Bivand [EMAIL PROTECTED] writes:
On Thu, 1 Sep 2005, Trond Hasle Amundsen wrote:
Hi,
I have built R 2.1.1 on Solaris and OSF1, and used the operating
system's native compiler set. This worked fine. However, when I try
using install.packages() it insists on using g77 etc., which
Sean Davis schrieb:
I think barplot2 in the gregmisc (gplots) bundle will do that.
Sean
Ok thank´s but I've got an error and do not found the solution:
ci.l
[1] 304.09677 202.49907 0.0 63.14547 0.0 0.0 0.0
[8] 0.0 0.0 0.0
ci.h
[1] 633.50323
On Thu, 01 Sep 2005 14:44:53 +0200 Uwe Ligges wrote:
Achim Zeileis wrote:
On Thu, 01 Sep 2005 08:25:14 -0400 Chuck Cleland wrote:
Should the Value: section of the help page read:
For an array (and hence in particular, for a matrix) dim retrieves
thedim attribute of the object. It
If I use the following command to plot points:
plot(c(1,2,2,3,3,3), type=p, pch=20, ylab=Y Label, xlab=X Label, xaxt=n)
there is a large amount of space between the label X Label and the
actual x-axis. If I change the xaxt=n to xaxt=s, the label X
Label don't move at all. Is there a way to get
Hi all,
Could anyone tell me if there is any package (or function) that generates
values from a multivariate skew normal distribution?
Thanks all,
Caio
__
[[alternative HTML version deleted]]
I forgot the error message ...
Error in barplot2.default(xrow, plot.ci = TRUE, ci.l = ci.l, ci.h = ci.h) :
confidence interval values are missing
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
On Thu, 1 Sep 2005, Trond Hasle Amundsen wrote:
Roger Bivand [EMAIL PROTECTED] writes:
On Thu, 1 Sep 2005, Trond Hasle Amundsen wrote:
Hi,
I have built R 2.1.1 on Solaris and OSF1, and used the operating
system's native compiler set. This worked fine. However, when I try
On Thu, 2005-09-01 at 15:00 +0200, Knut Krueger wrote:
Sean Davis schrieb:
I think barplot2 in the gregmisc (gplots) bundle will do that.
Sean
Ok thank´s but I've got an error and do not found the solution:
ci.l
[1] 304.09677 202.49907 0.0 63.14547 0.0
Trond Hasle Amundsen wrote:
Roger Bivand [EMAIL PROTECTED] writes:
On Thu, 1 Sep 2005, Trond Hasle Amundsen wrote:
Hi,
I have built R 2.1.1 on Solaris and OSF1, and used the operating
system's native compiler set. This worked fine. However, when I try
using install.packages() it insists
hy all,
When i plot under R it generates a 440x440px image, is it possible to modify
and increase this ?
thks all
guillaume.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Knut Krueger wrote:
Sean Davis schrieb:
I think barplot2 in the gregmisc (gplots) bundle will do that.
Sean
Ok thank´s but I've got an error and do not found the solution:
ci.l
[1] 304.09677 202.49907 0.0 63.14547 0.0 0.0 0.0
[8] 0.0
There is an error in your function call. The argument 'ci.h' is
incorrect, as it should be 'ci.u'. Thus, use:
... I think I will need glasses
Tank's a lot Knut
__
R-help@stat.math.ethz.ch mailing list
[EMAIL PROTECTED] wrote:
hy all,
When i plot under R it generates a 440x440px image, is it possible to modify
and increase this ?
What device are we talking about?
See the corresponding help file, maybe starting at ?Devices
Uwe Ligges
thks all
guillaume.
Jamieson Cobleigh wrote:
If I use the following command to plot points:
plot(c(1,2,2,3,3,3), type=p, pch=20, ylab=Y Label, xlab=X Label,
xaxt=n)
there is a large amount of space between the label X Label and the
actual x-axis. If I change the xaxt=n to xaxt=s, the label X
Label don't
Roger Bivand [EMAIL PROTECTED] writes:
On Thu, 1 Sep 2005, Trond Hasle Amundsen wrote:
Roger Bivand [EMAIL PROTECTED] writes:
On Thu, 1 Sep 2005, Trond Hasle Amundsen wrote:
Hi,
I have built R 2.1.1 on Solaris and OSF1, and used the operating
system's native compiler set.
Hi R-users,
I would like to draw a cube with a grid on it and labels along all
three axes. I have trouble printing the labels correctly. My
best attempt is described below.
Can somebody explain me how I can change the 0,20,40,80,100
along the x axis into character vectors like
no, light,
This technote explains the margin area (mar) and how to modify it to control
white space around a graphic:
http://research.stowers-institute.org/efg/R/Graphics/Basics/mar-oma/index.htm
When you have multiple figures on a graphic, you may also want to learn to
control the outer margin area (oma),
Hi Vincent,
I regularly have a look at your web site.
Very big thanks and bravo for this very useful work.
(Just a little remark, why not let the french html
version fully available on line, (findable by google, etc)
This would/could be useful to promote R for french peoples ?)
Many thanks.
That worked and gave me enough information so to make it look exactly
the way I want.
Thanks!
Jamie
On 9/1/05, Chuck Cleland [EMAIL PROTECTED] wrote:
How about this:
par(mar=c(2,4,1,1))
plot(c(1,2,2,3,3,3), type=p, pch=20, ylab=Y Label, xlab=, xaxt=n)
mtext(side=1, line=0.5, X Label)
Similar to my last question, I want to tighten up the spacing and
margins in a plot I am doing with lattice.
Here are the commands I'm using:
data - data.frame(x=c(1:3, 1:3), y=c(1:3, 1:3*2),
cat=c(foo,foo,foo,bar, bar,bar))
xyplot(panel=panel.superpose, y~x, data=data, groups=cat, type=b,
Earl F. Glynn [EMAIL PROTECTED] writes:
AFAIK, the only way to get the axis label closer to the axis is to
suppress the actual axis labels and use the mtext command to display
alternative text where you want it. For example, look at the blue text in
Figure 2B (at the above link) that is
Hello. Following pages 342-347 of Pinheiro Bates, I am trying to
write a self-starting nonlinear function (a non-rectagular hyperbola) to
be used in nonlinear least squares regression (and eventually for a
mixed model). When I use the getInitial function for my self-starting
function I get the
Hi,
I've got a quick question about what happens when indexing into matrices
with a single column. I was wondering if anyone can help ...
For example:
x - matrix(1:10)
y - cbind(x,x)
x[4:6,]
[1] 4 5 6
y[4:6,]
[,1] [,2]
[1,] 4 4
[2,] 5 5
[3,] 6 6
class(x[4:6,])
[1] integer
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] Behalf Of Caio Lucidius
Naberezny Azevedo
Sent: 01 September 2005 12:09
To: Help mailing list - R
Subject: [R] Multivariate Skew Normal distribution
Hi all,
Could anyone tell me if there is any package (or
There's an optional drop argument to the indexing operations, which is by
default TRUE. You just want
x[4:6,,drop=FALSE]
[,1]
[1,]4
[2,]5
[3,]6
Reid Huntsinger
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Crispin Miller
Sent:
On 01-Sep-05 Caio Lucidius Naberezny Azevedo wrote:
Hi all,
Could anyone tell me if there is any package (or function) that
generates values from a multivariate skew normal distribution?
Thanks all,
Caio
Hello, Caio
Please tell us what you mean by skew normal distribution.
Since
Crispin Miller wrote:
Hi,
I've got a quick question about what happens when indexing into matrices
with a single column. I was wondering if anyone can help ...
For example:
x - matrix(1:10)
y - cbind(x,x)
x[4:6,]
[1] 4 5 6
y[4:6,]
[,1] [,2]
[1,] 4 4
[2,] 5
Jamieson Cobleigh wrote:
Similar to my last question, I want to tighten up the spacing and
margins in a plot I am doing with lattice.
Here are the commands I'm using:
data - data.frame(x=c(1:3, 1:3), y=c(1:3, 1:3*2),
cat=c(foo,foo,foo,bar, bar,bar))
xyplot(panel=panel.superpose, y~x,
Thanks everyone!
Crispin
This email is confidential and intended solely for the use o...{{dropped}}
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
Chris Wallace [EMAIL PROTECTED] wrote in message
news:[EMAIL PROTECTED]
how about
plot(..., xlab=)
title(xlab=label text, line=2)
Yes, Chris, I like your idea, especially when I can fix both X and Y axes
at the same time:
plot(0, xlab=,ylab=)
title(xlab=X axis, ylab=Y axis, line=2)
I'd
On 9/1/05, Bill Shipley [EMAIL PROTECTED] wrote:
Hello. Following pages 342-347 of Pinheiro Bates, I am trying to
write a self-starting nonlinear function (a non-rectagular hyperbola) to
be used in nonlinear least squares regression (and eventually for a
mixed model). When I use the
On 9/1/05, Doran, Harold [EMAIL PROTECTED] wrote:
If you are indeed using lme and not lmer then the needed function is
VarCorr(). However, 2 recommendations. First, this is a busy list and
better emails subject headers get better attention. Second, I would
recommend using lmer as it is much
Dear list,
This is a simple question but I spent an hour on it without success. So I
sorry to bother you.
I fit a linear model lm and then did anova.
My idea is to run multiple pairwise comparision for several factor
variables (f1, f2, f3, f4 each with its own levels, say 0, 1, 3 etc
)at
hy,
I need to have the 0 on the bottom left corner of the graph being joined , not
with this little hole between x axis and y axis...
I've saw this question with the answer one time but i'm unable to find it
again..
thks.
guillaume
__
You are correct, VarCorr IS compatible with lmer. It must be what you state
below. I should have saved the message but VarCorr() complained that I was not
dealing with an lme() object. I had multiple packages (nlme, Matrix, mlmRev,
among others) in the search path.
On 9/1/05, Doran, Harold
On Thu, 2005-09-01 at 16:58 +0100, [EMAIL PROTECTED] wrote:
hy,
I need to have the 0 on the bottom left corner of the graph being
joined , not with this little hole between x axis and y axis...
I've saw this question with the answer one time but i'm unable to find
it again..
thks.
OK, my journey to make lab machines automagically install update all
desirable R packages is nearing an end! The only question I have now is
this: How can I control where the system prints the *.Rout file that is
created automatically when the R batch program runs. In man R I don't
find any
In R, sys.on.exit() retrieves the expression stored for use by 'on.exit' in
the function currently being evaluated. In S, you get the list of expressions
for the current frame and its parents. Is there any way I can do the same in
R? If not, can this be added? Maybe 'sys.on.exit' could take a
Paul Johnson wrote:
OK, my journey to make lab machines automagically install update all
desirable R packages is nearing an end! The only question I have now is
this: How can I control where the system prints the *.Rout file that is
created automatically when the R batch program runs.
That did the trick.
Thanks!
Jamie
On 9/1/05, Sundar Dorai-Raj [EMAIL PROTECTED] wrote:
Jamieson Cobleigh wrote:
Similar to my last question, I want to tighten up the spacing and
margins in a plot I am doing with lattice.
Here are the commands I'm using:
data -
Hi I'm working with the function Kest in the package SpatStat (under LINUX
with R 2.1.0). In order to evaluate the statistical significance of my
point pattern I'm doing 999 Montecarlo replications. The script that use
the Kest function runs OK for most of the different point patterns that I
have
Greetings,
I am attempting to build a R package, however the build fails and the
following message is ouput:
C:\ConstellaGroup\EPA\RAGG\packageR CMD BUILD --binary --force RAGG
* checking for file 'RAGG/DESCRIPTION' ... OK
* preparing 'RAGG':
* checking DESCRIPTION meta-information ... ERROR
Zajd, John wrote:
Greetings,
I am attempting to build a R package, however the build fails and the
following message is ouput:
C:\ConstellaGroup\EPA\RAGG\packageR CMD BUILD --binary --force RAGG
* checking for file 'RAGG/DESCRIPTION' ... OK
* preparing 'RAGG':
* checking DESCRIPTION
On 8/31/05, Maciej Kalisiak [EMAIL PROTECTED] wrote:
I'm doing bwplot(x ~ y | z, ...) with lattice, but would like the
z-labels to appear to the *side* of each bwplot, rather than on top...
is this possible?
Not currently (it's on my TODO list).
Deepayan
Most powerful in what way? Quite a lot depends on the jobs you're going to run.
- To run CPU-bound jobs, more CPUs is better. (Even though R doesn't
do threading, you can manually split some CPU-bound jobs in several
parts and run them simultaneously.) Apart from multiple CPUs and
scaling-4
xywidth-480
resolution-150
png(filename = c:/r/anschluss/plots/4.png, width = xywidth*scaling,
height = xywidth*scaling,pointsize = 12, bg = white, res =
resolution*scaling)
..
barplot(xrow,col = barcolors,cex.axis=scaling, ylab=mean time till attachment
in
On Thu, 2005-09-01 at 21:51 +0200, Knut Krueger wrote:
scaling-4
xywidth-480
resolution-150
png(filename = c:/r/anschluss/plots/4.png, width = xywidth*scaling,
height = xywidth*scaling,pointsize = 12, bg = white, res =
resolution*scaling)
..
barplot(xrow,col =
dear R wizards:
I am getting to play more and more with fun fonts (irony warning). I
now know that I can safely use my TeXtext encoding with the postscript
device, but not with the pdf device. Unfortunately, I believe that the
postscript device does not support translucent colors---or is
Platform: Windows
R version: 2.1.1
This is a request that in the savePlot function, type=eps be included
in the next release. This would be an alias for type=ps but with a
different file extension.
Rationale: The current version of R outputs excellent EPS files. It is
strictly correct to
Jia-Shing,
I missed your original message, but would like to reiterate Bogdan's
comments and suggestions.
In a former life, a colleague of mine led the way for us to construct a
small farm of Opteron servers that all had 2 AMD64 CPU's, SUSE Enterprise
Server OS, and the ability to have up to 16GB
Can someone provide help to resolve this build failure?
Thank you,
John Zajd
C:\ConstellaGroup\EPA\RAGG\packagercmd check RAGG
* checking for working latex ... OK
* using log directory 'C:/ConstellaGroup/EPA/RAGG/package/RAGG.Rcheck'
* using R version 2.1.0, 2005-04-18
* checking for file
On Thu, 1 Sep 2005 [EMAIL PROTECTED] wrote:
On 01-Sep-05 Caio Lucidius Naberezny Azevedo wrote:
Hi all,
Could anyone tell me if there is any package (or function) that
generates values from a multivariate skew normal distribution?
Thanks all,
Caio
Hello, Caio
Please tell us what you
Folks:
In answer to a query, Andy Liaw recently submitted some code to construct a
block diagonal matrix. For what seemed a fairly straightforward task, the
code seemed a little overweight to me (that's an American stock analyst's
term, btw), so I came up with a slightly cleaner version (with
y8011Koh?=
Date: Fri, 02 Sep 2005 08:24:49 +0800
MIME-Version: 1.0
Content-Type: multipart/alternative;
boundary==_NextPart_001_0018_01C39816.2C871AA0
X-Priority: 3
X-MSMail-Priority: Normal
X-Mailer: Microsoft Outlook Express 6.00.2800.1106
X-MimeOLE: Produced By Microsoft MimeOLE
I haven't seen a reply to this question, so I will attempt a few
remarks in spite of some confusion about what you are asking.
1. The function to use for parameter estimation depends on ths
structure of the data. My all-around preference for many purposes is
for optim,
1. I just got 10 hits from 'RSiteSearch(convolution of exponential
and normal)', at least two of which look like they might interest you.
2. I similarly got 27 hits from 'RSiteSearch(convolution
distribution)', several of which mention the distr package, which you
might
my last post was filtered,so I post it again with another title.
If I do not make a mistake,the partial association model is an
extension of log-linear model.I read a papers which gives an example
of it.(Sloane and Morgan,1996,An Introduction to Categorical Data
Analysis,Annual Review of
Probably a better first question is, why are you using a bitmapped
graphics format if you need the image to be re-scaled?
I need a 1000 dpi tif file in a size of appr. 10 to 10 cm for applied
animal behaviour science:
http://authors.elsevier.com/GuideForAuthors.html?PubID=503301dc=GFA
89 matches
Mail list logo