Re: [R] Incompatibility with VGAM

2005-05-30 Thread Uwe Ligges
Thomas Yee wrote: Hello, It seems that if glm used a namespace then the conflict would be avoided? No. I already descibed why *not*, and I descibed how to work with stats' family functions even if VGAM has been loaded. Please reread my former message and tell me what was unclear. You might

Re: [R] joining files after canonical correlation

2005-05-30 Thread Uwe Ligges
Brett Stansfield wrote: Dear R, I recently did a canonical correlation analysis on two subsets of data (location and weather). So I now have canonical scores for location and weather. but I'd now like to do a scatterplot matrix using the pairs statement. Is there a way to somehow join

[R] how to invert the matrix with quite small eigenvalues

2005-05-30 Thread huang min
Dear all, I encounter some covariance matrix with quite small eigenvalues (around 1e-18), which are smaller than the machine precision. The dimension of my matrix is 17. Here I just fake some small matrix for illustration. a-diag(c(rep(3,4),1e-18)) # a matrix with small eigenvalues

[R] Menu Language

2005-05-30 Thread Hafner, Reinhold (Risklab)
Hello, In my windows xp system the regional settings are set to German. However, I still would like to use the English version of R with English menu items and texts. RWinEdt tells me this is necessary in order to use the mdi framework. How can this be achieved with the 2.1 version. In 2.0.1 this

RE: [R] how to invert the matrix with quite small eigenvalues

2005-05-30 Thread Ted Harding
On 30-May-05 huang min wrote: Dear all, I encounter some covariance matrix with quite small eigenvalues (around 1e-18), which are smaller than the machine precision. The dimension of my matrix is 17. Here I just fake some small matrix for illustration. a-diag(c(rep(3,4),1e-18)) # a

[R] how to singlify entries

2005-05-30 Thread zhihua li
hi netters I have a rather simple question. I have a data frame with two variables X and Y, both of which are factors. X has 100 levels while Y has 10 levels only. The data frame has 100 rows in all, so for X the values are unique, and Y has many replicate values. Now I wanna reduce the

Re: [R] Menu Language

2005-05-30 Thread Uwe Ligges
Hafner, Reinhold (Risklab) wrote: Hello, In my windows xp system the regional settings are set to German. However, I still would like to use the English version of R with English menu items and texts. RWinEdt tells me this is necessary in order to use the mdi framework. How can this be achieved

Re: [R] how to invert the matrix with quite small eigenvalues

2005-05-30 Thread huang min
Maybe I should state more clear that I define b to get the orthogonal matrix bb$vectors. We also can define diag(b)-diag(b)+100, which will make the eigenvalues of b much bigger to make sure the orthogonal matrix is reliable. My intention is to invert the covariance matrix to perform some

Re: [R] how to singlify entries

2005-05-30 Thread Petr Pikal
Hallo On 30 May 2005 at 9:15, zhihua li wrote: hi netters I have a rather simple question. I have a data frame with two Well, I do not understand you simple question fully. You have something like that dat-data.frame(X=1:100, Y=sample(1:10,10)) dat$X-factor(dat$X) dat$Y-factor(dat$Y)

Re: [R] Chars as numbers

2005-05-30 Thread Martin Maechler
Package 'sfsmisc' has a function AsciiToInt() and a few useful related R-code only functions such as chars8bit(); see the help pages once you've installed and attached the package. But do note that these things do depend on the encoding, as Uwe Ligges has already told you. Things work fine for

Re: [R] joining files after canonical correlation

2005-05-30 Thread Duncan Murdoch
Brett Stansfield wrote: Dear R, I recently did a canonical correlation analysis on two subsets of data (location and weather). So I now have canonical scores for location and weather. but I'd now like to do a scatterplot matrix using the pairs statement. Is there a way to somehow join

Re: [R] how to invert the matrix with quite small eigenvalues

2005-05-30 Thread Ted Harding
On 30-May-05 huang min wrote: Maybe I should state more clear that I define b to get the orthogonal matrix bb$vectors. OK. Certainly bbv-bb$vectors is close to orthogonal: bbv%*%bbv differs from the unit matrix only in that the off-diagonal terms are O(10^(-16)). We also can define

Re: [R] sapply following using by with a list of factors

2005-05-30 Thread Frank E Harrell Jr
McClatchie, Sam (PIRSA-SARDI) wrote: Background: OS: Linux Mandrake 10.1 release: R 2.0.0 editor: GNU Emacs 21.3.2 front-end: ESS 5.2.3 - Colleagues I am having some trouble extracting results from the function by, used to average variables in a data.frame first

