Re: [R] RWinEdt installation problems with Vista
Hi, I was trying to install the packageRWinEdt in my computer with Vista OS. O Use Tinn-R it works with Vista although has some minor issues (which are probably system specific). https://sourceforge.net/projects/tinn-r Stefan -=-=- ... Time is an illusion, lunchtime doubly so. (Ford Prefect) __ 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, reproducible code.
Re: [R] kendall test
On Thursday 06 September 2007 09:48:22 elyakhlifi mustapha wrote: em I thougth that there is a function which does the kendall test in R, em I writed on the console apropos(kendall) and I didn't found anything em can you tell me how could I do to use the kendall test? ?cor.test btw.: rseek.org is a very good help for such questions Stefan __ 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, reproducible code.
Re: [R] R-2.5.1 RedHat EL5 compilation failed
Original Message Subject: [R] R-2.5.1 RedHat EL5 compilation failed From: Wang Chengbin [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: 26.08.2007 15:22 I can't get R-2.5.1 compiled under RedHat EL5 with gcc 4.1.1. Configure failed at the following: You don't need to compile, you could also use the Fedora Core 6 Extras repository package(s) of R (current: is R-2.5.1-2.fc6.i386.rpm) to install the necessary rpm packages from there. (Best is to use the smart package manager, there you can easily activate channels which are repositories.) As far as I understood FC6 is the base of RHEL 5. Stefan -=-=- ... Time is an illusion, lunchtime doubly so. (Ford Prefect) __ 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, reproducible code.
Re: [R] combine 2 data frames with missing values
I have 2 data frames as follows: abc 1 NA 1 2 NA 2 NA 3 3 So a, b are the input values and c is the output which I am interested in. NA - Missing values. I used rbind, but its not working. Let me know if anyone can help me What exactly is your problem? To create such a data frame or what? What did you do with rbind? Please be more explicit if you ask questions, sometimes its hard to guess what the people want... Stefan __ 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, reproducible code.
Re: [R] Newbie
I'm a bit new to the world of R so forgive my ignorance. That has nothing to do with the knowledge of R but with the model. Age has only 2 different values: 0 and 1 and if it is 0 there is no scars, so what exactly have you expected from the model? I would say if you just want to prove that older deer have more scars try the Mann Whitney non parametric test... Stefan I'm trying to do a zero-inflated negative binomial regression and have received an error message and i'm not sure what it means. I'm running R 2.5.1 on XP. I have just tried a really simple version of the model to see if it would run before I put all the variables in. I have attached all the variables to the object alan. Here is the table: Date Location Deer Code Sex Age Pel Per Lac Scars Emb Lar Nym Bee Mi Mass 1 12/06/2007 Ballysallagh0 B1 FALSE 1 A 1 1 4 0 0 0 0 0 21.59 2 15/06/2007 Ballysallagh0 B10 FALSE 1 A 1 1 5 0 0 0 0 2 26.19 3 15/06/2007 Ballysallagh0 B12 FALSE 1 A 1 1 0 5 0 0 0 1 28.55 4 15/06/2007 Ballysallagh0 B13 FALSE 1 A 0 1 5 0 0 0 4 45 23.93 5 15/06/2007 Ballysallagh0 B16 FALSE 1 A 1 1 0 4 0 0 0 1 34.19 6 15/06/2007 Ballysallagh0 B17 FALSE 1 A 1 1 5 0 0 0 0 0 25.02 7 15/06/2007 Ballysallagh0 B18 FALSE 1 A 1 1 5 0 0 0 0 7 33.06 8 12/06/2007 Ballysallagh0 B5 FALSE 1 A 1 1 4 0 0 0 0 2 23.55 9 12/06/2007 Ballysallagh0 B6 FALSE 1 A 1 1 4 0 0 0 0 2 22.67 10 12/06/2007 Ballysallagh0 B7 FALSE 1 A 1 1 5 0 0 0 0 2 22.57 11 12/06/2007 Ballysallagh0 B8 FALSE 1 A 1 1 5 0 0 0 0 10 24.01 12 24/06/2007 Caledon1 C12 FALSE 1 A 1 0 0 0 1 0 0 3 21.68 13 25/06/2007 Caledon1 C17 FALSE 1 A 1 0 4 0 1 0 0 23 20.56 14 25/06/2007 Caledon1 C18 FALSE 1 A 1 0 0 0 0 0 2 2 17.47 15 25/06/2007 Caledon1 C20 FALSE 1 A 1 0 8 0 0 0 0 3 19.97 16 25/06/2007 Caledon1 C21 FALSE 1 A 1 0 3 0 2 0 0 20 20.58 17 25/06/2007 Caledon1 C24 FALSE 1 A 1 0 0 4 2 0 0 7 17.37 18 25/06/2007 Caledon1 C27 FALSE 1 A 1 0 0 0 3 0 0 6 24.14 19 25/06/2007 Caledon1 C28 FALSE 1 A 1 0 0 5 1 0 0 6 21.58 20 25/06/2007 Caledon1 C33 FALSE 1 A 0 0 0 0 1 0 0 0 19.06 21 25/06/2007 Caledon1 C35 FALSE 1 A 1 0 3 0 2 0 0 20 24.55 22 24/06/2007 Caledon1 C4 FALSE 1 A 1 1 4 3 0 0 0 5 22.65 23 24/06/2007 Caledon1 C8 FALSE 1 A 0 0 0 0 3 0 0 3 23.08 24 01/06/2007 Hillsborough0 H11 FALSE 1 A 1 1 4 0 0 0 0 0 23.01 25 01/06/2007 Hillsborough0 H16 FALSE 1 A 1 1 5 0 0 0 0 0 18.18 26 01/06/2007 Hillsborough0 H17 FALSE 1 A 1 1 4 0 0 0 0 0 24.45 27 01/06/2007 Hillsborough0 H19 FALSE 1 A 0 1 0 0 0 0 0 3 21.01 28 01/06/2007 Hillsborough0 H2 FALSE 1 A 1 1 5 0 0 0 0 1 24.86 29 01/06/2007 Hillsborough0 H21 FALSE 1 A 1 1 5 0 0 0 0 0 20.48 30 01/06/2007 Hillsborough0 H25 FALSE 1 A 1 1 5 0 0 0 0 0 22.09 31 01/06/2007 Hillsborough0 H26 FALSE 1 A 1 1 5 0 0 0 0 0 20.18 32 01/06/2007 Hillsborough0 H3 FALSE 1 A 0 1 5 0 0 0 0 0 25.82 33 01/06/2007 Hillsborough0 H34 FALSE 1 A 1 0 0 4 0 0 0 1 18.78 34 01/06/2007 Hillsborough0 H35 FALSE 1 A 0 1 4 0 0 0 0 0 20.02 35 01/06/2007 Hillsborough0 H5 FALSE 1 A 0 1 5 0 0 0 0 0 19.58 36 01/06/2007 Hillsborough0 H6 FALSE 1 A 1 1 5 0 0 0 0 6 28.84 37 16/06/2007 LoughNavar1 LN1 FALSE 1 A 1 1 5 0 3 0 0 1 22.71 38 16/06/2007 LoughNavar1 LN10 FALSE 1 A 1 1 5 0 0 0 0 0 23.53 39 16/06/2007 LoughNavar1 LN11 FALSE 1 A 1 1 4 0 0 0 0 2 33.24 40 16/06/2007 LoughNavar1 LN13 FALSE 1 A 1 1 5 0 0 0 0 0 20.20 41 16/06/2007 LoughNavar1 LN14 FALSE 1 A 1 1 0 7 0 0 1 1 31.35 42 16/06/2007 LoughNavar1 LN16 FALSE 1 A 1 1 0 5 0 0 0 2 19.43 43 16/06/2007 LoughNavar1 LN18 FALSE 1 A 1 0 6 0 0 0 0 12 24.30 44 16/06/2007 LoughNavar1 LN23 FALSE 1 A 1 1 5 0 0 0 0 1 24.34 45 16/06/2007 LoughNavar1 LN7
Re: [R] Newbie
I would say if you just want to prove that older deer have more scars try the Mann Whitney non parametric test... Forgive me but even that does not really make sense since the values are all 0 so it is to obvious... __ 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, reproducible code.
Re: [R] ECDF, distribution of Pareto, distribution of Normal
Original Message Subject: [R] ECDF, distribution of Pareto, distribution of Normal From: livia [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: Tue Jul 10 2007 18:35:04 GMT+0200 Hello all, I would like to plot the emperical CDF, normal CDF and pareto CDF in the same graph and I amusing the following codes. z is a vector and I just need the part when z between 1.6 and 3. plot(ecdf(z), do.points=FALSE, verticals=TRUE, xlim=c(1.6,3),ylim=c(1-sum(z1.6)/length(z), 1)) x - seq(1.6, 3, 0.1) lines(x,pgpd(x, 1.544,0.4373,-0.2398), col=red) There is something wrong with your pgpd function, see ?pgpd for help and parameters... (I wonder how you got something plotted here...) y - seq(1.6, 3, 0.1) lines(y,pnorm(y, mean(z),sqrt(var(z))), col=blue) The emperical CDF and normal CDF look rather resonable, but the pareto CDF looks quite odd. I am not sure whether I plot the pareto CDF correctly e.g. in the right yaxs or any other mistake? At the same time, let t represents the vector whose values are larger than 1.6(the part we want). If I implement the following codes and plot the emperical CDF and pareto CDF, the pareto CDF seems fit. plot(ecdf(t), do.points=FALSE, verticals=TRUE) x - seq(1.6, 3, 0.1) lines(x,pgpd(x, 1.544,0.4373,-0.2398), col=red) Could anyone give me some advice on this? Many thanks. __ 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, reproducible code.
Re: [R] Help with write.foreign (exporting data to Stata)
I am not sure what you are doing there but what you need is library(foreign) and write.dta() see ?write.dta once you have loaded the foreign package Stefan Original Message Subject: [R] Help with write.foreign (exporting data to Stata) From: kdestler [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: Tue Jul 10 2007 19:37:54 GMT+0200 Hi. I'm trying to export a dataframe from R into Stata to use a statistical function I have there. I attached library write.foreign and renamed my variables to get them to match Stata's required format, and now have the following error: file /tmp/Rtmps7rmrM/file1c06dac8.raw not found Other than typing write.foreign, do I need to do something in R to get it to save the file on my hard drive? When I search for the file name on my computer nothing comes up. I'm using a Mac in case that makes a difference. Thanks, Kate __ 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, reproducible code.
Re: [R] help on fisher.test(stats)?
Original Message Subject: [R] help on fisher.test(stats)? From: zhijie zhang [EMAIL PROTECTED] To: R-help@stat.math.ethz.ch Date: 09.07.2007 09:03 Dear friends, My dataset have many zeros, so i must use fisher exact test . Unfortunately, the fisher.test(stats) function fail to do it. Anybody knows how to do the fisher exact test with many zeros in the dataset? My dataset is: a-matrix(c(0,1,0,0,0,0,1,0,1,0,0,0,0,1,0,1,1,0,2,1,5,1,1,6,4,4,1,17,2,8,5,7,1,1,24,3,6,1,1,3,2,16,7,4,0,2,4,0,17,0,1,0,0,0,1,2),nrow=8,byrow=TRUE) data.frame(a) b-a[,-7] as.matrix(b) c-as.matrix(b) c [,1] [,2] [,3] [,4] [,5] [,6] [1,]010000 [2,]010000 [3,]011021 [4,]116441 [5,]285711 [6,]361132 [7,]740240 [8,]010001 fisher.test(c,workspace=20) ŽíÎóÓÚfisher.test(c, workspace = 2e+17) : ÍâœÓº¯Êýµ÷ÓÃʱ²»ÄÜÓÐNA(arg10) ŽËÍâ: Warning message: Ç¿Öƞıä¹ý³ÌÖвúÉúÁËNA Any suggestion or help are greatly appreciated. Your workspace is by far to large. I have done it with fisher.test(c,workspace=4000) Fisher's Exact Test for Count Data data: c p-value = 0.01548 alternative hypothesis: two.sided (btw. it took half an hour...) Simulation would also be an alternative approach: fisher.test(c,simulate=T) Fisher's Exact Test for Count Data with simulated p-value (based on 2000 replicates) data: c p-value = 0.01349 alternative hypothesis: two.sided As you see the p-value is not that different, you could use more replications: fisher.test(c,simulate=T,B=100) Fisher's Exact Test for Count Data with simulated p-value (based on 1e+06 replicates) data: c p-value = 0.01514 alternative hypothesis: two.sided and it is still much faster... Stefan -=-=- ... The most incomprehensible thing about the world is that it is comprehensible. (A. Einstein) __ 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, reproducible code.
Re: [R] Windows Binary for ncdf package
Original Message Subject: [R] Windows Binary for ncdf package From: amna khan [EMAIL PROTECTED] To: R-help@stat.math.ethz.ch, R-help@stat.math.ethz.ch Date: 08.07.2007 19:52 Dear Sir There is no window binary version of package ncdf in the latest release of R 2.5.1. i dont have any information about the old versions. Please guid in thie regard Thank you It is there! install.packages(ncdf) trying URL 'http://www.stats.ox.ac.uk/pub/RWin/bin/windows/contrib/2.5/ncdf_1.6.zip' Content type 'application/zip' length 231140 bytes opened URL downloaded 225Kb package 'ncdf' successfully unpacked and MD5 sums checked Maybe you should try another mirror. Stefan version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 5.1 year 2007 month 06 day27 svn rev42083 language R version.string R version 2.5.1 (2007-06-27) -=-=- ... Money: There's nothing in the world so demoralizing as money. (Sophocles) __ 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, reproducible code.
Re: [R] How to install R 2.5 with Synaptic in Ubuntu?
to end of the entry making it: deb http://my.favorite.cran.mirror/bin/linux/ubuntu feisty main However after this it still complains that it can't find packages.gz Just a guess: have you replaced the my.favorite.cran.mirror by a mirror which is close to you? If you're in UK it would be for example deb http://www.stats.bris.ac.uk/R/bin/linux/ubuntu feisty main ;o) Stefan It appears to be looking in http://my.favorite.cran.mirror/bin/linux/ubuntu/distsfeisty which isn't the directory structure of the cran repository, but I can see anyway to modify this behaviour. Every other Ubuntu repositoy I have looked at contains the dists directory. Any suggestions for modifying this behaviour are gratefully recieved. Many thanks Mike Smith -=-=- ... The simple truth is that interstellar distances will not fit into the human imagination - (The Hitchhiker's Guide to the Galaxy) __ 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, reproducible code.
Re: [R] How to install R 2.5 with Synaptic in Ubuntu?
msmith schrieb: Hi, Thanks for the suggestion and I wish the solution was that obvious, but I have changed it to really point at my favourite mirror. Using your example Synaptic reports the following error when I try to update the repositories: http://www.stats.bris.ac.uk/R/bin/linux/ubuntu/dists/feisty/main/binary-i386/Packages.gz: 404 Not Found This is understandable since that location doesn't exist, but it makes me think that the directory structure of the R mirrors is not compatible with Ubuntu and Synaptic, since it automatically seeks /dists/feisty/ rather than just /feisty/ as it is on the CRAN mirrors. Hm. I have Fedora 7 so I cannot really check what my entry would look like. I recently installed Kubuntu on a friends notebook so there the readme from http://www.stats.bris.ac.uk/R/bin/linux/ubuntu/README.html did actually work. I just had a second look at your mail so could it be that there is a missing slash after feisty? If that does not work try another mirror like deb http://stat.ethz.ch/CRAN/bin/linux/ubuntu/ feisty/ there is also no main stated there in the readme... sorry I copied that from your mail so thats obviously not following the readme... I hope this works now. Stefan -=-=- ... Satisfaction does not come with achievement, but with effort. Full effort is full victory. (M.Gandhi) __ 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, reproducible code.
Re: [R] How to install R 2.5 with Synaptic in Ubuntu?
Stefan Grosse schrieb: deb http://www.stats.bris.ac.uk/R/bin/linux/ubuntu feisty main Sorry, I copied a mistake there, it should be: deb http://www.stats.bris.ac.uk/R/bin/linux/ubuntu feisty/ Stefan -=-=- ... The simple truth is that interstellar distances will not fit into the human imagination - (The Hitchhiker's Guide to the Galaxy) __ 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, reproducible code.
Re: [R] A More efficient method?
Cat=c('a','a','a','b','b','b','a','a','b')# Categorical variable C1=vector(length=length(Cat)) # New vector for numeric values # Cycle through each column and set C1 to corresponding value of Cat. for(i in 1:length(C1)){ if(Cat[i]=='a') C1[i]=-1 else C1[i]=1 } C1 [1] -1 -1 -1 1 1 1 -1 -1 1 Cat [1] a a a b b b a a b how about: Cat-c('a','a','a','b','b','b','a','a','b') c1- -2*(Cat==a)+1 -=-=- ... Time is an illusion, lunchtime doubly so. (Ford Prefect) __ 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, reproducible code.
Re: [R] A More efficient method?
Gabor Grothendieck wrote: set.seed(1) C - sample(c(a, b), 10, replace = TRUE) system.time(s1 - ifelse(C == a, 1, -1)) user system elapsed 0.370.010.38 system.time(s2 - 2 * (C == a) - 1) user system elapsed 0.020.000.02 system.time(s1 - ifelse(C == a, 1, -1)) user system elapsed 0.040.010.08 system.time(s2 - 2 * (C == a) - 1) user system elapsed 0 0 0 I am just wondering: how comes the time does add up to 0.05 while elapsed states 0.08 on my system? (Vista+R2.5.1) Stefan -=-=- ... Time is an illusion, lunchtime doubly so. (Ford Prefect) __ 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, reproducible code.
Re: [R] help again
hi i try to use arima and holtwinter to predict drought from 1895-2006 but i cannot read whole period of time and i try to do the exponent fitting, but it comes out as the coordinate x-y error i send the source code and data to take a look if anyone can help me, i am really new in R PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. This makes it much easier to help! (how do your numbers look like e.g.) You don't need to send the whole dataset but a few lines would be nice plus what commands you've done to receive those errors... Stefan -=-=- ... Lotteries are a tax on ignorance. (A. Smith - attributed) __ 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, reproducible code.
Re: [R] fractional calculations
is there a function in R that allows me to work with fractions without transforming them to floats (or whatever) in between? You could use the ryacas (CAS) package: http://code.google.com/p/ryacas/ and/or http://cran.r-project.org/doc/vignettes/Ryacas/Ryacas.pdf e.g: library(Ryacas) yacas(expression((1/2 + 1/8) * 1/2 )) [1] Starting Yacas! expression(5/16) Stefan Something that would calculate something like: (1/2 + 1/8) * 1/2 = 5/16 without ever transforming to 0.5 and 0.125? __ 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, reproducible code.
Re: [R] R-excel
http://cran.r-project.org/doc/manuals/R-data.html#Reading-Excel-spreadsheets plus there is a package xlsReadWrite that might be of your interest. Stefan Original Message Subject: [R] R-excel From: Erika Frigo [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: 25.06.2007 14:47 Good morning to everybody, I have a problem : how can I import excel files in R??? thank you very much Dr.sa. Erika Frigo Università degli Studi di Milano Facoltà di Medicina Veterinaria Dipartimento di Scienze e Tecnologie Veterinarie per la Sicurezza Alimentare (VSA) Via Grasselli, 7 20137 Milano Tel. 02/50318515 Fax 02/50318501 [[alternative HTML version deleted]] __ 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, reproducible code. No virus found in this incoming message. Checked by AVG Free Edition. Version: 7.5.472 / Virus Database: 269.9.6/865 - Release Date: 24.06.2007 08:33 __ 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, reproducible code.
Re: [R] transposing data.frames
Original Message Subject: [R] transposing data.frames From: Christoph Heibl [EMAIL PROTECTED] To: R-help@stat.math.ethz.ch Date: 25.06.2007 15:13 Hello, This must be simple... Thanks a lot - Christoph # Imagine you have a list, e.g: K - list(1:10, 2:11, 9:18) K # Transforming to dataframe... KK - as.data.frame(K) # ... one obtaines the list elements as column. KK # But I need the list elements as rows # How can I achieve this? Is there a simple way to transpose data.frames? yes: t(as.data.frame(K)) __ 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, reproducible code. __ 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, reproducible code.
Re: [R] Speed up R
Original Message Subject: Re:[R] Speed up R From: Matthew Keller [EMAIL PROTECTED] To: Robert McFadden [EMAIL PROTECTED] Date: 19.06.2007 23:09 but I think that windows has problems addressing that much RAM (surely the 64bit Vista is OK with it though... surely). Linux or Apple (the powermac) might be better bets if you're wanting to work with programs that use a lot of RAM. BTW, Intel does make 64 bit chips now. They use them in macs. The Core 2 Duo has 64bit processing and is not only used in Macs... The thing with the 64bit Vista is that there is much less software driver support at the moment. So if you will not use more than 4GB and want to use it for other purposes than you might prefer the 32bit Vista... My personal subjective feeling is that Vista is slower then XP. So I would say thats not the only reason to go for Linux. Stefan __ 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, reproducible code.
Re: [R] xtable results doesn't correspond to data.frame
It would be helpfull if you would state what exactly you did expect. As improvements you could reduce the digits, see ?xtable Original Message Subject: [R] xtable results doesn't correspond to data.frame From: Vumani Dlamini [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: 20.06.2007 11:20 Dear useRs, Am trying to use xtable on the following data.frame and I don't get what I expect: example.table - data.frame(rbind( c(Gender, , , ), cbind(rep( ,2),c(Male,Female),c(3.0,4.0),c(3/7,4/7)) )) colnames(example.table) - c( , ,number of patients,%) example.xtable - xtable(example.table) print.xtable(example.xtable, include.rownames=FALSE) I can seem to get latex output which corresponds to the data.frame which is, \begin{table}[ht] \begin{center} \begin{tabular}{} \hline number of patients \% \\ \hline Gender\\ Male 3 0.428571428571429 \\ Female 4 0.571428571428571 \\ \hline \end{tabular} \end{center} \end{table} __ 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, reproducible code. __ 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, reproducible code.
Re: [R] xtable results doesn't correspond to data.frame
Original Message Subject: Re:[R] xtable results doesn't correspond to data.frame From: Stefan Grosse [EMAIL PROTECTED] To: Vumani Dlamini [EMAIL PROTECTED] Date: 20.06.2007 12:20 Sorry, there is a line missing, I was asking what exactly was not as expected. As improvements you could reduce the digits, see ?xtable Original Message Subject: [R] xtable results doesn't correspond to data.frame From: Vumani Dlamini [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: 20.06.2007 11:20 __ 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, reproducible code.
Re: [R] xtable results doesn't correspond to data.frame
That table below is not the table from your first mail. However I would do something like: library(xtable) testdata-data.frame(gender=c(rep(male,3),rep(female,4)) testtab-as.data.frame(table(testdata)) pct-function(x){ x/sum(x)*100 } testtab$percent-pct(testtab$Freq) rownames(testtab)-c(1,2) colnames(testtab)-c(gender,Freq,%) #it looks like: testtab #and latex: xtable(testtab,digits=2) % latex table generated in R 2.5.0 by xtable 1.4-6 package % Wed Jun 20 18:25:22 2007 \begin{table}[ht] \begin{center} \begin{tabular}{rlrr} \hline gender Freq \% \\ \hline 1 female3 42.86 \\ 2 male4 57.14 \\ \hline \end{tabular} \end{center} \end{table} I hope this is now how expected. Cheers Stefan Original Message Subject: Re:[R] xtable results doesn't correspond to data.frame From: Vumani Dlamini [EMAIL PROTECTED] To: [EMAIL PROTECTED] Date: Wed Jun 20 2007 13:22:59 GMT+0200 This is what is not expected. Having Gender repeated twice and the labels for gender missing. Will format the digits and align once everything is right. Thanks. \begin{table}[ht] \begin{center} \begin{tabular}{llrr} \hline number of patients \% \\ \hline Gender Gender \\ 3 0.428571428571429 \\ 4 0.571428571428571 \\ \hline \end{tabular} \end{center} \ __ 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, reproducible code.
Re: [R] to read table
Why are you not using read.csv or read.csv2 when you are reading a csv file? Original Message Subject: [R] to read table From: elyakhlifi mustapha [EMAIL PROTECTED] To: R-help@stat.math.ethz.ch Date: 18.06.2007 09:23 Hello, I have a problem to read a csv table. To read it I used this syntax donParCara - read.table(C:/Documents and Settings/melyakhlifi/Mes documents/feuilles excel/calcul2.csv,header=TRUE,sep=;,quote=,dec=,) I don't understand my errors Erreur dans scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : la ligne 16 n'avait pas 21 éléments Can you help me please? __ 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, reproducible code.
Re: [R] JGR, Java and Kubuntu 7.04 ...
You could download the latest R (2.5.0) directly from CRAN, this just as a side remark, You find instructions here: http://cran.au.r-project.org/bin/linux/ubuntu/README With java you should make sure that you have the Java 5 JDK installed. I dont know how it is with Ubuntu but on fedora I have to set the java alternative (sun java installation does not change the default java there). Which java is installed you can check with java -version there you see it it is properly installed and a R CMD javareconf also indicates if your java is configured correctly. Additionally on Fedora I had to set JAVA_HOME to the JDK directory manually. Stefan Original Message Subject: [R] JGR, Java and Kubuntu 7.04 ... From: John Logsdon [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: 14.06.2007 19:16 R-ists Yet again Java rears its ugly head. I have Kubuntu 7.04 running the Kubuntu-repository version of R 2.4.1-1. Yes it isn't the very latest version but this is not the issue here. I want a Windows-like environment and everyone is talking about JGR. I downloaded it and installed it along with rJava. Both compile and install satisfactorily. But when I come to run it: library('JGR') Loading required package: rJava Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library '/usr/local/lib/R/site-library/rJava/libs/rJava.so': /usr/local/lib/R/site-library/rJava/libs/rJava.so: undefined symbol: JNI_GetCreatedJavaVMs Error: .onLoad failed in 'loadNamespace' for 'rJava' Error: package 'rJava' could not be loaded When I first tried to run it without Java being installed, I got a message saying that JDK wasn't installed but mentioned 1.4.2. The version of Java actually installed as the latest from the Ubuntu repository is Sun 1.5.0.11. I don't see the point in installing old versions of Java just for one application because the language, or at least the writing, should be backwards compatible. In all aspects I have seen Kubuntu is a very impressive in checking compatibility. Unfortunately this is frequently not the case with Java. I steer clear of Java as much as possible. Can anyone suggest what I should do? Use Windows perhaps? Run Windows in a kvm virtual machine just to run R? Put my head in a bucket of cold water? Is there an alternative IDE? Is there a later JGR somewhere that is not yet on CRAN? TIA __ 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, reproducible code.
Re: [R] Panel data
Panel data is that you have e.g. several subjects (=individuals,firms,households) being tracked over a certain time (meaning there are repeated observations). One could roughly say that it is a collection of (quite similiar) time series. This makes the analysis better than pooling the whole set but it needs certain econometric tools. An example would be estimating wage elasticities. Time series analysis usually just involves one subject (interest rates, share prices etc.) over some time. Standard books on panel data analysis are also Baltagi: econometric analysis of panel data and Wooldridge: econometric analysis of cross sectional and panel data. Both are not introductionary texts since you do not know the difference between time series and panel data maybe you should start with an introductionary econometrics book like Wooldridge introductory econometrics or Verbeek: modern econometrics. Stefan on Thursday 14 June 2007 07:33:05 Megh Dal wrote: MD Then what is the difference between panel data and time series data? You said panel data is data on same subject being tracked over time. But time series data also do the same. Please correct me if I am wrong. MD __ 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, reproducible code.
Re: [R] Linux equivalent to windows menus and script window
JGR is pretty close to the windows R-GUI (with the menus) but even better since it is system independent and e.g. has syntax highlightning in the editor. You could have a look at the screenshots at http://rosuda.org/JGR/index.htm What you need is a Java 5 SDK installed at the system and probably R has to be reconfigured to find java. Stefan John, on Tuesday 12 June 2007 23:53:07 you wrote: JS Stefan, JS Am I correct that neither option will provide the menu options provided under windows, i.e. the MISC, PACKANGE, or WINDOWS menus? JS Thanks, JS John __ 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, reproducible code.
Re: [R] update packages with R on Vista: error
Dear all, I figured out that the best solution for me is to install R in the /User/me/Documents/R/ directory instead of the program files directory since I am the only user on my notebook and R has then full write access and I do not need to change any environment variables... Thanks for the answers,- Stefan Original Message Subject: [R] update packages with R on Vista: error From: Stefan Grosse [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: 07.06.2007 12:00 Dear R-list, I have encountered the following error message trying to update R packages: update.packages(ask='graphics') Warning in install.packages(update[instlib == l, Package], l, contriburl = contriburl, : 'lib' is not writable Error in install.packages(update[instlib == l, Package], l, contriburl = contriburl, : unable to install packages I remember did not have the problem on the last update where R installed the files then in the Documents/R folder on my user account. Any ideas how to handle this? I made the directories completely writable so I do not know where the problem is now (especially since update worked before...) Stefan PS: Tinn-R 1.19.2.3 + R 2.5.0 on Vista Business __ 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, reproducible code. __ 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, reproducible code.
Re: [R] Problems with Vista, R 2.5.0 and function save
Hi Maja, if you are working on your own computer an alternative could be to install R in your user directory: c:\Users\your_loginname_\R or c:\Users\your_loginname_\Documents\R I also had a problem with the writing permissions when I tried updating R packages. Btw. when I work with files I set the path directly (since I have several different paths) so this is also a work around ( save(x,y,file=c:/Users/Maja/Documents/masterfile.Rda) just note that the windows backslash has been replaced). Stefan Original Message Subject: Re:[R] Problems with Vista, R 2.5.0 and function save From: Prof Brian Ripley [EMAIL PROTECTED] To: Maja Schröter [EMAIL PROTECTED] Date: 12.06.2007 13:48 So you are running R somewhere in your file system where you do not have permission to write. Did you create a shortcut with a working directory set as the rw-FAQ advised you to? On Tue, 12 Jun 2007, Maja Schröter wrote: Hi everyone, I want to make use of the save function but it did not work. I'm using vista and R 2.5.0, winzip is installed too. Here's the code (from example ?save): x - runif(20) y - list(a = 1, b = TRUE, c = oops) save(x, y, file = xy.Rdata) Fehler in gzfile(file, wb) : kann Verbindung nicht öffnen Zusätzlich: Warning message: kann komprimierte Datei 'xy.Rdata' nicht öffnen Thank you so much for your help. Background: I want to crate a variable masterfile that I can start with data(masterfile) or attach(masterfile). I.g. Town-c(London,Miami,Rio,Lansing) Pollution-c(34,32,50,17) masterfile-data.frame(Town,Pollution) save(masterfile,file=masterfile.Rda) Kindly regards, Maja Schröter __ 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, reproducible code. No virus found in this incoming message. Checked by AVG Free Edition. Version: 7.5.472 / Virus Database: 269.8.15/847 - Release Date: 12.06.2007 21:42 __ 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, reproducible code.
Re: [R] Linux equivalent to windows menus and script window
There is, two candidates are ESS-emacs and JGR as a gui. Have a look at: http://www.r-project.org/GUI , on Tuesday 12 June 2007 19:42:33 John Sorkin wrote: JS R 2.5 under windows XP JS R 2.5 under Linux FC 6 JS JS I currently use R under windows and when I do so, I get to work in a nice windows environment complete with a very useful menu bar that includes a FILE, EDIT, MISC, PACKAGE, WINDOWS, and HELP menus. I have begun to use R under LINUX. When I run under windows, I get a standard terminal window, but I don't get the menus specific to R. Under windows I also have the ability to open a SCRIPT window, edit R code in the window and submit the code from the window. Is there any way to get a window with R specific menus and a script window? JS Confidentiality Statement: JS This email message, including any attachments, is for the so...{{dropped}} funny... __ 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, reproducible code.
Re: [R] update packages with R on Vista: error
Thanks for pointing at this. But you know, the user is writable. R is installing Packages in /Documents/R/win-library which works fine so I find it absolutely naturally that update should work as well. Especially since when I install the packages it gets the latest version, library loads this latest version but update still does want to update this latest package with the package I already installed and fails ... In my opinion the update on windows is simply buggy. I think one should definitely not turn UAC off. ( its a good security feature). Btw. MikTeX 2.6 is able to deal with UAC - I can update my latex packages without any problems even though they are in the Program File directory (and also on-the-fly installation does work) ... Stefan Original Message Subject: Re:[R] update packages with R on Vista: error From: R. Villegas [EMAIL PROTECTED] To: Stefan Grosse [EMAIL PROTECTED] Date: 07.06.2007 23:13 If R is installed within Program Files, one of Vista's security settings may interfere with the -update- process. The setting may be disabled globally by choosing: Windows (Start) menu, Control Panels, User Accounts and Family Safety (green title), User Accounts (green title), and Turn User Account Control on or off (very bottom). You will be prompted for permission to continue; click continue. On the screen you will see a checkbox titled Use User Account Control (UAC) to help protect your computer. Uncheck this and click the OK button to save the changes. Windows Vista will now allow programs, including R, to update files in Program Files. Rod. __ 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, reproducible code.
Re: [R] update packages with R on Vista: error
I was pointed at that my message might be considered as impolite. It was not intended so. I was just trying to formulate that there should be some improvement since the solutions offered were either not optimal for me (disabling security features) or where not working (FAQ). I apologize for any possible inconvenience caused by my frustration spiced up possibly with my inabilities. Stefan Original Message Subject: Re:[R] update packages with R on Vista: error From: Stefan Grosse [EMAIL PROTECTED] To: R. Villegas [EMAIL PROTECTED] Date: 08.06.2007 11:07 Thanks for pointing at this. But you know, the user is writable. R is installing Packages in /Documents/R/win-library which works fine so I find it absolutely naturally that update should work as well. Especially since when I install the packages it gets the latest version, library loads this latest version but update still does want to update this latest package with the package I already installed and fails ... In my opinion the update on windows is simply buggy. I think one should definitely not turn UAC off. ( its a good security feature). Btw. MikTeX 2.6 is able to deal with UAC - I can update my latex packages without any problems even though they are in the Program File directory (and also on-the-fly installation does work) ... Stefan __ 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, reproducible code.
[R] update packages with R on Vista: error
Dear R-list, I have encountered the following error message trying to update R packages: update.packages(ask='graphics') Warning in install.packages(update[instlib == l, Package], l, contriburl = contriburl, : 'lib' is not writable Error in install.packages(update[instlib == l, Package], l, contriburl = contriburl, : unable to install packages I remember did not have the problem on the last update where R installed the files then in the Documents/R folder on my user account. Any ideas how to handle this? I made the directories completely writable so I do not know where the problem is now (especially since update worked before...) Stefan PS: Tinn-R 1.19.2.3 + R 2.5.0 on Vista Business __ 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, reproducible code.
Re: [R] update packages with R on Vista: error
Actually the packages R wants to update are: VR, cluster, lattice, mgcv, nlme and rcompgen. I did how described in the R-Win-FAQ create a .Renviron File containing the path to the win-library that R already created (R_LIBS=C: ... ). I also tried to add R_LIBS= as Rgui parameter from within Tinn-R. Additionally I tried to leave a file named Renviron.site in the etc library. Nothing worked thus far. Interestingly installing packages does work fine even without specifying the R_LIBS path manually with any of the above mentioned methods. Even more puzzling is that even when I install eg. nlme manually via install.packages(nlme) it works but R still wants to update it. Even though e.g. library(nlme), ?nlme shows that the latest version is installed. I would guess there is some problem with the library path variable in the update program... Stefan Original Message Subject: Re:[R] update packages with R on Vista: error From: Prof Brian Ripley [EMAIL PROTECTED] To: Stefan Grosse [EMAIL PROTECTED] Date: 07.06.2007 13:01 See the rw-FAQ, which describes this in detail. Almost certainly you are trying to update the package 'cluster' which is in the main library. But as you used the GUI, we can't see that. On Thu, 7 Jun 2007, Stefan Grosse wrote: Dear R-list, I have encountered the following error message trying to update R packages: update.packages(ask='graphics') Warning in install.packages(update[instlib == l, Package], l, contriburl = contriburl, : 'lib' is not writable Error in install.packages(update[instlib == l, Package], l, contriburl = contriburl, : unable to install packages I remember did not have the problem on the last update where R installed the files then in the Documents/R folder on my user account. Any ideas how to handle this? I made the directories completely writable so I do not know where the problem is now (especially since update worked before...) Stefan PS: Tinn-R 1.19.2.3 + R 2.5.0 on Vista Business __ 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, reproducible code. __ 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, reproducible code.
Re: [R] rJava installation under linux: configuration failed
Which java have you installed? ( java -version ) You need a Sun Java and better is a 1.5 series Java, JGR seems to have some problem with the new 1.6 series. And you need the JDK, not the JRE. You can google how to do those installations on suse linux. Stefan Original Message Subject: [R] rJava installation under linux: configuration failed From: zhihua li [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: 05.06.2007 19:55 Hi netter, Recently I was trying to install rJava. The operating system is suse 10.0, and the R versionis 2.5.0. Following the instructions of R Wiki for rJava, I did configuration first: R CMD javareconf and then it showed a series of information, from what it seems that java is in the system and the configuration succeeded. Then I tried to install rJava: install.packages(rJava) following which rJava was downloaded and being installed, but during the last test step it said: can't complie a simple JNL program. Configuration Failed! Did I do something wrong? Or there's something I should do that I didn't? Thanks a lot! Sincerely Yours: Zhihua Li __ 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, reproducible code. No virus found in this incoming message. Checked by AVG Free Edition. Version: 7.5.472 / Virus Database: 269.8.9/834 - Release Date: 05.06.2007 14:38 __ 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, reproducible code.
Re: [R] rJava installation under linux: configuration failed
Precautionary note: the '1.5 series Java' (aka Java 5) does not work with JNI (and hence rJava) on some platforms, including AMD64 Linux. You can google how to do those installations on suse linux. It can be tricky, especially if you need a later Java than your OS version supports. Since the Suse Version is 10.0 there is an installation repository provided by Suse: http://download.opensuse.org/distribution/SL-10.0-OSS/inst-source-java/ you can add that via Yast and this should then provide Java 5 JDK. (provided that you are on the correct platform ...) __ 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, reproducible code.
Re: [R] opening vignetten
?vignette Original Message Subject: [R] opening vignetten From: Schmitt, Corinna [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: 06.06.2007 10:41 Dear R-Users, I have a quite stupid question. I load the GO-package with the command require(GO). Now I want to read the corresponding vignetten but I forgot the command for opening it. Please help me, Corinna __ 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, reproducible code. __ 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, reproducible code.
Re: [R] GUI for R running under Linux
Hi John, I use JGR when I am on Fedora. This is a good java gui which is similiar (and better) compared to the R windows gui . You can install it via install.packages(JGR,dep=T) and start with library(JGR). Sometimes Java is not configured correctly so you might end up with an error. In that case do R CMD javareconf on a root shell. Stefan Original Message Subject: [R] GUI for R running under Linux From: John Sorkin [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: 05.06.2007 02:15 Colleagues, I have been using R under windows for some time, and have become accustomed to the GUI that accompanies R. I have recently begun using R under Linux (Fedora), but can not find an equivalent GUI interface; at a terminal window I enter R Am I missing something? Can I get a GUI for R (version 2.5) running under Linux? Thanks, John John Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics University of Maryland School of Medicine Division of Gerontology Baltimore VA Medical Center 10 North Greene Street GRECC (BT/18/GR) Baltimore, MD 21201-1524 (Phone) 410-605-7119 (Fax) 410-605-7913 (Please call phone number above prior to faxing) Confidentiality Statement: This email message, including any attachments, is for the so...{{dropped}} __ 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, reproducible code. __ 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, reproducible code.
Re: [R] read table
Original Message Subject: [R] read table From: jiqiang yao [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Date: 05.06.2007 18:17 Hi, I'm a novice of R. I want to read the following table into R: names mpgcyl disp hp drat Mazda RX4 21.0 6160.0 110 3.90 Mazda RX4 Wag 21.0 6160.0 110 3.90 The command I used is: test - read.table(file.choose(),header=T) Does: test - read.table(file.choose(),header=T,sep=\t) work? __ 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, reproducible code.
Re: [R] BATCH
it states that it does not find the object donParEssai which is not there. It is nor there since you have given the object the name donParCara from what I see... elyakhlifi mustapha wrote: hello, I tried to run programs in BATCH like you told me but to read the results it's a lil hard first to read the results it's important to write down an outfile but when I do this I've got stil the same answer R version 2.4.1 (2006-12-18) Copyright (C) 2006 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R est un logiciel libre livré sans AUCUNE GARANTIE. Vous pouvez le redistribuer sous certaines conditions. Tapez 'license()' ou 'licence()' pour plus de détails. R est un projet collaboratif avec de nombreux contributeurs. Tapez 'contributors()' pour plus d'information et 'citation()' pour la façon de le citer dans les publications. Tapez 'demo()' pour des démonstrations, 'help()' pour l'aide en ligne ou 'help.start()' pour obtenir l'aide au format HTML. Tapez 'q()' pour quitter R. [Sauvegarde de la session précédente restaurée] #programme pour calculer les caractères donParCara - read.table(C:/Documents and Settings/melyakhlifi/Mes documents/feuilles excel/copi_donnees6.csv,header=TRUE,sep=;,quote=,dec=,) #print(subset(donParCara, select = c(Id_Cara,Form_C,X))) #valeurs observées pour les caractères observés C103 - as.numeric(as.character(subset(donParEssai, Id_Essai == 1006961 Id_Cara == 103, c(Val_O,Surf_O))[,1])) Erreur dans subset(donParEssai, Id_Essai == 1006961 Id_Cara == 103, c(Val_O, : objet donParEssai non trouvé Exécution arrêtée I thnik that it's a problem in the options but I don't know how to do correct the errors. _ Ne gardez plus qu'une seule adresse mail ! Copiez vos mails vers Yahoo! Mail [[alternative HTML version deleted]] __ 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, reproducible code. No virus found in this incoming message. Checked by AVG Free Edition. Version: 7.5.467 / Virus Database: 269.7.3/809 - Release Date: 17.05.2007 17:18 __ 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, reproducible code.
Re: [R] Problem with lme4
the lme4 function you want is probably lmer() type ?lmer btw. your R is very old we are at 2.5.0 ... Stefan Amelie LESCROEL wrote: Hi - I'm having a problem trying to use the function GLMM() from lme4. Here is what happens: library(lme4) Loading required package: Matrix Loading required package: lattice f1 - GLMM(success~yearF, data=quality, random=~1|bandnumb, family=binomial, method=PQL) Error: couldn't find function GLMM I remember having used lme4 before, without any problem. Can someone help me understanding what is wrong? I'm using R 2.2.1, could that be a problem? I also tried installing lme4 again, and updating the package, with no success. Thanks for your help, Amelie [[alternative HTML version deleted]] __ 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, reproducible code. __ 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, reproducible code.
Re: [R] A simple question regarding plot of nls objects
The example works fine for me on R 2.5.0 and WinXP... Mielke Hans wrote: Hi, I was trying to run the example of Indomethacin kinetics from the book: ## From Pinheiro/Bates, Mixed-Effects-Models in S and S-Plus, ## Springer, Second Printing 2001, Section 6.2 library(nlme) plot(Indometh) fm1Indom.nls - nls(conc~SSbiexp(time,A1,lrc1,A2,lrc2), data=Indometh) summary(fm1Indom.nls) plot(fm1Indom.nls,Subject~resid(.),abline=0) ## the last plot command gives me the error message: Subject not found in data What point am I missing? Thanks in advance, Hans __ 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, reproducible code.
Re: [R] Fehlermeldung
First: the language of this mailing list is english (not to bother you but potentially you adress more people then who could possibly help) Second: on such errors specifying your R version and your System (windows/Linux) would help Third: a minimal example is also helpfull: hopw have used pmvnorm... I suppose you use windows... I cannot replicate your error with R 2.5.0 using the example(pmvnorm) Stefan [EMAIL PROTECTED] wrote: Sehr geehrte Damen und Herren, um es kurz zu machen: Ich habe mir kürzlich das Package mvtnorm heruntergelden. Leider kommt beim Aufruf des Befehles pmvnorm folgende Fehlermeldung: .Fortran(mvtdst, N = as.integer(n), NU=as.integer(df), lower = as.double(lower), : Fortran Symbolname mvtdst nicht in der DLL für Paket mvtnorm Ich habe überhaupt keine Ahnung, was das bedeutet. Können Sie mir da weiterhelfen? Herzlichen Dank für Ihre Hilfe. Liebe Grüße, Andreas Faller __ 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, reproducible code. __ 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, reproducible code.
Re: [R] Fehlermeldung
I do not get this error message with that example you have given. Hm. Dunno. Have you installed R-2.5.0 over an old version? Maybe there could be your problem (although I doubt that) If so, try |update.packages(checkBuilt=T) this installs new packages made for 2.5.0 over old made for 2.4.1 if the version number is the same. (Or make a clean install of R-2.5.0- Remove everything and do a new install). Stefan | Original Message Subject: Re:[R] Fehlermeldung From: [EMAIL PROTECTED] To: Stefan Grosse [EMAIL PROTECTED] Date: Wed May 09 2007 15:02:10 GMT+0200 Dear Stefan, My operating system is Windows XP, my version of R is the latest (R-2.5.0). Recently I have downloaded the package mvtnorm and a problem with the command pmvnorm occured. Trying to enter the lines ... A - diag(3) A[1,2] -0.5 A[1,3] - 0.25 A[2,3] - 0.5 pvmnorm(lower=c(-Inf,-Inf,-Inf), upper=c(2,2,2),mean = c(0,0,0), corr=A) I got the following error message: .Fortran(mvtdst, N = as.integer(n), NU=as.integer(df), lower = as.double(lower), : Fortran Symbolname mvtdst nicht in der DLL für Paket mvtnorm As I can make no sense of that whatsoever, I would like to ask you for some qualified help. Thank you very much indeed. Best Regards, Andreas Faller ___ SMS schreiben mit WEB.DE FreeMail - einfach, schnell und kostenguenstig. Jetzt gleich testen! http://f.web.de/?mc=021192 __ 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, reproducible code.
Re: [R] ? R 2.5.0 alpha bug
The version 2.5.0 has left Alpha status long time ago and its final version has been released so please try the new version. Inman, Brant A. M.D. wrote: This email is intended to highlight 2 problems that I encountered running R 2.5.0 alpha on a Windows XP machine. #1 - Open script error If I click the Open folder icon on the toolbar, R opens my script files perfectly. However, when I select File Open Script MyFileLocation, I get a fatal error that causes R to close immediately. This error was reproduced on 3 consecutive occasions but has been intermittent thereafter. One of these fatal errors resulted in a typical error reporting box being generated which I sent off. I was not able to verify if this error has been reported and corrected in subsequent versions of 2.5. #2 - Bug reporting link on CRAN website broken I tried to report the bug listed above on the CRAN website but when I clicked on the bug reporting link on the left-hand side panel of the main site (http://bugs.r-project.org/cgi-bin/R) , I get an error page with the following message: The system encountered a fatal error cannot open config file /home/sfe/r-bugs/jitterbug/R : No such file or directory The last error code was: No such file or directory uid/gid=30/8 This has been submitted to r-devel. Brant Inman __ 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, reproducible code. __ 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, reproducible code.
Re: [R] Query about finding correlations
I am sorry if I misinterpreted this but this questions looks/looked very like a very basic one about statistical inference... (have a look at the Manuals/contributed documentation e.g. Faraway) Lalitha Viswanath wrote: Hi This is not a homework assignment :) Me and my manager are trying to understand the problem better. In the meanwhile, we thought we would post the problem on this forum to seek some input from statisticians who possibly do this kind of analyses everyday and hence are possibly more proficient with R and/or any recommended methodologies. Lalitha On 5/2/07, *Stefan Grosse* [EMAIL PROTECTED] mailto:[EMAIL PROTECTED] wrote: How about making your homeworks yourselfes? lalitha viswanath wrote: Hi I have a dataframe which has 3 columns of numeric data A,B,C each of which has been obtained independent of the other. We are trying to find out, which of A or B cause C i.e. We are hypothesising that C is the effect and either A or B, not both is the cause. i.e. A causes C and this cause-effect relationship explains B. The data for A contains more noise than that for B. We are working with around 1000 points. I would greatly appreciate any inputs on the best statistcal approach to tackle this problem. I am thinking that we can find correlation coefficients between A and C, and between B and C, but I am not sure this answers the question. Also we do not know whether the correlation between them is linear or non linear. Thanks Lalitha __ R-help@stat.math.ethz.ch mailto: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 http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. No virus found in this incoming message. Checked by AVG Free Edition. Version: 7.5.467 / Virus Database: 269.6.2/784 - Release Date: 01.05.2007 14:57 __ 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, reproducible code.
Re: [R] Query about finding correlations
How about making your homeworks yourselfes? lalitha viswanath wrote: Hi I have a dataframe which has 3 columns of numeric data A,B,C each of which has been obtained independent of the other. We are trying to find out, which of A or B cause C i.e. We are hypothesising that C is the effect and either A or B, not both is the cause. i.e. A causes C and this cause-effect relationship explains B. The data for A contains more noise than that for B. We are working with around 1000 points. I would greatly appreciate any inputs on the best statistcal approach to tackle this problem. I am thinking that we can find correlation coefficients between A and C, and between B and C, but I am not sure this answers the question. Also we do not know whether the correlation between them is linear or non linear. Thanks Lalitha __ 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, reproducible code. __ 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, reproducible code.
Re: [R] tcltk support....
tcltk must be installed on your system, it is not a R package. Which Linux are you using? W Eryk Wolski wrote: Rusers, library(Rcmdr) Loading required package: tcltk Error in firstlib(which.lib.loc, package) : Tcl/Tk support is not available on this system While installing R I run configure with configure --prefix=/bla/bla --with-tcltk make make install configure did not returned any errors. Hence, I assume that the the tcl/tk libraries have been found. installin.packages(tcltk) returns an error telling me that there is no package tcltk on CRAN. And indeed there no such package on CRAN. Whats wrong? Eryk [[alternative HTML version deleted]] __ 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, reproducible code. __ 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, reproducible code.
Re: [R] tcltk support....
To my excuse I can say that I would have expected configure to issue an error or warning. I am absolutely certain that there was a remark during configure... you probably just have overseen that. I did run it configure with --with-tcltk. But there was not a singel message which would point me to the fact that some packages on my system are missing. __ 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, reproducible code.
Re: [R] tcltk support....
I just read another post so I know youre on (K)ubuntu: you need to install tcl8.4-dev and tk84-dev packages. If there is no tcl/tk R installs nevertheless but is stating that at the configuration that there is no tcl/tk support... so thats where you have to look at. Maybe it is better you install the deb package of 2.5.0 instead of compiling yourself... Stefan W Eryk Wolski wrote: Rusers, library(Rcmdr) Loading required package: tcltk Error in firstlib(which.lib.loc, package) : Tcl/Tk support is not available on this system While installing R I run configure with configure --prefix=/bla/bla --with-tcltk make make install configure did not returned any errors. Hence, I assume that the the tcl/tk libraries have been found. installin.packages(tcltk) returns an error telling me that there is no package tcltk on CRAN. And indeed there no such package on CRAN. Whats wrong? Eryk [[alternative HTML version deleted]] __ 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, reproducible code. __ 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, reproducible code.
Re: [R] Graphic plots in pdf (position)
you can place the figures with print command, e.g. pdf(output.pdf) print(text or textplot,position=c(0,0,.5,.5),more=T) print(persp(zPERSP_1),position=c(.5,0,1,0),more=T) print(persp(zPERSP_2),position=c(0,0.5,0.5,1),more=T) print(persp(zPERSP_3),position=c(.5,.5,1,1)) dev.off() pdf has options where you can adjust the size. I use something like: pdf(filename, width=13, height=13, paper=special) lokk ?pdf for options probably you want paper=a4 Felix Wave wrote: Hello, I have got three graphics. pdf(output.pdf) persp(zPERSP_1) persp(zPERSP_2) persp(zPERSP_3) dev.off() I will write all those on one page of a pdf file. On the left side I will write some lines text. Is this possible, positioning graphic and text elements in pdf? I has already thought I could write a HTML file and convert it anyway to pdf. Maybe anybody does know a toturial or has an idea. Thanks a lot. Felix My aim (pdf-file): - text text - text text | | text text | | | | | | - - - | | | | | | | | | | | | | | | | - - __ 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, reproducible code. __ 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, reproducible code.
Re: [R] from R to html
Is it so much effort to use the search on r-project.org? elyakhlifi mustapha wrote: hello, If I wanna export data.frame from R to html have you got some ideas to do this? ___ [[alternative HTML version deleted]] __ 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, reproducible code. __ 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, reproducible code.
Re: [R] data frame
how about: ?str ever considered reading an introductory text? find some here: http://cran.r-project.org/other-docs.html Stefan elyakhlifi mustapha wrote: hello, I wanna print something like this Class Levels Values Id_TrT1 1 2 Id_Geno764208 64209 64210 64211 64212 64213 64214 Id_Rep 2 12 Is it possible? I have some problem I think taht I should use data.frame with matrix but I'm not sure and perhaps it's false __ 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, reproducible code.
Re: [R] subset
What format does your date have? This is essential here. However it must be something like subset(yourdata, year %in% 2004) how to extract the year from your date you must find out yourself... (depending on the dates format...) ever considered reading an introductory text? find some here: http://cran.r-project.org/other-docs.html Stefan elyakhlifi mustapha wrote: hi, ok I understand how to use the subset function but sometimes I need to use it to extract data by date and its format it isn't so easy (like below) for example like in using SQL I thougth that it was possible to write %/2004 but it doesn't run. Can you help me please about this? subset(don, Date_O in %/2004, select = c(Annee_O, Date_O)) ___ Découvrez une nouvelle façon d'obtenir des réponses à toutes vos questions ! Profitez des connaissances, des opinions et des expériences des internautes sur Yahoo! Questions/Réponses [[alternative HTML version deleted]] __ 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, reproducible code. __ 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, reproducible code.
Re: [R] matlab style plotting in R
There are many ways if I understood what you indend to do: example data: toplot-data.frame(x=c(1:5,1:5),y=rnorm(10),sbs=c(rep(a,5),rep(b,5))) 1st with plot itself: plot(y~x,data=toplot,type=n) lines(y~x,data=subset(toplot,sbs==a),type=b,pch=4,col=blue) lines(y~x,data=subset(toplot,sbs==b),type=b,pch=2,col=red) 2nd with xyplot out of the lattice package: library(lattice) xyplot(y~x,groups=sbs,data=toplot,pch=c(2,4),main=Title) to see some variants. My advice is to look at a good documentation at the homepage (contributed documentation) or book ... Maria Vatapitakapha wrote: Hello I was wondering how I can achieve matlab style plotting in R, in the sense that matlab allows you to plot multiple sets of variables within the same x-y axes. plot in R does not seem to cater for this. I tried 'overplot' from the gplots package but this assumes different y axes for the variables. any suggestions would be very appreciated Maria [[alternative HTML version deleted]] __ 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, reproducible code. __ 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, reproducible code.
Re: [R] problems installing R on Linux
Thank you both Pr. Ripley and Dr. Skyler for taking the time to help. I got your further responses in my email here at work. I really think that Dr. Ripley is correct : The symptoms look like a very old Linux distro and a missing/damaged/incompatible libg2c. # gcc --version 2.96 # g77 --version GNU Fortran 0.5.26 2731 (Red Hat Linux 7.3 2.96-110) (embarrassingly old?) My collegues and I are going to try to update it to the most recent free version (9.something). Yes it looks terribly old. You should try Fedora Core 6, a re-spin (all updates till January) is available here: http://fedoraunity.org/news-archives/core-6-re-spin-20070111 It is Red Hat based (follow up?), up to date, and a R rpm is provided by the r-project so you do not need to compile by yourself. Get it from here: http://stat.ethz.ch/CRAN/bin/linux/redhat/fc6/ Stefan __ 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, reproducible code.
Re: [R] R in Industry
I just ran on my Windows PC the benchmark from http://www.sciviews.org/benchmark/index.html which is pretty old now. Thats probably the reason for the errors which I did not correct. As you see R has some advantages but Matlab has also some advantages. However the differences are not to big. There is only one big difference which indeed includes loops (Creation of a 220x220 Toeplitz matrix) where Matlab is much faster. But maybe a simple change in the programmation can change that... Has someone in the list an updated script? Stefan Grosse The benchmarks: R 2.4.1 R Benchmark 2.3 === I. Matrix calculation - Creation, transp., deformation of a 1500x1500 matrix (sec): 0.865 800x800 normal distributed random matrix ^1000__ (sec): 0.136 Sorting of 2,000,000 random values__ (sec): 0.615 700x700 cross-product matrix (b = a' * a)___ (sec): 0.557 Linear regression over a 600x600 matrix (c = a \ b') (sec): 0 #ERROR II. Matrix functions FFT over 800,000 random values__ (sec): 0.557 Eigenvalues of a 320x320 random matrix__ (sec): 0.495 Determinant of a 650x650 random matrix__ (sec): 0.276 Cholesky decomposition of a 900x900 matrix__ (sec): 0 #ERROR Inverse of a 400x400 random matrix__ (sec): 0 #ERROR III. Programmation -- 750,000 Fibonacci numbers calculation (vector calc)_ (sec): 0.469 Creation of a 2250x2250 Hilbert matrix (matrix calc) (sec): 1.016667 Grand common divisors of 70,000 pairs (recursion)___ (sec): 0.3966671 Creation of a 220x220 Toeplitz matrix (loops)___ (sec): 0.552 Escoufier's method on a 37x37 matrix (mixed) (sec): 2.66 --- End of test --- Matlab 7.0.4 Matlab Benchmark 2 == Number of times each test is run__: 3 I. Matrix calculation - Creation, transp., deformation of a 1500x1500 matrix (sec): 0.29047 800x800 normal distributed random matrix ^1000__ (sec): 0.42967 Sorting of 2,000,000 random values__ (sec): 0.71432 700x700 cross-product matrix (b = a' * a)___ (sec): 0.14748 Linear regression over a 600x600 matrix (c = a \ b') (sec): 0.12831 -- Trimmed geom. mean (2 extremes eliminated): 0.26403 II. Matrix functions FFT over 800,000 random values__ (sec): 0.24591 Eigenvalues of a 320x320 random matrix__ (sec): 0.38507 Determinant of a 650x650 random matrix__ (sec): 0.091612 Cholesky decomposition of a 900x900 matrix__ (sec): 0.11059 Inverse of a 400x400 random matrix__ (sec): 0.069414 -- Trimmed geom. mean (2 extremes eliminated): 0.13556 III. Programmation -- 750,000 Fibonacci numbers calculation (vector calc)_ (sec): 1.2386 Creation of a 2250x2250 Hilbert matrix (matrix calc) (sec): 3.0541 Grand common divisors of 70,000 pairs (recursion)___ (sec): 1.7637 Creation of a 220x220 Toeplitz matrix (loops)___ (sec): 0.0045972 Escoufier's method on a 37x37 matrix (mixed) (sec): 0.50481 -- Trimmed geom. mean (2 extremes eliminated): 1.0331 Total time for all 15 tests_ (sec): 9.1786 Overall mean (sum of I, II and III trimmed means/3)_ (sec): 0.33316 --- End of test --- -- --- lic. oec. Stefan Grosse University of Erfurt Microeconomics Nordhäuser Str. 63 99089 Erfurt Germany phone +49-361 - 737 45 23 fax+49-361 - 737 45 29 mobile +49-1609- 760 33 01 web http://www.uni-erfurt.de/mikrooekonomie mail [EMAIL PROTECTED] __ 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, reproducible code.
Re: [R] LyX 1.4.3-5 problem generating PDF documents
You are posting in the wrong mailing list, mail your question to the lyx users list: http://www.lyx.org/internet/mailing.php John Kane wrote: I just changed machines and, at first, could not get LyX 1.4.3-4 to even load. I tried 1.5 which loaded nicely but would not give me PDF output. Removed 1.5 tried using 1.4.3-5.exe. Result : On loading LyX : lyx: Disabling socket cont class std::basic_stringchar struct std: char::char_traits ... Removed and installed from LyXWin143Complete-2-9.exe. Loads okay but no PDF generation. My PDF reader, Foxit Reader loads but nothing is produced. I have managed to produce a postscript doc and convert it using Ghostscript. Could this just be a Foxit problem? As soon as I figure out how to reset the Adobe to default pdf reader I'll try it out. __ 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, reproducible code. __ 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, reproducible code.
Re: [R] CPU or memory
64bit does not make anything faster. It is only of use if you want to use more then 4 GB of RAM of if you need a higher precision of your variables The dual core question: dual core is faster if programs are able to use that. What is sure that R cannot make (until now) use of the two cores if you are stuck on Windows. It works excellent if you use Linux. So if you want dual core you should work with linux (and then its faster of course). The Core 2 duo is the fastest processor at the moment however. (the E6600 has a good price/performance ration) What I already told Taka is that it is probably always a good idea to improve your code for which purpose you could ask in this mailing list... (And I am very sure that you have there a lot of potential). Another speeding up possibility is e.g. using the atlas library... (where I am not sure if you already use it) Stefan John C Frain schrieb: *Can I extend Taka's question?* ** *Many of my programs in (mainly simulations in R which are cpu bound) on a year old PC ( Intel(R) Pentium(R) M processor 1.73GHz or Dell GX380 with 2.8Gh Pentium) are taking hours and perhaps days to complete on a one year old PC. I am looking at an upgrade but the variety of cpu's available is confusing at least. Does any one know of comparisons of the Pentium 9x0, Pentium(r) Extreme/Core 2 Duo, AMD(r) Athlon(r) 64 , AMD(r) Athlon(r) 64 FX/Dual Core AM2 and similar chips when used for this kind of work. Does anyone have any advice on (1) the use of a single core or dual core cpu or (2) on the use of 32 bit and 64 bit cpu. This question is now much more difficult as the numbers on the various chips do not necessarily refer to the relative speed of the chips. * *John * On 06/11/06, Taka Matzmoto [EMAIL PROTECTED] wrote: Hi R users Having both a faster CPU and more memory will boost computing power. I was wondering if only adding more memory (1GB - 2GB) will significantly reduce R computation time? Taka, _ Get FREE company branded e-mail accounts and business Web site from Microsoft Office Live __ 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, reproducible code. __ 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, reproducible code.
[R] qplot of ggplot package how to plot different size according to the values and not to the weights?
I have the following Data structure $ step45 : Factor w/ 2 levels $ obserror : num 6.2 6.2 5.6 6.6 6.6 ... $ Mon: num 2.2 2.0 1.0 3.2 2.0 ... $ inc.comp : num 4 5 2 5 5 5 5 5 4 4 ... all I wanted to do is plotting Mon against obserror, the colors should be by step45 and the size of the symbol should be according to inc.comp so I did this: qplot(obserror,Mon,data=obscomp,col=inc.comp,col=step45) unfortunately the size of the is something I do not want it to be, the legend for inc.comp says: 4, 2.25, 1 ,0.25 , 0 I suppose this are weights? I got closer to what I want specify a different symbol instead of a different size: qplot(obserror,Mon,data=obscomp,shape=as.factor(inc.comp),col=step45) I needed the as.factor otherwise the plot is again different. But that did not work with size. OK I thought this is a second best solution but unfortunately the symbols are to small. But setting a size = 3 which I thout adds a constant size everywhere (somehow expectedly) did add a third legend 3 with 3.61, 3.42, ... (The color worked fine) Unfortunately the documentation of qplot is rather short (FIXME)... So is there a kind of option to use for the size just the values and not the weights? Or a possibility to switch of 3 for legend plotting? Another thing what I like to complain about is that qplot is using only half of the windows graphics device space (especially wen adding this size =3) which is not only ugly but a waste of space... Stefan R-2.4.0 Windows XP Prof. __ 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, reproducible code.
Re: [R] qplot of ggplot package how to plot different size according to the values and not to the weights?
here an example with the car data qplot(drat,wt,data=mtcars,shape=as.factor(carb),col=am,size=3) it gets a bit better if one takes away the as.factor (but that of course changes other things) The picture is saved via metafile, screenshot did not work (mostly black whyever) I know its deleted for the list... Stefan __ 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, reproducible code.
Re: [R] qplot of ggplot package how to plot different size according to the values and not to the weights?
hadley wickham schrieb: I think you mean qplot(obserror,Mon,data=obscomp, size=inc.comp,col=step45) Yes sorry, that was a typo (from playing around with the options) unfortunately the size of the is something I do not want it to be, the legend for inc.comp says: 4, 2.25, 1 ,0.25 , 0 I suppose this are weights? I think there are two things here that are confusing you (both my fault) * the size aesthetic actually modifies the area (ie. sqrt(size)) of the points, because that is better perceptually * the legend isn't very good But it is mapping the inc.comp to the size of the variable. OK So I was confused by the legend values. qplot currently doesn't have any way to set an aesthetic to a fixed value (it can't tell that you're specifying an aesthetic value not a data value). You can do this in two steps though p - qplot(wt, mpg, data=mtcars, type=NULL, shape=cyl) ggpoint(p, size=3) Thanks that looks better. And when I do the as.factor I have the different values (and as.factor in the legend but that does not matter). Another thing what I like to complain about is that qplot is using only half of the windows graphics device space (especially wen adding this size =3) which is not only ugly but a waste of space... That sounds like a bug. Could you please provide a reproducible example? already sent to you and the list. Thanks a lot. Stefan __ 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, reproducible code.
Re: [R] pdf/ convert/ font problem
Sam, Have you tried to produce a ps/eps figure with the postscript device? Does this work? you could still convert it to pdf with ghostscript or include eps figures... Stefan __ 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, reproducible code.
Re: [R] Redhat compilers and lme4 with R-2.4.0
Prof Brian Ripley schrieb: Not surprising: you have not compiled R with -std=gnu99 to allow C99 features, as recommended in the R-admin manual, and selected by default. (Comment to Doug Bates: it is probably better still to write in C90 if you can.) This seems to be a problem with the flags used to build the FC5 RPMs. I think you can circumvent it by adding PKG_CFLAGS=-std=gnu99 to lme4/src/Makevars. I am having the same problem on Fedora Core 5 with the 2.4.0 rpm from CRAN. Update from lme4 0.9975-3 to lme4 0.9975-6 failed due to exactly the same error message. I would have expected the cran rpm is compiled according to the R-admin guide...? message (the german version ;-) ): * Installing *source* package 'lme4' ... ** libs gcc -I/usr/lib/R/include -I/usr/lib/R/include -I/usr/local/include -I/usr/lib/R/library/Matrix/include -fpic -O2 -g -pipe -Wall -Wp,-D_FORTIFY_SOURCE=2 -fexceptions -fstack-protector --param=ssp-buffer-size=4 -m32 -march=i386 -mtune=generic -fasynchronous-unwind-tables -c glmer.c -o glmer.o In file included from lmer.h:4, from glmer.h:4, from glmer.c:1: lme4_utils.h: In Funktion »internal_symmetrize«: lme4_utils.h:82: Fehler: Anfangsdeklaration in »for«-Schleife außerhalb C99-Modus verwendet lme4_utils.h:83: Fehler: Anfangsdeklaration in »for«-Schleife außerhalb C99-Modus verwendet lme4_utils.h: In Funktion »internal_make_named«: lme4_utils.h:105: Fehler: Anfangsdeklaration in »for«-Schleife außerhalb C99-Modus verwendet lme4_utils.h: In Funktion »internal_getElement«: lme4_utils.h:124: Fehler: Anfangsdeklaration in »for«-Schleife außerhalb C99-Modus verwendet glmer.c: In Funktion »random_effects_deviance«: glmer.c:433: Fehler: Anfangsdeklaration in »for«-Schleife außerhalb C99-Modus verwendet glmer.c: In Funktion »glmer_init«: glmer.c:583: Fehler: Anfangsdeklaration in »for«-Schleife außerhalb C99-Modus verwendet make: *** [glmer.o] Fehler 1 ERROR: compilation failed for package 'lme4' __ 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, reproducible code.
Re: [R] Multiple histograms in one plot
Have a look at the lattice package. It is the best choice for such a task imho. Johann Hibschman schrieb: Hi all, I'm trying to plot multiple histograms in one plot (cross-validation values of model parameters), but I cannot seem to reduce the margins enough to fit as many of them in as I would like. I'm using split.screen to divide the window into a 5x4 grid, then plotting with hist. I've tried explicitly reducing the margins with par(mar=c(1,1,1,1)), but it doesn't seem to have any effect. Visually, there is a lot of whitespace and very little histogram in my results. Can anyone suggest either a better method to visualize these results, a better way to plot histograms, or a way to actually reduce the margins used? The intent is to give a sense of how well-constrained the various model parameters are. Thanks, Johann __ 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, reproducible code. __ 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, reproducible code.
Re: [R] Latex editor recommendations
For LaTeX I use the winedt editor: http://www.winedt.com which is extremely powerful and there exists an add-on for the use with R (but I prefer using Tinn-R for R). Winedt is shareware, you have a trial period of 30 days. It has an affordable price and discount for academics if I remember correctly. A very good freeware is TeXnicCenter: http://www.texniccenter.org/ There are many other editors being capable of LaTeX editing like Jedit (which has a LaTeX plugin), xemacs plus auctex and so on but the two mentioned are the LaTeX specialists and in my opinion easier to use. (PS if you look for a nice reference-library for the use with LaTeX Jabref is very good http://jabref.sourceforge.net/ free and works excellent together with winedt (inserts citation into winedt by mouse-click), has a connection to citeseer and medline for receiving citations, can import endnote libraries and so on and so forth) Tom Backer Johnsen schrieb: This question is not oriented towards R, but is posted here because I have the impression that there are at least some Latex users among the contributors. The question is: What editors for Latex are to be recommended? I have located one: http://www.latexeditor.org/ What other alternatives are there? I am (for the most part) using Windows XP. Tom __ 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, reproducible code. __ 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, reproducible code.
Re: [R] Redhat compilers and lme4 with R-2.4.0
I had no problem on Fedora Core 5. Patrick Connolly schrieb: I've tried unsuccessfully to install lme4 on two different Redhat systems but I've not noticed anything on this list to indicate anyone else has had similar problems. __ 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, reproducible code.
Re: [R] Latex editor recommendations
There is also a rather old version running natively under Windows 1.0.5 something. But I think TeXmacs is somewhat hard to use and would always recommend LyX unless someone likes to embed calculations in a LaTeX like (which is I admit a really nice feature of TeXmacs. especially with maxima ) I even recommend LyX for Windows since you don't need the cygwin stuff - it runs nicely and native under Windows. So especially if someone is a newbie the installation and use of LyX is really much easier especially compared to any cygwin approach if one does not like the direct programming of LaTeX code in a nice editor environment like that of WinEdt or TeXnicCenter or whatever Niels Vestergaard Jensen schrieb: In a similar vein there's TeXmacs, which is somewhat WYSIWYG. (My apologies if it's already been mentioned). It's originally a linux program, but it's possible to run it under cygwin: http://www.fmi.uni-passau.de/~seidl/texmacs/tmtour.html#2 best Niels __ 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, reproducible code.
Re: [R] Mixed effect model in R
Interesting packages for you might be the nlme and lme4 packages and as a book Pinheiro/Bates, Mixed-Effects Models in S and S-Plus Lina Jansen schrieb: Hi, I am analysing an experiment that has one fixed (6 conditions) and two random factors (11 subjects, 24 images in the conditions). I read somewhere else that you can also see such a design as a nested experiment with the hierarchy: subjects - condition - image. For some analysis I have one respond variable and for others I have more. The response variables are non-normally distributed. Now the question: Is there a package that can deal with such a design? I would like to use a generalized linear model. Are there glms that are extended to do multivariate analysis (for the 2 random + 1 fixed variable design)? And how do you call such a design? Last question: Can you suggest me some literature about such a problem? I am quite unsure concerning the analysis. Thanks for any advice lisra [[alternative HTML version deleted]] __ 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, reproducible code. __ 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, reproducible code.
Re: [R] Mixed effect model in R
Please always reply to the list as well as there always might be someone faster/better answering (or it could be that I am wrong, so someone might correct me) Indeed Pinheiro/Bates assume gaussian error terms... but I am not really sure whether you meant that with non normally distributed respond variable resp. with non-normal data however: / Mixed-effects models: / The recommended nlme http://cran.r-project.org/src/contrib/Descriptions/nlme.html package, associated with Pinheiro and Bates, / Mixed-Effects Models in S and S-PLUS / (Springer, 2000), fits linear and nonlinear mixed-effects models, commonly used in the social sciences for hierarchical and longitudinal data. Generalized linear mixed-effects models may be fit by the glmmPQL function in the MASS package, and by the lmer function in the Matrix http://cran.r-project.org/src/contrib/Descriptions/Matrix.html package (related to the lme4 http://cran.r-project.org/src/contrib/Descriptions/lme4.html package, which largely supersedes nlme http://cran.r-project.org/src/contrib/Descriptions/nlme.html for / linear / mixed models). Also see the lmeSplines http://cran.r-project.org/src/contrib/Descriptions/lmeSplines.html and lmm http://cran.r-project.org/src/contrib/Descriptions/lmm.html packages. [ http://cran.r-project.org/src/contrib/Views/SocialSciences.html ] Lina Jansen schrieb: 2006/10/17, Stefan Grosse [EMAIL PROTECTED] mailto:[EMAIL PROTECTED]: Interesting packages for you might be the nlme and lme4 packages and as a book Pinheiro/Bates, Mixed-Effects Models in S and S-Plus Thank you for the answer. I am always unsure concerning the non-normality. Can I use the nlme and lme4 with non-normal data? First, I thought they would work like an ANOVA but with random and fixed effects. __ 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, reproducible code.
Re: [R] Barplot legend position
e.g.: barplot(x,col=c(lightgrey,darkgrey),besid=T) legend(topleft,c(left,right),fill=c(lightgrey,darkgrey)) try: ?legend and example(legend) for documentation!!! Ingmar Visser schrieb: Thanks, this could work! However, the legend does not reproduce the color/shading used in the original barplot, are those available somehow? Best, Ingmar From: David Hajage [EMAIL PROTECTED] Date: Fri, 13 Oct 2006 14:11:21 +0200 To: Ingmar Visser [EMAIL PROTECTED] Cc: R-help@stat.math.ethz.ch Subject: Re: [R] Barplot legend position For example : x=matrix(1:10,2,5) barplot(x,besid=T) legend(topleft, c(left,right), density= c(0,1000)) __ 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, reproducible code.
Re: [R] Need help with tables
I suggest you look at one of the guides: http://cran.r-project.org/other-docs.html before answering questions like this to the mailing list... please read the posting guide! laba diena schrieb: I have a data file with 3 columns and I need to take only 2, how to do that ? [[alternative HTML version deleted]] __ 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, reproducible code. __ 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, reproducible code.
Re: [R] 2 linux/R environment questions
Joe Byers schrieb: JGR, when I manually compiled it for my system by specifying my location of java 1.5.0_06, worked fairly nicely. The autoinstall does not find my installation of this java version since Redhat EL4 does not support it. Redhat EL5 will in the near future. Maybe JGR will autoinstall then. Can you tell how you worked that out (the manual compilation)? I failed on Fedora Core 5 with JRE 1.5.0_08 with install.packages and dep=true. I tried a hint from the JGR mailing list but failed... Stefan Grosse __ 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, reproducible code.
Re: [R] 2 linux/R environment questions
Prof Brian Ripley schrieb: Which arch and how did you install R? For me, it works with R compiled from the sources (but not installed from RPM) on FC5 i686, but not x86_64. The issue is a run-time one: I have i686, no x86_64 ... OK, so I have to compile R from source ... Thanks for the reply. Stefan Grosse __ 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, reproducible code.
Re: [R] regarding bootstrapping
Although it does not directly answer your question you might be interested in: Michalis Ioannides, A comparison of yield curve estimation techniques using UK data, Journal of Banking Finance, Volume 27, Issue 1, , January 2003, Pages 1-26. Abstract: I compare different methods of estimating the term structure of interest rates on a daily UK treasury bill and gilt data that spans the period from January 1995 to January 1999. In-sample and out-of-sample statistics reveal the superior pricing ability of certain methods characterised by an exponential functional form. In addition to these standard goodness of fit statistics, model performance is judged in terms of two trading strategies based on model residuals. Both strategies reveal that parsimonious representations of the term structure perform better than their spline counterparts characterised by a linear functional form. This is valid even when abnormal returns are adjusted for market movements. Linear splines overfit the data and are likely to give misleading results. Author Keywords: Bonds; Trading strategy; Term structure of interest rates; Pricing; Splines [EMAIL PROTECTED] schrieb: Hi All, Thanks for all your discussions which I read offline and enhance my knowledge. I just want to know how to make yield curve by bootstrapping using R. Please give me some subject-matter links and code links. Especially the Spot Yield Curve and YTM Curve. Sayonara With Smile With Warm Regards :-) G a u r a v Y a d a v Senior Executive Officer, Economic Research Surveillance Department, Clearing Corporation Of India Limited. Address: 5th, 6th, 7th Floor, Trade Wing 'C', Kamala City, S.B. Marg, Mumbai - 400 013 Telephone(Office): - +91 022 6663 9398 , Mobile(Personal) (0)9821286118 Email(Office) :- [EMAIL PROTECTED] , Email(Personal) :- [EMAIL PROTECTED] DISCLAIMER AND CONFIDENTIALITY CAUTION:\ \ This message and ...{{dropped}} __ 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, reproducible code. DISCLAIMER AND CONFIDENTIALITY CAUTION:\ \ This message and ...{{dropped}} __ 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, reproducible code. __ 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, reproducible code.
Re: [R] R Graphics: Saving PDF and other formats from Windows Graphic Device for LaTeX
pdf(sgr6100.pdf, horizontal=FALSE, onefile=FALSE, height=3, width=3, pointsize=6) Reducing point-size below 6 does not seem to make any difference to the size of text and symbols. Any suggestions to get smaller font sizes? I have never used the pointsize option. Increasing the height and the width decreases the relative size of the font so if you scale down the figure in LaTeX you would get smaller fonts. I am using WinEdt with MikTeX set-up. as well do I ... Latex complains about the [scale...] part for any scale in \includegraphics[scale=1]{} A scale of one does not make any sense since it is 100% or with other words you could leave that out. If you want to scale the figure to the text width you could write \begin{figure} \includegraphic[width=\textwidth]{sgr6100.pdf } \end{figure} btw. I use the postscript device to produce an eps file and dvipdfm to convert the dvi to pdf postscript(sgr6100.eps, width = 6, height = 6, horizontal = FALSE, onefile =T, paper=special) and then in LaTeX: \begin{figure} \includegraphic[width=\textwidth]{sgr6100} \end{figure} and pauses. Suggestions will be appreciated about what is the best way to scale R graphics for inclusion in LaTeX. see above. Have a look at a good LaTeX guide like lshort http://tobi.oetiker.ch/lshort/lshort.pdf (e.g. p.73/73) Stefan Grosse __ 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, reproducible code.
Re: [R] Gnuplot epslatex format also in R?
I did an eps/ps file with e.g.: postscript(c:/Temp/test.eps, width = 8.0, height = 6.0, horizontal = FALSE, onefile = FALSE, paper = special) curve(x^2) dev.off() but I am not sure what you mean with the second file with the latex commands. Paul Smith schrieb: Dear All Is there some way of exporting R plots to epslatex, i.e., to a file with the eps file and another one with the LaTeX commands (representing the text in the plots), likewise Gnuplot does? If so, could you please indicate it to me? Thanks in advance, Paul __ 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, reproducible code. __ 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, reproducible code.
Re: [R] converting decimal - hexadecimal
how about learning to use help.search? (In may you already asked a similiar question...) help.search(hexadecimal) which would lead you to format.hexmode ?format.hexmode Romain Lorrilliere schrieb: Hi, do you know, a method to convert an decimal value (integer) to the corresponding hexadecimal value ? thinks for help. Romain __ 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, reproducible code.
Re: [R] screen resolution effects on graphics
how about producing pdf output- should be the same on every PC... Charles Annis, P.E. schrieb: Greetings, R-Citizens: I have the good fortune of working with a 19 1280 X 1024 pixel monitor. My R-code produces nice-looking graphics on this machine but the same code results in crowded plots on an older machine with 800 X 600 resolution. In hindsight this seems obvious, but I didn't anticipate it. My code will be used on machines with varying graphics (and memory) capacity. Is there a way I can check the native resolution of the machine so that I can make adjustments to my code for the possible limitations of the machine running it? Thanks. Charles Annis, P.E. [EMAIL PROTECTED] phone: 561-352-9699 eFax: 614-455-3265 http://www.StatisticalEngineering.com __ 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, reproducible code. __ 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, reproducible code.
Re: [R] Authoring a book
//people.ee.ethz.ch/~oetiker/lshort/lshort.pdf It really was a not too short intro. I'll have a look at it. Yes definitly not too short. But it states in LaTeX in 133 min ... Seems to be for Linux only. My server is Windows, even if I have the rest of the components. Hm at the projects' homepage its stated that it is OS independent: http://sourceforge.net/projects/uniwakka/ I also dont see a hint that it is Linux dependent in the installation notes. Maybe you think its Linux stuff because of the *.tar.gz ending? You can unzip those files also on Windows, its just a packing file format. For unpacking on a windows machine use either 7zip: http://www.7-zip.org/ or the total commander: http://www.ghisler.com . The first is a free packing program and the second is the best ever windows file manager (if you want to open it there just double click at the archive). Stefan Grosse The unpacked instl. notes: Wakka Installation Not much to it (as long as it works, ahem). Unpack/upload the distribution files into a directory that can be accessed via the web. Then go to the corresponding URL. A web-based installer will walk you through the rest. Example: If your website, say, http://www.mysite.com, is mapped to the directory /home/jdoe/www/, and you place the Wakka distribution files into /home/jdoe/www/wakka/, you should go to http://www.mysite.com/wakka/. Note that Wakka distributions normally unpack into directories that include the version in their name; you'll probably want to rename those to just wakka -- or, if you're on a unixoid system, set up a symbolic link. During first installs, the installer will try to create a file called wakka.config.php in your Wakka directory. In order to do this, you will need to either make the Wakka directory writable by the web server, or create a new (empty) file called wakka.config.php which is writable by the web server. If the installer still fails to create the file, it will dump the file's contents which you can then upload manually. IMPORTANT: for installing or upgrading Wakka, do NOT access any of the files contained in the setup/ subdirectory. They're used by the web-based installer/updater, but you should really just access the Wakka directory itself, and it will (or at least should) work perfectly. Detailed instructions are available at http://www.wakkawiki.com/WakkaInstallation. - Hendrik Mans [EMAIL PROTECTED] __ 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, reproducible code.
[R] exact Wilcoxon signed rank test with ties and the no longer under development exactRanksumTests package
Dear List, after updating the exactRanksumTests package I receive a warning that the package is not developed any further and that one should consider the coin package. I don't find the signed rank test in the coin package, only the Wilcoxon Mann Whitney U-Test. I only found a signed rank test in the stats package (wilcox.test) which is able to calculate the exact pvalues but unfortunately the procedure cannot calculate the exact values with ties. Is there any other package that is providing a similiar test? Or is there an easy work out how to take the ties into account? (Or a chance that the correction is taken into account for the stats package?) Stefan Grosse Take the following example from Bortz/Lienert/Boehnke: x1-c(9,14,8,11,14,10,8,14,12,14,13,9,15,12,9) x2-c(13,15,9,12,16,10,8,13,12,16,9,10,16,12,9) # exactRankTests package: wilcox.exact(x1,x2,paired=TRUE) Exact Wilcoxon signed rank test data: x1 and x2 V = 13, p-value = 0.1367 alternative hypothesis: true mu is not equal to 0 # wilcox.test by stats package: wilcox.test(x1,x2,paired=TRUE,exact=TRUE) Wilcoxon signed rank test with continuity correction data: x1 and x2 V = 13, p-value = 0.1436 alternative hypothesis: true mu is not equal to 0 Warning messages: 1: cannot compute exact p-value with ties in: wilcox.test.default(x1, x2, paired = TRUE, exact = TRUE) 2: cannot compute exact p-value with zeroes in: wilcox.test.default(x1, x2, paired = TRUE, exact = TRUE) __ 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, reproducible code.
[R] xyplot with different symbols and colors?
Dear List, I try to make a xyplot with different colors and symbols, I came this far: library(DAAG) xyplot(csoa~it|sex*agegp,data=tinting,groups=target,pch=list(1,2),auto.key=list(space = right)) this produces a plot with different colors and symbols but unfortunately the legend does only follow the color scheme and not the different symbols. Any suggestions what to change? And how do I change or switch of that background colors in each plots title? Stefan Grosse __ 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, reproducible code.
Re: [R] xyplot with different symbols and colors?
R trellis.par.set(superpose.symbol=list(pch=c(1, 2))) Warning message: Note: The default device has been opened to honour attempt to modify trellis settings in: trellis.par.set(superpose.symbol = list(pch = c(1, 2))) R xyplot(csoa~it|sex*agegp,data=tinting,groups=target,auto.key=list(space=right)) Thanks that did it! __ 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, reproducible code.
Re: [R] Authoring a book
I think Peter Dalgaard is right. Since you are able to use R I believe you will be very fast in learning LaTeX. I think it needs less then a week to learn the most common LaTeX commands. And setting up a wiki and trying then to convert this into a printable document format plus learning the wiki syntax is probably more time consuming. Beside this R is able to work perfectly together with LaTeX, it creates LaTeX output and is doing excellent graphics in the EPS/PS format. The best introduction for LaTeX is the not so short introduction: http://people.ee.ethz.ch/~oetiker/lshort/lshort.pdf If you still are not convinced have a look at UniWakkaWiki: http://uniwakka.sourceforge.net/HomePage It is a Wiki for Science and University purposes and claims to be able to export to Openoffice as well as to LaTeX. Stefan Grosse I have my own server, Windows, based on Apache, PhP, and MySQL. SVN and LaTeX would be my tools of choice. A very different approach. SVN is not something I am aquainted with, but should be worth looking into. As to LaTex, the closest I have worked with is Lyx. The problem is, there are two other authors I have to persuade to learn new tools. So, it might be too complex. Tom __ 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, reproducible code.
Re: [R] installing the x86_64 R Binary on Fedora Core 5
Hi Roger, I dunno what exactly might be the source of that mistake but I would strongly recommend to install R while you are online. Often other packages must be installed for dependencies. (And then I recommend using the smart package manager ( http://labix.org/smart ) which is a great tool and handling dependencies better.) You could use rpm at the command line level. open a shell, type su and give your password, change to the folder where the rpm is downloaded to, type rpm -ivh R-2.3.1-1.fc5.x86_64.rpm will try to install and give you information on whats missing... Stefan Grosse roger bos schrieb: I am looking for help install the x86_64 R Binary onto my FC5 machine. At the risk of subjecting myself to tons of criticism, I must confess that I don't know anything about Linux and I have never compiled R from source. Therefore, I choose FC5 because I see that a 64-bit binary is already available. Here is what I tried: I installed FC5 with all options (productivity, software development, and web server). FC5 boots up fine. I downloaded all the R binary files in that FC5 directory to my USB drive and copied them onto my linux machine (where I don't yet have internet access). I created a folder in my rbos's Home directory and copied the files there. I clicked on the R-2.3.1-1.fc5.x86_64.rpm and right-clicked to choose 'open with install software' It asked me for my root password. I put in the same root password I choose when I installed FC5. I get a installing packages screen that shows the R filename and I click Apply. It then give me an Error: Unable to retrieve software information. Can anyone tell me what steps I am missing? The R install guide states that binary installs are platform specific so it only considers building from the sources. I look forward to learning a lot about Linux and using more than just the GUI, but to get started, I just want to learn how to install a binary of R. Thanks so much, Roger [[alternative HTML version deleted]] __ 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, reproducible code. __ 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, reproducible code.
Re: [R] How to interrupt running computation?
Lothar, Which system do you use? Windows, Linux, Mac? Stefan Lothar Schmid schrieb: If I start a computation in R, how can I interrupt it? I' using R 2.1.0. Thanks for your help, Lothar __ 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, reproducible code. __ 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, reproducible code.
Re: [R] matching pairs in a Dataframe?
I found a solution: pairs-merge(subset(data.frame,cond==2,select=c(group,x),subset(data.frame,cond==3,select=c(group,x)),by=c(group)) That yields something like group cond.x cond.y A 4 6.5 D 5 6 E 4 6 now I just need to figure out how I embed this into a function together with the test, so that I don't have to write to much Stefan Grosse Dear list, I want to extract pairs of values out of a dataframe where one criteria/condition does match. I have an experiment with 3 conditions which were not always applied: e.g.: group cond x A 1 2 A 2 4 A 3 6.5 B 1 3 B 2 4.5 C 1 2.5 C 3 4 D 2 5 D 3 6 E 1 1 E 2 4 E 3 6 Now I wanted to extract the x of those groups where condition 2 and condition 3 do both exist. In this example that would be groups A, D and E and the extracted pairs e.g.: cond2 cond3 4 6.5 5 6 4 6 __ 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, reproducible code.
Re: [R] graphic output file format
What speaks against using the jpeg (or png) device? type: ?jpeg for help Have also a look at the wiki: http://wiki.r-project.org/rwiki/doku.php?id=tips:graphics-base:0savegraphs Stefan Conan Phelan schrieb: I would like to save graphics I produce in jpeg or gif formats. The GDD package sounds like it should let me to do this, but I cannot get it to install (error: Can't find gd.h!). Anyone know what's up? Thnx, C __ 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, reproducible code. __ 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, reproducible code.
Re: [R] graphic output file format
Please mail always also to the list, since it may be that there is someone who could reply better and/or faster. Especially in this case since I am only experienced in M$ Windows and Linux so I am unable to guess whats wrong with OSX. does not even the postscript device work? With postscript I get the best quality for printing. e.g.: postscript(test.eps, width = 8.0, height = 6.0, horizontal = FALSE, onefile =FALSE, paper = special) curve(x^2) dev.off() Conan Phelan schrieb: Thanks for the feedback and link Stefan, Perhaps if I provide a little more background. (Pls forgive both my ignorance and my long windedness.) I am on Mac OSX and am only able to load the Quartz graphic device (is that normal?). From quartz, I've only been able to generate pdf files, which do not match the quartz display to my satisfaction. I've been attempting to find a way to generate other output files for comparison. (I've managed to generate a bmp file that looked awful, after acquiring Gostscript). Anyway, I figure R must able to produce high quality graphic outputs, so I'm trying different things. To this end I was hoping that the GDD package would allow me to load the x11 module but it seems to be a no-go (is there some system requirement I'm missing?). Conan __ 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, reproducible code.