[R] Multithreading

Dear list, is R capable of multithreading? At the moment it uses only 1 CPU out of 4 CPUs on my system. Is there a compilation option or start option so that R uses multiple CPUs at the same time to accomplish a calculation? Some information on my system: R.version _ platform

[R] Question about R

Hi, I would like to ask how can I do the problem as follows by R. Thank you very much. Q: Create a list of twenty distinct first names, ten male and ten female, and store them into an array. Write a routine for selecting a simple random sample of five names and counts the number of

[R] Deming regression

Dear R helpers, I would like to know whether there is any package in R for Deming regression? Many thanks Nguyen __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] problem for use R

hello, we are french student, we have problem for export the test on Excel. we didn't succed for save and name the test so we can't export it... can you help us please thanks _ mobile comme sur PC !

Re: [R] Multithreading

1) Please update your R: 2.5.0 is imminent. 2) Please read the 'R Installation and Adminstration manual'. R itself in not multithreaded but can make use of multithreaded BLAS. 3) A search of the list archives would have produced many hits, as does RSiteSearch(multithreaded), with more

[R] Odp: Question about R

Hi isnt it a homework? FYI see ?list ?vector ?data.frame ?sample ?lapply and some basic stuff like Introduction manual provided with your R instalation to see how some simple data manipulations can be done. Regards Petr Pikal [EMAIL PROTECTED] [EMAIL PROTECTED] napsal dne 05.04.2007

[R] Odp: problem for use R

Hi It is difficult from your question to derive what you exactly want. For reading Excel files there are several ways - see archives. I personally select some data in Excel press Ctrl-C and then in R read.delim(clipboard), however sometimes can be prefered a way through saving txt or csv

Re: [R] Generate a serie of new vars that correlate with existingvar

Hello, list why not add the smart proposal by Greg Snow as a built-in function in {stats}, just changing the x234 and newc lines to allow for more distributions to be generated ? Or do I miss an already existing function to do that ? Regards. Olivier # slight modification of the original code

[R] or, you could consider R-spatial.

http://r-spatial.sourceforge.net Linda Smith wrote: I have a netcdf gridded file with LCC projection. I can easily use image.plot to visualize it. However, as the axises are in X,Y, not Lat and Lon, I could not add state or country maps onto it (or lat lon information). I do have a

[R] reading of a matrix

Dear R-experts, I still have problems with the reading of a matrix. Input: matrixData6.txt A-Paar B-Paar C-Paar D-Paar E-Paar A 1 3 5 7 9 B 2 4 6 8 10 R-commands: y=read.table(file=Z:/Software/R-Programme/matrixData6.txt) y Result: A.Paar B.Paar C.Paar D.Paar E.Paar A 1 3

Re: [R] Problems loading package GeoXp

On Wed, 4 Apr 2007, Carlos Guerra wrote: Dear useR's, I wanted to load the package GeoXP in my computer but I couldn't because this warning message appeared: require(GeoXp) Carregando pacotes exigidos: GeoXp Carregando pacotes exigidos: tcltk Loading Tcl/Tk interface ... Erro em

Re: [R] gpclib library

