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!)

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

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

[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

[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

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

[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] 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. JAROSLAW.W.TUSZYNSKI at 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).

[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] 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 graphically

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 the

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] 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

[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
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

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 from the

[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

[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 -

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 so in R I do

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

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

[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 -

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

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() objects

[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

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 panel

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 ... ... 300 2 NA 6

[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.

[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

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

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, say

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

2005-02-15 Thread Gabor Grothendieck
Omar Lakkis uofiowa at 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

[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,

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 being

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 en) and use = etc

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

2005-02-15 Thread Gabor Grothendieck
Duncan Murdoch murdoch at stats.uwo.ca writes: : : On Thu, 10 Feb 2005 21:27:03 + (UTC), Gabor Grothendieck : ggrothendieck at myway.com wrote : : : Tuszynski, Jaroslaw W. JAROSLAW.W.TUSZYNSKI at saic.com writes: : : : : : Hi, : : : : I run into a problem with R CMD Rd2dvi command: it

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

2005-02-15 Thread Gabor Grothendieck
Gabor Grothendieck ggrothendieck at 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,

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

[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

[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

Re: [R] Making a Package

2005-02-15 Thread Gabor Grothendieck
KKThird at Yahoo.Com kkthird at 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

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

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 I

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

2005-02-15 Thread Gabor Grothendieck
RenE J.V. Bertin rjvbertin at 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? :

[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] 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

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 Tue, 15

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

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

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] Could anyone answer for the following question

2005-02-15 Thread Gabor Grothendieck
Cuichang Zhao cuiczhao at 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

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 the

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