[R] Help in installing rggobi in ubuntu linux
Hi R users. I am experimenting with ubuntu 7.04 Feisty. I install the ggobi package with apt-get. I got almost all the packages, but when I try to obtain rggobi, I got this message: - install.packages("rggobi") Aviso en install.packages("rggobi") : argument 'lib' is missing: using '/usr/local/lib/R/site-library' --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk interface ... done probando la URL 'http://cran.at.r-project.org/src/contrib/rggobi_2.1.4-4.tar.gz' Content type 'application/x-gzip' length 401451 bytes URL abierta == downloaded 392Kb * Installing *source* package 'rggobi' ... checking for pkg-config... /usr/bin/pkg-config checking pkg-config is at least version 0.9.0... yes checking for GGOBI... configure: creating ./config.status config.status: creating src/Makevars ** libs gcc -std=gnu99 -I/usr/share/R/include -I/usr/share/R/include -g -DUSE_EXT_PTR=1 -D_R_=1 -fpic -g -O2 -c brush.c -o brush.o En el fichero incluÃdo de brush.c:1: RSGGobi.h:5:22: error: GGobiAPI.h: No existe el fichero ó directorio In file included from RSGGobi.h:6, from brush.c:1: conversion.h:174: error: expected ‘=’, ‘,’, ‘;’, ‘asm’ or ‘__attribute__’ before ‘asCLogical’ conversion.h:176: error: expected ‘=’, ‘,’, ‘;’, ‘asm’ or ‘__attribute__’ before ‘asCRaw’ --- snip --- brush.c:124: error: ‘t’ no se declaró aquà (primer uso en esta función) brush.c:124: error: ‘s’ no se declaró aquà (primer uso en esta función) brush.c:124: error: el objeto ‘GGOBI()’ llamado no es una función brush.c: En el nivel principal: brush.c:135: error: expected ‘)’ before ‘cid’ make: *** [brush.o] Error 1 chmod: no se puede acceder a `/usr/local/lib/R/site-library/rggobi/libs/*': No existe el fichero ó directorio ERROR: compilation failed for package 'rggobi' ** Removing '/usr/local/lib/R/site-library/rggobi' The downloaded packages are in /tmp/RtmpVCacJd/downloaded_packages Warning message: installation of package 'rggobi' had non-zero exit status in: install.packages("rggobi") --- What am I doing wrong? Thank you for your help. -- Kenneth Roy Cabrera Torres Cel 315 504 9339 __ 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 reproduce the same sampled units?
?set.seed On Wed, 02 May 2007 18:25:18 -0500, Santanu Pramanik <[EMAIL PROTECTED]> wrote: > Hi all, > Is it possible to generate the same sample number of times in R? In > SAS, using the option "seed" it is possible to reproduce exactly the > same sample. Is there any such feature in R which I can use? > For further clarity, > for (i in 1:2) > { > samp = sample(1:1000,100,replace = FALSE) > print(samp) > } > For the above simulation, is it possible to generate the same sampled > units twice? > Any suggestion would be gratefully acknowledged. > Thanks, > Santanu > Santanu Pramanik > JPSM,1218J Lefrak Hall > University of Maryland,College Park > Phone no.-301-314-9916 > > [[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. -- Kenneth Roy Cabrera Torres Cel 315 504 9339 __ 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] A problem with svIDE in Tinn-R?
Hi R users: I want to know if any one of you had the problem with Tinn-R, when you call the library svIDE on the new R 2.5.0, (because in the old R 2.4.1 works with out any problem). I got this message: > library(svIDE) Loading required package: tcltk Loading Tcl/Tk interface ... done Warning messages: 1: '\A' is an unrecognized escape in a character string 2: no se reconoce el valor 'save' Is it a critical warning? What would work and what not? What is it new in the new version that is not in the old version that makes this message appears? Thank you for your help. -- Kenneth Roy Cabrera Torres Cel 315 504 9339 __ 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] Cluster in large dataset with nominal variables.
Hi R users: Is it any R implementation of a cluster procedure for large data sets (clara()) but with dissimilary that can handle continuous, categorical and nominal variables (daisy()) like CLARANS (Clustering Large Applications based up on RANdomized Search), por example? Thank you for your help. -- Kenneth Roy Cabrera Torres Cel 315 504 9339 __ 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] A problem with "copy()" in the svIO package
Hi R-users: When I type the following code in R, using the svIO package I got a problem only with the "latex" type option.(With "raw", "ascii" and "html" options, it works fine). x<-1:50 copy('x', type='latex', objname='x') Erro en file.info(fn <- c(...)) : argumento del nombre de archivo inválido Thank you for your help. Kenneth -- __ 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 with plotmath
Hi R users: How can I have several subscript number with a comma in a plot. I would like to have the LaTeX equivalent of x_{i,j}. I try: plot(1:10,1:10,type="n") text(5,5,expression(x[i,j])) but it doesn´t work. Thank you for your help. __ 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] Problem installing local additional packages of Rcmdr (WinXP)
Hi R users: I download the .zip files and also PACKAGES and PACKAGES.gz files in a local directory, but when I call library(Rcmdr) It shows me a tcl/tk windows telling me to install the additional packages, the it takes me to another window where I must specify the path where the .zip files are located. I do it, but I obtain the following message: Loading required package: tcltk Loading Tcl/Tk interface ... done Loading required package: car Erro en gzfile(file, "r") : no fue posible abrir la conexión Versión del Rcmdr 1.1-7 And the Rcmdr frame opens but without installing the additional packages. I am sure I download the required packages and the PACKAGES and PACKAGES.gz files. What am I doing wrong? If I make the same procedure with the internet conexion it works fine, but I need to install all the packages in computers without internet conexion. Thank you for your help. Kenneth $platform [1] "i386-pc-mingw32" $arch [1] "i386" $os [1] "mingw32" $system [1] "i386, mingw32" $status [1] "" $major [1] "2" $minor [1] "3.1" $year [1] "2006" $month [1] "06" $day [1] "01" $`svn rev` [1] "38247" $language [1] "R" $version.string [1] "Version 2.3.1 (2006-06-01)" -- __ 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
Re: [R] Repeat dataframe
As I said, a very more elegant solution! Thank you! On Mon, 01 May 2006 11:52:44 -0500, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > my_df[rep(1:nrow(my_df), times=n), ] -- Using Opera's revolutionary e-mail client: http://www.opera.com/mail/ __ 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
Re: [R] Repeat dataframe
Well, I find a solution! If DFe is a data frame and n is an integer then DFr<-data.frame(t(matrix(rep(t(DFe),n),dim(DFe)[2],dim(DFe)[1]*n))) names(DFr)<-names(DFe) Will work!! Maybe somebody has a more elegant solution. Again, thank you for your help. On Mon, 01 May 2006 11:23:14 -0500, Kenneth Cabrera <[EMAIL PROTECTED]> wrote: > Hi R list: > > How can I "repeat" a data frame n times (with n>1000), > and obtain a new data frame where all the n data frames > are binded by rows? > > Thank you for your help > > Kenneth > -- Using Opera's revolutionary e-mail client: http://www.opera.com/mail/ __ 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
[R] Repeat dataframe
Hi R list: How can I "repeat" a data frame n times (with n>1000), and obtain a new data frame where all the n data frames are binded by rows? Thank you for your help Kenneth -- Using Opera's revolutionary e-mail client: http://www.opera.com/mail/ __ 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
Re: [R] VARIANCE across each ROW
?apply apply(data,1,var) On Sat, 01 Apr 2006 21:27:56 -0500, mark salsburg <[EMAIL PROTECTED]> wrote: > I have a very large matrix. I would like to display the variance across > each > row. > > In other words, I want to output a vector containing the values of > variance > across row. > > When I use the function var(), it seems to give me the variability of > each > column. > > Any ideas?? > > [[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 -- Using Opera's revolutionary e-mail client: http://www.opera.com/mail/ __ 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
Re: [R] Still problems with "step()" function
Ok Here is the complete command secuence -- testData<-structure(list(MCHNV = c(215, 325, 435, 327, 150, 949, 417, 528, 572, 1255, 517, 221, 306, 274, 307, 203, 258, 126, 88, 85, 75, 64, 77, 49, 133, 192, 32, 73, 139, 52, 17, 307, 127, 182, 24, 51, 141, 40, 76, 40, 100, 42, 18, 20, 18, 75, 26, 21, 216, 81, 90, 69, 138, 79, 69, 52, 59, 69, 173, 29, 118, 23, 1456, 156, 26, 208, 294, 84, 57, 84, 121, 112, 25, 9281, 36, 71, 812, 135, 215, 94, 302, 148, 42, 84, 114, 54, 309, 148, 87, 75, 275, 171, 210, 209, 150, 69, 157, 110, 93, 98, 40, 132, 155, 152, 62, 52, 50, 124, 248, 52, 181, 76, 118, 73, 66, 65, 75, 33, 71, 84, 223, 229), IA = c(26.9863956381346, 21.443816698054, 24.3172115328382, 22.9861973271014, 18.1768608561301, 10.6563785744498, 15.0910419778194, 13.6340158642317, 15.2055383362367, 16.9860295444068, 17.9909267660402, 26.714962883463, 20.1507989217486, 27.4221304422647, 11.806433704469, 25.7458563535912, 27.2816742408892, 26.0150722854506, 21.6368286445013, 25.1843952723718, 18.6270125223614, 29.9515311742675, 26.2046204620462, 12.6096949435855, 18.7464907355418, 10.5534246575342, 17.2654236774291, 17.0988006853227, 15.353091109001, 19.3008705114255, 6.27750073335289, 26.0375981018434, 17.2681438238026, 15.443279313632, 7.7125328659071, 30.1453193646502, 14.8701064405557, 24.7907034861927, 17.3296398891967, 10.7364199781261, 19.3800588668138, 10.3085483055134, 8.80185302168878, 16.6951761888471, 4.50797141286421, 6.84546862896979, 8.21845174973489, 6.19588744588745, 16.9889711436773, 5.59826387611726, 8.68694955964772, 3.65186966980101, 16.1123827476325, 9.67841682127396, 13.5423542354235, 5.6497484139138, 8.7686432740111, 11.9602529358627, 6.86294126054023, 19.2194992008524, 4.64579055441478, 9.5090599349543, 2.79913286437223, 10.2515243902439, 8.66738894907909, 3.6576768917803, 12.4953778763831, 15.2370770972205, 21.2257484260568, 7.8413206434768, 12.5950805097303, 6.23108480731724, 4.72859966550099, 3.11622388590868, 15.9319508257792, 6.69567319313512, 2.47749261819931, 24.8112304583644, 3.94306597422581, 2.71509876129896, 1.79760050991181, 4.80243967117475, 18.7818696883853, 5.24955116696589, 19.0728713579217, 12.6099706744868, 3.58964845561636, 11.3911041683795, 7.48350816481021, 15.0857882809971, 19.6855832458419, 12.2966073388367, 4.08228296275315, 21.4270750408514, 9.83262770545642, 6.28415300546448, 4.2754914877677, 10.4337099078042, 13.9745206313318, 14.0491684743312, 8.79131920024113, 15.2282010300124, 19.1379750063759, 6.89338235294118, 5.51527398074409, 18.3089033659066, 23.5368956743003, 10.3041872432253, 13.7563263088668, 12.0249520153551, 12.008763169, 9.25797503467406, 11.4912313643021, 14.2689601250977, 15.2254249815225, 14.0718562874251, 9.21820303383897, 9.9353894916494, 8.6011885743498, 14.7106815614409, 8.6682462336269, 10.9155102240965), IM = c(1.07456622007468, 1.05477685097601, 1.08970363580813, 1.09309079792113, 0.983952398124775, 0.97821952889272, 1.02464131763495, 0.913249100021176, 0.999321512574634, 1.00064160401003, 1.00311566921978, 1.03642655495856, 1.05350982751705, 1.03894736842105, 1.02694423605901, 1.0255795363709, 1.08604552865165, 0.9960092095165, 1.10034379028792, 1.05835010060362, 1.09069658719028, 1.08689655172414, 1.11247966534976, 1.06872703695699, 1.09098914000587, 0.947705442902882, 1.0462158808933, 1.06583294006607, 1.03178310316816, 1.02785822645152, 0.95582329317269, 1.03257159815996, 1.05732223283745, 1.05083088954057, 1.0375, 1.18054532056006, 1.08226897069872, 1.08357198646186, 1.11061739943873, 1.09309423884014, 1.14100039385585, 1.00222807372899, 1.02085106382979, 1.11658218682114, 1.08840413318025, 0.940623825629463, 1.07633027522936, 1.24954351795496, 1.06553284443751, 1.02449921811006, 1.03136617100372, 0.890711135611907, 1.01692767267131, 1.05692478931468, 1.00617965129111, 1.01988510826337, 1.02820400409177, 1.00325732899023, 0.969987048820317, 1.0319350473613, 0.948, 0.976430976430976, 0.912491945614184, 1.03805825242718, 1.01016333938294, 0.934502240405221, 0.936918076139056, 1.03298389158783, 1.06829657188414, 1.04918881601657, 1.0427807486631, 0.99369804237061, 1.00763125763126, 0.859808006713988, 1.06752246469833, 0.976588628762542, 0.900774538520513, 1.014353041988, 0.92939648374078, 0.934334930708457, 0.863916626564004, 1.00074278438031, 1.02873563218391, 0.882901764586573, 1.05623621661358, 1.09948152948801, 0.971704768125408, 1.0102470870176, 0.949198988195616, 0.999676323029616, 0.99394006170119, 1.02639721523883, 0.963297650493535, 1.02262711864407, 0.990626969124134, 1.02701607456437, 0.962433460076046, 1.05496997998666, 0.995055517002082, 0.99971081550029, 1.06194323596437, 1.0324851109908, 1.06694781233527, 1.01801801801802, 1.01774793140664, 1.01807723911257, 1.03578996939016, 1.00624427608026, 1.01527025469632, 0.990068754774637, 0.967759900990099, 0.954370905805286, 1.03321371306406, 1.01163888015099, 1.1162669447341, 0.998290598290598, 0.96477699
[R] Still problems with "step()" function
Hi R users: I don´t know if anybody have had the same problem with 'step()' funtion. If I type the commands --- testData<-read.table("testData.dat",header=T) model1J<-glm(MCHNV~offset(Offset1),data=testData,family="poisson") step(model1J,direction="forward", scope=(list(upper=~IA+IM+altura+Region+Region:IA+Region:IM+Region:altura+offset(Offset1), lower=~offset(Offset1 --- I got the folowing error message: --- . . Df Deviance AIC + IA:Region 8 540.24 1377.73 + IM 1 586.84 1410.32 + altura:Region 8 574.21 1411.70 600.74 1422.23 Step: AIC= 1377.73 MCHNV ~ Region + IA + altura + Region:IA + offset(Offset1) Error in factor.scope(ffac, list(add = fadd, drop = fdrop)) : upper scope does not include model --- What am I doing wrong? Thank you for your help. I put here the database I use, so you can replicate the problem. (I am using R221 on WinXP) --- testdata.dat "MCHNV" "IA" "IM" "altura" "Region" "Offset1" 215 26.9863956381346 1.07456622007468 2 "URABA" -0.185125484126689 325 21.443816698054 1.05477685097601 4 "URABA" 0.118671529717499 435 24.3172115328382 1.08970363580813 8 "URABA" 0.363948427905231 327 22.9861973271014 1.09309079792113 200 "URABA" 0.241376319575269 150 18.1768608561301 0.983952398124775 30 "BAJO CAUCA" -0.683196849706777 949 10.6563785744498 0.97821952889272 25 "URABA" 1.30128098679293 417 15.0910419778194 1.02464131763495 28 "URABA" 0.314810739840034 528 13.6340158642317 0.913249100021176 50 "BAJO CAUCA" 0.986562784859487 572 15.2055383362367 0.999321512574634 34 "URABA" 0.714908672341458 1255 16.9860295444068 1.00064160401003 2 "URABA" 1.54158724552593 517 17.9909267660402 1.00311566921978 50 "BAJO CAUCA" 0.67141268841394 221 26.714962883463 1.03642655495856 100 "BAJO CAUCA" -0.122167633974208 306 20.1507989217486 1.05350982751705 50 "BAJO CAUCA" 0.0411419433311752 274 27.4221304422647 1.03894736842105 125 "BAJO CAUCA" 0.0344014267173323 307 11.8064337044690 1.02694423605901 650 "NORDESTE" 0.242161557149972 203 25.7458563535912 1.02557953637090 75 "URABA" -0.608806032126194 258 27.2816742408892 1.08604552865165 1550 "NORTE" 0.340748793388473 126 26.0150722854506 0.9960092095165 1165 "NORTE" -0.551647618286246 88 21.6368286445013 1.10034379028792 1200 "NORTE" -0.901402119380404 85 25.1843952723718 1.05835010060362 1700 "NORTE" -0.763569644856491 75 18.6270125223614 1.09069658719028 1850 "NORTE" -1.03282454813011 64 29.9515311742675 1.08689655172414 1200 "OCCIDENTE" -1.03845836584836 77 26.2046204620462 1.11247966534976 650 "OCCIDENTE" -1.10563690360507 49 12.6096949435855 1.06872703695699 1475 "NORTE" -1.06421086195078 133 18.7464907355418 1.09098914000587 1535 "NORDESTE" -0.616186139423817 192 10.5534246575342 0.947705442902882 2300 "NORTE" 0.508021696433256 32 17.2654236774291 1.0462158808933 1875 "NORTE" -1.18417017702976 73 17.0988006853227 1.06583294006607 850 "OCCIDENTE" -1.04128722204884 139 15.3530911090010 1.03178310316816 1550 "NORDESTE" -0.245900538436826 52 19.3008705114255 1.02785822645152 1675 "NORTE" -0.473208760194684 17 6.27750073335289 0.95582329317269 2550 "NORTE" -1.69826912614072 307 26.0375981018434 1.03257159815996 450 "OCCIDENTE" 0.131028262406404 127 17.2681438238026 1.05732223283745 980 "NORDESTE" -0.520875959619492 182 15.4432793136320 1.05083088954057 700 "NORDESTE" -0.16960278438618 24 7.7125328659071 1.0375 1800 "NORTE" -1.76609172247948 51 30.1453193646502 1.18054532056006 25 "URABA" -2.21640739675299 141 14.8701064405557 1.08226897069872 1300 "OCCIDENTE" -0.138113302129634 40 24.7907034861927 1.08357198646186 1625 "OCCIDENTE" -1.18090753139494 76 17.3296398891967 1.11061739943873 1250 "NORDESTE" -1.11474167059799 40 10.7364199781261 1.09309423884014 700 "OCCIDENTE" -0.711311151187616 100 19.3800588668138 1.14100039385585 75 "MAGDALENA MEDIO" -0.86988435906 42 10.3085483055134 1.00222807372899 1800 "NORTE" -0.898942093539542 18 8.80185302168878 1.02085106382979 1925 "OCCIDENTE" -1.48722027970985 20 16.6951761888471 1.11658218682114 500 "OCCIDENTE" -1.99510039324608 18 4.50797141286421 1.08840413318025 2300 "NORTE" -1.22417551164346 75 6.84546862896979 0.940623825629463 1050 "NORDESTE" -0.648173814917214 26 8.21845174973489 1.07633027522936 2550 "NORTE" -1.48722027970985 21 6.19588744588745 1.24954351795496 1920 "OCCIDENTE" -1.9241486572738 216 16.9889711436773 1.06553284443751 1350 "OCCIDENTE" 0.188137942115395 81 5.59826387611726 1.02449921811006 2550 "NORTE" 0.162968828278140 90 8.68694955964772 1.03136617100372 1450 "NORDESTE" -0.233193887167711 69 3.65186966980101 0.89
[R] Help with "step()" function(2)
Hi R users: Sorry if the attached file didn't arrive. Im sending you in a zip format, hope it arrives. The data base is not too large so I copy here in ascii format: Thank you for your help Kenneth testData.dat "MCHNV" "IA" "IM" "altura" "Region" "Offset1" 215 26.9863956381346 1.07456622007468 2 "URABA" -0.185125484126689 325 21.443816698054 1.05477685097601 4 "URABA" 0.118671529717499 435 24.3172115328382 1.08970363580813 8 "URABA" 0.363948427905231 327 22.9861973271014 1.09309079792113 200 "URABA" 0.241376319575269 150 18.1768608561301 0.983952398124775 30 "BAJO CAUCA" -0.683196849706777 949 10.6563785744498 0.97821952889272 25 "URABA" 1.30128098679293 417 15.0910419778194 1.02464131763495 28 "URABA" 0.314810739840034 528 13.6340158642317 0.913249100021176 50 "BAJO CAUCA" 0.986562784859487 572 15.2055383362367 0.999321512574634 34 "URABA" 0.714908672341458 1255 16.9860295444068 1.00064160401003 2 "URABA" 1.54158724552593 517 17.9909267660402 1.00311566921978 50 "BAJO CAUCA" 0.67141268841394 221 26.714962883463 1.03642655495856 100 "BAJO CAUCA" -0.122167633974208 306 20.1507989217486 1.05350982751705 50 "BAJO CAUCA" 0.0411419433311752 274 27.4221304422647 1.03894736842105 125 "BAJO CAUCA" 0.0344014267173323 307 11.8064337044690 1.02694423605901 650 "NORDESTE" 0.242161557149972 203 25.7458563535912 1.02557953637090 75 "URABA" -0.608806032126194 258 27.2816742408892 1.08604552865165 1550 "NORTE" 0.340748793388473 126 26.0150722854506 0.9960092095165 1165 "NORTE" -0.551647618286246 88 21.6368286445013 1.10034379028792 1200 "NORTE" -0.901402119380404 85 25.1843952723718 1.05835010060362 1700 "NORTE" -0.763569644856491 75 18.6270125223614 1.09069658719028 1850 "NORTE" -1.03282454813011 64 29.9515311742675 1.08689655172414 1200 "OCCIDENTE" -1.03845836584836 77 26.2046204620462 1.11247966534976 650 "OCCIDENTE" -1.10563690360507 49 12.6096949435855 1.06872703695699 1475 "NORTE" -1.06421086195078 133 18.7464907355418 1.09098914000587 1535 "NORDESTE" -0.616186139423817 192 10.5534246575342 0.947705442902882 2300 "NORTE" 0.508021696433256 32 17.2654236774291 1.0462158808933 1875 "NORTE" -1.18417017702976 73 17.0988006853227 1.06583294006607 850 "OCCIDENTE" -1.04128722204884 139 15.3530911090010 1.03178310316816 1550 "NORDESTE" -0.245900538436826 52 19.3008705114255 1.02785822645152 1675 "NORTE" -0.473208760194684 17 6.27750073335289 0.95582329317269 2550 "NORTE" -1.69826912614072 307 26.0375981018434 1.03257159815996 450 "OCCIDENTE" 0.131028262406404 127 17.2681438238026 1.05732223283745 980 "NORDESTE" -0.520875959619492 182 15.4432793136320 1.05083088954057 700 "NORDESTE" -0.16960278438618 24 7.7125328659071 1.0375 1800 "NORTE" -1.76609172247948 51 30.1453193646502 1.18054532056006 25 "URABA" -2.21640739675299 141 14.8701064405557 1.08226897069872 1300 "OCCIDENTE" -0.138113302129634 40 24.7907034861927 1.08357198646186 1625 "OCCIDENTE" -1.18090753139494 76 17.3296398891967 1.11061739943873 1250 "NORDESTE" -1.11474167059799 40 10.7364199781261 1.09309423884014 700 "OCCIDENTE" -0.711311151187616 100 19.3800588668138 1.14100039385585 75 "MAGDALENA MEDIO" -0.86988435906 42 10.3085483055134 1.00222807372899 1800 "NORTE" -0.898942093539542 18 8.80185302168878 1.02085106382979 1925 "OCCIDENTE" -1.48722027970985 20 16.6951761888471 1.11658218682114 500 "OCCIDENTE" -1.99510039324608 18 4.50797141286421 1.08840413318025 2300 "NORTE" -1.22417551164346 75 6.84546862896979 0.940623825629463 1050 "NORDESTE" -0.648173814917214 26 8.21845174973489 1.07633027522936 2550 "NORTE" -1.48722027970985 21 6.19588744588745 1.24954351795496 1920 "OCCIDENTE" -1.9241486572738 216 16.9889711436773 1.06553284443751 1350 "OCCIDENTE" 0.188137942115395 81 5.59826387611726 1.02449921811006 2550 "NORTE" 0.162968828278140 90 8.68694955964772 1.03136617100372 1450 "NORDESTE" -0.233193887167711 69 3.65186966980101 0.890711135611907 2200 "NORTE" -0.415515443961666 138 16.1123827476325 1.01692767267131 550 "OCCIDENTE" -0.139262067333508 79 9.67841682127396 1.05692478931468 750 "OCCIDENTE" -0.638658995275876 69 13.5423542354235 1.00617965129111 950 "MAGDALENA MEDIO" -0.97021907389971 52 5.6497484139138 1.01988510826337 2475 "NORTE" -0.136965855073157 59 8.7686432740111 1.02820400409177 1975 "NORDESTE" -0.512493680866688 69 11.9602529358627 1.00325732899023 780 "OCCIDENTE" -0.807436326962073 173 6.86294126054023 0.969987048820317 1300 "VALLE DE ABURRA" 0.436963775167535 29 19.2194992008524 1.0319350473613 1800 "SUROESTE" -1.26584820804402 118 4.64579055441478 0.948 1425 "VALLE DE ABURRA" 0.392717535285662 23 9.50905993495431 0.976430976430976 1875 "ORIENTE" -1.21739582465808 1456 2.79913286437223 0.912491945614184 1450 "VALLE DE ABURRA" 2.50894864426824 156 10.2515243902439 1.03805825242718 1475 "NORDESTE" -0.0693500781347932 26 8.66738894907909 1.01016333938294 1650 "ORIENTE" -1.53247687129797 208 3.6576768917803 0.934502240405221 1425 "VALLE DE ABURRA" 0.718814927308523 294 12.4953778763831 0.936918076139056 125 "MAGDALENA MEDIO" 0.3
[R] Help with "step()" function
Hi R users: When I use this commands: - testData<-read.table("testData.dat",header=T) model1J<-glm(MCHNV~offset(Offset1),data=testData,family="poisson") step(model1J,direction="forward", scope=(list(upper=~IA+IM+altura+Region+Region:IA+Region:IM+Region:altura+offset(Offset1 - (the "testData" (testData.dat) is attached to the e-mail to reproduce the problem): I got this error: - . . . MCHNV ~ Region + IA + altura + offset(Offset1) Df Deviance AIC + IA:Region 8 540.24 1377.73 + IM 1 586.84 1410.32 + altura:Region 8 574.21 1411.70 600.74 1422.23 Step: AIC= 1377.73 MCHNV ~ Region + IA + altura + Region:IA + offset(Offset1) Error in factor.scope(ffac, list(add = fadd, drop = fdrop)) : upper scope does not include model - What am I doing wrong? I try also with: step(model1J,direction="forward", scope=(list(upper=~IA+IM+altura+Region+IA:Region+Region:IM+Region:altura+offset(Offset1 and with step(model1J,direction="forward", scope=(list(upper=~IA+IM+altura+Region+IA:Region+Region:IA+Region:IM+Region:altura+offset(Offset1 I got the same message Thank you for your help. Kenneth __ 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
Re: [R] Linux Distribution Choice
Here we use Scientific Linux Distribution 4.2 (https://www.scientificlinux.org/ ), that is based on RHEL4. In this platform we both have R(2.2.1.) and GRASS(6.0.2) working. Graham Smith wrote: I am making some tentative steps into using Linux (Mandriva at the moment) and notice that not all the Linux binaries on CRAN are the latest release. As R (plus Grass) will be key programs for me on Linux, is there a preferred Linux distribution that people in the R communiuty use? Many thanks, Graham [[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 __ 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
Re: [R] Update problems of Rcmdr?
Hi R users: I am not sure how do I fix it. But I erase an old R library installation (2.2.0), and also I erase all "file" directories on the R.2.2.1 library directory, and then I try library(Rcmdr) again and it works!. Any way, how could I trace or debug, so I get a hint whats going on with a similar problem in the future? Thank you very much for your help. Kenneth John Fox wrote: Dear Kenneth, After just running update.packages() myself, I can't duplicate this error (using R 2.2.1 under Win XP). Can you provide a little more information? Are any but the standard packages loaded when you issue the library(Rcmdr) command? Have you made any modifications to a startup file, such as Rprofile.site? I'm sorry that I can't offer any additional suggestions at this point. John John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Kenneth Cabrera Sent: Monday, February 13, 2006 7:54 PM To: r-help@stat.math.ethz.ch Subject: [R] Update problems of Rcmdr? Importance: High Hi R users: I got this messages when I try to use the Rcmdr library, after I make an update with update.packages() command. I am using R221 in windows environment. > library(Rcmdr) Loading required package: tcltk Loading Tcl/Tk interface ... done Loading required package: car Error in parse(file, n, text, prompt) : syntax error in "*" Error: .onAttach failed in 'attachNamespace' Error: package/namespace load failed for 'Rcmdr' Thank you for your help Kenneth __ 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 __ 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
[R] Update problems of Rcmdr?
Hi R users: I got this messages when I try to use the Rcmdr library, after I make an update with update.packages() command. I am using R221 in windows environment. > library(Rcmdr) Loading required package: tcltk Loading Tcl/Tk interface ... done Loading required package: car Error in parse(file, n, text, prompt) : syntax error in "*" Error: .onAttach failed in 'attachNamespace' Error: package/namespace load failed for 'Rcmdr' Thank you for your help Kenneth __ 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
[R] Data frame index?
Hi, R users: I have a data.frame (not a matrix), I got a vector with the same length as the number of records (rows) of the data frame, and each element of that vector is the column number (in a specific range of columns) of the corresponding record that I must set to zero. How can I do this without a "for" loop? Thank you for your help. Kenneth __ 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
[R] Which is the best hardware?
Hi R users: In your opinion and experience, which hardware configuration is the best to run R over LINUX ? With "best" I mean best performance, and also cheapest. (about U$ 2.000 the whole basic system: mother board+CPUs+RAM+HD) By the way, which LINUX distribution is the best to run R with high computing technics (simulation, bayesian, etc) and huge data base? and in combination with what kind of (cheap) hardware? Thank you very much for your help. Kenneth Happy Holydays!!! __ 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
Re: [R] How to sort data sets?
If you want to sort a data frame according a variable, the use it like datfram[order(datfram$var),] to sort all the data.frame, not just a variable Spencer Graves wrote: ?order Martin Lam wrote: Hi, I was wondering if someone know how to sort a data set by column. I've tried sort() but without luck. I would think there should be a function for it somewhere. An example with the iris data set would be appreciated. Thanks, Martin __ 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 __ 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
Re: [R] Memory question
Thank you Dr. Spencer Graves for your answer. What kind of matrices? They come form an image of about 3000x5000, and I need to generate arround 1024 matrices of the same size, they are not sparse matrices. What function can I use to, once generated one matrix, I could save into disk and then use the same space for the following, and so on. Thank you very much for your help Kenneth Spencer Graves wrote: What kinds of matrices? There are facilities in the Matrix and SparseM packages that might help for sparse matrices. If they are N x k where N is large and k is not, can you compute something like the QR decomposition and get away with keeping only the R part for most of your matrices? One could potentially define a class of matrices that are only kept in memory only when needed; I think S-Plus may do that. It would take a lot of work to make that work generally, but you might be able to accomplish what you need with a much smaller effort. spencer graves Kenneth Roy Cabrera Torres wrote: Hi R users and developers: I want to know how can I save memory in R for example: - saving on disk a matrix. - using again the matrix (changing their values) - saving again the matrix on disk in a different file. The idea is that I have a process that generate several matrices, but if I keep them all in memory it will overflow. How can I save them in different files, so I use the same amount of memory for each processed matrix? Thank you for your help. __ 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
[R] Problem installing packages.
Hello, R users: I got the following problem when I try to install some packages using the 2.1.0 patched version on Window$. This is only an example. The same happens when I try to install ade4. (it download the packages and their dependency, but when it finish to download it shows the same error) > utils:::menuInstallPkgs() trying URL 'http://cran.at.r-project.org/bin/windows/contrib/2.1/CoCo_0.1.6.5.zip' Content type 'application/zip' length 1895592 bytes opened URL downloaded 1851Kb package 'CoCo' successfully unpacked and MD5 sums checked bundle 'CoCo' successfully unpacked and MD5 sums checked Error in sprintf(gettext("unable to move temp installation '%d' to '%s'"), : use format %s for character objects What am I doing wrong? Thank you for your help. Kenneth __ 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
[R] Big databases
Hi R users: I am using R 2.01. over a LINUX(Scientific Linux CERN) platform and I got the following problem: It takes too much time (more than 6 hours, because I stop the process) to read a 270MB database. I am using read.table() function. Is there any workaround to read faster a big data base? Thank you for your help. Kenneth __ 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
[R] Legal issues
Hi R Users: I have some legal questions about R development. R is under GNU Licence, version 2. How the core team deal with their employer (most of the cases it Universities, I don't know if public or private) about their contribution to R development? Do you have to make all the work outside the institution? Is it part of the work at the institution, the R development? Do the institution give resources to the team, for example, the computer server where R is allocated? Thank you very much for your help! -- Kenneth Roy Cabrera Torres Universidad Nacional de Colombia Sede MedellÃn Cel 315 504 9339 __ 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
[R] Problem with graphic window in LINUX
Hi R users: I got a problem when I make graphics in LINUX, the graphic window is static, I can resize it, but I can not move it (there is no frame), * what should I configure in LINUX to make them work again? * Is it posible to have a history list of the graphics like in W2K, so I can see them again? Thank you for your help Kenneth __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] How about a mascot for R?
I vote for R-madillo! Erin Hodgess wrote: How about an R-madillo? Erin from Texas __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html -- Kenneth Roy Cabrera Torres Celular +57 (315) 504 9339 __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] Standard error of a sum
If you supose that each mean is independent then s.err_total_chem <- sqrt(sum(s.err_event^2)) If not you must estimate de var-covariance matrix to find this s.err. This estimation will be difficult given the scarce data available. Hope it will helps. Kenneth Jose A. Hernandez wrote: Folks, This is so simple is driving me crazy. It's not really an R question is more an Statistics question. I applied a chemical in 3 different events during a growing season. I took 3 samples in each event to estimate mean applied chemical. Thus, I can easily estimate means and st.err by event. I can then estimate the total chemical applied for the whole season by adding the 3 means for each event. How can I estimate the standard error of this total applied chemical. Thanks in advance, Jose event <- c(1,1,1,2,2,2,3,3,3) chem <- c(131.75,125.92,130.09,42.67,48.26,54.56,85.05,91.17,90.86) d.1 <- data.frame(event,chem) attach(d.1) s.err <- function(x) sd(x)/(sqrt(length(x))) print(mean_event <- tapply(d.1$chem,factor(data.1$event),mean)) print(s.err_event <- tapply(d.1$chem,factor(data.1$event),s.err)) print(total_chem <- sum(mean_event)) s.err_total_chem <- -- Kenneth Roy Cabrera Torres Celular +57 (315) 504 9339 __ [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] Problem with Rcmr in R 1.9.0patched
Dear R users: I install the new R version (1.9.1 patched) in a W2K platform in the E:\rw1091patch path. Then I reinstall ALL the packages from CRAN, using the "Packages > Install package(s) from CRAN" menu option. But when I call the library library(Rcmdr) It doen't do any thing. What am I missing? Thank you for your help. PS: Also when I look for the help files using help.start(), it answers me: "Could not open file" file://localhost/E:/rw1091/library/Rcmdr/html/00Index.html -- Kenneth Cabrera Universidad Nacional de Colombia Tel Of 430 9351 Cel 315 504 9339 Medellín __ [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] for() to lapply()
Hi dear R-users: I have the following problem: I have a list of data.frames (12 variables and 6 rows, each) I have to merge from an specific point of the list to the end of the list, I am doing so with a for() loop but it is too inefficient and it exhausts memory. How can I convert this for() loop in a function and then use lapply? - LIST: list of data frames. m:point of start merging. result<-merge(LIST[[m]],LIST[[m+1]],by="key",all.x=T) for (i in (m+2):length(LIST)) { result<-merge(result,LIST[[i]],by="key",all.x=T) } The problem is the acumulative: "result<-merge(result,...)" I can think in a function like this: mergeListelem<-function(i,LS,m) { merge(LS[[m]],LS[[i]],by="key",all.x=T) } ind<-m:length(LIST) lapply(ind,mergeListelem,LIST,m) But I just obtain a list of merged data.frames, and I donÂt know how to join all the elements of the list in just one final data.frame. Thank you for your help -- Kenneth Cabrera Universidad Nacional de Colombia Tel Of 430 9351 Cel 315 504 9339 MedellÃn __ [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] model.tables (se=T) with interaction
Hi R users. Why if I have a two-way or multifactor analysis of variance and I call the function model.tables with the se=T options there is not estimation when I have an interaction on the model? Thank you very much for your help -- Kenneth Cabrera Universidad Nacional de Colombia Tel Of 430 9351 Cel 315 504 9339 Medellín __ [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
[Fwd: Re: [R] normal scores test]
-- Kenneth Roy Cabrera Torres Celular +57 (315) 504 9339 --- Begin Message --- see in package quantreg Hope this helps. Padmanabhan, Sudharsha wrote: Hello, I need help in performing a Van_der_Waerden normal scores test in R. I have two arrays of scores(final on therapy scores from drug and placebo) and want to use the normal scores procdeure to test for significance. (observations are unequal in number - due to dropouts). Could you please help me out with the coding or let me know if there is a package that can be used (for example, suppdist etc.) Thanks in advance Regards ~S __ [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 -- Kenneth Roy Cabrera Torres Celular +57 (315) 504 9339 --- End Message --- __ [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] citation() doesn´t work
Hi R users: I want to know if you have the same problem with the citation() function, it doesnÂt work! Thank you R in W2K, version 1.8.1. -- Kenneth Cabrera Universidad Nacional de Colombia Tel 430 9351 Cel 315 504 9339 MedellÃn __ [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] Reading formated databases
Hi R-users: How can I read an ascii database that is controled by the column number? For example: 7349593Luis Miguel Ariza Gutierrez 32342123 9394583X XX 34234930 39483 CCC CC39203230 3484932YY YYZZ ZZZ39402343 39203 WWW V 342343 There are 4 variables, ID, Name, Last Name, Numeric Variable. 1 Variable column 1-8 2 Variable column 9-25 3 Variable column 26-51 4 Variable column 51-59 Thank you for your help Kenneth -- Kenneth Roy Cabrera Torres Celular +57 (315) 405 9339 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] ts format for daily time serie
Prof Brian Ripley <[EMAIL PROTECTED]> wrote: It's probably best to handle this an an x-y plot (?plot.POSIXct) or as an irregular time series (see package its on CRAN). A daily time serie is not a regulary spaced time serie? The ts function is designed for time series that you want to do time-series analysis on. That is what I want, make time-series analysis, but at the same time I want to plot it directly. Achim Zeileis <[EMAIL PROTECTED]> wrote: ts() only allows regulary spaced time series. Try the package its or the irts() function in package tseries for irregularly spaced time series. hth, Z A daily time serie is not a regulary spaced time serie? Excerp of the database: FECHA TRM 101/01/2000 1873.77 202/01/2000 1873.77 303/01/2000 1873.77 404/01/2000 1874.35 505/01/2000 1895.97 . . . 1397 10/11/2003 2843.82 1398 11/11/2003 2840.41 1399 12/11/2003 2840.41 1400 13/11/2003 2845.69 1401 14/11/2003 2850.24 1402 18/11/2003 2842.53 1403 19/11/2003 2831.97 1404 20/11/2003 2826.60 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] ts format for daily time serie
Hi R-users: How can I format a daily time serie with ts function so the plot of the time shows the date right (dd/mm/yy) or yy. ? Excerp of the database: FECHA TRM 101/01/2000 1873.77 202/01/2000 1873.77 303/01/2000 1873.77 404/01/2000 1874.35 505/01/2000 1895.97 . . . 1397 10/11/2003 2843.82 1398 11/11/2003 2840.41 1399 12/11/2003 2840.41 1400 13/11/2003 2845.69 1401 14/11/2003 2850.24 1402 18/11/2003 2842.53 1403 19/11/2003 2831.97 1404 20/11/2003 2826.60 Thank you for your help. -- __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] file not found?
Hi R maintainers: when I use update.packages() And I try to update the package "maps" maps : Version 2.0-9 in D:/rw1080/library Version 2.0-10 on CRAN Update (y/N)? y The following message appears: trying URL `http://cran.r-project.org/bin/windows/contrib/1.8/maps_2.0-10.zip' Error in download.file(url, destfile, method, mode = "wb") : cannot open URL `http://cran.r-project.org/bin/windows/contrib/1.8/maps_2.0-10.zip' In addition: Warning message: cannot open: HTTP status was `404 Not Found' All the rest of packages work fine. I use R 1.8.0 on W2K platform. Thank you for your help Kenneth -- Kenneth Roy Cabrera Torres Celular +57 (315) 405 9339 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Handling memory
Dear R-Users: Which is the best way to fee memory inside a function to make room for some other process that uses most of the RAM available and latter retrive this data again to finish the function? Thank you very much for your help Kenneth Cabrera -- __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] geoR vector length error
Dear Dr Thomas What should I do in this case? How can I estimate the variogram using the full data set? Thank you very much for your help. Kenneth On Tue, 25 Mar 2003, Kris Nackaerts wrote: Dear, we have a problem with geoR. We try to read an ASCII table (x,y,z) with 4 lines. With read.geodata we get the error: Error in vector("double", length) : cannot allocate vector of length 79998 We can read the file without any problem with read.table, but trying to convert it to the geodata class gets the same error. Looking at the source of as.geodata() it seems that it uses dist() to check for coincident points, and this will involve creating a vector of length approximately 8x10^8. You can't do that in R. -thomas __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help -- __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Indicator Kriging
Hi R-Users Is there any package in R that work Indicator Kriging? Or somebody is working in a link between GSLIB library in FORTRAN for R? Thank you very much for your help. -- __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] Trasparent graphs?
Yes, Dr. Blackwell: What I mean is the second part: If you need to give the appearance of color mixing, you will have to identify which points coincide, and plot those again using a third color. I don't use lattice graphics, but I presume the situation is much the same. Is it possible with pixmap library or image options? Thank you again for your help. Kenneth __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Trasparent graphs?
Hi R-users and R-experts: Is it posible to make transparent graphics on R? I mean, one graphic with colors over another one but both visible with transparent options? Thank you for your help Kenneth -- __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] convert a Character-string to a number
Try as.numeric(gsub(",",".",a)) Hope it helps!! Best regards Kenneth On Fri, 12 Sep 2003 14:29:13 +0200, Michael Kirschbaum <[EMAIL PROTECTED]> wrote: Hi Everyone. I have a simple problem but don't know, how to get along. how can I convert the vector a<-c("0,01","1,00") in a vector b<-c(0.01,1.00) Thank you for suggestions M.Kirschbaum [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help -- __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] identify with image
Hi, Dr Halvorsen Try identify(expand.grid(1:512,1:512)) But, be careful with the index of the matrix (image) On Tue, 02 Sep 2003 15:31:45 -0400, kjetil brinchmann halvorsen <[EMAIL PROTECTED]> wrote: Hola! I will want to identify pixels in an image with the mouse, for so getting the image data from the matrix(es), for use in subsequent discriminant analysis. But the following bombs R: (windows XP, rw1071) str(baboon) list() - attr(*, "size")= int [1:2] 512 512 - attr(*, "cellres")= num [1:2] 1 1 - attr(*, "bbox")= num [1:4] 0 0 512 512 - attr(*, "channels")= chr "grey" - attr(*, "bbcent")= logi FALSE - attr(*, "class")= chr "pixmapGrey" - attr(*, "grey")= num [1:512, 1:512] 0.537 0.510 0.345 0.259 0.322 ... class(baboon) [1] "pixmapGrey" library(pixmap) plot(baboon) identify(baboon, n=1) ... and then R bombs! What to do? Kjetil Halvorsen __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help -- __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] Gram-Schmidt orthonormal factorization
See ?qr and http://ikpe1101.ikp.kfa-juelich.de/briefbook_data_analysis/node224.html Note the compact form of how QR decomposition in R works. On Mon, 01 Sep 2003 08:54:55 -0300, CENDOYA, Gabriela <[EMAIL PROTECTED]> wrote: Hi: Does R have a function as gsorth is SAS, that perform a the Gram-Schmidt orthonormal factorization of the m ×n matrix A, where m is greater than or equal to n? That is, the GSORTH subroutine in SAS computes the column-orthonormal m ×n matrix P and the upper triangular n ×n matrix T such that A = P*T. or any other version of Gram-Schmidt orthonormal factorization? I search the help, but I didn’t find any thing about. Thanks, Gabriela PD: I am using R 1.7.1, for Windows 98. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help -- __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Re: [S] subscripts error
Hi Gianluca: try : x2<-x[1:(length(x)-2)]-x[3:length(x)] OR x2<-diff(x,lag=2) Hope it helps. Gianluca wrote: Hi all, I'm back to use S-Plus after a lot programming with SAS, so it's hard for me to face many simple problems Let's come to what is going to make me crazy. I have a simple 1000 obs time series, called "x". If I want to compute second order differences just typing: x2<-x[1:length(x)-2]-x[3:length(x)] but I get the following error message: Error in x[1:length(x) - 2]: Only 0's may be mixed with negative subscripts Much more interesting is that with first order differences it works perfectly Thanks in advance to all. Gianluca -- Kenneth Roy Cabrera Torres Celular +57 (315) 405 9339 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Intervention/Impact analysis in time series
Hi R users: Does any one knows about a R library for deal with intervention/impact analysis in time series (eg. Box-Tiao et. al. theory?). Thank you for your help -- Kenneth Roy Cabrera Torres Celular +57 (315) 405 9339 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Three graphs
Hi R-users: How can I obtain 3 graphics in a page but, one on the top and centered, and the other two in the bottom ? Thank you for your help -- Kenneth Roy Cabrera Torres Celular +57 (315) 405 9339 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] Windows and Linux version
If you are asking about .Rdata files, the answer is yes... I do that several times and it works without any problem. The .Rhistory file is an ASCII file, so, no problem... Best Regards Saludos desde Colombia Kenneth [EMAIL PROTECTED] wrote: Hi, I would like to know if my files saved on R windows version will run normally on Linux version? Thanks Cecília _ Voce quer um iGMail protegido contra vírus e spams? Clique aqui: http://www.igmailseguro.ig.com.br __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help -- Kenneth Roy Cabrera Torres Celular +57 (315) 405 9339 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Sorry for the mistake!
DDear R Users! SSorry for the mistake!!! II send you a .zip file of some database. SI ask you for pardon! KKenneth -- Using M2, Opera's revolutionary e-mail client: http://www.opera.com/m2/ __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] fft help
Hi R-users: I want to know if there is an easy way to obtain a Fourier Transform form a vector or an array (just like fft does), but with a more density base. I mean, if I have a vector of 512 of length, I want the Fourier Transform to be 1024, or 2048, etc, in length (de u domain). Or should I modify the fft C code to do that? If I want to modify the precision of the fft function, which parameter of the .Machine option do I modify, and how? Thank you for your help. Kenneth Cabrera __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R]
use the option sep="\t" on write.table. see ?write.table On Mon, 10 Mar 2003 11:55:46 -0800, Francisco J Molina <[EMAIL PROTECTED]> wrote: Subject: separator = TAB in write.table X-Mailer: VM 7.00 under 21.4 (patch 6) "Common Lisp" XEmacs Lucid Reply-To: [EMAIL PROTECTED] FCC: /home/f/.xemacs/mail/sent How can I choose my separator to be TAB in functions like write.table ()? __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help -- Using M2, Opera's revolutionary e-mail client: http://www.opera.com/m2/ __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] fptex link?
Hi R developers: I am collecting the elements to Build R on W2K. But I found that the www.fptex.org link, lead me to www.dante.de site. In this site I don't found where to download the fptex software. And worst: ich spreche kein Deutch!!! Danke sehr! __ [EMAIL PROTECTED] mailing list http://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] What´s wrong with update.packages()
Dear Dr. Ripley: [EMAIL PROTECTED] wrote: Please take a look at library/tseries/DESCRIPTION and see what version it says. My guess is that it says 0.9-6. Your are right! My guess is that you got a cached copy of the tseries.zip file. Please try again, perhaps using another CRAN mirror, or get the file yourself via ftp. I will do that! Where the file is cached, in the proxy server, maybe? Thank you again On Fri, 24 Jan 2003, Kenneth Cabrera wrote: Hi R users: I'm using R 1.6.2 version in W2K platform. When I use update.packages() to update from the CRAN site it shows me several outdated packages, for example tseries : Version 0.9-6 in C:/rw1062/library Version 0.9-7 on CRAN The system download them without any problem. When the program ask me: Delete downloaded files (y/N)? I answer "n" But when I use again the same update.packages() function, it appears the same packages that it was supposed to be updated, with the same messages. Or if I call a package, for example 'tseries' it appears: `tseries' version: 0.9-6 `tseries' is a package for time series analysis and computational finance. See `library (help=tseries)' for details. The old package version, not the 0.9-7 new version __ [EMAIL PROTECTED] mailing list http://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] What´s wrong with update.packages()
Hi R users: I'm using R 1.6.2 version in W2K platform. When I use update.packages() to update from the CRAN site it shows me several outdated packages, for example tseries : Version 0.9-6 in C:/rw1062/library Version 0.9-7 on CRAN The system download them without any problem. When the program ask me: Delete downloaded files (y/N)? I answer "n" But when I use again the same update.packages() function, it appears the same packages that it was supposed to be updated, with the same messages. Or if I call a package, for example 'tseries' it appears: `tseries' version: 0.9-6 `tseries' is a package for time series analysis and computational finance. See `library (help=tseries)' for details. The old package version, not the 0.9-7 new version. What am I doing wrong? Thank you for your help Kenneth Cabrera __ [EMAIL PROTECTED] mailing list http://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] ConTEXT .chl file for R?
Hi R-Users Does any body on the list have build a .chl file for R using the ConTEXT editor? Thank you for your help Kenneth __ [EMAIL PROTECTED] mailing list http://www.stat.math.ethz.ch/mailman/listinfo/r-help