On Thu, 5 Apr 2007, Leonardo Lami wrote: my OS is Linux Debian distribution, my R vers is 2.4.0 Patched (2006-11-25 r39997) I tried several repositories but the output is ever the same And you have *still* not had the courtesy to show us what it actually is (as asked for both by me and the

Re: [R] reading of a matrix

Thank you very much. It realizes the things I want. Thanks a lot, Corinna -Ursprüngliche Nachricht- Von: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Im Auftrag von Vladimir Eremeev Gesendet: Donnerstag, 5. April 2007 13:25 An: r-help@stat.math.ethz.ch Betreff: Re: [R] reading of a

Re: [R] reading of a matrix

Hi [EMAIL PROTECTED] napsal dne 05.04.2007 12:50:09: Dear R-experts, I still have problems with the reading of a matrix. Input: matrixData6.txt A-Paar B-Paar C-Paar D-Paar E-Paar A 1 3 5 7 9 B 2 4 6 8 10 R-commands:

Re: [R] gpclib library

On Thu, 5 Apr 2007, Leonardo Lami wrote: Hi my OS is Linux Debian distribution, my R vers is 2.4.0 Patched (2006-11-25 r39997) I tried several repositories but the output is ever the same Please do not write to the author of the original non-R-project C code - the package maintainer is

[R] clear console

Hi R Users, In Windows I can clear console using CTRL-L, but can I do this by certain command in my programs? e.g for (i in 1:10){ something clear console } Best, Rob [[alternative HTML version deleted]] __

Re: [R] reading of a matrix

Schmitt, Corinna wrote: I still have problems with the reading of a matrix. Input: matrixData6.txt A-Paar B-Paar C-Paar D-Paar E-Paar A 1 3 5 7 9 B 2 4 6 8 10 R-commands: y=read.table(file=Z:/Software/R-Programme/matrixData6.txt) y Result: A.Paar B.Paar C.Paar

Re: [R] clear console

Check out: http://tolstoy.newcastle.edu.au/R/help/06/02/21634.html On 4/5/07, Robert McFadden [EMAIL PROTECTED] wrote: Hi R Users, In Windows I can clear console using CTRL-L, but can I do this by certain command in my programs? e.g for (i in 1:10){ something clear console }

Re: [R] gpclib library

Hi my OS is Linux Debian distribution, my R vers is 2.4.0 Patched (2006-11-25 r39997) I tried several repositories but the output is ever the same Thank yor for your help All the best Leonardo Prof Brian Ripley ha scritto: It's a _package_ that you are looking for (there is something else

Re: [R] Question about R

?sample nams - c(rep(a,10), rep(b,10)) table(sample(nams, 5, replace = TRUE)) --- Chia-Yi Tseng [EMAIL PROTECTED] wrote: Hi, I would like to ask how can I do the problem as follows by R. Thank you very much. Q: Create a list of twenty distinct first names, ten male and ten

[R] R 'could not find any X11 fonts'

I'm having this problem - and I have seen earlier post on the same problem - but the solutions either do not work for me or (being fairly new to linux) I do not understand the suggested actions. So here goes: Whenever I try to add a 'main' title R responds with an error that it can't find any

Re: [R] reading of a matrix

On Apr 5, 2007, at 6:50 AM, Schmitt, Corinna wrote: Dear R-experts, I still have problems with the reading of a matrix. Input: matrixData6.txt A-Paar B-Paar C-Paar D-Paar E-Paar A 1 3 5 7 9 B 2 4 6 8 10 R-commands: y=read.table(file=Z:/Software/R-Programme/matrixData6.txt) y

[R] glmmPQL and wald test or score test

Dear All, I would like to know if I can use the Wald test or the Score test for testing a set of parameter in a glmm model estimated by the glmmPQL method. I am asking this question because I don't know whether the glmmPQL's estimates are maximum likelihood estimates or not, and also

[R] Cluster in large dataset with nominal variables.

Hi R users: Is it any R implementation of a cluster procedure for large data sets (clara()) but with dissimilary that can handle continuous, categorical and nominal variables (daisy()) like CLARANS (Clustering Large Applications based up on RANdomized Search), por example? Thank you for your

Re: [R] how to image.plot a XY grid file into a lat-lon map

If your country or state borders are polygons or polylines, you could convert them to desired projection using the function project from the package rgdal. Latitude-longitude grid also could be added by generating desired polylines in lat-lon and converting them to the desired projection using

[R] Likelihood returning inf values to optim(L-BFGS-B) other options?

Dear R-help list, I am working on an optimization with R by evaluating a likelihood function that contains lots of Gamma calculations (BGNBD: Hardie Fader Lee 2005 Management Science). Since I am forced to implement lower bounds for the four parameters included in the model, I chose the

Re: [R] clear console

This has been asked previously and I happened to have saved the final answer: cls - function() { require(rcom) wsh - comCreateObject(Wscript.Shell) comInvoke(wsh, SendKeys, \014) invisible(wsh) } -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL

Re: [R] clear console

Gabor thank you very much. It works wonderful. Rob -Original Message- From: Gabor Grothendieck [mailto:[EMAIL PROTECTED] Sent: Thursday, April 05, 2007 1:39 PM To: Robert McFadden Cc: r-help@stat.math.ethz.ch Subject: Re: [R] clear console Check out:

[R] read.spss (package foreign) and SPSS 15.0 files

Hello, does anyone have experience with reading SPSS Version 15.0 files into R (version 2.4.1, WinXP)? I have long been sucessfully reading SPSS files with read.spss from the wonderful foreign package, but somehow after upgrading from SPSS14 to SPSS15 I seem to have problems. Trying a simple

[R] R Graphics training course - Basel - 21.22.23 may 2007

Mango Solutions are pleased to announce the above 3 days R course in Basel as part of our schedule for Q2 2007. --- R Graphics -21th - 23rd May 2007 (Basel)

[R] 20% off selected R and S-plus books from Cambridge University Press

Get 20% off selected Cambridge books for the R-help community. Until the end of June, R-help users can get an exclusive discount on the following Cambridge University Press books that use R or S-plus! Simply visit http://www.cambridge.org/rbooks and click the appropriate link to get more

Re: [R] how to image.plot a XY grid file into a lat-lon map

On Thu, 5 Apr 2007, Vladimir Eremeev wrote: If your country or state borders are polygons or polylines, you could convert them to desired projection using the function project from the package rgdal. Latitude-longitude grid also could be added by generating desired polylines in lat-lon

[R] p value for coefficients in multinomial model

Dear all, 1)how can I easily get p value for the coefficients of factors in a multinomial model? 2)why the p values for type III test with Anova are not identical to that from SAS? for example: A,B and C are categorical variables,but the proportions of each level in each categorical

