RE: [R] matrix inverse in C

Hello! For calculating the Choleski factor, you should look at: LINPACK: dpofa, dchdc, (dppfa, although i think this one is not included with R) The R code of function chol() is a good starting point. It is located at %RHOME%/library/base/R/base (line 4115 on the windows version of R 1.8.1

Re: [R] Phase Plane

Erin Hodgess wrote: Dear R People: Is there a function available to for phase plane plotting, please? Thanks, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: [EMAIL PROTECTED] PS R 1.8.1 on Windows XP 2000

[R] Proportions again

Hello. I asked before and it was great, cause as a beginner I learned a lot. But, if I have this in R (1 and 2 are codes for sex): sex-c(1,2,2,1,1,2,2,2) sex [1] 1 2 2 1 1 2 2 2 I´d like to obtain the proportion according to sex.So I type: prop.table(sex) [1] 0.07692308 0.15384615

Re: [R] Proportions again

table(sex)/length(sex) Andrew On Sunday 29 February 2004 06:35, Carlos Mauricio Cardeal Mendes wrote: Hello. I asked before and it was great, cause as a beginner I learned a lot. But, if I have this in R (1 and 2 are codes for sex): sex-c(1,2,2,1,1,2,2,2) sex [1] 1 2 2 1 1 2 2 2 I´d

[R] graphics device problems

I'm using R 1.8.0 under windows XP. I can't get certain of the graphics devices set up. For example, when I copy this line directly from the postscript help screen, I get the error messages that follow it: postscript(foo.ps) Error in PS(file, old$paper, old$family, old$encoding, old$bg,

Re: [R] Proportions again

Andrew Robinson [EMAIL PROTECTED] writes: table(sex)/length(sex) or, as is the intended usage: prop.table(table(sex)) sex 1 2 0.375 0.625 I asked before and it was great, cause as a beginner I learned a lot. But, if I have this in R (1 and 2 are codes for sex):

Re: [R] graphics device problems

David Parkhurst wrote: I'm using R 1.8.0 under windows XP. I can't get certain of the graphics devices set up. For example, when I copy this line directly from the postscript help screen, I get the error messages that follow it: postscript(foo.ps) Error in PS(file, old$paper, old$family,

Re: [R] Proportions again

Carlos Mauricio Cardeal Mendes wrote: I asked before and it was great, cause as a beginner I learned a lot. But, if I have this in R (1 and 2 are codes for sex): sex-c(1,2,2,1,1,2,2,2) sex [1] 1 2 2 1 1 2 2 2 I´d like to obtain the proportion according to sex.So I type: prop.table(sex) [1]

RE: [R] Proportions again

Dear Carlos, prop.table() takes a table as its argument, so you could specify prop.table(table(sex)). See ?prop.table for more details. John -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Carlos Mauricio Cardeal Mendes Sent: Sunday, February 29,

[R] stripchart and axes

Hi, I'd like to remove the axes from a plot produced by stripchart(). However, when trying stripchart(..., axes = FALSE), I get the error meassage Error in stripchart(hypokvot1 ~ treatment, jitter, pch = 1, vert = TRUE, : unused argument(s) (axes ...) using R 1.8.1 on Windows. Can it

RE: [R] Proportions again

Several people have alrady answered you by this time and in addition to their answers you might also be interested in CrossTable in package gregmisc. --- Date: Sun, 29 Feb 2004 11:35:06 -0300 From: Carlos Mauricio Cardeal Mendes [EMAIL PROTECTED] To: [EMAIL PROTECTED] Subject: [R]

[R] RMySQL Not Loading

Hi, I am having a problem getting RMySQL to run under FreeBSD 4.9. I am using R 1.8.1 with latest patches, MySQL 4.0.17 and RMySQL_0.03.tar.gz. Everything does appear to compile properly. However when I access the library I get the following error: library(RMySQL) Error in dyn.load(x,

Re: [R] stripchart and axes

On Sun, 2004-02-29 at 12:32, Henric Nilsson wrote: Hi, I'd like to remove the axes from a plot produced by stripchart(). However, when trying stripchart(..., axes = FALSE), I get the error meassage Error in stripchart(hypokvot1 ~ treatment, jitter, pch = 1, vert = TRUE, : unused

RE: [R] Proportions again

On Sun, 2004-02-29 at 12:40, Gabor Grothendieck wrote: Several people have alrady answered you by this time and in addition to their answers you might also be interested in CrossTable in package gregmisc. Gabor, Thanks for pointing out CrossTable(). Just as a quick heads up/clarification for

Re: [R] Rcmd SHLIB

Thanks tons. If you get a moment and feel like answering another question: Regarding the readme.packages file excerpt below: How can I set the R_HOME path? Is it always the R install directory? Also, is creating libR.a and libRblas.a a neccessary step in creating ANY dll for use with R?

RE: [R] Proportions again

That's true; however, CrossTable(x,x) does give the desired counts and proportions in the margin line at the bottom. See the row labelled Column Total in the following example based on Carlos' vector: sex-c(1,2,2,1,1,2,2,2) CrossTable(sex,sex) Cell Contents |-| |

Re: [R] Rcmd SHLIB

On Sun, 29 Feb 2004 15:43:25 -0800, you wrote: Thanks tons. If you get a moment and feel like answering another question: Regarding the readme.packages file excerpt below: How can I set the R_HOME path? Is it always the R install directory? Also, is creating libR.a and libRblas.a a neccessary

RE: [R] Proportions again

On Sun, 2004-02-29 at 18:27, Gabor Grothendieck wrote: That's true; however, CrossTable(x,x) does give the desired counts and proportions in the margin line at the bottom. See the row labelled Column Total in the following example based on Carlos' vector: sex-c(1,2,2,1,1,2,2,2)

[R] glm logistic model, prediction intervals on impact af age 60 compared to age 30

Dear R-list. I have done a logistic glm using Age as explanatory variable for some allergic event. #the model model2d-glm(formula=AEorSAEInfecBac~Age,family=binomial(logit),data=emrisk) #predictions for age 30 and 60 preds-predict(model2d,data.frame(Age=c(30,60)),se.fit=TRUE) # prediction

RE: [R] Proportions again

I agree its overkill but prop.table generalizes to mulitple dimensions, in which case it becomes comparable to CrossTable, so I thought it was worth mentioning. By the way, as the author of CrossTable, perhaps you might might consider generalizing CrossTable to the 1D case too? --- Date:

[R] se.contrast ....too hard??? .... Too easy????? .....too trivial???? ...... Too boring.....too????????

Hi all, Regular and avid readers of this column will know that Don Driscoll and I have recently posted two messages requesting assistance concerning an apparent failure of se.contrast to produce an se for a contrast. So far, an ominous silence rings in our ears, but read on Gentle Reader, and see

RE: [R] Proportions again

On Sun, 2004-02-29 at 19:26, Gabor Grothendieck wrote: I agree its overkill but prop.table generalizes to mulitple dimensions, in which case it becomes comparable to CrossTable, so I thought it was worth mentioning. By the way, as the author of CrossTable, perhaps you might might consider

[R] Re: R-help Digest, Vol 12, Issue 27

Dear R helpers, I am interested in doing an overlapping moving block bootstrap (Davinson and Hinkley(1997)) on financial market data using R. I have daily stock returns and have tried ( unsuccessfully using sample and boot function) to create a sample containing consecutive trading days For

[R] Confused in simplest-possible function

I wrote the following code: --- oneindex - function(x) { summary(x) } A - read.table(try.data, col.names=c(date, lNifty)) summary(A) oneindex(A$lNifty)

[R] question about mixed effects model

Hello. I have some trouble with mixed effects in R, similar to problems that other people had with not nested models and lme, as I understand from the mailing list archive. Unfortunately, I could not understand the solutions that were proposed... I have a data set with response variable (y) and