Re: [R] effects in ANCOVA

2006-11-18 Thread Tomas Goicoa
Dear Chuck, thank you very much indeed. I was looking for that and I could not find it. Cheers Tomas __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Problems in plot.lm with option which=5

2006-11-18 Thread Prof Brian Ripley
[PLease don't send HTML mail: your message was badly formatted on arrival. See the posting guide and the double-spaced mess below.] This _is_ a bug: the author had tried to reorder the factor levels by mean fitted value, but forgot to reorder the residuals (etc) to match. You may have noticed

Re: [R] Gantt chart problem after upgrade to R 2.4.0

2006-11-18 Thread Jim Lemon
John Kane wrote: I am having a problem with a gantt chart since moving to R2.4.0. from 2.3.1 ... Okay, I think I have fixed the problem. I don't yet know why it worked on Windows and not Linux, but this should work on both. I have added another argument to the gantt.chart function and it

[R] Random sample from log-normal distribution

2006-11-18 Thread Megh Dal
Dear all R users, Please forgive me if my question is too trivial. Suppose I have two variables, (x,y) which is log-normally distributed with expected value (mu1, mu2) and some variance-covariance matrix. Now I want to draw a random sample of size 1000 from this distribution. Is there any

[R] (no subject)

2006-11-18 Thread Céline Henzelin
Hello, I need help to understand my error in this code... I would like to make a direct sampler... thanks Celine p-matrix(c(1,2,3,2,1,4),3,2,byrow=T) p-p/sum(p) ky-ncol(p) kx-nrow(p) p.vec-as.vector(t(p)) l-rmultinom(1,1,p) l-(t(l)) m-(l%%ky) ifelse (m==0, i-l/ky (j-ky), i-((l/ky)+1)

Re: [R] deriv when one term is indexed

2006-11-18 Thread Gabor Grothendieck
This works for me in terms of giving results without error messages except for the confint(dd.plin) which I assume you don't really need anyways. gg - model.matrix(~ Gun/GL - Gun, dd) dd.plin - nls(Lum ~ gg^gamm, dd, start = list(gamm = 2.4), +alg = plinear +

Re: [R] deriv when one term is indexed

2006-11-18 Thread Gabor Grothendieck
I mixed up examples. Here it is again. As with the last one confint(dd.plin) gives an error (which I assume is a problem with confint that needs to be fixed) but other than that it works without issuing errors and I assume you don't need the confint(dd.plin) in any case since dd.plin is just

Re: [R] Random sample from log-normal distribution

2006-11-18 Thread jim holtman
?rlnorm On 11/18/06, Megh Dal [EMAIL PROTECTED] wrote: Dear all R users, Please forgive me if my question is too trivial. Suppose I have two variables, (x,y) which is log-normally distributed with expected value (mu1, mu2) and some variance-covariance matrix. Now I want to draw a random

Re: [R] (no subject)

2006-11-18 Thread David Barron
I'm not sure what the last line is trying to achieve, but this might be what you want: i - ifelse(m==0, l/ky, l/ky + 1) j - ifelse(m==0,ky, l %% ky) On 18/11/06, Céline Henzelin [EMAIL PROTECTED] wrote: Hello, I need help to understand my error in this code... I would like to make a direct

Re: [R] Random sample from log-normal distribution

2006-11-18 Thread Prof Brian Ripley
On Sat, 18 Nov 2006, jim holtman wrote: ?rlnorm Not for a bivariate lognormal, nor specifying via mean and variance (which is not the normal way to specify a univariate lognormal). I think we need clarification of exactly what is meant, but the answer is most likely 'no'. On 11/18/06, Megh

[R] .First.lib

2006-11-18 Thread Charles Annis, P.E.
Greetings: I’m running R version 2.4.0 (2006-10-03) on WindowsXP in a 3 year old DELL box with 2 gig. I have successfully created a zipped package, mh1823_2.1.zip, using R CMD build --binary mh1823. R CMD build --binary automatically adds a file, mh1823, containing .First.lib and its default

Re: [R] Gantt chart problem after upgrade to R 2.4.0

2006-11-18 Thread John Kane
--- Jim Lemon [EMAIL PROTECTED] wrote: John Kane wrote: I am having a problem with a gantt chart since moving to R2.4.0. from 2.3.1 ... Okay, I think I have fixed the problem. I don't yet know why it worked on Windows and not Linux, but this should work on both. Looks pretty

Re: [R] Gantt chart problem after upgrade to R 2.4.0

2006-11-18 Thread John Kane
--- No No [EMAIL PROTECTED] wrote: I am on ubuntu linux with R 2.4.0 compiled and I get the same picture as with yours with R 2.4.0. Thanks. From another reply it looks like this is a 2.4.0 problem that is solved in 2.4.0-patched and not a Windows vs Linux problem. I upgraded to the patch

Re: [R] Random sample from log-normal distribution

2006-11-18 Thread Ted Harding
On 18-Nov-06 Megh Dal wrote: Dear all R users, Please forgive me if my question is too trivial. Suppose I have two variables, (x,y) which is log-normally distributed with expected value (mu1, mu2) and some variance-covariance matrix. Now I want to draw a random sample of size 1000 from

Re: [R] .First.lib

2006-11-18 Thread Prof Brian Ripley
On Sat, 18 Nov 2006, Charles Annis, P.E. wrote: Greetings: I’m running R version 2.4.0 (2006-10-03) on WindowsXP in a 3 year old DELL box with 2 gig. I have successfully created a zipped package, mh1823_2.1.zip, using R CMD build --binary mh1823. R CMD build --binary automatically adds a

[R] lattice layout multiple pages

2006-11-18 Thread Patrick Giraudoux
Dear all, Using lattice I would like to print a conditionnal plot of 32 panels in a limited number of panels (eg 4) in each of several pages: xyplot(pds~time|Idnid,data=croispond3,layout=c(4,4)) This works well in principle, but the pages are printed without any possibility of stopping the

[R] Questions regarding integrate function

2006-11-18 Thread Le Wang
Hi there. Thanks for your time in advance. I am using R 2.2.0 and OS: Windows XP. My final goal is to calculate 1/2*integral of (f1(x)^1/2-f2(x)^(1/2))^2dx (Latex codes: $\frac{1}{2}\int^{{\infty}}_{\infty} (\sqrt{f_1(x)}-\sqrt{f_2(x)})^2dx $.) where f1(x) and f2(x) are two marginal densities.

[R] sorry for repeated msgs

2006-11-18 Thread Le Wang
Hi there, Sorry for the repeated msgs. I set my option to receive my own email, but it didn't seem to work. So, I didn't realize my previous posts actually got through. Le - ~~~ Le Wang, Ph.D Population Center University of Minnesota

Re: [R] .First.lib

2006-11-18 Thread Charles Annis, P.E.
Thank you Professor Ripley. I don't have compiled code and do not have a namespace, which is required (I understand) to use onLoad. The package is of moderate size, however with about 100 R functions. I know that I could emulate the technique used by Rcmdr, but had hoped to use R CMD build

[R] Multinomial sampling WAS: Re: (no subject)

2006-11-18 Thread Charles C. Berry
Celine, The posting guide directs you to Use an informative subject line (not something like `question'). Doing so will increase the probability of getting an informative response. You need to spend a more time with the help pages and with An Introduction to R. -- Specifically for the

Re: [R] Random sample from log-normal distribution

2006-11-18 Thread Juan Pablo Lewinger
Dear all R users, Please forgive me if my question is too trivial. Suppose I have two variables, (x,y) which is log-normally distributed with expected value (mu1, mu2) and some variance-covariance matrix. Now I want to draw a random sample of size 1000 from this distribution. Is there any

Re: [R] .First.lib

2006-11-18 Thread Charles Annis, P.E.
Perseverance, I guess, is the answer. Thank you exceedingly Professor Ripley for your help. Things do look awfully simple on this side of success. At your suggestion (after several false starts) I defined a .First.lib function within a file named zzz.R and built with R CMD build --binary,

Re: [R] Questions regarding integrate function

2006-11-18 Thread Matthias Kohl
Hello, your integrand needs to be a function which accepts a numeric vector as first argument and returns a vector of the same length (see ?integrate). Your function does not fulfill this requirement. Hence, you have to rewrite your function or use sapply, apply or friends; something like

Re: [R] Gantt chart problem after upgrade to R 2.4.0

2006-11-18 Thread No No
It is working on linux fine now. __ 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 commented, minimal, self-contained,

Re: [R] Trellis Plot Labels

2006-11-18 Thread Deepayan Sarkar
On 11/17/06, Deepayan Sarkar [EMAIL PROTECTED] wrote: On 11/17/06, Turgut Durduran [EMAIL PROTECTED] wrote: On 11/17/06, Turgut Durduran [EMAIL PROTECTED] wrote: Hello everyone, I am ploting a groupeddata object with formula: formula(mydatausegroup) BF ~ HO | ID/Infar/Day

Re: [R] Gantt chart problem after upgrade to R 2.4.0

2006-11-18 Thread John Kane
--- [EMAIL PROTECTED] wrote: This is interesting. I am on Win2000 running R version 2.4.0 Patched (2006-11-15 r39915) and plotrix 2.1.5. I get the plot that scaled a little differently that the 2.3.1 example but the dates look fine. My fault there. I simply copied the metafile and

[R] deriv when one term is indexed

2006-11-18 Thread Ken Knoblauch
I don't know why I get a different result than you do, though I'll play around with it some and check the code to see if I can figure out what is going on, on my end. I will add that by transforming GL so that it spans the interval [0, 1] instead of [0, 255], the correlations among the parameters

Re: [R] lattice layout multiple pages

2006-11-18 Thread Deepayan Sarkar
On 11/18/06, Patrick Giraudoux [EMAIL PROTECTED] wrote: Dear all, Using lattice I would like to print a conditionnal plot of 32 panels in a limited number of panels (eg 4) in each of several pages: xyplot(pds~time|Idnid,data=croispond3,layout=c(4,4)) This works well in principle, but the

Re: [R] gzfile with multiple entries in the archive

2006-11-18 Thread Duncan Temple Lang
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 Apologies for entering this late, but last week was extremely busy. I hadn't realized that you would implement something so quickly and I could have saved you some time. I was in the process of adding facilities for gzipped tar files in the

[R] Ryacas not working properly

2006-11-18 Thread Paul Smith
Dear All I have just installed the package Ryacas, but getting the following: library(Ryacas) Loading required package: XML yacas(1/2 * 3/4) [1] Starting Yacas! Error in socketConnection(host = 127.0.0.1, port = 9734, server = FALSE, : unable to open connection In addition: Warning

[R] questions on adding reference line?

2006-11-18 Thread zhijie zhang
Dear Ruser, I use abline() function to add the reference line successfully, but i can't display the values corresponding to the reference line on the x/y axis, anybody knows how to display it? *My simulated programs:* y-rnorm(50) plot(x,y) abline(v=0.5) *#my question is how to display x=0.5

Re: [R] Ryacas not working properly

2006-11-18 Thread Gabor Grothendieck
Checkout the INSTALLATION - UNIX and TROUBLESHOOTING sections of the home page: http://code.google.com/p/ryacas/ On 11/18/06, Paul Smith [EMAIL PROTECTED] wrote: Dear All I have just installed the package Ryacas, but getting the following: library(Ryacas) Loading required package: XML

Re: [R] questions on adding reference line?

2006-11-18 Thread Gabor Grothendieck
Try: axis(1, 0.5, tick = TRUE) On 11/18/06, zhijie zhang [EMAIL PROTECTED] wrote: Dear Ruser, I use abline() function to add the reference line successfully, but i can't display the values corresponding to the reference line on the x/y axis, anybody knows how to display it? *My simulated