[R] Transform distance constraints into a 2D-map

Dear all, Could you please point me out to a function or set of functions that would allow me to optimize positions on the 2D map of points for which I are distance contrainted? for example, if I have: a b c a 0.0 0.5 0.7 b 0.5 0.0 0.3 c 0.7 0.3 0.0 I would like to obtain (x,y)

[R] package for Matlab

Hallo, does a package for Matlab exist in R? If yes, where can I find it and how can I install it under R? Thanks, Corinna __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] Likelihood returning inf values to optim(L-BFGS-B) other options?

Hi, In your code, the variables x (which I assume is the observed data), Tvec, and flag are not passed to the function as arguments. This could be a potential problem. Another problem could be that you have to use negative log-likelihood function as input to optim, since by default it minimizes

[R] permutations of a binary matrix with fixed margins

Dear R Users, How can I randomize a matrix (with contains only 0 and 1) with the constraint that margin totals (columns and rows sums) are kept constant? Some have called this swap permutation or something like that. The principle is to find bloc of 10 01 and change it for 01 10 there can be

Re: [R] package for Matlab

Hallo, here again my modified information: does a package for Matlab exist in R? I'm using Windows. I know that there exist one for UNIX. Does one exist for Windows? If yes, where can I find it and how can I install it under R? Thanks, Corinna __

Re: [R] package for Matlab

[EMAIL PROTECTED] napsal dne 05.04.2007 15:24:57: Hallo, does a package for Matlab exist in R? You probably could get quicker answer trying to use some search posibilities provided by CRAN. In first few hits from Rseek there is a package R.matlab. I hope it can be installed by a standard

Re: [R] Deming regression

Looking in the archives for Deming regression brings up a previous reply of mine: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/57384.html , which may shed some light. David L. Reiner Rho Trading Securities, LLC Chicago IL 60605 312-362-4963 -Original Message- From: [EMAIL PROTECTED]

[R] Extent of time zone vulerability for POSIX date and time classes

