Re: [R] MCD distance
There first two *are* in the lqs package: you just need to use mahalanobis() with cob.rob(). I've never heard of BACON. On Fri, 6 Jun 2003, osama wrote: I would like to compare robust distance measures,e.g. MVE, MCD, BACON, as measures of leverage and to detect outliers in X-space . I could not find it in LQS or MASS backages. Why look in only two packages? -- 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
[R] XML package for R
Hi, Does anyone have a binary of XML package for 1.7.0 on Windows? I have searched it at the RSXML pages (http://www.omegahat.org/RSXML/ ) but it seems it is no longer support. Any pointers for the compiled version will be appreciated. Regards, mushtaq [[alternate HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] XML package for R
thanks :) -mushtaq - Original Message - From: Prof Brian Ripley [EMAIL PROTECTED] To: Dr Mushtaq Ahmed [EMAIL PROTECTED] Cc: [EMAIL PROTECTED] Sent: Monday, June 09, 2003 12:44 PM Subject: Re: [R] XML package for R See the ReadMe in the appropriate area on CRAN (XML is also a CRAN packages). http://cran.r-project.org/bin/windows/contrib/1.7/ReadMe which points you to http://www.stats.ox.ac.uk/pub/RWin On Mon, 9 Jun 2003, Dr Mushtaq Ahmed wrote: Does anyone have a binary of XML package for 1.7.0 on Windows? I have searched it at the RSXML pages (http://www.omegahat.org/RSXML/ ) but it seems it is no longer support. Any pointers for the compiled version will be appreciated. You could always compile it yourself! -- 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
Re: [R] looking for Prof Bates' file
Viet Nguyen,,, [EMAIL PROTECTED] writes: Hello I'm reading up on fitting truncated Weibull distribution to data. There are posts in 2002 that point to this presentation by Prof Bates: http://www.stat.wisc.edu/~bates/JSM2001.pdf but now the file is not there. I can't find it anywhere else, Google doesn't have a cached copy for it. Could someone please give me a copy of this file, if they have it? Thanks and regards, viet. Looks like www.stat.wisc.edu is no longer identical to Doug's desktop machine. The file is still there, but I'm not too sure whether he would want us to tell the world how to get to it, so perhaps we should wait for him to reply... -- 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
Re: [R] Beginner questions
Michele Grassi wrote: Hi, when i create a new variable e.g var-c(1,2,3,...), how can i delete one or more osservation? See An Introduction to R, Section 2.7. When i draw my scatterplot,i want to identify, in the graph,some point whit x,y coordinates. How can i do it? See An Introduction to R, Section 12.3. Thanks. Michele. Please read the manuals! Uwe Ligges __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] installing XML in linux
Hi, I was trying to install XML package in a linux (dabian) machine and got the following error message. I am not sure whether there is any error in my linux installion. I appreciate suggestions to install XML properly in this machine. Thanks in advance. Mahbub. ##Errors $ R CMD INSTALL -l /usr/local/lib/R/library/ XML_0.93-4.tar.gz * Installing *source* package 'XML' ... creating cache ./config.cache checking for gcc... gcc checking whether the C compiler (gcc ) works... yes checking whether the C compiler (gcc ) is a cross-compiler... no checking whether we are using GNU C... yes checking whether gcc accepts -g... yes checking how to run the C preprocessor... gcc -E checking for xml-config... no checking for libxml/parser.h... no checking for gnome-xml/parser.h... no checking for libxml/parser.h... (cached) no checking for gnome-xml/parser.h... (cached) no checking for libxml/parser.h... (cached) no Cannot find parser.h. Set the value of the environment variable LIBXML_INCDIR to point to where it can be found. ERROR: configuration failed for package 'XML' version _ platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major1 minor7.0 year 2003 month04 day 16 language R __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] executable R scripts
On 20030608 (Sun) at 1735:14 -0700, John Zedlewski wrote: Hi, I'm a newbie trying to make an R program executable on UNIX, just like one would write an executable perl script by putting #!/usr/bin/perl in the first line, and so on. It seems, though, that this would only work if I use the BATCH command to tell R to execute the program in its first argument. This would have the unfortunately side-effect of dumping all output to a file rather than stdout. Additionally, I'd want to see only the results of print statements on stdout, not all off R's output, just as when you source a script with echo=FALSE. I've seen the other replies, but thought this might be of interest too: I hacked a hashbang wrapper so you can start an R script with #! /bin/sh /usr/bin/setR and then invoke it with command-line arguments, which get passed to the script in a character vector called argv. See http://www.bgl.nu/~glouis/setR.html if you're interested. I think the version displayed has a couple of ampersand-lt; that need to be changed to if you use 'save as' rather than cutting and pasting. -- | G r e g L o u i s | gpg public key: finger | | http://www.bgl.nu/~glouis | [EMAIL PROTECTED] | | http://wecanstopspam.org in signatures fights junk email | __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] looking for Prof Bates' file
Peter Dalgaard BSA [EMAIL PROTECTED] writes: Viet Nguyen,,, [EMAIL PROTECTED] writes: Hello I'm reading up on fitting truncated Weibull distribution to data. There are posts in 2002 that point to this presentation by Prof Bates: http://www.stat.wisc.edu/~bates/JSM2001.pdf but now the file is not there. I can't find it anywhere else, Google doesn't have a cached copy for it. Could someone please give me a copy of this file, if they have it? Thanks and regards, viet. Looks like www.stat.wisc.edu is no longer identical to Doug's desktop machine. The file is still there, but I'm not too sure whether he would want us to tell the world how to get to it, so perhaps we should wait for him to reply... Thanks for your interest in the file, Nguyen. As Peter guessed, I changed a configuration so that www.stat.wisc.edu/~bates is served from a departmental machine instead of a private machine. I have copied the file JSM2001.pdf to the departmental server so the URL http://www.stat.wisc.edu/~bates/JSM2001.pdf is valid again. By the way, that PDF file opens in 'full-screen' mode, which sometimes confuses browsers. I recommend that you download the file then open it in a PDF viewer application like acrobat. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Re: [BioC] installing XML in linux
You will need to install the Debian package libxml-dev before you can install the R XML package. Mahbub Latif [EMAIL PROTECTED] writes: I was trying to install XML package in a linux (dabian) machine and got the following error message. I am not sure whether there is any error in my linux installion. I appreciate suggestions to install XML properly in this machine. Thanks in advance. Mahbub. ##Errors $ R CMD INSTALL -l /usr/local/lib/R/library/ XML_0.93-4.tar.gz * Installing *source* package 'XML' ... creating cache ./config.cache checking for gcc... gcc checking whether the C compiler (gcc ) works... yes checking whether the C compiler (gcc ) is a cross-compiler... no checking whether we are using GNU C... yes checking whether gcc accepts -g... yes checking how to run the C preprocessor... gcc -E checking for xml-config... no checking for libxml/parser.h... no checking for gnome-xml/parser.h... no checking for libxml/parser.h... (cached) no checking for gnome-xml/parser.h... (cached) no checking for libxml/parser.h... (cached) no Cannot find parser.h. Set the value of the environment variable LIBXML_INCDIR to point to where it can be found. ERROR: configuration failed for package 'XML' __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
RE: [R] Basic question on applying a function to each row of ada taframe
Another neat way is: with(DF, foo(x, w)) HTH, Andy -Original Message- From: peter leonard [mailto:[EMAIL PROTECTED] Sent: Sunday, June 08, 2003 4:35 PM To: [EMAIL PROTECTED] Subject: [R] Basic question on applying a function to each row of a dataframe Hi, I have a function foo(x,y) and a dataframe, DF, comprised of two vectors, x w, as follows : x w 1 1 1 2 2 1 3 3 1 4 4 1 etc I would like to apply the function foo to each 'pair' within DF e.g foo(1,1), foo(2,1), foo(3,1) etc I have tried apply(DF,foo) apply(DF[,],foo) apply(DF[DF$x,DF$w],foo) However, none of the above worked. Can anyone help ? Thanks in advance, Peter __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help -- Notice: This e-mail message, together with any attachments, cont... {{dropped}} __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Nested anovas
Hello, Is anyone out there doing nested anovas? I'm familiar with the notation to designate the nesting in SAS but can't seem to find any for R. How can I specify the nesting levels? Thanks, Suzanne __ McAfee VirusScan Online from the Netscape Network. Comprehensive protection for your entire computer. Get your free trial today! http://channels.netscape.com/ns/computing/mcafee/index.jsp?promo=393397 Get AOL Instant Messenger 5.1 free of charge. Download Now! __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Re: [BioC] installing XML in linux
Thanks a lot Prof Bates. One more thing... What Debian packages should I install to install R packages Rgraphviz and rhdf5? Mahbub. --- Douglas Bates [EMAIL PROTECTED] wrote: You will need to install the Debian package libxml-dev before you can install the R XML package. Mahbub Latif [EMAIL PROTECTED] writes: I was trying to install XML package in a linux (dabian) machine and got the following error message. I am not sure whether there is any error in my linux installion. I appreciate suggestions to install XML properly in this machine. Thanks in advance. Mahbub. ##Errors $ R CMD INSTALL -l /usr/local/lib/R/library/ XML_0.93-4.tar.gz * Installing *source* package 'XML' ... creating cache ./config.cache checking for gcc... gcc checking whether the C compiler (gcc ) works... yes checking whether the C compiler (gcc ) is a cross-compiler... no checking whether we are using GNU C... yes checking whether gcc accepts -g... yes checking how to run the C preprocessor... gcc -E checking for xml-config... no checking for libxml/parser.h... no checking for gnome-xml/parser.h... no checking for libxml/parser.h... (cached) no checking for gnome-xml/parser.h... (cached) no checking for libxml/parser.h... (cached) no Cannot find parser.h. Set the value of the environment variable LIBXML_INCDIR to point to where it can be found. ERROR: configuration failed for package 'XML' __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
RE: [R] Basic question on applying a function to each row of a dataframe
Another another neat way is: DF - data.frame(x=1:4, y=rep(1,4)) foo - function(z,x,y)z[x]+z[y] apply(DF,1,foo,x=x,y=y) 1 2 3 4 2 3 4 5 or DF - data.frame(x=1:4, y=rep(1,4)) foo - function(z)z-x+y foo(DF) [1] 2 3 4 5 Liaw, Andy [EMAIL PROTECTED] wrote: Another neat way is: with(DF, foo(x, w)) HTH, Andy -Original Message- From: peter leonard [mailto:[EMAIL PROTECTED] Sent: Sunday, June 08, 2003 4:35 PM To: [EMAIL PROTECTED] Subject: [R] Basic question on applying a function to each row of a dataframe Hi, I have a function foo(x,y) and a dataframe, DF, comprised of two vectors, x w, as follows : x w 1 1 1 2 2 1 3 3 1 4 4 1 etc I would like to apply the function foo to each 'pair' within DF e.g foo(1,1), foo(2,1), foo(3,1) etc I have tried apply(DF,foo) apply(DF[,],foo) apply(DF[DF$x,DF$w],foo) However, none of the above worked. Can anyone help ? Thanks in advance, Peter __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help -- Notice: This e-mail message, together with any attachments, cont... {{dropped}} __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help - [[alternate HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Loops question
Hello Mr Ripley, i was waiting you to ask and if you don't mind if there is a fast way in R to do these loops made with R and that takes a week because my matrix is with thousand of rows. I got another method that is more fast and is utilised by many packages : interfacing with other languages, but it is only a question for information foo-function() {for (i in 1:(n-1)){ for (k in (i+1):n){ b_0 for (j in 1:m){ if(bin[i,j]==TRUEbin[k,j]==TRUE) b_b+1} print (b) } } } Thanks a lot Ramzi - [[alternate HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] Basic question on applying a function to each row of a dataframe
Yes, but very slow for this example when the data.frame gets large R DF - data.frame(x = rnorm(4), y = rnorm(4)) R foo - function(z, x, y) z[x] + z[y] R bar - function(x, y) x + y R system.time(x - apply(DF, 1, foo, x = x, y = y)) [1] 6.94 0.04 7.37 NA NA R system.time(y - with(DF, bar(x, y))) # from A. Liaw [1] 0.01 0.00 0.01 NA NA R all(x == y) [1] TRUE R system.time(z - mapply(bar, DF[,1], DF[,2])) Timing stopped at: 145.49 0.26 156.74 NA NA R # not sure why mapply doesn't work here... Best, Sundar Ramzi Feghali wrote: Another another neat way is: DF - data.frame(x=1:4, y=rep(1,4)) foo - function(z,x,y)z[x]+z[y] apply(DF,1,foo,x=x,y=y) 1 2 3 4 2 3 4 5 or DF - data.frame(x=1:4, y=rep(1,4)) foo - function(z)z-x+y foo(DF) [1] 2 3 4 5 Liaw, Andy [EMAIL PROTECTED] wrote: Another neat way is: with(DF, foo(x, w)) HTH, Andy -Original Message- From: peter leonard [mailto:[EMAIL PROTECTED] Sent: Sunday, June 08, 2003 4:35 PM To: [EMAIL PROTECTED] Subject: [R] Basic question on applying a function to each row of a dataframe Hi, I have a function foo(x,y) and a dataframe, DF, comprised of two vectors, x w, as follows : x w 1 1 1 2 2 1 3 3 1 4 4 1 etc I would like to apply the function foo to each 'pair' within DF e.g foo(1,1), foo(2,1), foo(3,1) etc I have tried apply(DF,foo) apply(DF[,],foo) apply(DF[DF$x,DF$w],foo) However, none of the above worked. Can anyone help ? Thanks in advance, Peter __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help -- Notice: This e-mail message, together with any attachments, cont... {{dropped}} __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help - [[alternate 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] Loops question
Ramzi Feghali wrote: Hello Mr Ripley, Whom is this mail adressed to? The mailing list R-help or Professor Ripley? i was waiting you to ask and if you don't mind if there is a fast way in R to do these loops made with R and that takes a week because my matrix is with thousand of rows. I got another method that is more fast and is utilised by many packages : interfacing with other languages, but it is only a question for information So you are going to optimize the following code? foo-function() {for (i in 1:(n-1)){ for (k in (i+1):n){ b_0 for (j in 1:m){ if(bin[i,j]==TRUEbin[k,j]==TRUE) b_b+1} print (b) } } } Thanks a lot Ramzi There seems to be a logical n x m matrix called bin. And the sollution of your problem should be solved within seconds by bin %*% t(bin) where the upper triangular matrix of the result consists of the results for all i and k in your loops. Uwe Ligges __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] questions about nls
Dear R users, I am new in R and I want to use the nls package to analyze some experimental data. The data is in the attached file data. It is the response Sav measured at different C0. Basically, the C0 is a function of C1, K2, and r, and the Sav is a function of C0, C1, K2, and r. The math equations are shown in the attached fileequations. The parameters K2 and r are the physical properties I want to get from the non-linear regression. The R codes I wrote is in the attached fileRcode. Basically, I wrote two functions. First , I calculated the C1 for different C0 at the estimated K2 and r using the binary search method and implemented it in the function fn1. Then for the calculated C1 and estimated K2 and r, I used the function fn2 to get the estimated Sav for different C0. nls was used to minimize the differences between the calculated Sav and observed Sav. When I run my R script, I got the error message step factor reduced below minFactor . If I changed the minFactor to zero, the nls continued but did not converge (exceed the maxiteration). I changed the tolerance to higher value, nls finished but the fitting is bad. From the published results, the best fitted values for K2 and r are 0.2237 and 1.296*10^7, respectively. I can't get these numbers using my R script. I know there are a lot math and R experts on the R mailing-list. I will appreciate it if anyone can tell me what is wrong in my R script or in the methods I used to get these parameters. Mei-Chu Lo C0 Sav 0.6766 6.0875 1.6165 6.4249 1.8796 6.5374 2.4436 7.025 4.8872 7.5125 5.3759 7.625 5.7518 7.8499 7.218 8.0749 9.2105 8.1125 12.7067 9.4624 12.5939 10.025 16.203 11.7125 17.2932 12.0124 18.1578 12.5749 21.5413 12.6875 21.5413 13.0625 fn1-function(C0,k2,r){ m-matrix(nrow=26,ncol=length(C0)) Ct-vector(numeric) C1-vector(numeric) dC-vector(numeric) j=1 repeat{ C1[j]-C0[j]/2;dC[j]-C0[j]/4 repeat{ m[1,j]=C1[j];m[2,j]=k2*C1[j]^2 m[26,j]=26*(k2/2)^25*r*C1[j]^26 i=3 repeat{ m[i,j]-i*(k2/2)^(i-1)*C1[j]^i i-i+1 if(i=26) break} Ct[j]-sum(m[,j]) if(abs(Ct[j]-C0[j])10^-6*C0[j]) break if(Ct[j]C0[j]) { C1[j]=C1[j]-dC[j];dC[j]-dC[j]/2} else{ C1[j]=C1[j]+dC[j];dC[j]-dC[j]/2} } j=j+1 if (jlength(C0)) break } C1 } fn2-function(C1,C0,k2,r){ m-matrix(nrow=26,ncol=length(C0)) m[1,]-5.8*(1-0.018*C0)*C1; m[2,]-2^(2/3)*5.8*(1-0.018*C0)*k2*C1^2 m[26,]-42*(1-0.019*C0)*26*(k2/2)^25*r*C1^26 j=3 repeat{ m[j,]-j^(5/3)*5.8*(1-0.018*C0)*(k2/2)^(j-1)*C1^j j=j+1 if(j=26)break } p-vector(numeric) k=1 repeat{ p[k]-sum(m[,k])/C0[k] k=k+1 if (klength(C0))break } p } # load data ass-read.table(data.txt,header=T) plot(ass$C0,ass$Sav,cex=1.5, xlab=C0(mg/ml),ylab=Sav(S)) # give initial guess of k2 and r k2=0.3;r=10^6 fit-nls(Sav~fn2(fn1(ass$C0,k2fit,rfit),ass$C0,k2fit,rfit),data=ass,start=list(k2fit=k2,rfit=r)) __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] Re: [BioC] installing XML in linux
On 9 Jun 2003, Douglas Bates wrote: You will need to install the Debian package libxml-dev before you can install the R XML package. and note that the package is really looking for 'libxml/parser.h'. On my machine, for example, parser.h is in /usr/include/libxml2/libxml/, so I have LIBXML_INCDIR=/usr/include/libxml2. -- Robert Burrows, PhD New England Biometrics [EMAIL PROTECTED] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] Re: [BioC] installing XML in linux
On Mon, Jun 09, 2003 at 06:03:55AM -0700, Mahbub Latif wrote: Thanks a lot Prof Bates. One more thing... What Debian packages should I install to install R packages Rgraphviz and rhdf5? For any given 'foo' I usually start by $ apt-cache search foo # may return lots $ apt-cache search libfoo # more focussed for libraries In this case, assumimg you want standard hdf5 (and not the mpich and lam variants that are also available) and $ apt-get install libhdf5-serial-dev graphviz but I do not know if the graphviz package has all you need. One day we may have suitable .deb packages for this to make installation of complex CRAN / bioC packages a litte easier. Dirk -- Don't drink and derive. Alcohol and analysis don't mix. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] Nested anovas
Did you try http://www.r-project.org/; - Search - R Site Search? I entered variants of nest and formula for nested model and got many hits, some of which mentioned |. I have not done that for a while, but I'm confident that some variant of this will work. hth. spencer graves Suzanne E. Blatt wrote: Hello, Is anyone out there doing nested anovas? I'm familiar with the notation to designate the nesting in SAS but can't seem to find any for R. How can I specify the nesting levels? Thanks, Suzanne __ McAfee VirusScan Online from the Netscape Network. Comprehensive protection for your entire computer. Get your free trial today! http://channels.netscape.com/ns/computing/mcafee/index.jsp?promo=393397 Get AOL Instant Messenger 5.1 free of charge. Download Now! __ [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] questions about nls
There is too much here for me to parse it right now. Did you walk through the code line by line looking at the output anc checking what it did at each step? If you don't get a satisfactory reply from someone else and you can't figure it out from a line-by-line analysis of your code, please simplify the question: Make up a toy problem that illustrates in only a few lines what you don't understand. You are more likely to get the answer you want if it takes seconds rather than minutes to understand your question. hth. spencer graves Mei-Chu Lo wrote: Dear R users, I am new in R and I want to use the nls package to analyze some experimental data. The data is in the attached file data. It is the response Sav measured at different C0. Basically, the C0 is a function of C1, K2, and r, and the Sav is a function of C0, C1, K2, and r. The math equations are shown in the attached fileequations. The parameters K2 and r are the physical properties I want to get from the non-linear regression. The R codes I wrote is in the attached fileRcode. Basically, I wrote two functions. First , I calculated the C1 for different C0 at the estimated K2 and r using the binary search method and implemented it in the function fn1. Then for the calculated C1 and estimated K2 and r, I used the function fn2 to get the estimated Sav for different C0. nls was used to minimize the differences between the calculated Sav and observed Sav. When I run my R script, I got the error message step factor reduced below minFactor . If I changed the minFactor to zero, the nls continued but did not converge (exceed the maxiteration). I changed the tolerance to higher value, nls finished but the fitting is bad. From the published results, the best fitted values for K2 and r are 0.2237 and 1.296*10^7, respectively. I can't get these numbers using my R script. I know there are a lot math and R experts on the R mailing-list. I will appreciate it if anyone can tell me what is wrong in my R script or in the methods I used to get these parameters. Mei-Chu Lo C0 Sav 0.6766 6.0875 1.6165 6.4249 1.8796 6.5374 2.4436 7.025 4.8872 7.5125 5.3759 7.625 5.7518 7.8499 7.218 8.0749 9.2105 8.1125 12.7067 9.4624 12.5939 10.025 16.203 11.7125 17.2932 12.0124 18.1578 12.5749 21.5413 12.6875 21.5413 13.0625 fn1-function(C0,k2,r){ m-matrix(nrow=26,ncol=length(C0)) Ct-vector(numeric) C1-vector(numeric) dC-vector(numeric) j=1 repeat{ C1[j]-C0[j]/2;dC[j]-C0[j]/4 repeat{ m[1,j]=C1[j];m[2,j]=k2*C1[j]^2 m[26,j]=26*(k2/2)^25*r*C1[j]^26 i=3 repeat{ m[i,j]-i*(k2/2)^(i-1)*C1[j]^i i-i+1 if(i=26) break} Ct[j]-sum(m[,j]) if(abs(Ct[j]-C0[j])10^-6*C0[j]) break if(Ct[j]C0[j]) { C1[j]=C1[j]-dC[j];dC[j]-dC[j]/2} else{ C1[j]=C1[j]+dC[j];dC[j]-dC[j]/2} } j=j+1 if (jlength(C0)) break } C1 } fn2-function(C1,C0,k2,r){ m-matrix(nrow=26,ncol=length(C0)) m[1,]-5.8*(1-0.018*C0)*C1; m[2,]-2^(2/3)*5.8*(1-0.018*C0)*k2*C1^2 m[26,]-42*(1-0.019*C0)*26*(k2/2)^25*r*C1^26 j=3 repeat{ m[j,]-j^(5/3)*5.8*(1-0.018*C0)*(k2/2)^(j-1)*C1^j j=j+1 if(j=26)break } p-vector(numeric) k=1 repeat{ p[k]-sum(m[,k])/C0[k] k=k+1 if (klength(C0))break } p } # load data ass-read.table(data.txt,header=T) plot(ass$C0,ass$Sav,cex=1.5, xlab=C0(mg/ml),ylab=Sav(S)) # give initial guess of k2 and r k2=0.3;r=10^6 fit-nls(Sav~fn2(fn1(ass$C0,k2fit,rfit),ass$C0,k2fit,rfit),data=ass,start=list(k2fit=k2,rfit=r)) __ [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] Questions for package ts prediction
On Mon, 2003-06-09 at 03:56, Prof Brian Ripley wrote: You have a scoping problem: the predict method needs to find the data: please supply an explicit newdata argument. Your first example works because `y' happens to be globally visible and be the right object. Thanks. That makes sense. On 9 Jun 2003, zhu wang wrote: I am trying to write a function to return prediction values using package ts. I have written three different versions since I am not sure what's wrong with my func2. func and func1 return the same results.But func1 and func2 don't. In particular, the only difference between func1 and func2 is the function variable name being y and data, respectively. But running the last line of the following script will give the message: Error in ts(x): object is not a matrix. I am confused. Also, could somebody kindly let me what's the answer if any for the following sunspot example from the package help: data(sunspot) (sunspot.ar - ar(sunspot.year)) # why not just sunspot.ar - ar(sunspot.year) ? Have you tried it? Please do so, and you will learn the difference! This idiom is widely used in the R help files. I have tried it before but now I know the trick here. predict(sunspot.ar, n.ahead=25) -- zhu wang [EMAIL PROTECTED] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] Re: [BioC] installing XML in linux
Robert Burrows wrote: On 9 Jun 2003, Douglas Bates wrote: You will need to install the Debian package libxml-dev before you can install the R XML package. and note that the package is really looking for 'libxml/parser.h'. On my machine, for example, parser.h is in /usr/include/libxml2/libxml/, so I have LIBXML_INCDIR=/usr/include/libxml2. Thanks for pointing out the fact that one can set LIBXML_INCDIR and it is definitely a good thing. In this particular case, however, I hope that the configuration should work automatically without any additional manual settings. The configure script should be using xml2-config --cflags and that should return -I/usr/include/libxml2 and libxml/parser.h is then duly found within that directory. If this isn't the case, I'd be happy to find out and change things. Unfortunately, conflicts between libxml and libxml2, use of libxml/ and gnome-xml/, etc... make this a complex configuration. -- Robert Burrows, PhD New England Biometrics [EMAIL PROTECTED] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help -- ___ Duncan Temple Lang[EMAIL PROTECTED] Bell Labs, Lucent Technologiesoffice: (908)582-3217 700 Mountain Avenue, Room 2C-259 fax:(908)582-3340 Murray Hill, NJ 07974-2070 http://cm.bell-labs.com/stat/duncan __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
RE: [R] Introductory Resources
[EMAIL PROTECTED] wrote ... 1) I want a test suite for R. I noted in the messages (Date: Mon Feb 24 2003 - 22:18:03 EST) that Prof Ripley wrote Well, R itself has lots of tests in its test suite (see directory tests in the sources) packages... but I was too stupid to find them. Q1: Can someone provide directions to this test suite that even an idiot can follow? This is an area where R really shines. I still test my code with Splus (3.3) but I switched to maintaining it in R because of the test (quality assurance) facilities for packages. The tests indicated above are in $(RHOME)/tests and can be run when you build from source (in Linux or Unix) with ./configure make make check I expect what you really want is your own test suite. To do this put your code in an R package (see Writing R Extensions) and put the tests in the tests directory. If the test code executes stop() or ends abnormally for another reason, then checking or building the package will fail with an error message indicating a problem. If you are really interested in this I will expand sometime. Paul Gilbert Head Statistician/Statisticien en chef, Department of Monetary and Financial Analysis, /Département des Études monétaires et financiers, Bank of Canada/Banque du Canada 234 Wellington St., Ottawa, Canada K1A 0G9 (613) 782-7346 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] questions abtou nls
The attached file equation in my previous e-mail does not go through (It was created by Microsoft word). I re-type the equations in Notepad and attach it now. Mei-Chu Lo C0 = C1 + K2*C1^2+ [Sum(i=3 to 25)(i*(K2/2)^(i-1) *C1^i]+ 26 * (K2/2)^25*r*C1^26 Sav = {5.8*(1-0.0018*C0)*C1+2^(2/3)*5.8*(1-0.018*C0)*K2*C1^2 +[sum(i=3 to 25)(i^(5/3)*5.8*(1-0.018*C0)*(K2/2)^(i-1)*C1^i] + 42* (1-0.019*C0)*26*(K2/2)^25*r*C1^26}/C0 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] OT: Xemacs config help
Marc, In reply to your post of /Thu Feb 6 17:34:05 MET 2003/ re: Xemacs config help you wrote, in part: I'd like to ps pretty print R code to study after I've written it. I can do this with no problem from xemacs under Windows, but in Linux, I keep getting the same error: ps-postscript-code-directory isn't set properly I've scoured the net, have posted to the xemacs group and all I can find are others with the same problem, but no solution. Your query came to my attention while I was searching for a solution to this problem. My search also came up dry. I am also using xemacs 21.4 on RedHat Linux 8.0 and made the following fix: Where the symbol ~ represents your home directory-- add the next line to the file ~/.xemacs/init.el *(setq ps-postscript-code-directory /usr/share/emacs/21.2/etc)* Best Regards, Craig Bourne [[alternate HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Re: [BioC] installing XML in linux
apt-get install libhdf5-serial-dev To install libhdf5 what I am getting $ apt-get install libhdf5-serial-dev libhdf5-serial-dev: Depends: libhdf5-serial (= 1.4.5-2) but it is not going to be installed then I try... $ apt-get install libhdf5-serial libhdf5-serial: Depends: libc6 (= 2.3.1-1) but 2.2.5-11.5 is to be installed That means I need libc6 2.3.1-1 or higher. How can I get that? Mahbub. unless you really want the MPI or PVM versions. best, -tony -- A.J. Rossini / [EMAIL PROTECTED] / [EMAIL PROTECTED] Biomedical/Health Informatics and Biostatistics, University of Washington. Biostatistics, HVTN/SCHARP, Fred Hutchinson Cancer Research Center. FHCRC: 206-667-7025 (fax=4812)|Voicemail is pretty sketchy/use Email CONFIDENTIALITY NOTICE: This e-mail message and any attachments may be confidential and privileged. If you received this message in error, please destroy it and notify the sender. Thank you. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] Ordering long vectors
On Mon, 9 Jun 2003, Prof Brian Ripley wrote: On Sun, 8 Jun 2003, Göran Broström wrote: On Sun, 8 Jun 2003, Thomas Lumley wrote: On Sat, 7 Jun 2003, [ISO-8859-1] Göran Broström wrote: I need to order a long vector of integers with rather few unique values. This is very slow: I think the culprit is src/main/sort.c: orderVector1 /* Shell sort isn't stable, but it proves to be somewhat faster to run a final insertion sort to re-order runs of ties when comparison is cheap. */ This also explains: aa-sample(rep(1:10,5)) system.time( order(aa, 1:length(aa))) [1] 3.67 0.01 3.68 0.00 0.00 system.time( order(aa)) ^C Timing stopped at: 49.33 0.01 49.34 0 0 which is perhaps the simplest work-around :). There is a simple and very much faster solution if you don't care about ordering of ties: system.time(ord4 - sort(x, method=quick, index.return = TRUE)$ix) Thanks, it is indeed much faster; about 7 times on my example. That is in the help, and I am very surprised you failed to find it. Well, I found it, but I was cooled off by the somewhat faster than Shellsort and poor performance in the worst case, together with the fact that I only needed the order and not the sorted vector. Maybe a more positive description is in order :-). Especially on the help page of 'order'. G. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R]parameters calculated by lda and qda
Hi, I am not sure if I remember correctly a mail last week to this messageboard where Prof Ripley said Fisher's discriminants are not used in lda within R and that Rao's are used instead, is this true? I looked in MASS (1996) but I couldn't find the reference to this in chapter 12 (multivariate analysis), the (help)lda doesn't tell me either. Could anyone tell me where can I find out this information as I am also interested in the method for qda? Thanks Anne __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R]parameters calculated by lda and qda
Power, Anne Marie wrote: Hi, I am not sure if I remember correctly a mail last week to this messageboard where Prof Ripley said Fisher's discriminants are not used in lda within R and that Rao's are used instead, is this true? I looked in MASS (1996) but I couldn't find the reference to this in chapter 12 (multivariate analysis), Chapter 12 Classification in the 4th edition (2002) of MASS has a description and gives references. the (help)lda doesn't tell me either. Could anyone tell me where can I find out this information as I am also interested in the method for qda? See the recent edition of MASS for qda() as well. Uwe Ligges Thanks Anne __ [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] converting by to a data.frame?
Hi, Don: Thanks for your suggestion to use do.call in my get.Index. I discovered that your version actually produces cosmetically different answers in R 1.6.3 and S-Plus 6.1 for Windows. Fortunately, in the context, this difference was unimportant. Since yours is faster, it is clearly superior. To check my understanding, I generalized my toy example as follows: by.df - data.frame(A=rep(c(A1, A2), each=3), + B=rep(c(B1, B2), each=3), C=rep(c(C1, C2), each=3), + x=1:6, y=rep(0:1, length=6)) With this, your get.Index produced the following in R 1.6.2: get.Index - function(x, INDICES) do.call('paste',c(x[INDICES],sep=':')) get.Index(by.df, c(A, B, C)) [1] A1:B1:C1 A1:B1:C1 A1:B1:C1 A2:B2:C2 A2:B2:C2 A2:B2:C2 In S-Plus 6.1 for Windows, I got the following: get.Index - function(x, INDICES) do.call(paste, c(x[INDICES], sep = :)) get.Index(by.df, c(A, B, C)) [1] 1:1:1 1:1:1 1:1:1 2:2:2 2:2:2 2:2:2 Fortunately, this difference is unimportant in this context, as by.to.data.frame produces the same answer in both cases. Moreover, your answer converts to a single call to paste, which means that it should be faster. For someone who understands do.call, your version is also easier to read. Thanks again for your help. Spencer Graves ## Don MacQueen wrote: Glad to hear it was helpful. You can also use the do.call trick for the paste indices business. Try get.Index - function(x, INDICES) do.call('paste',c(x[INDICES],sep=':')) This works because a data frame is actually a list, albeit a special kind of list, and do.call() wants a list for its second arg. -Don ### Thanks to Thomas Lumley, Sundar Dorai-Raj, and Don McQueen for their suggestions. I need the INDICES as part of the output data.frame, which McQueen's solution provided. I generalized his method as follows: by.to.data.frame - function(x, INDICES, FUN){ # Split data.frame x on x[,INDICES] # and lapply FUN to each data.frame subset, # returning a data.frame # # Internal functions get.Index - function(x, INDICES){ Ind - as.character(x[,INDICES[1]]) k - length(INDICES) if(k 1) Ind - paste(Ind, get.Index(x, INDICES[-1]), sep=:) Ind } FUN2 - function(data., INDICES, FUN){ vec - FUN(data.) Vec - matrix(vec, nrow=1) dimnames(Vec) - list(NULL, names(vec)) cbind(data.[1,INDICES], Vec) } # Combine INDICES Ind - get.Index(x, INDICES) # Apply ...: Do the work. Split - split(x, Ind) byFits - lapply(Split, FUN2, INDICES, FUN) # Convert to a data.frame do.call('rbind',byFits) } Applying this to my toy problem produces the following: by.df - data.frame(A=rep(c(A1, A2), each=3), + B=rep(c(B1, B2), each=3), x=1:6, y=rep(0:1, length=6)) by.to.data.frame(by.df, c(A, B), function(data.)coef(lm(y~x, data.))) A B (Intercept) x A1:B1 A1 B1 0.333 -1.517960e-16 A2:B2 A2 B2 0.667 3.282015e-16 Thanks for the assistance. I can now tackle the real problem that generated this question. Best Wishes, Spencer Graves Don MacQueen wrote: Since I don't have your by.df to test with I may not have it exactly right, but something along these lines should work: byFits - lapply(split(by.df,paste(by.df$A,by.df$B)), FUN=function(data.) { tmp - coef(lm(y~x,data.)) data.frame(A=unique(data.$A), B=unique(data.$B), intercept=tmp[1], slope=tmp[2]) }) byFitsDF - do.call('rbind',byFits) That's assuming I've got all the closing parantheses in the right places, since my email software (Eudora) doesn't do R syntax checking! This approach can get rather slow if by.df is big, or when the computations in FUN are extensive (or both). If by.df$A has mode character (as opposed to being a factor), then replacing A=unique(data.$A) with A=I(unique(data.$A)) might improve performance. You want to avoid character to factor conversions when using an approach like this. -Don At 2:54 PM -0700 6/5/03, Spencer Graves wrote: Dear R-Help: I want to (a) subset a data.frame by several columns, (b) fit a model to each subset, and (c) store a vector of results from the fit in the columns of a data.frame. In the past, I've used for loops do do this. Is there a way to use by? Consider the following example: byFits - by(by.df, list(A=by.df$A, B=by.df$B), + function(data.)coef(lm(y~x, data.))) byFits A: A1 B: B1 (Intercept) x 3.33e-01 -1.517960e-16 A: A2 B: B1 NULL A: A1 B: B2 NULL
Re: [R] Basic question on applying a function to each row of a dataframe
On Mon, 9 Jun 2003, Sundar Dorai-Raj wrote: Yes, but very slow for this example when the data.frame gets large Indeed. However, it works when the function is not already vectorised, and if the function is already vectorised it is unnecessary. -thomas __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] estimate the number of clusters
Dear All, I am using Silhouette to estimate the number of clusters in a microarray dataset. Initially, I used the iris data to test my piece of code as follows: library(cluster) data(iris) mydata-iris[,1:4] maxk-15# at most 15 clusters myindex-rep(0,maxk) # hold the si values for each k clusters mdist-1-cor(t(mydata)) #dissimlarity mdist-as.dist(mdist) for(k in 2:maxk) { hc-diana(mdist,diss =TRUE, stand = FALSE) si-silhouette.default(cutree(as.hclust(hc),k=k),mdist) myindex-summary(si)$avg.width } myk-rev(order(myindex))[1] #select the number of k clusters with the #largest si value I met the following problems: for(k in 2:maxk) + { + hc-diana(mdist,diss =TRUE, stand = FALSE) + si-silhouette.default(cutree(as.hclust(hc),k=k),mdist) + myindex-summary(si)$avg.width + } Error in [-(*tmp*, iC, sil_width, value = s.i) : number of items to replace is not a multiple of replacement length In addition: Warning messages: 1: longer object length is not a multiple of shorter object length in: b.i - a.i 2: number of rows of result is not a multiple of vector length (arg 2) in: cbind(mmm, as.vector(each)) Could any one help me how I can solve the problems??? Your kind help is highly appreciated!! ping - Post your free ad now! Yahoo! Canada Personals [[alternate HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Appending elements to an array
Hi all, I am having a bit of trouble with the array structure of R. What I want to do is dynamically add/remove elements to an array. For example: Let's say I have created an array: myArray - array(c(3,8), dim=c(1,2)) myArray [,1] [,2] [1,]38 And I now want to, for example, push an element (5,6) on to this array so it will read: [,1] [,2] [1,]38 [2,]56 And then pop the first element of the array so the array now reads: [,1] [,2] [1,]56 How would I do this? So far I've only read how to create an array if you know the dimensions beforehand, but I've been unable to find how to dynamically add/remove elements to/from an array. I've figured out how to do this with lists and vectors, but not with arrays. For this particular structure I need to work with arrays. Any help on how to do something like this would be much appreciated. Thanks, Jonck __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] Appending elements to an array
On Mon, 9 Jun 2003, Jonck van der Kogel wrote: Hi all, I am having a bit of trouble with the array structure of R. What I want to do is dynamically add/remove elements to an array. For example: Let's say I have created an array: myArray - array(c(3,8), dim=c(1,2)) myArray [,1] [,2] [1,]38 And I now want to, for example, push an element (5,6) on to this array so it will read: [,1] [,2] [1,]38 [2,]56 And then pop the first element of the array so the array now reads: [,1] [,2] [1,]56 How would I do this? So far I've only read how to create an array if you know the dimensions beforehand, but I've been unable to find how to dynamically add/remove elements to/from an array. I've figured out how to do this with lists and vectors, but not with arrays. For this particular structure I need to work with arrays. myArray - array(c(3,8), dim=c(1,2)) #add to the bottom myArray - cbind(myArray, c(5,6)) # display the top myArray # remove the top myArray - myArray[-1,] -thomas __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] understanding eigen(): getting non-normalized eigenvectors
Hi, dear R pros I try to understand eigen(). I have seen, that eigen() gives the eigenvectors normalized to unit length. What shall I do to get the eigenvectors not normalized to unit length? E.g. take the example: A [,1] [,2] V1 0.7714286 -0.2571429 V2 -0.4224490 0.1408163 Calculating eigen(A) by hand gives the eigenvectors (example from Backhaus, multivariate analysis): 0.77143 and 0.25714 -0.42245 0.14082 but even eigen(solve(Derror)%*%Dtreat, symmetric = FALSE, EISPACK =TRUE) which according to ?eigen should not necessarily give the normalized eigenvectors give the vectors (such as eigen()): $vectors [,1] [,2] [1,] 0.8770963 0.3162278 [2,] -0.4803146 0.9486833 - how can I replicate the result we get by hand (I ask because for students it is nice to see the same results with R as the results written in textbooks, derived manually? Thanks a lot Christoph -- Christoph Lehmann [EMAIL PROTECTED] University Hospital of Clinical Psychiatry __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] ESRI shapefiles and EMME/2 packages
I just uploaded two packages to CRAN. shapefiles_0.1.tar.gz - functions to read and write ESRI shapefiles (including dbfs) emme2_0.1.tar.gz - functions to read binary data from an EMME/2 databank data (EMME/2 is a transportation modeling program) Please let me know if you find any bugs or have some suggestions. Thanks. Regards, Benjamin Stabler Transportation Planning Analysis Unit Oregon Department of Transportation 555 13th Street NE, Suite 2 Salem, OR 97301 Ph: 503-986-4104 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] looking for Prof Bates' file
Actually the file is there when I check today. The server was down yesterday. Thanks, viet Message: 31 Date: 09 Jun 2003 11:44:25 +0200 From: Peter Dalgaard BSA [EMAIL PROTECTED] Subject: Re: [R] looking for Prof Bates' file To: Viet Nguyen,,, [EMAIL PROTECTED] Cc: [EMAIL PROTECTED] Message-ID: [EMAIL PROTECTED] Content-Type: text/plain; charset=us-ascii Viet Nguyen,,, [EMAIL PROTECTED] writes: Hello I'm reading up on fitting truncated Weibull distribution to data. There are posts in 2002 that point to this presentation by Prof Bates: http://www.stat.wisc.edu/~bates/JSM2001.pdf but now the file is not there. I can't find it anywhere else, Google doesn't have a cached copy for it. Could someone please give me a copy of this file, if they have it? Thanks and regards, viet. Looks like www.stat.wisc.edu is no longer identical to Doug's desktop machine. The file is still there, but I'm not too sure whether he would want us to tell the world how to get to it, so perhaps we should wait for him to reply... -- 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
Re: [R] Re: [BioC] installing XML in linux
Mahbub Latif [EMAIL PROTECTED] wrote: To install libhdf5 what I am getting $ apt-get install libhdf5-serial-dev libhdf5-serial-dev: Depends: libhdf5-serial (= 1.4.5-2) but it is not going to be installed then I try... $ apt-get install libhdf5-serial libhdf5-serial: Depends: libc6 (= 2.3.1-1) but 2.2.5-11.5 is to be installed That means I need libc6 2.3.1-1 or higher. How can I get that? libc6-2.3.1 is in testing (sarge), while libc6-2.2.5 is still in the stable (woody) distribution of debian. What confuses me a bit from your message is that libhdf5 is _not_ in woody (libhdf4 is), and that the version that you are trying to install is not even in sarge. On a sarge system, the current package libhdf5-serial-dev has version number 1.4.4-0.2 and depends on libc6-2.2.5-13, not the newer libc6-2.3.1. In other words, it looks like your /etc/apt/sources.list points towards the unstable (sid) distribution, while your debian system has not been upgraded. My 1¢ piece of advice would be to steer away from sid until you master the intricacies of the debian packaging system. Aside from fixing this mess, you may just retrieve the sources and compile them against your current libc. -- Philippe __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] R hits the New Zealand Herald
The following story about R appeared this morning in NZ's largest circulation newspaper, the New Zealand Herald. http://www.nzherald.co.nz/storydisplay.cfm?storyid=3506567 Chris Wild Department of Statistics The University of Auckland * BIRTHPLACE OF R !! * [[alternate HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help