[R] Manova and specifying the model
Hi, I would like to conduct a MANOVA. I know that there 's the manova() funciton and the summary.manova() function to get the appropriate summary of test statistics. I just don't manage to specify my model in the manova() call. How to specify a model with multiple responses and one explanatory factor? If I type: pcor.manova-manova(isol+hcom+habarea+inclin+windprot+shrubcov+herbh+baregr+flowcov~pcor, data=pcor.df) I always get error messages like: Error in manova(isol + hcom + habarea + inclin + windprot + shrubcov + : need multiple response Can someone help me? Thanks Steffi Stefanie von Felten Käppelistrasse 24 4600 Olten Tel: 062/296 13 14 e-mail: [EMAIL PROTECTED] [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] discriminant analysis
Hi, I have done different discriminant function analysis of multivariat data. With the CV=True option I was not able to perform the predict() call. What do I have to do? Or is there no possibility at all? You also need the predicted values to produce a plot of the analysis, as far as I know. Here my code: pcor.lda2-lda(pcor~habarea+hcom+isol+flowcov+herbh+inclin+windprot+shrubcov+baregr, data=pcor.df, CV=T) table2-table(pcor.df$pcor, pcor.lda2$class) table2 #doesn't work, becoause CV=True? pcor.ld2-predict(pcor.lda2, dimen=1)$x plot(pcor.ld2) plot(pcor.lda2, type=density, dimen=1) #kernel density estimates I am happy if I get an answer from somebody! Stefanie von Felten Institut für Umweltwissenschaften Universität Zürich Winterthurerstrasse 190 CH-8057 Zürich e-mail: [EMAIL PROTECTED] telefon: +41 (0)1 635 61 23 felefax: +41 (0)1 635 57 11 [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Manova and specifying the model
See the examples in ?summary.manova, please. Y - cbind(tear, gloss, opacity) shows how to make a multiple response. On Mon, 24 May 2004, Stefanie von Felten wrote: I would like to conduct a MANOVA. I know that there 's the manova() funciton and the summary.manova() function to get the appropriate summary of test statistics. I just don't manage to specify my model in the manova() call. How to specify a model with multiple responses and one explanatory factor? If I type: pcor.manova-manova(isol+hcom+habarea+inclin+windprot+shrubcov+herbh+baregr+flowcov~pcor, data=pcor.df) I always get error messages like: Error in manova(isol + hcom + habarea + inclin + windprot + shrubcov + : need multiple response -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] discriminant analysis
On Mon, 24 May 2004, Stefanie von Felten wrote: I have done different discriminant function analysis of multivariat Using lda in contributed package MASS, uncredited. data. With the CV=True option I was not able to perform the predict() call. What do I have to do? Or is there no possibility at all? You also It makes no sense. You ask for LOO cross-validation, and that is n separate fits, not a single fit from which you can predict. need the predicted values to produce a plot of the analysis, as far as I know. Here my code: pcor.lda2-lda(pcor~habarea+hcom+isol+flowcov+herbh+inclin+windprot+shrubcov+baregr, data=pcor.df, CV=T) table2-table(pcor.df$pcor, pcor.lda2$class) table2 #doesn't work, becoause CV=True? pcor.ld2-predict(pcor.lda2, dimen=1)$x plot(pcor.ld2) plot(pcor.lda2, type=density, dimen=1) #kernel density estimates I am happy if I get an answer from somebody! -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Cannot call R's ISNAN() from a C code in 1.7 versions.
First, there are no versions 1.7, 1.8 and 1.9. Was your version of R compiled against MSVC++ 6.0? The binary on CRAN was not, and binaries for different versions of R were compiled with different versions of MinGW. The entry point isnan is part of the statically linked runtime on modern MinGW. MSVC++ 6.0 does supply _isnan (as it really should as it is part of the C99 ISO standard), and you need to link against it appropriately. Hint: it may have an extra underline, since it seems it is known to C as _isnan. You may need to add #undef ISNAN #define ISNAN(x) _isnan(x) We don't support adding extensions to R using a different compiler to the one used to build R. Changes already made for future releases of R will make this less likely to work in R 2.0.x. On Mon, 24 May 2004, Paul Y. Peng wrote: Dear R users, Have you experienced any difficulty in calling R's ISNAN() from a C code? I have C codes including ISNAN() calls and they worked well until I upgraded my R from 1.7 to later versions. When I tried to compile the codes in the version 1.8 and 1.9, I got error messages like this: test.obj : error LNK2001: unresolved external symbol _isnan .\testR.dll : fatal error LNK1120: 1 unresolved externals NMAKE : fatal error U1077: 'link.exe' : return code '0x460' Stop. I checked Writing R Extensions and did not find any changes in the API entry point for ISNAN in the later versions. Could any one enlighten me on why _isnan cannot be resolved? The compiler I used is MSVC++ 6.0 and the platform is WinXP. Many thanks. -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] bug in cor (..., use= ...)?
Dear R users, I have not found anything on this in the archives. Does anyone know whehther the parameter use= is not functioning in cor or enlighten me what it is supposed to do? My R version is R version 1.8.1, 2003-11-21 on Windows 2000. I am hoping to be able to update to 1.9.1 as soon as it has appeared (we are not allowed here to install software on our own and thus I am trying to be able to have the .1 versions installed ...). Test code: x - 1:10 y - 2:11 x [1] - NA y [10] - 12 cor (x, y, use= 'all.obs', method= 'kendall') cor (x, y, use= 'complete.obs', method= 'kendall') cor (x, y, use= 'pairwise.complete.obs', method= 'kendall') As I understand, the first one of this should result in an error which it does not. All the results are the same and seemingly treat the NA as if it was 0. Any ideas are appreciated. Thanks and regards, Lorenz - Lorenz Gygax, Dr. sc. nat. Tel: +41 (0)52 368 33 84 / [EMAIL PROTECTED] Tag der offenen Tür, 11./12. Juni 2004: http://www.fat.ch/2004 Center for proper housing of ruminants and pigs Swiss Veterinary Office agroscope FAT Tänikon, CH-8356 Ettenhausen / Switzerland Fax : +41 (0)52 365 11 90 / Tel: +41 (0)52 368 31 31 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] bug in cor (..., use= ...)?
From the NEWS file for 1.9.0 o The cor() function did not remove missing values in the non-Pearson case. Your example works correctly there. (I am fairly sure this has been discussed on the mailing lists.) On Mon, 24 May 2004 [EMAIL PROTECTED] wrote: Dear R users, I have not found anything on this in the archives. Does anyone know whehther the parameter use= is not functioning in cor or enlighten me what it is supposed to do? My R version is R version 1.8.1, 2003-11-21 on Windows 2000. I am hoping to be able to update to 1.9.1 as soon as it has appeared (we are not allowed here to install software on our own and thus I am trying to be able to have the .1 versions installed ...). Test code: x - 1:10 y - 2:11 x [1] - NA y [10] - 12 cor (x, y, use= 'all.obs', method= 'kendall') cor (x, y, use= 'complete.obs', method= 'kendall') cor (x, y, use= 'pairwise.complete.obs', method= 'kendall') As I understand, the first one of this should result in an error which it does not. All the results are the same and seemingly treat the NA as if it was 0. -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Tramo-seats
Working - among other things- in the field of (short long term) electricity forecast, we are now using too many too expensive pieces of licensed software: SAS, SPSS, EViews. This sedimentation is due to the fact that my predecessors in the past used different consultant companies to manage each procedure. Having attended the useR2004! Conference with the aim of assessing if R ALONE could glue all those fragmented and isolated procedures, I'm almost convinced now that YES it could do the job! Now - to start with - a first problem to solve: we have to comply with the Tramo-seats closed-source procedure (http://www.bde.es/informes/be/docs/dt0014e.pdf) to deal with seasonality of electricity monthly time-series, in line with the methodology officially adopted by our National Bureau of Statistics. Searching in R-help mailing list I didn't find anything about a tramo-seats R version. Does anyone know of unofficial R translation of tramo-seats? OR What do you suggest? Vittorio __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] problems with starting R
When I try to start R from my desktop, an information-window pops up: Fatal error: Invalid HOMEDRIVE. I already deleted and reinstalled the program, but the problem lasts. How can I solve this problem? I would be very glad, if you can give me a hint, what´s going wrong. Best regards Sabine Bader -- Sabine Bader Hamburger Str. 83 44135 Dortmund 0231/7950683 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
RE: [R] problems with starting R
Hi Sabine, -Original Message- From: Sabine Bader [SMTP:[EMAIL PROTECTED] Sent: Monday, May 24, 2004 2:22 PM To: [EMAIL PROTECTED] Subject: [R] problems with starting R When I try to start R from my desktop, an information-window pops up: Fatal error: Invalid HOMEDRIVE. This problem appeared dozens of times during the last few weeks. Best thing is probably to read one of the previous threads related to this topic. If you go, for example, to: http://maths.newcastle.edu.au/~rking/R/ there is a searchable R-help archive. Enter the keywords: INVALID HOMEDRIVE and you will get many many results (I just tried it and it gave 128 hits). Hope this helps, Roland + This mail has been sent through the MPI for Demographic Rese...{{dropped}} __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] problems with starting R
Hi! This problem was discussed on the list. For searching the mail archives you can use eg. http://maths.newcastle.edu.au/~rking/R/ The link to this page can be found at: http://cran.r-project.org/search.html Entering Fatal error: Invalid HOMEDRIVE. Will provide you many answers to your question. Eg.: http://tolstoy.newcastle.edu.au/R/help/04/04/1171.html Sincerely Eryk *** REPLY SEPARATOR *** On 5/24/2004 at 2:22 PM Sabine Bader wrote: When I try to start R from my desktop, an information-window pops up: Fatal error: Invalid HOMEDRIVE. I already deleted and reinstalled the program, but the problem lasts. How can I solve this problem? I would be very glad, if you can give me a hint, what´s going wrong. Best regards Sabine Bader -- Sabine Bader Hamburger Str. 83 44135 Dortmund 0231/7950683 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Dipl. bio-chem. Eryk Witold Wolski@MPI-Moleculare Genetic Ihnestrasse 63-73 14195 Berlin 'v' tel: 0049-30-83875219 / \ mail: [EMAIL PROTECTED]---W-Whttp://www.molgen.mpg.de/~wolski __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] non-hierarchical non-exclusive clustering of large data sets
Hi, I'm trying to use R to cluster words with related meanings. Does anyone know of a non-hierarchical clustering method in R that produces non-exclusive clusters? With non-exclusive, I mean that words should be allowed to be part of multiple clusters. So my data matrix would look something like: T1 T2 T3 CLOWN_N 0 1 0 BANK_N 3 0 2 RIVER_N 0 0 2 FLOW_V 0 0 3 MONEY_N 2 0 0 PAY_V 2 0 0 The first line indicates the noun clown occurred only once in my text collection, namely in text 2. Ideally, the clustering method would produce the clusters [bank_n,river_n,flow_v], [bank_n,money_n,pay_v] and [clown_n]. The data matrix I would use would be much bigger than the one above, its dimensions would be in the order of (10,10). Does anyone know if this would cause practical problems, perhaps very slow clustering? Best wishes, Murk Wuite, MA student Department of Language and Speech Katholieke Universiteit Nijmegen, The Netherlands __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] replacing backslashes with slashes using gsub
Hi! I am trying to replace backslashes with slashes using gsub (R1.9.0 on XP) gsub(,/,D:\Prog\R\rw1090\library\cluster\libs) [1] D:ProgR\rw1090libraryclusterlibs ? Sincerely Eryk __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Month names
This is how I get the month names from within R: mon - rep(strptime(01/01/1952, format = %d/%m/%Y), 12) mon$mon - mon$mon + 0:11 mnam - months(mon, abbreviate = F) mnam [1] januarfebruar marts april maj juni juli augustseptember [10] oktober november december Surely someone on the list can beat that in elegance? Bendix Carstensen -- Bendix Carstensen Senior Statistician Steno Diabetes Center Niels Steensens Vej 2 DK-2820 Gentofte Denmark tel: +45 44 43 87 38 mob: +45 30 75 87 38 fax: +45 44 43 07 06 [EMAIL PROTECTED] www.biostat.ku.dk/~bxc __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] replacing backslashes with slashes using gsub
gsub(,/,D:\Prog\R\rw1090\library\cluster\libs) [1] D:ProgR\rw1090libraryclusterlibs Probably not the best way, but what about escaping all the backslashes in the original string? gsub(,/,D:\\Prog\\R\\rw1090\\library\\cluster\\libs) __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] (no subject)
Hello! Please!How do I download R from the internet? http://cran.r-project.org/doc/FAQ/R-FAQ.html#How%20can%20R%20be%20obtained%3f __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] (no subject)
Whittaker, Carly wrote: Hello! Please!How do I download R from the internet? Several ways, one ist to right click and say something like save to ... in your browser when visiting CRAN at http://cran.r-project.org/ and having browsed to the sources or the binary version that is the right one for your OS. Uwe Ligges Thankyou for your time carly [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Stopping the process after a certain time
Greetings R-community, I'm running simulations within R that I wrote in C. The simulations require fitting that occasionally fails to finish. I was wondering if there is any kind of tool for process control in R, such that after e.g. 15 minutes I could kill the process, record the state for post-hoc analysis, and move to the next simulation? I'm running FreeBSD and could almost surely do something in Perl but I'd rather stay inside R if possible. Thanks, Andrew Andrew Robinson Ph: 208 885 7115 Department of Forest Resources Fa: 208 885 6226 University of Idaho E : [EMAIL PROTECTED] PO Box 441133W : http://www.uidaho.edu/~andrewr Moscow ID 83843 Or: http://www.biometrics.uidaho.edu No statement above necessarily represents my employer's opinion. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Month names
How about `month.name'? -roger BXC (Bendix Carstensen) wrote: This is how I get the month names from within R: mon - rep(strptime(01/01/1952, format = %d/%m/%Y), 12) mon$mon - mon$mon + 0:11 mnam - months(mon, abbreviate = F) mnam [1] januarfebruar marts april maj juni juli augustseptember [10] oktober november december Surely someone on the list can beat that in elegance? Bendix Carstensen -- Bendix Carstensen Senior Statistician Steno Diabetes Center Niels Steensens Vej 2 DK-2820 Gentofte Denmark tel: +45 44 43 87 38 mob: +45 30 75 87 38 fax: +45 44 43 07 06 [EMAIL PROTECTED] www.biostat.ku.dk/~bxc __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Metafiiles into Word R 1.9.0
Dear R-helpers, I recently upgraded to R 1.9.0 in my computer at work and at home: 1.-The computer at home has Windows XP and Office XP and it seems to work perfectly and I copy-paste graphics perfectly. 2.-The computer at work has Microsoft Windows 2000 (5.00.2195 Service Pack 2) and Word 2000 (9.0.4402 SR-1) I cannot copy-paste windows metafiles into Office applications in my computer at work. The resulting object is empty. Is there any way of overcoming this other than going back in R versions? (an indication to patch from Microsoft Office that can take new metafiles :-). best regards, Jesus P.S.: this mail is probably a similar topic as Patrick Giraudoux email at the beginning of the month, I apologise for the duplication. platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major1 minor9.0 year 2004 month04 day 12 language R -- Jesús María Frías Celayeta School of Food Sci. and Env. Health. Faculty of Tourism and Food Dublin Institute of Technology Cathal Brugha St., Dublin 1. Ireland t +353 1 4024459 f +353 1 4024495 w www.dit.ie/DIT/tourismfood/science/staff/frias.html -- -- This message has been scanned for content and viruses by the DIT Information Services MailScanner Service, and is believed to be clean. http://www.dit.ie __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Month names
Roger == Roger D Peng [EMAIL PROTECTED] on Mon, 24 May 2004 10:38:09 -0400 writes: Roger How about `month.name'? English only. Note that he got the names in Danish (I think) Martin Roger BXC (Bendix Carstensen) wrote: This is how I get the month names from within R: mon - rep(strptime(01/01/1952, format = %d/%m/%Y), 12) mon$mon - mon$mon + 0:11 mnam - months(mon, abbreviate = F) mnam [1] januarfebruar marts april maj juni juli augustseptember [10] oktober november december Surely someone on the list can beat that in elegance? Bendix Carstensen -- Bendix Carstensen Senior Statistician Steno Diabetes Center Niels Steensens Vej 2 DK-2820 Gentofte Denmark tel: +45 44 43 87 38 mob: +45 30 75 87 38 fax: +45 44 43 07 06 [EMAIL PROTECTED] www.biostat.ku.dk/~bxc __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] RMySQL problem
Hi, I'm using R 1.9.0 with RMySQL 0.5-4 and MySQL 3.23.55 on a suse 8.2 box. I have a simulation study and (as usual for newbies in simulation, I guess) I have a lot of data that I want to store in MySQL. I want to write an R script that reads data from RData files and writes it to a MySQL database. I read some R documents (R Data Import/Export and DSC papers) but I'm finding differences between the documents and the packages (RMySQL and DBI). I don't find the methods to write data like dbWriteTable... On the other hand I've tryied to build a sql statement to insert data but I'm stucked because of ... who knows, my ignorance probably. I want to take advantage of MySQl INSERT INTO statement that deals with several rows at once to insert a complete data.frame into a table. I've tryied to use paste to build the sql string but It works by column and I need it by row ... The sql systax should be something like: INSERT INTO TABLEA(COL1, COL2, COL3) VALUES (VAL11, VAL12, VAL13), (VAL21, VAL22, VAL23), ... (VALN1, VALN2, VALN3); and I have a data.frame with 3 columns corresponding to that table columns. How can I do this ? Thanks EJ __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] RPM post-install scripts to update R
chris albert [EMAIL PROTECTED] writes: Hi, Paul Johnson asked: While I'm on the RPM subject, can I ask an R RPM packaging question? I want the R RPM to install so that post install then R starts and runs update.packages() as well as install.packages(c(Design,Hmisc,lmtest,car,rgl,effects)) well, you get the idea. I want to add in more packages, of course. Can I ask what might be the best way to package this? You can add scripts to the %post macro section. ...etc... Erm, but would you really want the result of installing an RPM to depend on the install date and the previous configuration of the system? I can appreciate Paul's desire to automate, but when designing an RPM, you also have to consider that you'd be doing things that the user may not even know what mean. That sort of thing easily becomes a support nightmare. -- O__ Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Re: Windows versus Unix packages in CRAN ...
On Mon, 24 May 2004, Peter Dalgaard wrote: Whereas we cannot reasonably require people to test code on platforms that they haven't got, we can and should require them to adhere to reasonable standards and test procedures (which, mind you, other people have invested a serious amount of time in working out). Well, yes, for CRAN. That's presumably one reason why the packages in question aren't on CRAN. -thomas __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Stopping the process after a certain time
On Mon, 24 May 2004, Andrew Robinson wrote: Greetings R-community, I'm running simulations within R that I wrote in C. The simulations require fitting that occasionally fails to finish. I was wondering if there is any kind of tool for process control in R, such that after e.g. 15 minutes I could kill the process, record the state for post-hoc analysis, and move to the next simulation? I'm running FreeBSD and could almost surely do something in Perl but I'd rather stay inside R if possible. As one component of this, if you send an R process SIGUSR1 it will quit and save .RData (approximately the equivalent of CTRL-C, q(yes)). -thomas __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Month names
Martin Maechler [EMAIL PROTECTED] writes: Roger == Roger D Peng [EMAIL PROTECTED] on Mon, 24 May 2004 10:38:09 -0400 writes: Roger How about `month.name'? English only. Note that he got the names in Danish (I think) Martin Yep. Notice, however, that even with Brian's solution, you do need to set the locale first: format(ISOdate(2004,1:12,1),%B) [1] January February March April May June [7] July AugustSeptember October November December Sys.setlocale(LC_TIME,da_DK) ;format(ISOdate(2004,1:12,1),%B) [1] da_DK [1] januarfebruar marts april maj juni [7] juli augustseptember oktober november december -- O__ Peter Dalgaard Blegdamsvej 3 c/ /'_ --- Dept. of Biostatistics 2200 Cph. N (*) \(*) -- University of Copenhagen Denmark Ph: (+45) 35327918 ~~ - ([EMAIL PROTECTED]) FAX: (+45) 35327907 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] as.matrix.data.frame() in R 1.9.0 converts to character when it should (?) convert to numeric
I don't think a POSIXt element *is* numeric (that's a basic atomic vector), so the new behaviour seems right to me. The Warning is wrong, though, and will be fixed. On Mon, 24 May 2004, Don MacQueen wrote: Conversion of a data frame to a matrix using as.matrix() when a column of the data frame is POSIXt and all other columns are numeric has changed in R 1.9.0 from R 1.8.1. The new behavior issues a warning message and converts to a character matrix. In R 1.8.1, such an object was converted to a numeric matrix. Here is an example. R 1.9.0 foo - data.frame( x=1:3,dt=ISOdatetime(2003,1,1:3,0,0,0)) as.matrix(foo) x dt 1 1 2003-01-01 2 2 2003-01-02 3 3 2003-01-03 Warning message: longer object length is not a multiple of shorter object length in: cl == c(Date, POSIXct, POSIXlt) version _ platform sparc-sun-solaris2.8 arch sparc os solaris2.8 system sparc, solaris2.8 status Patched major1 minor9.0 year 2004 month04 day 30 language R ### R 1.8.1 foo - data.frame( x=1:3,dt=ISOdatetime(2003,1,1:3,0,0,0)) foo x dt 1 1 2003-01-01 2 2 2003-01-02 3 3 2003-01-03 as.matrix(foo) x dt 1 1 1041408000 2 2 1041494400 3 3 1041580800 version _ platform sparc-sun-solaris2.8 arch sparc os solaris2.8 system sparc, solaris2.8 status Patched major1 minor8.1 year 2003 month12 day 03 language R In both versions: class(foo$dt) [1] POSIXt POSIXct In R 1.8.1, as.matrix.data.frame() has these lines: if (length(levels(xj)) 0 || !(is.numeric(xj) || is.complex(xj)) || (!is.null(cl - attr(xj, class)) any(cl == c(POSIXct, POSIXlt In R 1.9.0 there is instead if (length(levels(xj)) 0 || !(is.numeric(xj) || is.complex(xj)) || (!is.null(cl - attr(xj, class)) any(cl == c(Date, POSIXct, POSIXlt And that, I think, explains the warning message. From ?as.matrix() in R 1.9.0: 'as.matrix' is a generic function. The method for data frames will convert any non-numeric/complex column into a character vector using 'format' and so return a character matrix, except that all-logical data frames will be coerced to a logical matrix. The POSIXt element is numeric, and so should be converted to numeric is.numeric(foo$dt) [1] TRUE I think this might qualify for bug status, either in and of itself or relative to documentation. But I'm not, as the posting guide says, completely and utterly sure. So I'm posting to r-help first...I will send a bug report if an R-core member asks me to. Thanks -Don -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] R 1.8.1 - 1.9.0 incompatability: Underscore in syntactically valid names
Peter == Peter Dalgaard [EMAIL PROTECTED] writes: Gah! I could swear we discussed that particular issue leading up to 1.9.x and had plans for a compatibility option. You might file a bug report at least for the docs, since the example is clearly wrong... Done. I tried to write some version dependencies into a sample code, but I'm stumped by the fact that _ is not allowed before 1.9.0. For example, suppose I have a data file, example.dat, like this: a b x some_factor 1 1 0.4 orange 2 1 0.3 blue 1 1 0.2 dog 2 1 0.1 orange 1 2 0.4 blue 2 2 0.3 dog 1 2 0.2 orange 2 2 0.1 blue To read and use this in a version independent way, I've tried this: df - read.table('example.dat',header=T) if ( version['minor'] == 9.0 ) { plot(x ~ some_factor, data=df) } else { plot(x ~ some.factor, data=df) } This fails in R 1.8.1, because some_factor throws a syntax error: df - read.table('example.dat',header=T) if ( version['minor'] == 9.0 ) { + if ( version['minor'] == 9.0 ) { + plot(x ~ some_factor, data=df) Error: syntax error Ick. Is there a known idiom for handling this sort of version dependency in R? I'm going to avoid R 1.9.x for now, and I encourage any authors of contributing packages to do their best to maintain backwards compatibility for those of us who cannot make the switch quickly. Mike -- Michael A. Miller [EMAIL PROTECTED] Imaging Sciences, Department of Radiology, IU School of Medicine __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Stopping the process after a certain time
Thomas, that's very interesting - thanks! That will be helpful. Andrew On Mon, 24 May 2004, Thomas Lumley wrote: On Mon, 24 May 2004, Andrew Robinson wrote: Greetings R-community, I'm running simulations within R that I wrote in C. The simulations require fitting that occasionally fails to finish. I was wondering if there is any kind of tool for process control in R, such that after e.g. 15 minutes I could kill the process, record the state for post-hoc analysis, and move to the next simulation? I'm running FreeBSD and could almost surely do something in Perl but I'd rather stay inside R if possible. As one component of this, if you send an R process SIGUSR1 it will quit and save .RData (approximately the equivalent of CTRL-C, q(yes)). -thomas __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Andrew Robinson Ph: 208 885 7115 Department of Forest Resources Fa: 208 885 6226 University of Idaho E : [EMAIL PROTECTED] PO Box 441133W : http://www.uidaho.edu/~andrewr Moscow ID 83843 Or: http://www.biometrics.uidaho.edu No statement above necessarily represents my employer's opinion. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Month names
As far as I know the original solution uses the current locale, and mine is just a much-simplified version of the same underlying call. On 24 May 2004, Peter Dalgaard wrote: Martin Maechler [EMAIL PROTECTED] writes: Roger == Roger D Peng [EMAIL PROTECTED] on Mon, 24 May 2004 10:38:09 -0400 writes: Roger How about `month.name'? English only. Note that he got the names in Danish (I think) Martin Yep. Notice, however, that even with Brian's solution, you do need to set the locale first: Just as you did with the original code, which gives English names on my machine. format(ISOdate(2004,1:12,1),%B) [1] January February March April May June [7] July AugustSeptember October November December Sys.setlocale(LC_TIME,da_DK) ;format(ISOdate(2004,1:12,1),%B) [1] da_DK [1] januarfebruar marts april maj juni [7] juli augustseptember oktober november december -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Mac OS X jpg
I think the message is very clear. Try ?BATCH Using 'R CMD BATCH' sets the GUI to 'none', so none of 'x11', 'jpeg' and 'png' are available. and ?jpeg R can be compiled without support for either or both of these devices: this will be reported if you attempt to use them on a system where they are not supported. They will not be available if R has been started with '--gui=none' (and will give a different error message), and they may not be usable unless the X11 display is available to the owner of the R process. The bitmap() device will work under R CMD BATCH. On Mon, 24 May 2004, Rafael Najmanovich wrote: I am using R 1.9.0 (2004-04-12) for Mac OS X. I am having trouble printing to a jpg file. This only happens when I use R in batch (with the BATCH option) mode. If I launch the R GUI or the command line version and run the script from either, I have no problems. Only when running in BATCH mode I get a file with the same name as my script where I find the following message at the end: Error in jpeg(filename = file, width = 1500, height = 2000, pointsize = 12, : R_X11 module cannot be loaded In addition: Warning message: X11 module is not available under this GUI Execution halted When I run the same script with the only difference that it prints to a postscript file instead, all works fine. Any suggestions as to what could the problem be? Failure to consult the documentation, I believe. -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Cannot call R's ISNAN() from a C code in 1.7 versions.
On Mon, 24 May 2004, Paul Y. Peng wrote: Prof Brian Ripley wrote: First, there are no versions 1.7, 1.8 and 1.9. Sorry for my misuse of the version numbers. Was your version of R compiled against MSVC++ 6.0? The binary on CRAN was not, and binaries for different versions of R were compiled with different versions of MinGW. The entry point isnan is part of the statically linked runtime on modern MinGW. I used the binary on CRAN. Then please use the compilers described in readme.packages which match it. MSVC++ 6.0 does supply _isnan (as it really should as it is part of the C99 ISO standard), and you need to link against it appropriately. Hint: it may have an extra underline, since it seems it is known to C as _isnan. You may need to add #undef ISNAN #define ISNAN(x) _isnan(x) Many thanks for this suggestion. It works, as always. The extra underline is required. We don't support adding extensions to R using a different compiler to the one used to build R. Changes already made for future releases of R will make this less likely to work in R 2.0.x. I wish that the R API entry points documented in Writing R Extensions be supported in the future versions of R, because it will make programs built with R more portable than directly using compiler-specific functions, such as testing and generating the IEEE 754 special values. The existing entry points have saved me trouble to chase these values whenever I moved to an environment with a different compiler. Who said they would not supported? What I said is that using a different compilers to compile different parts of R is not supported (and never has been). And that is even less likely to work in future releases. All we guarantee is that ISNAN and R_FINITE produce calls to functions that work under the system used to configure R. As on Windows R_FINITE now produces a call to a MinGW macro, that is not going to work under VC++ 6. We are not talking about using `compiler-specific functions' here, like _isnan. What R core is doing is using *standard* C99 functions where available, and a problem with VC++ 6.0 is that it is far from compliant with IEC60559 aka IEEE754. Given that the recommended compiler is freely available, it is unreasonable to expect any support for less capable compilers. -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UKFax: +44 1865 272595 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] adress the current index
How can I adress the current index of a vector? I want to work with time series and therefore give the n-th element of a vector some value dependent on the value of the n-1th element. Thank you in advance. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] what does numeric mean? {was ... as.matrix.data.frame() ...}
BDR == Prof Brian Ripley [EMAIL PROTECTED] on Mon, 24 May 2004 17:02:21 +0100 (BST) writes: BDR I don't think a POSIXt element *is* numeric (that's a BDR basic atomic vector), well, yes, in some strict sense, but at least is.atomic( ISOdatetime(2003,1,1:3,0,0,0)) [1] TRUE BDR so the new behaviour seems right to me. The Warning is BDR wrong, though, and will be fixed. Further note that data.matrix(.) still does convert to numeric, and AFAIK all good books on S / R tell you to rather use data.matrix(df) in situations you want a numeric matrix from a data frame df. BDR On Mon, 24 May 2004, Don MacQueen wrote: Conversion of a data frame to a matrix using as.matrix() when a column of the data frame is POSIXt and all other columns are numeric has changed in R 1.9.0 from R 1.8.1. The new behavior issues a warning message and converts to a character matrix. In R 1.8.1, such an object was converted to a numeric matrix. Here is an example. R 1.9.0 foo - data.frame( x=1:3,dt=ISOdatetime(2003,1,1:3,0,0,0)) as.matrix(foo) x dt 1 1 2003-01-01 2 2 2003-01-02 3 3 2003-01-03 Warning message: longer object length is not a multiple of shorter object length in: cl == c(Date, POSIXct, POSIXlt) From ?as.matrix() in R 1.9.0: 'as.matrix' is a generic function. The method for data frames will convert any non-numeric/complex column into a character vector using 'format' and so return a character matrix, except that all-logical data frames will be coerced to a logical matrix. The POSIXt element is numeric, and so should be converted to numeric is.numeric(foo$dt) [1] TRUE that is a point, particularly together with the above is.atomic(foo$dt) |- TRUE I think this might qualify for bug status, either in and of itself or relative to documentation. But I'm not, as the posting guide says, completely and utterly sure. So I'm posting to r-help first...I will send a bug report if an R-core member asks me to. Thanks -Don Note that with a factor (instead of POSIX*t), things are similar str(ffoo - data.frame(x=1:3, f=gl(3,1))) `data.frame': 3 obs. of 2 variables: $ x: int 1 2 3 $ f: Factor w/ 3 levels 1,2,3: 1 2 3 as.matrix(ffoo) x f 1 1 1 2 2 2 3 3 3 data.matrix(ffoo) x f 1 1 1 2 2 2 3 3 3 but in that case, is.numeric(ffoo$f) [1] FALSE hence, all according to the docs -- Maybe we should either be more specific in help(as.matrix) {and more other places ?} about what numeric means there, or consider -- and this may well be unfeasible because it potentially breaks current code -- to redefine 'is.numeric(.)' to be more restrictive: e.g., if (is.object(x)) === ! is.numeric(x) ? Martin __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] barplot
I´ve tried version 1.9.0 barplot with these (and others) example from the help page: tN - table(Ni - rpois(100, lambda=5)) r - barplot(tN, col='gray') I get : ...OLE_Obj... Same example with version 1.8.1 gives the following result: ...OLE_Obj... What is wrong with v.1.9.0? Thanks, Carlos Guevel __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] identify() in script
hallo, I want a script of mine to run identify(), and possibly switch to interactive mode at some point. Is it possible? Thanks in advance __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Applications of the Distance from a Point to a Curve
I would like to know if anyone knows about (real) situations where it is necessary to compute the euclidean distance from a point to a curve (or even the point where it is reached), you may not assume that the point is very close to the curve (unlike the usual fitting problems). I have special interest in applications to phisics. Thanks in advance, Dimas __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] barplot
Carlos Guevel wrote: I´ve tried version 1.9.0 barplot with these (and others) example from the help page: tN - table(Ni - rpois(100, lambda=5)) r - barplot(tN, col='gray') I get : ...OLE_Obj... Same example with version 1.8.1 gives the following result: ...OLE_Obj... What is wrong with v.1.9.0? Thanks, Carlos Guevel This has been reported a while back: http://r-bugs.biostat.ku.dk/cgi-bin/R/Graphics?id=6777;expression=barplot;user=guest which also suggests a workaround. tN - table(Ni - rpois(100, lambda=5)) r - barplot(as.vector(tN), col='gray') Alternatively, barplot2 in the gregmisc package should also work. library(gregmisc) tN - table(Ni - rpois(100, lambda=5)) r - barplot2(tN, col='gray') --sundar __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Tramo-seats
On Monday 24 May 2004 11:00, [EMAIL PROTECTED] wrote: Working - among other things- in the field of (short long term) electricity forecast, we are now using too many too expensive pieces of licensed software: SAS, SPSS, EViews. This sedimentation is due to the fact that my predecessors in the past used different consultant companies to manage each procedure. Having attended the useR2004! Conference with the aim of assessing if R ALONE could glue all those fragmented and isolated procedures, I'm almost convinced now that YES it could do the job! Now - to start with - a first problem to solve: we have to comply with the Tramo-seats closed-source procedure (http://www.bde.es/informes/be/docs/dt0014e.pdf) to deal with seasonality of electricity monthly time-series, in line with the methodology officially adopted by our National Bureau of Statistics. Searching in R-help mailing list I didn't find anything about a tramo-seats R version. Does anyone know of unofficial R translation of tramo-seats? OR What do you suggest? Vittorio __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Not an R response I'm afraid, but have you looked at Demetra http://forum.europa.eu.int/irc/dsis/eurosam/info/data/demetra.htm which appears to be an officially sponsored windows-only implementation of Tramo-seats and X-12 ARIMA. The licence is not clear from the website, but it looks free as in beer but probably not open source. Paul Bivand __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] installing R on Fedora Core 2 test 2
Hello: Please give more info, how did you proceed?, when you install/upgrade Fedora what kind on instalation did you made? A few days ago I upgrade my computer from FC-1 to FC-2 (workstation). After that, I install the binaries (from CRAN) and everything went just fine. On Mon, 24 May 2004, Benjamin Yakir wrote: I am new to Linux and just installed Fedora Core 2. Tried to install R from Fedora core 1 binaries and from the source files according to the manual but failed. Would appreciate any help. Benny Yakir Department of Statistics The Hebrew University [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- Ulises M. Alvarez LAB. DE ONDAS DE CHOQUE FISICA APLICADA Y TECNOLOGIA AVANZADA UNAM [EMAIL PROTECTED] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Seasonal ARIMA question - stat package (formerly ts)
To whom it may concern: I am trying to better understand the functionality of 'R' when making arima predictions to avoid any Black Box disadvantages. I'm fitting a seasonal arima model using the following command (having already loaded 'stat' package). arimaSeason - arima(Data,order=c(1,0,1),seasonal=list(order=c(1,0,1),period=12)) I can then generate subsequent predictions using the 'predict' function. However, I can't seem to duplicate these predictions in a separate program using the model coefficients. From duplicating simpler models, I understand the input variables are adjusted by the intercept term. (e.g. for an arima(1,0,0) the prediction equation is y(t) = beta1 * ( y(t-1) - beta0 ) + beta0 where beta0 is the intercept) Currently, I've expected the prediction to follow the equation below: y(t) = beta0 + [beta1*( y(t-1) - beta0 )] + [beta2 * epsilon(t-1)] + [beta3 * (y(t-12) - beta0)] + [beta4 * epsilon(t-12)] This has proved unsuccessful. What equation underlies this arima prediction? Is there something different that happens Many thanks for your help, Keith Campbell Keith Campbell Researcher ROTELLA CAPITAL MANAGEMENT |180 N. Stetson Suite 5100 I Chicago, IL 60601 | t: 312.706.0442 | f: 312.861.0468 | _ This message is for the designated recipient only and may contain privileged, proprietary, or otherwise private information. If you are not the intended recipient or otherwise believe that you have received this message in error, please notify the sender immediately and delete the original. Any other use of this message by you is prohibited. [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Metafiiles into Word R 1.9.0
2.-The computer at work has Microsoft Windows 2000 (5.00.2195 Service Pack 2) and Word 2000 (9.0.4402 SR-1) I cannot copy-paste windows metafiles into Office applications in my computer at work. The resulting object is empty. Is it genuinely empty, or can you see the graphs in print layout view, or print preview? I've seen that particular Word bug before. It happens with Visio drawings too (Visio is a M$ product). Cheers Jason __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] installing R on Fedora Core 2 test 2
Dear Prof Ripley Thank you very much. Half the night is over but running the patch worked. Benny On Tue, 2004-05-25 at 00:02, Prof Brian Ripley wrote: What was the failure? For the sources, there is a bug in the X11 headers it ships. Please use the R-patched tarball from ftp://ftp.stat.math.ethz.ch/Software/R/, which has a workaround. On Mon, 24 May 2004, Benjamin Yakir wrote: I am new to Linux and just installed Fedora Core 2. Tried to install R from Fedora core 1 binaries and from the source files according to the manual but failed. Would appreciate any help. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Re: identify() in script
On Mon, 2004-05-24 at 19:43, Diogo FC Patrao wrote: On Mon, 2004-05-24 at 17:17, Diogo FC Patrao wrote: hallo, I want a script of mine to run identify(), and possibly switch to interactive mode at some point. Is it possible? Thanks in advance I think I didn`t explained my point properly. Consider I have a file named, say, test.R, which contains x=rnorm(100) y=rnorm(100) plot( x,y ) identify(x,y,1:100) and want to invoke it as a batch, running (in linux) $ R test.R just by doing this, no plot appear. That's my doubt: is it possible to do this? or, there is a command which halts the execution at some point and allow me to enter commands, as I've invoked R in interactive mode? thanks for your patience, dfcp __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] yeardays
Dear All, Is there a R (Windows) equivalent of yeardays in Splus 6.0 for Windows? I am using XP. ** Steve Su ([EMAIL PROTECTED]) PhD student. School of Accountancy Queensland University of Technology Postal Address: Steve Su, School of Accountancy, QUT, PO Box 2434, Brisbane, Queensland, Australia, 4001. Phone: +61 7 3864 2017 Fax:+61 7 3864 1812 Mobile: 0421 840 586 . _--_|\ / QUT \_.--._/ v ** [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Re: Windows versus Unix packages in CRAN ...
On 24 May 2004 17:18:52 +0200, Peter Dalgaard [EMAIL PROTECTED] wrote: Duncan Murdoch [EMAIL PROTECTED] writes: On Fri, 21 May 2004 15:41:18 +0200, Tamas Papp [EMAIL PROTECTED] wrote: You are probably right in saying that they _could_ have done better, but I would not use should in this context. AFAIK the package is free (as in beer) software, which means that you are not paying for it. The maintainers probably do not need a Linux version (yet), so it was not easy to use it under Linux. Feel free to contribute. Yes, indeed! Well, yes and no. Yes, people should feel free to help out with the maintenance, but no, it is not reasonable to leave cross-platform issues unaddressed. There's a list of packages on CRAN in /bin/windows/contrib/1.9/@ReadMe that fail to build on Windows for one reason or another. Should we critcize the authors of those packages for not addressing cross-platform issues? I don't think so. Why should it be any different when the situation is reversed? Duncan Murdoch __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] loess
Hi, When I was running the function loess(y~x, span=0.0020), I got a warning message k-d tree limited by memory. ncmax= 4231 Does that mean the function has not been computed correctly ? If it has not, is there any way to adjust it so that it will do correctly ? Thanks. Rupen. *** If you fail to plan, you are planning to fail. _ Download music tracks from 95c here: __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Scheire-Ray-Hare extension to Kruskal-Wallis
Hi, I was wondering if it was possible to do the Scheirer-Ray-Hare extension of the Kruskal-Wallis test using R, or if it was possible to use some other non-parametric test for the following situation: I have continuous data that is grouped by 2 factor variables. I would like a non-parametric test because the group variances are not equal. I would like to look for the effects of either factor variable and their interaction. Thanks for any help, Sarah [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] e1071, R1.9.0, Solaris 2.9, should I be worried?
In R 1.9.0 running under Solaris 2.9 on a SunBlade 100, with Sun WorkShop 6 update 2 C++ 5.3 2001/05/15 as the C++ compiler, I just did install.packages(e1071) The output includes these lines, which I have wrapped to fit nicely in mail: ** libs cc -I/users/local/lib/R/include -I/usr/local/include -KPIC -xlibmil \ -dalign -xO4 -c cmeans.c -o cmeans.o cc -I/users/local/lib/R/include -I/usr/local/include -KPIC -xlibmil \ -dalign -xO4 -c cshell.c -o cshell.o cc -I/users/local/lib/R/include -I/usr/local/include -KPIC -xlibmil \ -dalign -xO4 -c floyd.c -o floyd.o cc -I/users/local/lib/R/include -I/usr/local/include -KPIC -xlibmil \ -dalign -xO4 -c Rsvm.c -o Rsvm.o CC -I/users/local/lib/R/include -I/usr/local/include -KPIC -xlibmil \ -dalign -xO4 -c svm.cpp -o svm.o svm.cpp, line 444: Warning: l hides Solver::l. svm.cpp, line 444: Warning: Q hides Solver::Q. svm.cpp, line 444: Warning: Cp hides Solver::Cp. svm.cpp, line 444: Warning: Cn hides Solver::Cn. svm.cpp, line 444: Warning: eps hides Solver::eps. svm.cpp, line 507: Warning: String literal converted to char* in formal argument fmt in call to info(char*, ...). svm.cpp, line 517: Warning: String literal converted to char* in formal argument fmt in call to info(char*, ...). svm.cpp, line 690: Warning: String literal converted to char* in formal argument fmt in call to info(char*, ...). svm.cpp, line 881: Warning: l hides Solver::l. svm.cpp, line 881: Warning: Q hides Solver::Q. svm.cpp, line 881: Warning: b hides Solver::b. svm.cpp, line 881: Warning: y hides Solver::y. svm.cpp, line 881: Warning: alpha hides Solver::alpha. svm.cpp, line 881: Warning: Cp hides Solver::Cp. svm.cpp, line 881: Warning: Cn hides Solver::Cn. svm.cpp, line 881: Warning: eps hides Solver::eps. svm.cpp, line 881: Warning: si hides Solver_NU::si. svm.cpp, line 1278: Warning: String literal converted to char* in formal argument fmt in call to info(char*, ...). svm.cpp, line 1328: Warning: String literal converted to char* in formal argument fmt in call to info(char*, ...). svm.cpp, line 1404: Warning: String literal converted to char* in formal argument fmt in call to info(char*, ...). svm.cpp, line 1439: Warning: String literal converted to char* in formal argument fmt in call to info(char*, ...). svm.cpp, line 1483: Warning: String literal converted to char* in formal argument fmt in call to info(char*, ...). svm.cpp, line 1507: Warning: String literal converted to char* in formal argument fmt in call to info(char*, ...). svm.cpp, line 1719: Warning: String literal converted to char* in formal argument fmt in call to info(char*, ...). 24 Warning(s) detected. How worried should I be? I guess the Warning: x hides Solver::x warnings related to a deliberate style choice, but what about the String literal converted to char* ones? __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html