Hi. I frequently convert date and time data to and from character representations. I'm frustrated with chron, because 'seconds' are required to create a time object (my input data never has seconds). More importantly, I cannot make chron print the format 12/30/2006 (which my output data

Re: [R] Likelihood returning inf values to optim(L-BFGS-B) other options?

Hi, I got this message sometimes too. Try to change the initial value and lower bound value and see if it works. Jin [EMAIL PROTECTED] wrote: Dear R-help list, I am working on an optimization with R by evaluating a likelihood function that contains lots of Gamma calculations

[R] Plotting multiple curves with lattice graphs

Hi List, I would like to plot multiple curves (parametric density curves) in one plot. For example: # parameters for three normal density curves parms = data.frame(ID=c(1,2,3),mu=c(50,55,60),sigma=c(10,12,15)) # I can easily draw three normal density curves using curve():

[R] coxme residuals

Hi all, I am trying to use the coxme function to fit a random effect cox model with correlated frailty. The function is working well but I would like to know how to get predicted values and residuals. Thanks for your help, Youenn Drouet, Student in biostatistics from Lyon1 university.

[R] ISO Sims test

Does anyone know of a package that includes the Modified Sims test [Sims, 1972]? This test is used in econometrics and is a kind of alternative to the Granger test [Granger, 1969], which is in the package lmtest. Thanks in advance, chris Refernces: Granger, C.W.J. (1969), Investigating Causal

Re: [R] Likelihood returning inf values to optim(L-BFGS-B) other options?

On Thu, 5 Apr 2007, Ravi Varadhan wrote: In your code, the variables x (which I assume is the observed data), Tvec, and flag are not passed to the function as arguments. This could be a potential problem. I think scoping will probably find them. Another problem could be that you have to

Re: [R] package for Matlab

Hi Corinna! I don't if this is , but look that: http://cran-r.c3sl.ufpr.br/src/contrib/Descriptions/matlab.html http://cran-r.c3sl.ufpr.br/src/contrib/Descriptions/R.matlab.html -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22

Re: [R] package for Matlab

Petr PIKAL wrote: You probably could get quicker answer trying to use some search posibilities provided by CRAN. OTOH, if we try to google for R, we get 1.880.000.000 hits :-) Alberto Monteiro __ R-help@stat.math.ethz.ch mailing list

Re: [R] package for Matlab

Schmitt, Corinna wrote: Hallo, does a package for Matlab exist in R? To read and write MAT files, there is the R.matlab package: http://cran.r-project.org/src/contrib/Descriptions/R.matlab.html This package also enables bidirectional communication between R and Matlab. An alternative

Re: [R] changing existing graphs optimally

Using traditional and lattice(trellis) graphics, the usual approach is to have your plotting commands in a script file (under windows you can use the one with the GUI, in linux/unix you can have a text editor open with the script). Then either run the commands by copy/paste or using the source

[R] Densitymap