Re: [R] Re: Vector Manipulation

2005-05-30 Thread Gabor Grothendieck
x[cumsum(x!=0)!=0] or x[!!cumsum(!!x)] will also do it. On 5/30/05, ManojW [EMAIL PROTECTED] wrote: OK...x[min(which(x!=0)):length(x)] does the trick! I guess the coffee is slowly but surely working! . Manoj - Original Message - From: ManojW To: R-help Sent: Monday, May

Re: [R] how to singlify entries

2005-05-30 Thread Charles Plessy
On Mon, May 30, 2005 at 09:15:32AM +, zhihua li wrote : hi netters I have a rather simple question. I have a data frame with two variables X and Y, both of which are factors. X has 100 levels while Y has 10 levels only. The data frame has 100 rows in all, so for X the values are

[R] You can appreciate the reduced prices on rxdrugs and check our weekly specials.

2005-05-30 Thread michale baillet
Check our cybershop and you will uncover how others cut their expenses on rxmeds significantly. Check our cybershop for quality rxdrugs on male tissue dysfunction, pain, affliction, relaxants for muscles, man's care, over-wt. and other discomforts. Our company provides customers quick and

Re: [R] how to singlify entries

2005-05-30 Thread Gabor Grothendieck
On 5/30/05, Charles Plessy [EMAIL PROTECTED] wrote: On Mon, May 30, 2005 at 09:15:32AM +, zhihua li wrote : hi netters I have a rather simple question. I have a data frame with two variables X and Y, both of which are factors. X has 100 levels while Y has 10 levels only. The data

[R] values of bars in barplot

2005-05-30 Thread luc tardieu
Hi, I couldn't find how to have the values written on the top of each bar in a barplot. When using hist(), it is possible to use labels=T, but this option does not seem to exist for barplot(). Is there a trick I could use to do that ? Thanks to all Luc

Re: [R] how to singlify entries

2005-05-30 Thread Petr Pikal
On 30 May 2005 at 21:56, Charles Plessy wrote: On Mon, May 30, 2005 at 09:15:32AM +, zhihua li wrote : hi netters I have a rather simple question. I have a data frame with two variables X and Y, both of which are factors. X has 100 levels while Y has 10 levels only. The data

Re: [R] values of bars in barplot

2005-05-30 Thread Gabor Grothendieck
On 5/30/05, luc tardieu [EMAIL PROTECTED] wrote: Hi, I couldn't find how to have the values written on the top of each bar in a barplot. When using hist(), it is possible to use labels=T, but this option does not seem to exist for barplot(). Is there a trick I could use to do that ?

Re: [R] values of bars in barplot

2005-05-30 Thread joerg van den hoff
luc tardieu wrote: Hi, I couldn't find how to have the values written on the top of each bar in a barplot. When using hist(), it is possible to use labels=T, but this option does not seem to exist for barplot(). Is there a trick I could use to do that ? Thanks to all Luc

Re: [R] values of bars in barplot

2005-05-30 Thread Petr Pikal
On 30 May 2005 at 15:19, luc tardieu wrote: Hi, I couldn't find how to have the values written on the top of each bar in a barplot. When using hist(), it is possible to use labels=T, but this option does not seem to exist for barplot(). Hallo From help page Value: A numeric

[R] rbind wastes memory

2005-05-30 Thread lutz . thieme
Hello everybody, if I try to (r)bind a number of large dataframes I run out of memory because R wastes memory and seems to forget to release memory. For example I have 10 files. Each file contains a large dataframe ds (3500 cols by 800 rows) which needs ~20 MB RAM if it is loaded as the only

Re: [R] rbind wastes memory

2005-05-30 Thread Douglas Bates
[EMAIL PROTECTED] wrote: Hello everybody, if I try to (r)bind a number of large dataframes I run out of memory because R wastes memory and seems to forget to release memory. For example I have 10 files. Each file contains a large dataframe ds (3500 cols by 800 rows) which needs ~20 MB

Re: [R] rbind wastes memory

2005-05-30 Thread Duncan Murdoch
[EMAIL PROTECTED] wrote: Hello everybody, if I try to (r)bind a number of large dataframes I run out of memory because R wastes memory and seems to forget to release memory. For example I have 10 files. Each file contains a large dataframe ds (3500 cols by 800 rows) which needs ~20 MB RAM if

[R] Egarch

