Re: [R] Making a Package

2005-02-15 Thread Gabor Grothendieck
michael watson (IAH-C bbsrc.ac.uk> writes: : : I follow the guide here and it never fails: : : http://www.biostat.jhsph.edu/~kbroman/Rintro/Rwinpack.html : : So if I have the MyPackage directory, that was created with package.skeleton and subsequently edited, I : will cd to the directory th

Re: [R] lattice multiple plots per page

2005-02-15 Thread Paul Murrell
Hi Deepayan Sarkar wrote: On Monday 14 February 2005 18:12, Jeff Jorgensen wrote: Dear R-sters, I was wondering if anyone has encountered the following issues. I've figured out how to get multiple levelplots [library(lattice)] on a single plot. However, when I add text (adding axis labels for the

RE: [R] reading a 40kb csv file in R

2005-02-15 Thread Mulholland, Tom
I think you may need to check the list for posts relating to "memory", "large datasets", etc http://cran.r-project.org/search.html It would help those on the list, if you gave more details about the machine you are using. Given the vintage of R you areusing is considered quite old, there is th

Re: [R] Could anyone answer for the following question

2005-02-15 Thread Gabor Grothendieck
Cuichang Zhao yahoo.com> writes: : I am using R 2.0.1 under Windows XP. I want to write a function that : makes four graphs and stores each of them in graphics history. When the : function finishes, in other words, I want its graphical output to be : stored in a way that I can look at it usin

RE: [R] Could anyone answer for the following question

2005-02-15 Thread Berton Gunter
options(graphics.record=TRUE) turns on graphics recording. This can also be done from the History menu of the graphics window. All subsequent plots to the graphics device will be saved. -- Bert Gunter Genentech Non-Clinical Statistics South San Francisco, CA "The business of the statistician is

Re: [R] Problem with

2005-02-15 Thread Gabor Grothendieck
I am just in the middle of a conversation on the tex list but just in case you are about to post anything to your web site here is a slight simplification. It consists of updating the MiKTeX file database using initexmf instead of the GUI so you don't have to move back and forth during the ins

Re: [R] Tcl/tk

2005-02-15 Thread Peter Dalgaard
[EMAIL PROTECTED] writes: > I am new to R so I apologize if this is a simple question > > I have created a tkgrid with a button that, once pushed, will extract the > information from the grid, but I can not figure out the command to save the > entries from my grid. > Any thoughts? Several, all s

Re: [R] summary(aov(...)) into a string?

2005-02-15 Thread Peter Dalgaard
Prof Brian Ripley <[EMAIL PROTECTED]> writes: > It doesn't print anything: the summary.aov (or summary.aovlist) print > method does. > > ?summary.aov tells you the structure of the objects they return. Yes. but wouldn't capture.output() and par(family="mono") be closer to the mark? > On Tu

Re: [R] special symobol / character

2005-02-15 Thread Paul Murrell
Hi FISCHER, Matthew wrote: Hi Eric, Thanks for your speedy reply, I should have noted that I'm using a Linux machine. However, when I copy the symbol from Windows to Linux (using R/emacs) via an x-win 32 window it replaces the per mille symbol with a /211. R then produces the character (not a

[R] Could anyone answer for the following question

2005-02-15 Thread Cuichang Zhao
Hello , could anyone answer for the following question for me: I am using R 2.0.1 under Windows XP. I want to write a function that makes four graphs and stores each of them in graphics history. When the function finishes, in other words, I want its graphical output to be stored in a way

RE: [R] Making a Package

2005-02-15 Thread michael watson \(IAH-C\)
I follow the guide here and it never fails: http://www.biostat.jhsph.edu/~kbroman/Rintro/Rwinpack.html So if I have the MyPackage directory, that was created with package.skeleton and subsequently edited, I will cd to the directory that contains the MyPackage directory and: R CMD INSTALL --

[R] Tcl/tk

2005-02-15 Thread Alexander . O . BETTINARDI
I am new to R so I apologize if this is a simple question I have created a tkgrid with a button that, once pushed, will extract the information from the grid, but I can not figure out the command to save the entries from my grid. Any thoughts? Thanks Alex Bettinardi Transportation Analyst TDD -

Re: [R] summary(aov(...)) into a string?

2005-02-15 Thread Gabor Grothendieck
RenE J.V. Bertin hotmail.com> writes: : : Hello, : : I'd like to annotate a plot with the output of summary(aov(model)), ideally just with the significant : effects. I don't find a means to redirect what that command prints into a string. : Is this possible, and if so, how? : ?capture.outpu

Re: [R] Trellis barchart, column display order issue

2005-02-15 Thread Deepayan Sarkar
On Tuesday 15 February 2005 14:57, Bock, Michael wrote: > I have searched the archives and instructions high and low but have > not found what I need. > > I have a dataframe named Data with columns: > AdjResND0 - the numeric result > Parameter - the classification (chemical name) > Sort - the order

Re: [R] summary(aov(...)) into a string?

2005-02-15 Thread Prof Brian Ripley
It doesn't print anything: the summary.aov (or summary.aovlist) print method does. ?summary.aov tells you the structure of the objects they return. On Tue, 15 Feb 2005, RenE J.V. Bertin wrote: I'd like to annotate a plot with the output of summary(aov(model)), ideally just with the significant e

Re: [R] Making a Package

2005-02-15 Thread Gabor Grothendieck
KKThird Yahoo.Com yahoo.com> writes: : : Hello. : I have what I know to be a simple question, but never having done anything like this it is : pretty tough. : : I'm trying to write an R package. I have a collection of functions that I loaded into R and : then used package.skeleton(). After e

[R] Making a Package

2005-02-15 Thread [EMAIL PROTECTED]
Hello. I have what I know to be a simple question, but never having done anything like this it is pretty tough. I'm trying to write an R package. I have a collection of functions that I loaded into R and then used package.skeleton(). After editing everything in the resulting folder, call it Ne

[R] Trellis barchart, column display order issue

2005-02-15 Thread Bock, Michael
I have searched the archives and instructions high and low but have not found what I need. I have a dataframe named Data with columns: AdjResND0 - the numeric result Parameter - the classification (chemical name) Sort - the order I want the chemical names to appear, with leading zeros so all are t

Re: [R] Problem with "R CMD Rd2dvi": Rd.sty not found

2005-02-15 Thread Duncan Murdoch
On Tue, 15 Feb 2005 19:44:51 + (UTC), Gabor Grothendieck <[EMAIL PROTECTED]> wrote : >I tried playing around with MiKTeX and got some advice on the TeX >list and now am using this config. This example assumes you have >put your localtexmf file in \localtexmf and are running >R 2.0.1 from \P

Re: [R] Problem with "R CMD Rd2dvi": Rd.sty not found

2005-02-15 Thread Gabor Grothendieck
Gabor Grothendieck myway.com> writes: > I tried playing around with MiKTeX and got some advice on the TeX > list and now am using this config. This example assumes you have > put your localtexmf file in \localtexmf and are running > R 2.0.1 from \Program Files\R\rw2001; otherwise, make the obv

Re: [R] Problem with "R CMD Rd2dvi": Rd.sty not found

2005-02-15 Thread Gabor Grothendieck
Duncan Murdoch stats.uwo.ca> writes: : : On Thu, 10 Feb 2005 21:27:03 + (UTC), Gabor Grothendieck : myway.com> wrote : : : >Tuszynski, Jaroslaw W. saic.com> writes: : > : >: : >: Hi, : >: : >: I run into a problem with "R CMD Rd2dvi" command: it gives me "File `Rd.sty' : >: not found"

Re: [R] how many 7th of the month is there between two dates

2005-02-15 Thread Prof Brian Ripley
It is really easy using R's date and datetime classes: st <- as.Date("1998-12-17") en <- as.Date("2000-1-7") st0 <- st # Get the 7th of the month before st0 <- as.POSIXlt(st); st0$mday <- 7; st0$mon <- st0$mon - 1 st0 <- as.Date(st0) ll <- seq.Date(st0, en, by="month") sum(ll > as.Date(st) & ll < e

Re: Off topic -- large data sets. Was RE: [R] 64 Bit R Background Question

2005-02-15 Thread Prof Brian Ripley
On Tue, 15 Feb 2005, Graham Jones wrote: In message <[EMAIL PROTECTED]>, r-help- [EMAIL PROTECTED] writes [Actually quoting Bert Gunter, BTW] Can comeone give me an example (perhaps in a private response, since I'm off topic here) where one actually needs all cases in a large data set ("large" bein

[R] nlsList vs nlme with lsoda function

2005-02-15 Thread Manuel Morales
Hello list members, I'm trying to get nlme to work with lsoda from odesolve. Currently, I can use nlsList to fit a simple logistic growth model to some simulated data, but I get the following error message with nlme: "Error in model.frame(formula, rownames, variables, varnames, extras, extraname

Re: [R] how many 7th of the month is there between two dates

2005-02-15 Thread Gabor Grothendieck
Omar Lakkis gmail.com> writes: : I want to caculate how many an exact day of the month there is between : two dates. : : For example; How many 7th of the month is there between "1998/12/17" : and "2000/1/7". To make the problem simple, the day of the month (7) : is the day in the 2nd date. d1

Off topic -- large data sets. Was RE: [R] 64 Bit R Background Question

2005-02-15 Thread Graham Jones
In message <[EMAIL PROTECTED]>, r-help- [EMAIL PROTECTED] writes >Can comeone give me an example (perhaps in a private response, since I'm off >topic here) where one actually needs all cases in a large data set ("large" >being > 1e6, say) to do a STATISTICAL analysis? By "statistical" I exclude, >

Re: [R] shrinkage estimates in lme

2005-02-15 Thread Dimitris Rizopoulos
What exactly do you mean by slope estimates? Marginal or subject-specific? Presuming that you mean the fixed-effects (since you refer to OLS) then these are unbiased since they are Weighted Least Squares estimators based on the marginal model. They are unbiased even in the case where you missp

[R] reading a 40kb csv file in R

2005-02-15 Thread Fairouz Makhlouf
I am trying to read a matrix of 5000*2000 from a csv file of real numbers in R using the scan function. I am getting an error saying that R can not read a vector > 32KB and it can't allocate a memory of 191 MB etc. What should I do. I am using R1.7.1. Thanks = Thanks Fairouz Makhlouf

[R] how many 7th of the month is there between two dates

2005-02-15 Thread Omar Lakkis
This is a eaeir way to ask my prior question: I want to caculate how many an exact day of the month there is between two dates. For example; How many 7th of the month is there between "1998/12/17" and "2000/1/7". To make the problem simple, the day of the month (7) is the day in the 2nd date. __

RE: [R] memory problem with package mix

2005-02-15 Thread Ted Harding
On 15-Feb-05 [EMAIL PROTECTED] wrote: > Hello, > I think we have a memory problem with em.mix. > > We have done: > >>library(mix) >>Manq <- read.table("C:/.../file.txt") >>attach(Manq) >>Manq >>V1 V2 V3 V4 .V27 >> 1 1 1 1 1... >> 2 1 NA 3 6 >> 3 1 2 6 2 >> ...

Re: [R] lattice multiple plots per page

2005-02-15 Thread Deepayan Sarkar
On Monday 14 February 2005 18:12, Jeff Jorgensen wrote: > Dear R-sters, > > I was wondering if anyone has encountered the following issues. I've > figured out how to get multiple levelplots [library(lattice)] on a > single plot. However, when I add text (adding axis labels for the > entire four p

[R] shrinkage estimates in lme

2005-02-15 Thread Bill Shipley
Hello. Slope estimates in lme are shrinkage estimates which pull the OLS slope estimates towards the population estimates, the degree of which depends on the group sample size and the distance between the group-based estimate and the overall population estimate. Although these shrinkage estimates

Re: [R] Correct effect plots from lme() objects

2005-02-15 Thread Deepayan Sarkar
On Tuesday 15 February 2005 08:15, CG Pettersson wrote: > Hello all! > > R2.0.1, W2k > > I posted this question to the list last Sunday without getting any > replies on the list. I got two off the list though, suggesting me to > plot "manually" as a second step, from estimable() or intervals() > ob

RE: [R] convolution of functions

2005-02-15 Thread Huntsinger, Reid
If you mean the convolution (f*g)(x) = integral f(x-y)g(y) dy for integrable functions f and g on R^n, then I think using the fact that the Fourier transform of the convolution is the product of the Fourier transforms of f and g is the most efficient approach, unless f or g have some special st

[R] D(eval(g)) problem, since "Function `eval' is not in the derivatives table"

2005-02-15 Thread Christoph Lehmann
thanks Andy and Dimitris for your reply to my expression/eval - problem starting with the resulting expression g I need g's derivative as expression, but I get: "Function `eval' is not in the derivatives table": #function (double gamama) h <- expression((t/d1)^a1*exp(-(t-d1)/b1)-c*(t/d2)^a2*exp(-

Re: [R] matlab norm(h) command in R: sqrt(sum(h^2)) - use in anexpression

2005-02-15 Thread Dimitris Rizopoulos
you are trying to square an expression not an number! try this: h.fun <- function(tt) ((tt/d1)^a1) * exp(-(tt-d1)/b1) - cc*((tt/d2)^a2) * exp(-(tt-d2)/b2) h <- h.fun(tt) h <- h/sqrt(sum(h*h)) I hope it helps. Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Pu

RE: [R] Tests on contingency tables

2005-02-15 Thread Jacques VESLOT
Thanks a lot for your help ! Right ! According to tables, most factors look indeed highly dependent... but, because of strange p-values and warning messages, as I tried to test it with Chisquare test, and because Fisher's Exact Test function doesn't work on my data, I wondered whether there were o

RE: [R] matlab norm(h) command in R: sqrt(sum(h^2)) - use in an e xpression

2005-02-15 Thread Liaw, Andy
> From: Christoph Lehmann > > Hi > > in matlab I defined a function (double gamma, parameters at > the end of > this mail) as >h(i)=((t/d1)^a1)*exp(-(t-d1)/b1)-c*((t/d2)^a2)*exp(-(t-d2)/b2); >h=h/norm(h); > > I do know that norm() in matlab is equal to: > >sqrt(sum(x^2)) > in R >

[R] matlab norm(h) command in R: sqrt(sum(h^2)) - use in an expression

2005-02-15 Thread Christoph Lehmann
Hi in matlab I defined a function (double gamma, parameters at the end of this mail) as h(i)=((t/d1)^a1)*exp(-(t-d1)/b1)-c*((t/d2)^a2)*exp(-(t-d2)/b2); h=h/norm(h); I do know that norm() in matlab is equal to: sqrt(sum(x^2)) in R so in R I do it like: #function (double gamama) h <- expressi

[R] memory problem with package mix

2005-02-15 Thread Delphine . Gille
Hello, I think we have a memory problem with em.mix. We have done: >library(mix) >Manq <- read.table("C:/.../file.txt") >attach(Manq) >Manq >V1 V2 V3 V4 .V27 > 1 1 1 1 1... > 2 1 NA 3 6 > 3 1 2 6 2 > ... > ... > 300 2 NA 6 2... > Essaimanq <-prel

RE: [R] R, MacOSX G5 bi-proc

2005-02-15 Thread Liaw, Andy
You can use the snow package on CRAN. I don't know if RPVM/RMPI work on the Macs, but they are also on CRAN. Andy > From: Naji > > Hi all, > > > I've downloaded binary files from CRAN to install R and it's > working quite > fine..I'm wandering if I can take advantage of two processors > fro

Re: [R] Tests on contingency tables

2005-02-15 Thread Prof Brian Ripley
I forgot to mention a crucial statistical point: As you are doing this pairwise, remember Simpson's paradox. or if you don't know about it, Google it (it is not really due to Simpson, an instance of Stigler's Law of Eponymy). On Tue, 15 Feb 2005, Prof Brian Ripley wrote: You can test inde

[R] R, MacOSX G5 bi-proc

2005-02-15 Thread Naji
Hi all, I've downloaded binary files from CRAN to install R and it's working quite fine..I'm wandering if I can take advantage of two processors from the Mac. Where can I find lecture/help for this? Thanks in advance Naji __ R-help@stat.math.ethz.ch

[R] Correct effect plots from lme() objects

2005-02-15 Thread CG Pettersson
Hello all! R2.0.1, W2k I posted this question to the list last Sunday without getting any replies on the list. I got two off the list though, suggesting me to plot "manually" as a second step, from estimable() or intervals() objects respectively. As this was not really what I wished for, I take

Re: [R] Tests on contingency tables

2005-02-15 Thread Prof Brian Ripley
You can test independence via a log-linear model. More importantly, you can model that dependence and learn something useful about the data. I don't see your point here: the two factors are clearly highly dependent: who cares what the exact p value is? Did you do e.g. a mosaicplot as I suspe

RE: [R] using poly in a linear regression in the presence of NA f ails (despite subsetting them out)

2005-02-15 Thread Liaw, Andy
My apologies: It's another case of me not thinking statistically... It may also help those of us whose brains run at slow clock speeds to have ?poly, ?bs and ?ns mention how they react to NAs. Best, Andy > From: Prof Brian Ripley > > Andy, > > I don't think it is a bug. The problem is that

Re: [R] update.packages - delete downloaded - N

2005-02-15 Thread Heinz Tuechler
Thank you for the hint regarding the destdir parameter of update.packages. Heinz Tüchler At 12:15 15.02.2005 +0100, Uwe Ligges wrote: >Heinz Tuechler wrote: > >> Hello, >> >> Using Update Packages from CRAN in R Version 2.0.1 Patched (2005-01-15) >> under Windows 98, I found it tricky to save th

Re: [R] R + MacOSX + Emacs(XEmacs) + ESS

2005-02-15 Thread Jari Oksanen
On Tue, 2005-02-15 at 10:34 -0200, Ronaldo Reis-Jr. wrote: > Hi, > > I try to use Emacs or XEmacs with R in a MacOS X Panter without X11. > > Anybody can make this work? Did you try googling for "macos X emacs"? That's the way you get it. I have found two different versions, both work graphicall

[R] R + MacOSX + Emacs(XEmacs) + ESS

2005-02-15 Thread Ronaldo Reis-Jr.
Hi, I try to use Emacs or XEmacs with R in a MacOS X Panter without X11. Anybody can make this work? How? Thanks Ronaldo -- Um mapa-mundi que nao inclua a utopia, nao e digno de consulta, pois deixa de fora as terras a que a Humanidade esta sempre apontando. -- Oscar Wilde -- |> // |

Re: [R] Problem with "R CMD Rd2dvi": Rd.sty not found

2005-02-15 Thread Duncan Murdoch
On Thu, 10 Feb 2005 21:27:03 + (UTC), Gabor Grothendieck <[EMAIL PROTECTED]> wrote : >Tuszynski, Jaroslaw W. saic.com> writes: > >: >: Hi, >: >: I run into a problem with "R CMD Rd2dvi" command: it gives me "File `Rd.sty' >: not found" error (See the output message on the bottom). I get the

[R] Tests on contingency tables

2005-02-15 Thread Jacques VESLOT
Dear all, I have a dataset with qualitative variables (factors) and I want to test the null hypothesis of independance between two variables for each pair by using appropriate tests on contingency tables. I first applied chisq.test and obtained dependance in almost all cases with extremely small

Re: [R] update.packages - delete downloaded - N

2005-02-15 Thread Uwe Ligges
Heinz Tuechler wrote: Hello, Using Update Packages from CRAN in R Version 2.0.1 Patched (2005-01-15) under Windows 98, I found it tricky to save the downloaded files. Even if I answer "N" to the question, if the downloaded files should be deleted, they are deleted after R is quitted. I understood t

[R] update.packages - delete downloaded - N

2005-02-15 Thread Heinz Tuechler
Hello, Using Update Packages from CRAN in R Version 2.0.1 Patched (2005-01-15) under Windows 98, I found it tricky to save the downloaded files. Even if I answer "N" to the question, if the downloaded files should be deleted, they are deleted after R is quitted. I understood that I have to copy th

[R] convolution of functions

2005-02-15 Thread Mari Luz Gamiz Perez
Dear sir, we would like to know if there exist any R package containing the computational performance of the n-fold Stieljes' convolution of functions. We look forward to hearing from you. Thank you in advance. M.Luz Gámiz Pérez Dpt. Estadística e

RE: [R] special symobol / character

2005-02-15 Thread Prof Brian Ripley
On Tue, 15 Feb 2005, FISCHER, Matthew wrote: Thanks for your speedy reply, I should have noted that I'm using a Linux machine. However, when I copy the symbol from Windows to Linux (using R/emacs) via an x-win 32 window it replaces the per mille symbol with a /211. \211 I hope. (That is octal

Re :Off topic -- large data sets. Was RE: [R] 64 Bit R Background Question

2005-02-15 Thread Naji
Hi, Also I agree those cases are relatively rare in STATISTICAL analysis, you can encounter them for simulation topics (natural catalysm a 5 meter in the topographics can change all the simulations) Two ideas (in addition to loading several sections) is 1- to search for duplicate cases and estim

RE: [R] special symobol / character

2005-02-15 Thread FISCHER, Matthew
Hi Eric, Thanks for your speedy reply, I should have noted that I'm using a Linux machine. However, when I copy the symbol from Windows to Linux (using R/emacs) via an x-win 32 window it replaces the per mille symbol with a /211. R then produces the character (not a per mille symbol!) that