Re: [R] postgres

2003-02-12 Thread ripley
RPgSQL etc seem now to be unsupported. As for ODBC, there is some information in the README *but* setting up ODBC is not an R question. Once you have a working connection (tested by say isql -v foo), then odbcConnect(foo) will work too. There is lots of documentation on the unixODBC site.

Re: [R] rpart v. lda classification.

2003-02-12 Thread ripley
On Tue, 11 Feb 2003, Rolf Turner wrote: I've been groping my way through a classification/discrimination problem, from a consulting client. There are 26 observations, with 4 possible categories and 24 (!!!) potential predictor variables. I tried using lda() on the first 7 predictor

[R] multinomial conditional logit models

2003-02-12 Thread John Hendrickx
A little while ago, I asked for some help in estimating multinomial conditional logit models. If you restructure your data, you can estimate a multinomial logit as a conditional logit model. This gives flexibility in imposing constraints on your dependent variable. One application is to estimate a

[R] Receiving Rnews by mail

2003-02-12 Thread Ronaldo Reis Jr.
Hi all, I possible to receive Rnews automatically in my E-mail? Thanks Ronaldo -- O cerebro e um orgao maravilhoso. Comeca a funcionar assim que voce se levanta da cama e nao p ra ate voce chegar ao escritorio. -- Robert Frost -- | //|\\

Re[2]: [R] Interpolation

2003-02-12 Thread Remigijus Lapinskas
Many thanks to all who replied to my e-mail. My problem was that I had not known about the approx function. By the way, if I have x - c(1990,1994,1995,1997), is there an automated way to fill in the gaps, i.e., to get c(1991,1992,1993,1996)? Cheers, Remigijus Wednesday, February 12, 2003,

[R] models for square tables

2003-02-12 Thread John Hendrickx
I've posted a sample file for estimating loglinear models for square tables (mobility models) at http://www.xs4all.nl/~jhckx/mcl/R/ Comments and suggestions are welcome. John Hendrickx __ [EMAIL PROTECTED] mailing list

Re: [R] Interpolation

2003-02-12 Thread Sundar Dorai-Raj
Remigijus Lapinskas wrote: Many thanks to all who replied to my e-mail. My problem was that I had not known about the approx function. By the way, if I have x - c(1990,1994,1995,1997), is there an automated way to fill in the gaps, i.e., to get c(1991,1992,1993,1996)? Try this: R x -

[R] GLMMGibbs crashes R

2003-02-12 Thread Stefano Calza
Hi everybody. I'm trying to use the GLMMGibbs package (R 1.6.2, Linux/Debian 3.0) with the data scottish.lip.cancer, as described in the paper by J.Miles and D. Clayton. The problem is that the code at pag 18 crasches R: ***sparse_rd***too few elements (column 32) /n Process R exited

[R] rbind.data.frame: character comverted to factor

2003-02-12 Thread Markus Jäntti
Dear All, on rbind:ing together a number of data.frames, I found that character variables are converted into factors. Since this occurred for a data identifier, it was a little inconvenient and, to me, unexpected. (The help page explains the general procedure used. I also found that on forming a

RE: [R] Dynamic Linear Models for Times Series - Implemented?

2003-02-12 Thread Thomas Lumley
On Wed, 12 Feb 2003, Gavin Simpson wrote: Hi, Following an off-list reply to my original post, I realised that I hadn't really provided very much information for you to work with. So here's a second attempt: Following West Harrison (1989) and Pole et al. (1994) a DLM is defined as:

[R] `UNPROTECT' and `return'

2003-02-12 Thread Timur Elzhov
Dear R experts, In all R functions written in C one must unprotect result before returning them: { ... UNPROTECT(1) /* unprotecting `ans' */ return ans; } Why does one shure that memory occupied by `ans' won't be used by R immediately after unprotecting? Ok, is the next

Re: [R] Dynamic Linear Models for Times Series - Implemented?

2003-02-12 Thread Paul Gilbert
Gavin I am not familiar with the two texts you mention that define DLM, but I think since at least twenty years prior to those texts, the term has been used to mean the state-space model you describe, and also ARMA models, transfer function models, and possibly some other representations of

Re: [R] `UNPROTECT' and `return'

2003-02-12 Thread ripley
This _is_ covered in `Writing R Extensions'. What is unsafe is to allow a gc to occur. On Wed, 12 Feb 2003, Timur Elzhov wrote: In all R functions written in C one must unprotect result before returning them: { ... UNPROTECT(1) /* unprotecting `ans' */ return

Re: [R] rbind.data.frame: character comverted to factor

2003-02-12 Thread Rolf Turner
Brian Ripley wrote: d1 - data.frame(id=I(letters[1:20]), x = runif(20)) d2 - data.frame(id=I(paste(letters[1:20],letters[1:20],sep=)), x=rexp(20)) d3 - rbind(d1, d2) which of course works! Why ``of course''? It seems to me that there is no ``of course'' about it. It is

[R] Poesia

2003-02-12 Thread kjetil brinchmann halvorsen
Hola! Sorry for going off-topic, but here are something I found on the web yesterday - an explanation of statistics in poetic form: First, you see your data for what they seem to be Then, you ask them for the truth - are you what you seem to me? You see with broad expanse

[R] Three questions - loading(sourcing a file??), summaries to bitmaps,and error messages

2003-02-12 Thread Tony Vargas
R helpers, I have three questions about R and was wondering if I could get some help. First, I am trying to send the output of a a summary command - summary(info) to a bitmap file instead of the default device. Anybody know how to send a summary to a bitmap device? Second - this question is

[R] Various Errors using Survey Package

2003-02-12 Thread Thompson, Trevor
Hi, I have been experimenting with the new Survey package. Specifically, I was trying to use some of the functions on the public-use survey data from NHIS (2000 Sample Adult file). Error 1): The first error I get is when I try to specify the complex survey design.