Hello, we have a question concerning densitymaps. Actually, we found in the R help for package Geo Xp that the densitymap is of the following form: densitymap(long, lat, var,kernel='triweight',listvar=NULL, listnomvar=NULL,carte=NULL, criteria=NULL,label=,

Re: [R] Plotting multiple curves with lattice graphs

Try this: x - 1:100 xyplot(dnorm(x, 50, 10) + dnorm(x, 55, 12) + dnorm(x, 60, 15) ~ x, type = l) By the way, you can save classic graphics, at least on windows, like this: ### your code parms - data.frame(ID=c(1,2,3),mu=c(50,55,60),sigma=c(10,12,15))

Re: [R] p value for coefficients in multinomial model

Dear Xingwang Ye, As it says on the help page for Anova(), Be very careful in formulating the model for type-III tests, or the hypotheses tested will not make sense. It looks as if you fit the model with the default contrast coding for the factors, contr.treatment. To get sensible Type-III sums

[R] Running MCMCpack

Hi there, I am running MCMCpack (MCMCirt1d model) on some files (26 items x about 800 – 1200 individuals). I have a problem. When I am working on “big” files, the R program crashes. More precisely I got the following Microsoft Warning: ***

Re: [R] Extent of time zone vulerability for POSIX date and time classes

If you use Sys.putenv(TZ = GMT) at the beginning of your session then local time zone and GMT time zone will be the same so you should not have a problem. This was not possible, at least on Windows, at the time the R News article was written. On 4/5/07, Tim Bergsma [EMAIL PROTECTED] wrote:

Re: [R] Generate a serie of new vars that correlate with existingvar

Oliver, I have thought of adding something like this to a package, but here is my current thinking on the issue. This question (or similar) has been asked a few times, so there is some demand for a general answer, I see three approaches: 1. Have an example of the necessary steps archived in a

[R] Reasons to Use R

Dear All, The institute I work for is organizing an internal workshop for High Performance Computing (HPC). I am planning to attend it and talk a bit about fluid dynamics, but there is also quite a lot of interest devoted to data post-processing and management of huge data sets. A lot of people

[R] omega squared calculation

Hi list, Is there a function in R to calculate the omega squared (an estimate of the dependent variance accounted for by the independent variable in the population for a fixed effects model), or does anyone have some code to let me calculate this fairly easy out of aov() or anova() output?

Re: [R] Transform distance constraints into a 2D-map

Look at the cmdscale function (there are variations in some packages, but that should get you started). -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL

Re: [R] package for Matlab

Dear Tobias, the first packet was succesfull installed. Please check the link for the last packet again. I get the following error message if I try to install it: CC - http://odin.mdacc.tmc.edu/~roebuck/R; install.packages(rwt, dependencies = TRUE, CRAN = MDACC) Error in

Re: [R] package for Matlab

Dear R-experts, Thanks for the help. I just istalled everything and will test it after Easter. Happy Easter, Corinna __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] package for Matlab

The first 2 links from doing: RSiteSearch('matlab') Look like they may answer your question. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL

Re: [R] Transform distance constraints into a 2D-map

On 05-Apr-07 13:05:29, Emmanuel Levy wrote: Dear all, Could you please point me out to a function or set of functions that would allow me to optimize positions on the 2D map of points for which I are distance contrainted? for example, if I have: a b c a 0.0 0.5 0.7 b 0.5 0.0

Re: [R] Running MCMCpack

On 4/5/2007 10:22 AM, Thomas Christin wrote: Hi there, I am running MCMCpack (MCMCirt1d model) on some files (26 items x about 800 – 1200 individuals). I have a problem. When I am working on “big” files, the R program crashes. More precisely I got the following Microsoft

[R] about systemfit

Hello. I am still a newbie in R. Excuse me if I am asking something obvious. My efforts to get an answer through browsing the mailing archives failed. I want to perform an augmented Dickey-Fuller test and to obtain AIC and BIC and to be able to impose some linear restrictions on the ADF

Re: [R] Reasons to Use R

Dear Mr. Isella, I just started my PhD Thesis. I need to work with R. Good sources are Bioconductor (www.bioconductor.org). It is a DB based on R-programming. Another institute which has good experiences with R is the HKI in Jena, Germany. Perhaps you can contact Mrs. Radke to get more

Re: [R] Transform distance constraints into a 2D-map

Dear All, Thank you for your replies. That will definitly get me started. Ted: I said optimize because I suppose that an exact solution does not necessarily exists sometimes. Best, Emmanuel On 4/5/07, Ted Harding [EMAIL PROTECTED] wrote: On 05-Apr-07 13:05:29, Emmanuel Levy wrote: Dear

Re: [R] Extent of time zone vulerability for POSIX date and time classes

Gabor, thanks for the feedback. I like the idea of coercing the system time zone, and sticking with strptime/strftime. The balance would probably tip in favor of chron, however, if I could get the format correct. Apparently the default is two-digit year: library(chron) chron(1) [1] 01/02/70

Re: [R] about systemfit

Dear Martin, I'll address only part of your question, which is how to get the code for linear.hypothesis() in the car package: linear.hypothesis() is an S3 generic function with several methods: library(car) methods(linear.hypothesis) [1] linear.hypothesis.default* linear.hypothesis.glm*