2005-05-30 Thread Friedrich, Andreas (dit)
Hallo, I got to implement parameter estimation for Egarch(p,q) processes. My question: has anyone done that so far or are there implementations for other asymmetric garch models? Andreas [[alternative HTML version deleted]] __

[R] Reference Card?

2005-05-30 Thread Martin Klaffenboeck
Hello! For LaTeX I found a reference Card at http://www.cs.ualberta.ca/~c603/LaTeX_docs/Symbol_Source/latex_symbols.pdf Is there something available for R? thanks, Martin __ R-help@stat.math.ethz.ch mailing list

Re: [R] Reference Card?

2005-05-30 Thread Kevin Wang
Hi, Have you tried looking under Documentation - Contributed, under CRAN? Kev Martin Klaffenboeck wrote: Hello! For LaTeX I found a reference Card at http://www.cs.ualberta.ca/~c603/LaTeX_docs/Symbol_Source/latex_symbols.pdf Is there something available for R? thanks, Martin

Re: [R] Reference Card?

2005-05-30 Thread TEMPL Matthias
Hello! For LaTeX I found a reference Card at http://www.cs.ualberta.ca/~c603/LaTeX_docs/Symbol_Source/latex _symbols.pdf Is there something available for R? Hello Martin, See the reference cards on http://cran.r-project.org/other-docs.html Best, Matthias PLEASE do read the posting

[R] Formatting numbers with a limited amount of digits consistently

2005-05-30 Thread Henrik Andersson
I have tried to get signif, round and format to display numbers like these consistently in a table, using e.g. signif(x,digits=3) 17.01 18.15 I want 17.0 18.2 Not 17 18.2 Why is the last digit stripped off in the case when it is zero! Is this a feature of R or did I miss something?

Re: [R] how to singlify entries

2005-05-30 Thread Charles Plessy
On Mon, May 30, 2005 at 09:09:27AM -0400, Gabor Grothendieck wrote : Try using reshape, e.g. if dd is your data frame: reshape(dd, dir = wide, idvar = F1, timevar = F2, varying = list(c(VX,VY))) Thank you very much, and to Petr Pikal too. Reshape is exactly what I had forgotten. Now

Re: [R] rbind wastes memory

2005-05-30 Thread Roger D. Peng
Rather than 'rbind' in a loop, try putting your dataframes in a list and then doing something like 'do.call(rbind, list.of.data.frames)'. -roger [EMAIL PROTECTED] wrote: Hello everybody, if I try to (r)bind a number of large dataframes I run out of memory because R wastes memory and seems to

[R] a question about read.marrayRaw

