On Thu, 20 Jul 2006, Carlos J. Gil Bellosta wrote:
Dear Rusers,
I have developed two packages for a client of mine. After new features
are added or bugs corrected, I upload them to my own web repository. I
create both source and binary versions.
binary Linux packages, it seems. The latter
On Wed, 19 Jul 2006 11:01:02 -0500,
Sundar Dorai-Raj (SD) wrote:
Dieter Menne wrote:
Dear R-Listeners,
as the Sweave faq says:
http://www.ci.tuwien.ac.at/~leisch/Sweave/FAQ.html
creating several figures from one figure chunk does not work, and for
standard
Hallo
Is there an elegant way to do the following:
Dataset consists of 2 variables: var1: some measurements, and var2: a grouping
variable with two values, 1 and 2.
There are (say) 10 measurements from group 1 and 15 measurements from group 2.
The idea is to study the permutation distribution
data1 - expand.grid(var1=1:15, var2=1:2)
test - replicate(1000, with(data.frame(var1=data1$var1,
var2=sample(data1$var2)),
diff(tapply(var1, var2, mean
hist(test)
---
Jacques VESLOT
CNRS UMR 8090
I.B.L (2ème étage)
1
Hi all,
I wish to draw 2 hist (or barplot) on the same graph.
I can do it simply by using par(new=TRUE) , but it overlap with the
first drawn, how to tell R to put in front of the graph the min value of
the two graph in order for it to be seen and don't hide each other.
I've been
SpG == Spencer Graves [EMAIL PROTECTED]
on Wed, 19 Jul 2006 15:23:58 -0700 writes:
SpG Hi, Thomas: Thanks very much. I haven't tried it yet,
SpG but it looks very useful. Best Wishes, Spencer Graves
Hmm, ?methods has been containing for a while
methods Note:
methods
methods
s1 - runif(10,0,10)
s1
[1] 8.328396 2.840870 7.401377 9.998165 5.045539 9.568728 5.372493 5.232439
[9] 5.774790 4.224103
s2 - runif(10,0,10)
s2
[1] 1.230750 3.855060 8.652698 7.846725 9.100171 7.309179 9.235562 1.581741
[9] 6.979521 1.918997
ss - cbind(s1,s2)
smin -
Thank you for your answers.
So, I have to derive variables, reclassify, binning, aggregate, or merge
variables (data management). After that I have to do a logistic regression
to do scoring.
I have tried a MySQL database with RODBC but I did a mistake, I have 33
millions rows table and not 3
On 7/19/06, hadley wickham [EMAIL PROTECTED] wrote:
Can you be a bit more excact? I a biologist and relatively new to R
In that case, I would _strongly_ advise that you get advice from a
local statistician.
I am afraid that, by comparison, I am the local statistican. I am also
the local
Dear all,
I would like to draw a dot chart on a log scale.
What is the syntax for this? A barchart may use
log=x, but trying this with dotchart() leads
to an error message.
Greetings
Johannes
__
R-help@stat.math.ethz.ch mailing list
Thank you very much, Gabor Grothendieck and Muhammad Subianto. Both of
these work perfectly. I think I was misunderstanding gp and gpar()
before. Again, thank you both.
-KZ
Quoting Muhammad Subianto [EMAIL PROTECTED]:
Maybe like this:
mosaic(allmorph, direction = v, pop = FALSE,
I wrote a function f whose argument is an array
If I set
x - c(5,1,2)
and I run f(x) everything works fine.
On the other hand, If I extract the same array
from a matrix by:
x - my.matrix[1,1:3]
such that the first three elements of the
first row of my.matrix are exactly 5 1 2,
something
I tried to make the following function:
function(x, y){
dates-intersect(x[,1],y[,1])
m-matrix(NA,length(dates),3)
m[,1]-dates
j-1
k-1
for(i in fdax[,1]){
if(is.element(i,dates)){
m[j,3]-as.numeric(fdax[k,2])
j-j+1
}
k-k+1
}
return(m)
}
Dear Rusers,
Well, then it seems that the problem is that I am building linux binary
packages. Since I do not have any compiled code within --just R code--, their
contents should --and, in fact, are-- directly installable on Windows
platforms
(which is what I intend to do).
If I understand
Sorry to disturb again the list subscribers,
but I fixed my problem, by simply using the c
operator:
x - c(my.matrix[1,1:3])
Alessandro
* Alessandro Antonucci [EMAIL PROTECTED] [200706, 12:12]:
I wrote a function f whose argument is an array
If I set
x - c(5,1,2)
and I run f(x)
Ivan Kalafatic [EMAIL PROTECTED] writes:
I tried to make the following function:
function(x, y){
dates-intersect(x[,1],y[,1])
m-matrix(NA,length(dates),3)
m[,1]-dates
j-1
k-1
for(i in fdax[,1]){
if(is.element(i,dates)){
m[j,3]-as.numeric(fdax[k,2])
I am afraid that, by comparison, I am the local statistican. I am also
the local R-guru, and neither is saying much - so please bear with
me.
Do you know of some functions (built in hopefully) that I can try?
I'm a bit leery of offering advice without really sitting down and
discussing your
On Thu, 20 Jul 2006, Carlos J. Gil Bellosta wrote:
Dear Rusers,
Well, then it seems that the problem is that I am building linux binary
packages. Since I do not have any compiled code within --just R code--, their
contents should --and, in fact, are-- directly installable on Windows
Ivan Kalafatic wrote:
I tried to make the following function:
function(x, y){
dates-intersect(x[,1],y[,1])
m-matrix(NA,length(dates),3)
m[,1]-dates
j-1
k-1
for(i in fdax[,1]){
if(is.element(i,dates)){
m[j,3]-as.numeric(fdax[k,2])
j-j+1
List:
Thank you for the replies to my post yesterday. Gabor and Phil also gave
useful replies on how to improve the function by relying on mapply
rather than the explicit for loop. In general, I try and use the family
of apply functions rather than the looping constructs such as for, while
etc as
I would like to draw a dot chart on a log scale.
What is the syntax for this? A barchart may use
log=x, but trying this with dotchart() leads
to an error message.
You can do this easily with ggplot:
install.packages(ggplot)
library(ggplot)
qplot(mpg, factor(cyl), data=mtcars, log=x)
Il giorno mar, 18/07/2006 alle 12.25 -0300, raul sanchez ha scritto:
how can I select the countries of Lat. Am. thank you
hint: add to the data.frame a column with the continental ID so this
could be usefull for other needs.
hint: ?subset
--
Daniele Medri
Note that if you use mapply in the way I suggested, which is not
the same as in your post, then its just as fast. (Also the version
of mapply in your post gives different numerical results than
the for loop whereas mine gives the same.) like.mat is the for
loop version, like.mat2 is your mapply
By far, the cheapest and easiest solution (and the very first to try)
is to add more memory. The cost depends on what kind you need, but
here's for example 2 GB you can buy for only $150:
http://www.newegg.com/Product/Product.asp?Item=N82E16820144157
Project constraints?! If they don't want to
On Thu, 20 Jul 2006, Jacob van Wyk wrote:
Hallo
Is there an elegant way to do the following:
Dataset consists of 2 variables: var1: some measurements, and var2: a
grouping variable with two values, 1 and 2.
There are (say) 10 measurements from group 1 and 15 measurements from group 2.
Gabor Grothendieck wrote:
Note that if you use mapply in the way I suggested, which is not
the same as in your post, then its just as fast. (Also the version
of mapply in your post gives different numerical results than
the for loop whereas mine gives the same.) like.mat is the for
loop
Dear all
I have Kubuntu linux and have updated to the latest version (6.06 dapper). I
do not know why but now I can not make no plots. For instance, when I type
hist(...)
this is the message I get:
can't find X11 font
Error in X11 (display, width, heigth..)
unable to start devide X11
Yes, that significantly improves the speed so that the for loop and the
mapply function both return the same CPU time.
Also, thank you for your diagnosis of the likelihood matrix.
Harold
-Original Message-
From: Gabor Grothendieck [mailto:[EMAIL PROTECTED]
Sent: Thursday, July 20,
On Thu, 20 Jul 2006, Kie Zuraw wrote:
Thank you very much, Gabor Grothendieck and Muhammad Subianto. Both of
these work perfectly. I think I was misunderstanding gp and gpar()
before. Again, thank you both.
Two further additions to this:
Maybe like this:
mosaic(allmorph, direction = v,
Dear all,
I am unsure about how to specify a model in R and I thought of asking
some advice to the list. I have two groups (Group= A, B) of
subjects, with each subject undertaking a test before and after a
certain treatment (Time= pre, post). Additionally, I want to enter
the age of the
Hello everybody,
i'm having some trouble performing correspondence analysis with R for Mac OS
X. Does anyone know about some useful package?
And also, if i had found coordinates of points representing row and column
profiles, how do i put them in a single figure with labels identifying each
Hi, Martin:
I looked for 'showMethods' before I posted this question. I see two
possible reasons others might not have mentioned it previously -- and
why my post didn't mention it:
HARD TO FIND
Now that you have mentioned it, I recall having stumbled across it
and used
Using R on windows
I have the following code
for(i in 1:1) {
draw some random weights
perfom a weighted least squares regression
some simple addition and multiplication
}
The code works fine but is slow.
I have mingw installed and can dyn.load, although I am
more used to doing
Without details it is hard to say, but if you do the WLS via lm.wfit,
probably not. If you are using lm() and taking the overhead 1 times,
that is the first thing to avoid.
As an example of this sort of thing, look at the C code for lqs{MASS}
which does thousands of regressions quite
Notice that in Martin's example he used
showMethods(class = lmer, where = package:Matrix)
The result is still not optimal but it is slightly more useful than
trying the default location. The approach that Thomas contributed
earlier in this discussion is probably the better route.
On 7/20/06,
On 7/20/2006 10:38 AM, Emanuele Mazzola wrote:
Hello everybody,
i'm having some trouble performing correspondence analysis with R for Mac OS
X. Does anyone know about some useful package?
And also, if i had found coordinates of points representing row and column
profiles, how do i put
On Thu, 20 Jul 2006, Douglas Bates wrote:
Notice that in Martin's example he used
showMethods(class = lmer, where = package:Matrix)
The result is still not optimal but it is slightly more useful than
trying the default location. The approach that Thomas contributed
earlier in this
Hi, Hi all,
i'm having some trouble performing correspondence analysis with R for Mac OS
X. Does anyone know about some useful package?
And also, if i had found coordinates of points representing row and column
profiles, how do i put them in a single figure with labels identifying each
one
When printing a table it is broken at some point (depending how long are
the associated names)
see example below.
Is there a way to control number of columns being printed for a given
chunk of the table?
Best regards,
Ryszard
z5
AAA BBB CCC DDD EEE FFF
This is controlled by options(width): characters not columns.
On Thu, 20 Jul 2006, [EMAIL PROTECTED] wrote:
When printing a table it is broken at some point (depending how long are
the associated names)
see example below.
Is there a way to control number of columns being printed for a
Hi William,
Thanks a lot for your response. I checked the package and found that what I
want to solve was the opposite, that is, from mean and sd to parameters shape
and scale. Could anyone give some hints please? Any suggestion would be
appreciated!
Leaf
- Original Message -
Dear all,
I apologize if this is a FAQ (seems a bit like one, but I didn't find
anything).
I'm looking for an easy way to cut one value out of a vector and shorten
the vector accordingly. Something like:
x - c(1, 1, 0, 6, 2)
throwaway(x[3])
which will return x = 1 1 6 2, with length(x) = 4. I
In some cases it may be sufficient to abbreviate the colnames:
library(MASS); data(survey)
head(survey)
Sex Wr.Hnd NW.Hnd W.HndFold PulseClap Exer Smoke Height M.I
1 Female 18.5 18.0 Right R on L92Left Some Never 173.00 Metric
2 Male 19.5 20.5 Left R on
John Wiedenhoeft wrote:
Dear all,
I apologize if this is a FAQ (seems a bit like one, but I didn't find
anything).
I'm looking for an easy way to cut one value out of a vector and shorten
the vector accordingly. Something like:
x - c(1, 1, 0, 6, 2)
throwaway(x[3])
which will return
John Wiedenhoeft wrote:
Dear all,
I apologize if this is a FAQ (seems a bit like one, but I didn't find
anything).
I'm looking for an easy way to cut one value out of a vector and shorten
the vector accordingly. Something like:
x - c(1, 1, 0, 6, 2)
throwaway(x[3])
which will return
Hi all,
thank you for your answers; i've tried both cca from vegan library, and
dudi.coa from ade4 library; one last question: my deal is mainly with
contingency tables, like the one i'm posting here
acciaieria-c(.41,.02,.44,.04,.09)
laminatoio-c(.34,.28,.26,.01,.11)
On Wed, 19 Jul 2006, Ulrik Stervbo wrote:
[much deleted]
I am measureing the amount of DNA in cells, and I need to know the
percentage of cells in a part of the cell cycle; that the percentage
of cells in the first peak, in the second peak and so on. I want to
integrate the area between to
Try to insert a space before the last bracket! R
sometimes does not accept a curly bracket as a first
character on a line of a script.
Regards,
Valentin
--- Ivan Kalafatic [EMAIL PROTECTED] wrote:
I tried to make the following function:
function(x, y){
dates-intersect(x[,1],y[,1])
I just got 145 hits for RSiteSearch(ecology, functions) and 28 for
RSiteSearch(Moran's I, functions). Might any of these help you?
Spencer Graves
Ivan Rubio Perez wrote:
Hi!
Dear All,
I want to implement an autocorrelation analysis and estimated the Moran's I
in a set of
Thanks for providing such a self-contained example by which 'nlme'
crashes R. Could you please also give us 'sessionInfo()'? I don't have
time to test it myself now, but perhaps if you identify your platform,
you might interest someone else in checking it.
I'm sorry I
On Thu, 20 Jul 2006, Jacob van Wyk wrote:
Hallo
Is there an elegant way to do the following:
Dataset consists of 2 variables: var1: some measurements, and var2: a
grouping variable with two values, 1 and 2.
There are (say) 10 measurements from group 1 and 15 measurements from group 2.
Easier than I thought. Thanks a lot!
Am Donnerstag, den 20.07.2006, 17:25 +0100 schrieben zig Leute ;-):
x[-3]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
Thanks Spencer. Here is my sessionInfo():
Version 2.3.1 (2006-06-01)
i386-pc-mingw32
attached base packages:
[1] methods stats graphics grDevices utils datasets
[7] base
other attached packages:
nlme
3.1-75
On Thu, 20 Jul 2006, Spencer Graves wrote:
Thanks
Dear R-Helpers,
Can anyone please help me to interpret warning messages from zeroinfl
(package pscl) while fitting a zero inflated negative binomial model?
The console reports convergence and the parameters seam reasonable, but
these
Warning messages:
1: algorithm did not converge in:
Bernat -
Just a guess, but maybe you need to install the packages
xfonts-100dpi
and
xfonts-75dpi
- Phil
On Thu, 20 Jul 2006, Bernat Claramunt i Lopez wrote:
Dear all
I have Kubuntu linux and have updated to the
Hi there,
i´m have vector of kernels. just like:
kernels = c('gauss','epan','rectangular')
i know there are density.default$kernels, but thats not my question
here. my own kernel functions are running and working.
my problem is the following is not working:
dev.off()
par(mfrow=c(3,3))
Running R.app on Mac OS X 10.4
version
_
platform powerpc-apple-darwin8.6.0
arch powerpc
os darwin8.6.0
system powerpc, darwin8.6.0
status
On 7/18/2006 2:28 PM, [EMAIL PROTECTED] wrote:
I am having trouble fitting correlation structures within nlme. I would like
to
fit corCAR1, corGaus and corExp correlation structures to my data. I either
get the error step halving reduced below minimum in pnls step or
alternatively R
On 7/20/2006 3:49 PM, Steven McKinney wrote:
Running R.app on Mac OS X 10.4
version
_
platform powerpc-apple-darwin8.6.0
arch powerpc
os darwin8.6.0
system powerpc, darwin8.6.0
On 7/20/2006 3:42 PM, bunny , lautloscrew.com wrote:
Hi there,
i´m have vector of kernels. just like:
kernels = c('gauss','epan','rectangular')
i know there are density.default$kernels, but thats not my question
here. my own kernel functions are running and working.
my problem is the
On 7/20/06, bunny , lautloscrew.com [EMAIL PROTECTED] wrote:
Hi there,
i´m have vector of kernels. just like:
kernels = c('gauss','epan','rectangular')
i know there are density.default$kernels, but thats not my question
here. my own kernel functions are running and working.
my problem is
@yahoo.ca
Subject: Re: [R] Weibull distribution
Message-ID: [EMAIL PROTECTED]
X-mailer: Foxmail 6, 3, 103, 21 [cn]
Mime-Version: 1.0
Content-Type: multipart/alternative;
boundary==003_Dragon527446281311_=
This is a multi-part message in MIME format.
To all:
Genentech has immediate openings for two MS/PhD statisticians in its
preclinical/nonclinical statistics group in South San Francisco, CA. As the
name indicates, this group provides statistical services to all of
Genentech's RD, manufacturing, and marketing activities **except** clinical
I´m working on an R-implementation of the simulation-based finite-sample
null-distribution of (R)LR-Test in Mixed Models (i.e. testing for
Var(RandomEffect)=0) derived by C. M. Crainiceanu and D. Ruppert.
I'm in the beginning stages of this project and while comparing quick and dirty
Duncan, i could not find one seed which caused both corGaus and corExp to crash
in R-patched. but, i found a seed that caused each to fail individually.
Thanks for your help.
# For corExp:
set.seed(26)
for(i in 1:length(CO2$conc)){
CO2$conc[i]-(CO2$conc[i]+rnorm(1))
}
fm1CO2.lis -
As others have mentioned its not really a good idea
to modify the namespace of a package and writing
a wrapper as Duncan suggested is much preferable.
An intermediate approach that
is not as good as the wrapper but better than modifying
the namespace is to copy the objects of interest
to your
R tries to use the maximum precision (64 bit mantissa) in the floating
point registers when it can. When it stores results to memory, they are
stored in double precision (53 bit mantissa). There's unlikely to be
anything specific about conversion to a list that lost the precision,
but I'm
On Thu, 20 Jul 2006, [EMAIL PROTECTED] wrote:
Duncan, i could not find one seed which caused both corGaus and corExp
to crash in R-patched. but, i found a seed that caused each to fail
individually. Thanks for your help.
Running these under Valgrind they both show the same problem
On 7/20/2006 7:16 PM, Thomas Lumley wrote:
On Thu, 20 Jul 2006, [EMAIL PROTECTED] wrote:
Duncan, i could not find one seed which caused both corGaus and corExp
to crash in R-patched. but, i found a seed that caused each to fail
individually. Thanks for your help.
Running these under
Hello,
I looked for R packages which focused on order-restricted statistical
inference, but I could find only the isoreg() function.
I would need to test whether the means in my (repeated measures) data follow
a given order, e.g. AB=CD.
I took a look at the monograph by Barlow et al. (1972) on
Dear R mailing group,
The second parameter for the function optim()is a function whose
parameters are to be optimized. The description of this function given in
the help file is the following:
fn: A function to be minimized (or maximized), with first
argument the vector of
Consider the following example (based on an example in Pat Altham's GLM
notes)
pyears - scan()
18793 52407 10673 43248 5710 28612 2585 12663 1462 5317
deaths - scan()
2 32 12 104 28 206 28 186 31 102
Smoke - gl(2,1,10,labels=c(No,Yes))
Age - gl(5,2,10,labels=c(35-44,45-54,55-64,65-74,75-84),
72 matches
Mail list logo