[R] Reasons to Use R

Lorenzo Isella writes: (4)finally, a list of the advantages for using R over commercial statistical packages. Here's my entry on the list, as this was a topic of conversation over lunch: it's better than the proprietary statistical software I use most of the time. By better I mean that the

[R] about systemfit

Thank you very much for your responsiveness. Here are the tests that show the same results, as they must: linear.hypothesis(model=adfResulm,hypothesis.matrix=Rrestr10,test=Chisq): Res.Df RSS Df Sum of Sq Chisq Pr(Chisq) 1127 7.3782 2137 7.6848 -10

Re: [R] Plotting multiple curves with lattice graphs

Thanks for two nice simple solutions. I did not know about recordPlot() or replayPlot(). Unfortunately these do not give me quite what I want because I need (for use in another function) an object of class lattice. I should have mentioned this part earlier! However, your first solution using

Re: [R] Extent of time zone vulerability for POSIX date and time classes

There are several ways: 1. the default number of year digits is set by the chron.year.abb option. The R News article mentions this in the Avoiding Errors section but does suggest you not set it. options(chron.year.abb = FALSE) chron(1) 2. Alternately, define: out.fmt - function(x) with(

Re: [R] package for Matlab

On 4/5/07, Tobias Verbeke [EMAIL PROTECTED] wrote: Schmitt, Corinna wrote: Hallo, does a package for Matlab exist in R? To read and write MAT files, there is the R.matlab package: http://cran.r-project.org/src/contrib/Descriptions/R.matlab.html This package also enables bidirectional

Re: [R] Extent of time zone vulerability for POSIX date and time classes

Gabor, The out.fmt approach seems optimal. Thank you! Tim. Gabor Grothendieck wrote: There are several ways: 1. the default number of year digits is set by the chron.year.abb option. The R News article mentions this in the Avoiding Errors section but does suggest you not set it.

Re: [R] changing existing graphs optimally

On 4/5/07, Greg Snow [EMAIL PROTECTED] wrote: Using traditional and lattice(trellis) graphics, the usual approach is to have your plotting commands in a script file (under windows you can use the one with the GUI, in linux/unix you can have a text editor open with the script). Then either run

Re: [R] about systemfit

Dear Martin, Note that the F-test provided by linear.hypothesis() in the car package appears to agree with the Wald test provided by both linear.hypothesis() and waldtest.systemfit(), while ftest.systemfit() gives a much smaller p-value and an F that's even larger than the Wald chisquare. This

[R] Indexing in anova summary output of the form: summary(aov(y ~ x1, Error = (x1/x2)))

Hi, there I'm trying to get the value of the Mean Square from the ANOVA model summary that comes from specifying the error term, and am wondering if one can actually do this ( I know it's possible when using anova(lm) objects and the like, but I'm having a tough time with it under this

Re: [R] Indexing in anova summary output of the form: summary(aov(y ~ x1, Error = (x1/x2)))

The output is of class summary.anovalist so (hint, hint) it is a list. Try [[ ]] to get list elements. Without a reproducible example (see the footer and the posting guide) it is not possible to show you working code. On Thu, 5 Apr 2007, Michael D. Rennie wrote: Hi, there I'm trying to

[R] creating a data frame from a list

Dear all, A few months ago, I asked for your help on the following problem: I have a list with three (named) numeric vectors: lst = list(a=c(A=1,B=8) , b=c(A=2,B=3,C=0), c=c(B=2,D=0) ) lst $a A B 1 8 $b A B C 2 3 0 $c B D 2 0 Now, I'd love to use this list to create the following data

Re: [R] Indexing in anova summary output of the form: summary(aov(y ~ x1, Error = (x1/x2)))

Hi, there Sorry- I excluded some example code in the haste of my original e-mail. Please find it included below. Note, however, that the elements in the list in the summary output appear to be the anova tables themselves- what I need to do is get inside the anova tables (execute the code

[R] MANOVA with repeated measurements

Hello, I have a question regarding performing manova. I have an experiment where I want to measure 10 output variables for 3 different measurement methods. Since each of the methods requires some user interaction, I would also like to include repeated measures for each of the output

[R] binom.test() query

Hi Folks, The recent correspondence about strange fisher.test result, and especially Peter Dalgaard's reply on Tue 03 April 2007 (which I want to investigate further) led me to take a close look at the code for binom.test(). I now have a query! The code for the two-sided case computes the

Re: [R] binom.test() query

On 4/5/2007 5:35 PM, (Ted Harding) wrote: Hi Folks, The recent correspondence about strange fisher.test result, and especially Peter Dalgaard's reply on Tue 03 April 2007 (which I want to investigate further) led me to take a close look at the code for binom.test(). I now have a query!

[R] Logistic/Cox regression: Parameter estimates directly from model matrix

Hi out there Is there a way to get the estimated coefficients in a logistic / Cox regression without having to specify a 'formula' but by only giving the model matrix? Example for Cox regression: ## predictors n - 50 q1 - rnorm(n) q2 - rgamma(n, 2, 2) Z - cbind(q1, q2) ## response ttf -

Re: [R] creating a data frame from a list

Dmitri: As you apparently have not received a reply IMHO, one of the glories of R is the ease with which you can create de novo solutions for little problems like this yourself. While there may be more efficient,robust, and elegant solutions already available, it can frequently be

Re: [R] Reasons to Use R

--- Lorenzo Isella [EMAIL PROTECTED] wrote: (4)finally, a list of the advantages for using R over commercial statistical packages. The money-saving in itself is not a reason good enough and some people are scared by the lack of professional support, though this mailing list is simply

[R] lm() intercept at the end, rather than at the beginning

Hi, I wonder if someone has already figured out a way of making summary(mylm) # where mylm is an object of the class lm() to print the (Intercept) at the last line, rather than the first line of the output. I don't know about, say, biostatistics, but in economics the intercept is usually the

Re: [R] read.spss (package foreign) and SPSS 15.0 files

--- RINNER Heinrich [EMAIL PROTECTED] wrote: Hello, does anyone have experience with reading SPSS Version 15.0 files into R (version 2.4.1, WinXP)? I have long been sucessfully reading SPSS files with read.spss from the wonderful foreign package, but somehow after upgrading from SPSS14

Re: [R] changing existing graphs optimally

--- hadley wickham [EMAIL PROTECTED] wrote: On 4/5/07, Greg Snow [EMAIL PROTECTED] wrote: Using traditional and lattice(trellis) graphics, the usual approach is to have your plotting commands in a script file (under windows you can use the one with the GUI, in linux/unix you can have

Re: [R] lm() intercept at the end, rather than at the beginning

Try this. It captures the output of x and sets idx to the line numbers of the coefficients, rearranging their order in the next line and printing them out in the line after that. my.print.summary.lm - function(x, ...) { out - capture.output(x) idx - seq(grep(Intercept, out),

Re: [R] MANOVA with repeated measurements

Dear Mark, Take a look at the Anova() function in the car package. There's an example of a repeated-measures MANOVA in ?Anova (the last example). Also see ?anova.mlm. I hope this helps, John John Fox Department of Sociology McMaster University Hamilton, Ontario

Re: [R] read.spss (package foreign) and SPSS 15.0 files

Not being the developer I cannot answer definitively but, as a frequent user of SPSS files I can give you my experience. 1) The unrecognized coding is perhaps due to the locale of the SPSS installation. I have had success reading in files from version 15 but often encounter that error when the

[R] translating sas proc mixed to lme()

Hi All I am trying to translate a proc mixed into a lme() syntax. It seems that I was able to do it for part of the model, but a few things are still different. It is a 2-level bivariate model (some call it a pseudo-3-level model). PROC MIXED DATA=psdata.bivar COVTEST METHOD = ml; CLASS