[R] Debugging in R

2003-02-12 Thread Charles Sugnet
Hello *, Is it possible in R to set breakpoints and conditional breakpoints in the browser mode? I've discovered debug() to debug functions, but I find stepping throw loops or having to go back and put in browse() calls a little tedious. I would really like to be able to set a breakpoint while

[R] rl_callback_read_char error on Solaris 7

2003-02-12 Thread Yuelin Li
This question is about compiling R-1.6.2 from source on a Sparcstation-20 machine running Solaris 7, gcc-2.95.3 and readline-4.2 (both gcc and readline from sunfreeware.com). I searched the r-help archive and the FAQ. The errors seem to be related to the binary readline libraries, and I

[R] Matrix formatting

2003-02-12 Thread Tim Sharac
Hi R-users: I have a data formatting question. I have a data set that looks something like this: foo.dat - cbind(c(NA, 1, 2, 3, 4, 5), c(NA, NA, 0, 10 ,20, 30)) What I have: [,1] [,2] [1,] NA NA [2,]1 NA [3,]20 [4,]3 10 [5,]4 20 [6,]5 30 I want to

Re: [R] Three questions - loading(sourcing a file??), summaries to bitmaps, and error messages

2003-02-12 Thread Peter Dalgaard BSA
Tony Vargas [EMAIL PROTECTED] writes: R helpers, I have three questions about R and was wondering if I could get some help. First, I am trying to send the output of a a summary command - summary(info) to a bitmap file instead of the default device. Anybody know how to send a summary to

Re: [R] Matrix formatting

2003-02-12 Thread Jerome Asselin
Would this do what you want? Cheers, Jerome foo.dat - cbind(c(NA, 1, 2, 3, 4, 5), c(NA, NA, 0, 10 ,20, 30)) apply(foo.dat,2,function(x) x[order(as.logical(x))]) [,1] [,2] [1,]10 [2,]2 10 [3,]3 20 [4,]4 30 [5,]5 NA [6,] NA NA On Wednesday 12 February

[R] R CMD check .

2003-02-12 Thread Robin Hankin
Hello everybody I'm writing a package and am trying to get it past R CMD check . It has no C or Fortan code, just R code. R CMD check . reports that the examples don't work: r:Davies% R CMD check . * checking for working latex ... OK * using log directory

Re: [R] Matrix formatting

2003-02-12 Thread Roger Peng
Is apply(foo.dat, 2, sort, na.last = TRUE) what you want? -roger ___ UCLA Department of Statistics [EMAIL PROTECTED] http://www.stat.ucla.edu/~rpeng On 12 Feb 2003, Tim Sharac wrote: Hi R-users: I have a data formatting question. I have a data set that looks

Re: [R] Matrix formatting