2005-05-30 Thread shizhu zang
Hi, everyone, I ran the test case in Introduction to the Bioconductor marrayInput package first. When I ran the session in (read.marrayRaw) I got the following error messages. Other sessions work well. The functions and errors are labled with blue and red respectively. mraw - read.Spot(path

Re: [R] Incompatibility with VGAM

2005-05-30 Thread Thomas Lumley
On Mon, 30 May 2005, Thomas Yee wrote: Hello, It seems that if glm used a namespace then the conflict would be avoided? No. glm does use a namespace, so this can't be true. Remember that R passes arguments by value, and consider glm(y~x, family=poisson()) The namespace ensures that

[R] How to access to sum of dissimilarities in CLARA

2005-05-30 Thread Amir Safari
Dear All , Since dissimilarity is one of quality measures in clustering , I'm trying to access to the sum of dissimilarity as a whole measure. But after running my data using CLARA I obtain : 1128 dissimilarities, summarized : Min. 1st Qu. Median Mean 3rd Qu. Max. 0.033155

[R] R GUI for Linux?

2005-05-30 Thread White, Charles E WRAIR-Wash DC
I feel your pain. grin I am a new Linux user who has spent most of the weekend trying to get a functional R setup. When I installed Fedora Core 3 (FC3) on my home computer, I thought using R in a terminal would be a snap. I installed R using the rpm packages and tried to use it with the FC3

Re: [R] How to access to sum of dissimilarities in CLARA

2005-05-30 Thread TEMPL Matthias
Hello, #Example: data(xclara) p - clara(xclara,3) names(p) p$diss Best, Matthias Dear All , Since dissimilarity is one of quality measures in clustering , I'm trying to access to the sum of dissimilarity as a whole measure. But after running my data using CLARA I obtain : 1128

Re: [R] how to invert the matrix with quite small eigenvalues

2005-05-30 Thread Thomas Lumley
On Mon, 30 May 2005, huang min wrote: My intention is to invert the covariance matrix to perform some algorithm which is common in the estimating equations like GEE. In that case there is no benefit in being able to invert very extreme covariance matrices. The asymptotic approximations to the

[R] problem in R2.2-dev

2005-05-30 Thread ronggui
when i install the R2.2-dev(customed install,the with the chinese translation msg ,etc.),i found the tanslated GUI can not display correctly.totally unlike Chinese characters. __ R-help@stat.math.ethz.ch mailing list

Re: [R] Formatting numbers with a limited amount of digits consistently

2005-05-30 Thread Duncan Murdoch
Henrik Andersson wrote: I have tried to get signif, round and format to display numbers like these consistently in a table, using e.g. signif(x,digits=3) 17.01 18.15 I want 17.0 18.2 Not 17 18.2 Why is the last digit stripped off in the case when it is zero! signif() changes the value;

Re: [R] problem in R2.2-dev

2005-05-30 Thread Duncan Murdoch
ronggui wrote: when i install the R2.2-dev(customed install,the with the chinese translation msg ,etc.),i found the tanslated GUI can not display correctly.totally unlike Chinese characters. I don't think you've given nearly enough information for anyone to act on this. Please describe what

Re: [R] Formatting numbers with a limited amount of digits consistently

2005-05-30 Thread Gabor Grothendieck
On 5/30/05, Duncan Murdoch [EMAIL PROTECTED] wrote: Henrik Andersson wrote: I have tried to get signif, round and format to display numbers like these consistently in a table, using e.g. signif(x,digits=3) 17.01 18.15 I want 17.0 18.2 Not 17 18.2 Why is the

Re: [R] R GUI for Linux?

2005-05-30 Thread Sander Oom
Hi Charles, Warm felt sympathies for your struggles. I consider myself a happy GUI user and have also struggled with the 'command line' history and lack of out-of-the-box functionality associated with Linux. However, Linux does have many, many advantages over other OS's, so I will stick to

Re: [R] Formatting numbers with a limited amount of digits consistently

2005-05-30 Thread Duncan Murdoch
Gabor Grothendieck wrote: On 5/30/05, Duncan Murdoch [EMAIL PROTECTED] wrote: Henrik Andersson wrote: I have tried to get signif, round and format to display numbers like these consistently in a table, using e.g. signif(x,digits=3) 17.01 18.15 I want 17.0 18.2 Not 17 18.2 Why is the

[R] Unique arrangements of a vector

2005-05-30 Thread Tarmo Remmel
Dear List, Running on a PC (Windows 2000) with 256 MB RAM, Version R1.9.1 I have a relatively simple problem, which I can solve for relatively small datasets, but run into difficulties with larger ones. I believe that my approach is a hack rather than something elegant and I was hoping that

[R] FANNY function in R package cluster

2005-05-30 Thread Aamir M
Dear All, I am attempting to use the FANNY fuzzy clustering function in R (Kaufman Rousseeuw, 1990), found in the cluster package. I have run into a variety of difficulties; the two most crucial difficulties are enumerated below. 1. Where is the 'm' parameter in FANNY? In _Finding Groups in

RE: [R] R GUI for Linux?

2005-05-30 Thread John Fox
Dear Charles, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of White, Charles E WRAIR-Wash DC Sent: Monday, May 30, 2005 10:52 AM To: r-help@stat.math.ethz.ch Cc: [EMAIL PROTECTED] Subject: [R] R GUI for Linux? . . . Rcmdr: There are all

[R] persp, add lines/highlights

2005-05-30 Thread jjorgensen
Hello R-sters, I'm trying to add several lines to a response surface that I've plotted using persp(). I've tried lines() using the trans3d function but I've been unsuccessful in getting it to work (R v2.0.1). Essentially, I'm trying to highlight one or more of the surface wireframe lines in a

Re: [R] R GUI for Linux?

2005-05-30 Thread Patrick Connolly
On Mon, 30-May-2005 at 08:29PM +0200, Sander Oom wrote: | .. I have | given up on any instructions that tell me to run 'make'. Linux | distributions are just to[o] idiosyncratic for it to be worth the effort. That might be true of Linux distributions in general, but installing R is

[R] Piecewise Linear Regression

2005-05-30 Thread Abhyuday Mandal
Hi, I need to fit a piecewise linear regression. x = c(6.25,6.25,12.50,12.50,18.75,25.00,25.00,25.00,31.25,31.25,37.50,37.50,50.00,50.00,62.50,62.50,75.00,75.00,75.00,100.00,100.00) y =

Re: [R] Piecewise Linear Regression

2005-05-30 Thread Sean Davis
Abhyuday, There are a number of answers in the archives: http://www.google.com/u/newcastlemaths?q=piecewise+linear+regressionsa=Google+Search Do any of those meet your needs? Sean - Original Message - From: Abhyuday Mandal [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Sent:

[R] labels on map

2005-05-30 Thread jose silva
dear all: Im trying to obtain maps on R, under mapdata library, but I cannot define the labels. Here is an example: library(mapdata) map(worldHires, c(portugal,spain),ylim=c(34,46),xlim=c(-14,3.5)) axis(1,at=seq(-12,3,3)) axis(2) when I try the parameter xlab or ylab in axis, i get: parameter

Re: [R] persp, add lines/highlights

2005-05-30 Thread Paul Murrell
Hi [EMAIL PROTECTED] wrote: Hello R-sters, I'm trying to add several lines to a response surface that I've plotted using persp(). I've tried lines() using the trans3d function but I've been unsuccessful in getting it to work (R v2.0.1). Essentially, I'm trying to highlight one or more of

Re: [R] Piecewise Linear Regression

2005-05-30 Thread roger koenker
It is conventional to fit piecewise linear models by assuming Gaussian error and using least squares methods, but one can argue that median regression provides a more robust approach to this problem. You might consider the following fit: x = c

Re: [R] labels on map

2005-05-30 Thread Duncan Murdoch
jose silva wrote: dear all: Im trying to obtain maps on R, under mapdata library, but I cannot define the labels. Here is an example: library(mapdata) map(worldHires, c(portugal,spain),ylim=c(34,46),xlim=c(-14,3.5)) axis(1,at=seq(-12,3,3)) axis(2) when I try the parameter xlab or ylab in

RE: [R] labels on map

2005-05-30 Thread Mulholland, Tom
What makes you think that there is a ylab parameter? args(map) function (database = world, regions = ., exact = FALSE, boundary = TRUE, interior = TRUE, projection = , parameters = NULL, orientation = NULL, fill = FALSE, col = 1, plot = TRUE, add = FALSE, namesonly = FALSE, xlim =

Re: Re: [R] labels on map

2005-05-30 Thread jose silva
thanks Duncan, thats what I was looking for :) p.s. xlab in map() gives the same result as in axis() j. silva jose silva wrote: dear all: Im trying to obtain maps on R, under mapdata library, but I cannot define the labels. Here is an example: library(mapdata) map(worldHires,

RE: [R] R GUI for Linux?

2005-05-30 Thread White, Charles E WRAIR-Wash DC
John: Thank you for your interest. After more investigation I see that my problem is with linking tcltk to R. tcl 8.4.7-2 and tk 8.4.7-2 are tuned to fc3 (Fedora Core 3) and part of the standard installation. However, required file locations are different than what is expected by R. I set the

[R] Problem going back to a viewport with gridBase

2005-05-30 Thread Gabor Grothendieck
I am setting up base plots -- one in viewport A and and one in B. This part works fine. But if I go back to A after having done B and add horizontal lines it seems to not use the correct coordinates. How do I tell it to resume using A's coordinates? I am already using par(fig = gridFIG()) but

Re: [R] Formatting numbers with a limited amount of digits consistently

2005-05-30 Thread Gabor Grothendieck
On 5/30/05, Duncan Murdoch [EMAIL PROTECTED] wrote: Gabor Grothendieck wrote: On 5/30/05, Duncan Murdoch [EMAIL PROTECTED] wrote: Henrik Andersson wrote: I have tried to get signif, round and format to display numbers like these consistently in a table, using e.g. signif(x,digits=3)

[R] Trying to write a linear regression using MLE and optim()

2005-05-30 Thread Ajay Shah
I wrote this: # Setup problem x - runif(100) y - 2 + 3*x + rnorm(100) X - cbind(1, x) # True OLS -- lm(y ~ x) # OLS likelihood function -- ols.lf - function(theta, K, y, X) { beta - theta[1:K] sigma - exp(theta[K+1]) e - (y - X%*%beta)/sigma logl - sum(log(dnorm(e))) return(logl) }

[R] Installing RGL on SuSE 9.3

2005-05-30 Thread Shige Song
I am running R 2.10 on SuSE 9.3. When I tried to install the package rgl, I got error message: ... In file included from pixmap.cpp:13: pngpixmap.h: In static member function `static void PNGPixmapFormat::Load::info_callback(png_struct*, png_info*)': pngpixmap.h:149: error: invalid conversion