2003-02-12 Thread Jerome Asselin
Sorry, this is probably better if you want to preserve the order of the numbers in each column. apply(foo.dat,2,function(x) x[order(is.na(x))]) Jerome On Wednesday 12 February 2003 13:05, Jerome Asselin wrote: Would this do what you want? Cheers, Jerome foo.dat - cbind(c(NA, 1, 2, 3,

Re: [R] R CMD check .

2003-02-12 Thread David Scott
On Thu, 13 Feb 2003, Robin Hankin wrote: Hello everybody I'm writing a package and am trying to get it past R CMD check . It has no C or Fortan code, just R code. R CMD check . reports that the examples don't work: r:Davies% R CMD check . * checking for working latex ... OK *

Re: [R] rl_callback_read_char error on Solaris 7

2003-02-12 Thread Don MacQueen
At 9:03 PM + 2/12/03, [EMAIL PROTECTED] wrote: On Wed, 12 Feb 2003, Yuelin Li wrote: This question is about compiling R-1.6.2 from source on a Sparcstation-20 machine running Solaris 7, gcc-2.95.3 and readline-4.2 (both gcc and readline from sunfreeware.com). This is not an R question,

Re: [R] R CMD check .

2003-02-12 Thread Jeff Gentry
* checking for CRLF line endings in C sources/headers ... OK * creating Davies-Ex.R ... OK * checking examples ... ERROR Running examples failed. r:Davies% How do I tell which .Rd file is the problem? (I suspect that the examples aren't actually being executed because most of them take a

Re: [R] R CMD check .

2003-02-12 Thread David Scott
Sorry, there was a typo in my answer: On Thu, 13 Feb 2003, David Scott wrote: . . . lots of stuff cut . . examples aren't actually being executed because most of them take a long time and the error is pretty much instantaneous). Rcmd check makes a directory called libname.Rcheck. Inside

Re: [R] rl_callback_read_char error on Solaris 7

2003-02-12 Thread Peter Dalgaard BSA
Yuelin Li [EMAIL PROTECTED] writes: This question is about compiling R-1.6.2 from source on a Sparcstation-20 machine running Solaris 7, gcc-2.95.3 and readline-4.2 (both gcc and readline from sunfreeware.com). I had some problems with the sunfreeware readline too when I tried 1.5 years

[R] legend

2003-02-12 Thread Wilkinson, Mark
I think I'm missing something tonight in the usage of 'legend': plot(0, type=n) legend(locator(1), month.abb[1:5], pch=15, col=1:5) gives me something similar to what I want. But legend(locator(1), month.abb[1:5], fill=T, col=1:5) gives me 5 black boxes. What am I doing wrong? Thank you,

Re: [R] R CMD check .

2003-02-12 Thread Robin Hankin
Hello everybody thanks for your help with my R CMD check problem. The file Davies-Ex.Rout does indeed contain a transcript of the examples. I've found the .Rd file that contains the problem. The relevant command executes perfectly when I cut-and-paste it onto the R command line, but not when

Re: [R] legend

2003-02-12 Thread Hedderik van Rijn
legend(locator(1), month.abb[1:5], fill=T, col=1:5) gives me 5 black boxes. Try: legend(locator(1), month.abb[1:5], fill=1:5) What am I doing wrong? The T is interpreted as (equal to?) a 1, so you're requesting the boxes to be filled with color 1, which is black. - Hedderik.

[R] sorting in lmList object

2003-02-12 Thread emb7
Hi all, Forgive me if this is an obvious one I want to make a plot of confidence intervals from an lmList object with a collection of simple linear models (lm(y~x)) using: plot(intervals(mylmList)) and sort the plot by increasing mean values for the intercept. Is there an easy way of doing

Re: [R] Various Errors using Survey Package

2003-02-12 Thread Thomas Lumley
On Wed, 12 Feb 2003, Thompson, Trevor wrote: Hi, I have been experimenting with the new Survey package. Specifically, I was trying to use some of the functions on the public-use survey data from NHIS (2000 Sample Adult file). Error 1): The first error I get is when I try to specify the

[R] multi-color plot

2003-02-12 Thread Saurav Pathak
hi all, i am trying to make multi-color plots. that is, i generally use, plot(x, y, type=n) text(x, y, labels=class) here, the vector class denotes the class of each point. there are usually 3-4 classes of points. how may i display the different classes in different colors?

Re: [R] multi-color plot

2003-02-12 Thread Sundar Dorai-Raj
Saurav Pathak wrote: hi all, i am trying to make multi-color plots. that is, i generally use, plot(x, y, type=n) text(x, y, labels=class) here, the vector class denotes the class of each point. there are usually 3-4 classes of points. how may i display the different classes in different

Re: [R] multi-color plot

2003-02-12 Thread Saurav Pathak
Thus spake Sundar Dorai-Raj: + R plot(x=1:26,y=rep(0:1,13),type=n) + R text(x=1:26,y=0.5,labels=letters,col=terrain.colors(26)) now i get it. what i am doing now is: R mycol - cls R mycol[mycol==0] - red R mycol[mycol==1] - blue R mycol[mycol==2] - green R plot(x, y, type=n) R text(x, y,

[R] Urgent Help on Cumulant computation on two R.V's

2003-02-12 Thread Feng Zhang
Hey Now I am going to check the independence of random variables using cumulant function. So if R has such package or functions to calculate the sample cumulant of a random vector? Thanks a lot. Fred __ [EMAIL PROTECTED] mailing list

Re: [R] k- means cluster analysis

2003-02-12 Thread ripley
On Wed, 12 Feb 2003, Sundar Dorai-Raj wrote: Ngayee J Law wrote: Hi all, I am trying to run the k-means cluster analysis using the function kmeans in the package cluster. I think it's the one in package mva. [...] However, choosing bad starting values could cause kmeans to crash, as

RE: [R] Interpolation

2003-02-12 Thread Adaikalavan Ramasamy
Here is another way of doing it: x - c(1990, 1994, 1995, 1997) all.x - seq(min(x), max(x)) complementary.x - setdiff(all.x, x) -Original Message- From: Sundar Dorai-Raj [mailto:[EMAIL PROTECTED]] Sent: Wednesday, February 12, 2003 9:39 PM To: Remigijus Lapinskas Cc: [EMAIL PROTECTED]