RE: [R] core dump during make check when building 64-bit R on Solaris8/9

2004-11-30 Thread Prof Brian Ripley
We know of at least two builds which successfully used libsunperf, one (mine) with Studio ONE 7 and one with 9. So this may well be a problem with your particular compiler suite. On Wed, 1 Dec 2004, CHAN Chee Seng wrote: Hi, After reading some of the posting in this list, I came across this pos

Re: [R] Relative subscripting

2004-11-30 Thread Uwe Ligges
(Ted Harding) wrote: On 30-Nov-04 Tobias Muhlhofer wrote: [...] I have monthly a dataset with, among other things, "marketcap", and "return". I want to multiply return by the marketcap of the previous month. How do I do this? In STATA (which I have used frequently in the past) I would simply use an

[Fwd: Re: [R] Kernel Fisher Discriminant in R?]

2004-11-30 Thread Gorden Jemwa
Ooops...I meant "formulations shouldn't be difficult" Original Message Subject: Re: [R] Kernel Fisher Discriminant in R? Date: Wed, 01 Dec 2004 08:05:38 +0200 From: Gorden Jemwa <[EMAIL PROTECTED]> To: [EMAIL PROTECTED] CC: [EMAIL PROTECTED] You could us the kernlab package as a

Re: [R] Kernel Fisher Discriminant in R?

2004-11-30 Thread Gorden Jemwa
You could us the kernlab package as a building block. Implementation of Kernel Fisher Discriminant Analysis (or Generalized Fisher Discriminant Analysis) formulation are fairly straightforward and should be difficult. Rgds, Gorden Date: Mon, 29 Nov 2004 16:57:28 -0500 From: "Huh, Seungho" <[EMAI

RE: [R] core dump during make check when building 64-bit R on Solaris8/9

2004-11-30 Thread CHAN Chee Seng
Hi, After reading some of the posting in this list, I came across this posting from: From: Peter Dalgaard Date: Fri 29 Oct 2004 - 08:02:40 EST Replying to Re: [R] Errors during make check He described a problem similar to mine, that build 64-bit R (I am building R version 2.0.1) with sunperf li

[R] tuning SVM's

2004-11-30 Thread stephenc
Hi I am doing this sort of thing: POLY: > > obj = best.tune(svm, similarity ~., data = training, kernel = "polynomial") > summary(obj) Call: best.tune(svm, similarity ~ ., data = training, kernel = "polynomial") Parameters: SVM-Type: eps-regression SVM-Kernel: polynomial

[R] Request for Gamma Frailty Simulation

2004-11-30 Thread rahgozar
Dear Sir/Madam, I need to use " Gamma Frailty Model " in a simulation study ,so I should replicate the program many times and get estimates of " THETA and its Variance " or "Variance-Covariance" matrix Or Variace of random effect in Gamma frailty model in Survival analysis at each time of

Re: [R] Combined variable names

2004-11-30 Thread Richard A. O'Keefe
Tobias Muhlhofer <[EMAIL PROTECTED]> wrote: I am trying to define a large number of variables through a loop construct. He wants to do for (i in 1:100) { assign(paste("v", i, sep=""), something or other...) } This is, of course, a FAQ. It's

[R] time data

2004-11-30 Thread Karla Meurk
I have a data set with date, time and rainfall (sample below) my questions are (a) is there a command like as.date for time? and (b) Can R smooth data 6 hourly with such unevenly spaced data? Thanks Carla 09/29/02 19:33 0 09/29/02 19:34 0 09/29/02 19:35 0 09/29/02

Re: [R] Creating a factor from a combination of vectors

2004-11-30 Thread Richard A. O'Keefe
Yves Brostaux <[EMAIL PROTECTED]> wrote: I want to produce a factor from a subset of the combination of two vectors. I have the vectors a et b in a data-frame : > df <- expand.grid(a=c(0, 5, 10, 25, 50), b=c(0, 25, 50, 100, 200)) ... and want to cr

RE: [R] predict.gam problem

2004-11-30 Thread Liaw, Andy
I thought the documentation you quoted is rather clear. The GAM is something like: g(Y) = sum f_j(x_j) so predict.gam() returns sum f_j(x_j) by default, and if type="response", g^-1(sum f_j(x_j)). This is similar to predict.glm, I believe. HTH, Andy > From: Jon Egil Strand > > > ---

Re: [R] New trellis settings

2004-11-30 Thread Deepayan Sarkar
On Tuesday 30 November 2004 17:01, David Hinds wrote: > On Tue, Nov 30, 2004 at 04:33:50PM -0600, Deepayan Sarkar wrote: > > > I've fiddled with these and can pretty much get what I want. But > > > without really understanding what I'm doing. It isn't clear to > > > me, for instance, why the "key

Re: [R] Combined variable names

2004-11-30 Thread Peter Dalgaard
Tobias Muhlhofer <[EMAIL PROTECTED]> writes: > I am trying to define a large number of variables through a loop construct. > > I have my loop variable i being cycled through 1:100 and I would like > the variables produced by this to be called > > vi (i.e. v1 v2 v3 etc) > > so, for example I'm g

Re: [R] Attn Heinz Tuechler: Re: problem with sp ecial characters ( ä,ö,ü)

2004-11-30 Thread Heinz Tuechler
At 10:27 30.11.2004 -0500, Duncan Murdoch wrote: >[I tried to send this message privately, but the return address >bounced.] > >I think this has been fixed in R-patched, but I doubt if the fix has >been tested in Win98. Could you please download a copy from >

Re: [R] New trellis settings

2004-11-30 Thread David Hinds
On Tue, Nov 30, 2004 at 04:33:50PM -0600, Deepayan Sarkar wrote: > > > > I've fiddled with these and can pretty much get what I want. But > > without really understanding what I'm doing. It isn't clear to me, > > for instance, why the "key" settings affect my plot when no key is > > drawn. > > I

Re: [R] nls() error: Object not found

2004-11-30 Thread Douglas Bates
Yang, Richard wrote: Dear R-helpers; Using nls() to fit a function,Rdum, defined below I stumbled on an error: "Error in eval(expr, envir, enclos): Object "s0" not found." The function Rdum is defined as Rdum <- deriv(~ h1 * (s0 + sl0*sl + sm0*sm + sp01*sp1 + sp02*sp2 + sp03*sp3+sp

Re: [R] Combined variable names

2004-11-30 Thread Jean Eid
see ?assign and ?get i,e instead ob below have something like assign(paste("v", i, sep=""), a[i:N]) and if you need to loop over calling them say get(paste("v", i, sep="")) of course typing v1 will call the variable... Jean On Tue, 30 Nov 2004, Tobias Muhlhofer wrote: > I am trying to defin

Re: [R] Combined variable names

2004-11-30 Thread Sean Davis
See the R-FAQ here: http://cran.r-project.org/doc/FAQ/R-FAQ.html Number 7.21. However, as the FAQ also points out, you might be better served using lists. a <- list() for (i in 1:100) { a[[i]] <- some stuff } Sean On Nov 30, 2004, at 5:42 PM, Tobias Muhlhofer wrote: I am trying to define a lar

[R] predict.gam problem

2004-11-30 Thread Jon Egil Strand
- Mac OS X 10.3 R 2.0.1 gam 0.9.2 - Greetings I am facing difficulties in prediction on using Trevor Hasties GAM package. How can one interpret the response from predict.gam? Documentation on the type-parameter: The default (link) produces predictions on the

Re: [R] Combined variable names

2004-11-30 Thread Thomas Lumley
On Tue, 30 Nov 2004, Tobias Muhlhofer wrote: I am trying to define a large number of variables through a loop construct. I have my loop variable i being cycled through 1:100 and I would like the variables produced by this to be called vi (i.e. v1 v2 v3 etc) Look at FAQ 7.21, which explains how to

[R] nls() error: Object not found

2004-11-30 Thread Yang, Richard
Dear R-helpers; Using nls() to fit a function,Rdum, defined below I stumbled on an error: "Error in eval(expr, envir, enclos): Object "s0" not found." The function Rdum is defined as Rdum <- deriv(~ h1 * (s0 + sl0*sl + sm0*sm + sp01*sp1 + sp02*sp2 + sp03*sp3+sp04*sp4) *

[R] Combined variable names

2004-11-30 Thread Tobias Muhlhofer
I am trying to define a large number of variables through a loop construct. I have my loop variable i being cycled through 1:100 and I would like the variables produced by this to be called vi (i.e. v1 v2 v3 etc) so, for example I'm going: for(i in 1:100) { <- a[i:N] # or whatever else you want

Re: [R] New trellis settings

2004-11-30 Thread Deepayan Sarkar
On Tuesday 30 November 2004 13:33, David Hinds wrote: > On Mon, Nov 29, 2004 at 11:23:35PM -0600, Deepayan Sarkar wrote: > > > Specifically, I have some code that in the past changed the > > > relative proportions of text versus the plot area using: > > > > > > trellis.par.set('fontsize', list(te

Re: [R] Reading a tab delimitted text

2004-11-30 Thread Douglas Bates
Srinivas Iyyer wrote: Dear group, I have a tab delimitted text file with 21 column and 1988 rows. When I read the file: D1 <- read.table(file="gist_data1.txt",sep="\t",header=T) dim(D1) [1] 1811 20 It reads only 1811 rows from that file. I saved this file as CSV file and I read it again using:

Re: [R] augPred with lme(...,subset=...)

2004-11-30 Thread Douglas Bates
Manuel Morales wrote: Hello, Is there a way to get augPred to work with lme if a subset statement is used? For example, if I modify the example from ?augPred.lme to include a subset statement, I get the following error: fm1 <- lme(Orthodont, random = ~1, subset=distance>19) augPred(fm1, length.out

RE: [R] Relative subscripting

2004-11-30 Thread Ted Harding
On 30-Nov-04 Tobias Muhlhofer wrote: > [...] > I have monthly a dataset with, among other things, "marketcap", > and "return". > > I want to multiply return by the marketcap of the previous month. > How do I do this? > > In STATA (which I have used frequently in the past) I would simply > use an

Re: [R] Random effects and non-linear models

2004-11-30 Thread Douglas Bates
Certo, Trevis wrote: I am trying to model growth over time for a number of companies (my unit of analysis). My dependent variable is a count (the number of acquisitions each firm made each year), and I'm having some trouble figuring out exactly how to model this. My assumption is that modeling the

Re: [R] Relative subscripting

2004-11-30 Thread Tobias Muhlhofer
OK, yeah. I did think of that in the meantime, but I was also wondering if there was some other convenience function to do it. But this will do in any case. Toby Jean Eid wrote: wouldn't it be return[2:NROW(return)]*marketcap[1:(NROW(return)-1)] (note that return is also a function in R so mayb

Re: [R] Relative subscripting

2004-11-30 Thread Kjetil Brinchmann Halvorsen
Tobias Muhlhofer wrote: Hi! I'm trying to do the following. I have monthly a dataset with, among other things, "marketcap", and "return". I want to multiply return by the marketcap of the previous month. How do I do this? In STATA (which I have used frequently in the past) I would simply use a

Re: [R] Relative subscripting

2004-11-30 Thread Jean Eid
wouldn't it be return[2:NROW(return)]*marketcap[1:(NROW(return)-1)] (note that return is also a function in R so maybe you should stay away from calling your variable return. Anyways if you have a data frame you can add another variable to it but attach an NA to the first element (since you are l

Re: [R] seriesMerge

2004-11-30 Thread Yasser El-Zein
I want to merge two timeSeries (union) and get a resulting timeSeries. using union or merge I got a data.frame and when I tried to convert it to timeSeries the data was corrupt. Is there a more straight forward way to take two timeSeries and merege them for a resulting timeSeries of the union of bo

Re: [R] bitmap (blank image) plot rendering without X11

2004-11-30 Thread Peter Dalgaard
"Director, Beracah Yankama" <[EMAIL PROTECTED]> writes: > Hi, > I am trying to generate plots on our unix (no X) server to be included > in our web-pages -> specifically pngs. > > Without X, png() obviously doesn't work, so I have been trying to > define the graphics print device bitmap as descri

Re: [R] Choice modelling (was:(no subject))

2004-11-30 Thread Kjetil Brinchmann Halvorsen
[EMAIL PROTECTED] wrote: Hello, I am trying to estimate a choice model with varying choice set for each individual. I would like to fit different kinds of model (logit ,nested logit, probit...). So far I have found that package *mnp* allows me to estimate a probit model with varying choice set. But

Re: [R] Package for multivariate binary logistic regression?

2004-11-30 Thread Kjetil Brinchmann Halvorsen
Lynne Baker wrote: I am trying to find out if someone has implemented a (McFadden-type) multivariate binary logistic regresssion package for R? From what I can tell, this is not available for R. Thank you, Lynne Baker [[alternative HTML version deleted]] __

[R] calling R from Java

2004-11-30 Thread Lana Schaffer
Greetings, I would like to call a small R function from Java (actually to Bioconductor). From a little reading, looks like I can only do this on a Unix system and not on a PC. Is this indeed the case? Does anyone have experience with calling a R function from Java? Lana Lana Schaffer Research S

[R] SJava

2004-11-30 Thread Yasser El-Zein
How can I increase the JVM's memory using teh SJava library? __ [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

[R] Relative subscripting

2004-11-30 Thread Tobias Muhlhofer
Hi! I'm trying to do the following. I have monthly a dataset with, among other things, "marketcap", and "return". I want to multiply return by the marketcap of the previous month. How do I do this? In STATA (which I have used frequently in the past) I would simply use an expression of the form

[R] bitmap (blank image) plot rendering without X11

2004-11-30 Thread Director, Beracah Yankama
Hi, I am trying to generate plots on our unix (no X) server to be included in our web-pages -> specifically pngs. Without X, png() obviously doesn't work, so I have been trying to define the graphics print device bitmap as described in the help: > bitmap(file="plot.png", type = "png256", height

[R] Reading a tab delimitted text

2004-11-30 Thread Srinivas Iyyer
Dear group, I have a tab delimitted text file with 21 column and 1988 rows. When I read the file: D1 <- read.table(file="gist_data1.txt",sep="\t",header=T) > dim(D1) [1] 1811 20 It reads only 1811 rows from that file. I saved this file as CSV file and I read it again using: > D1 <- read.csv(

[R] augPred with lme(...,subset=...)

2004-11-30 Thread Manuel Morales
Hello, Is there a way to get augPred to work with lme if a subset statement is used? For example, if I modify the example from ?augPred.lme to include a subset statement, I get the following error: fm1 <- lme(Orthodont, random = ~1, subset=distance>19) augPred(fm1, length.out = 2, level = c(0,1)

[R] Random effects and non-linear models

2004-11-30 Thread Certo, Trevis
I am trying to model growth over time for a number of companies (my unit of analysis). My dependent variable is a count (the number of acquisitions each firm made each year), and I'm having some trouble figuring out exactly how to model this. My assumption is that modeling the effect of time on thi

Re: [R] xy_plot

2004-11-30 Thread Mike Saunders
Sean: You need to provide more information on your problem for anyone to help you much. One of the "apply" functions working in tandem with "points" or "lines" (or others) can often do multiple item quickly, but again I don't know the context of your problem. Mike - Original Message -

Re: [R] impute missing values in correlated variables: transcan?

2004-11-30 Thread Jonathan Baron
On 11/30/04 13:21, Frank E Harrell Jr wrote: >Jonathan Baron wrote: >> I would like to impute missing data in a set of correlated >> variables (columns of a matrix). It looks like transcan() from >> Hmisc is roughly what I want. It says, "transcan automatically >> transforms continuous and catego

Re: [R] RecordPlot

2004-11-30 Thread Paul Murrell
Hi The main problem here is that you really shouldn't be poking around inside a "recordedplot" object. From the help page: WARNING: The format of recorded plots may change between R versions. Recorded plots should *not* be used as a permanent storage format for R plots. R wil

Re: [R] A basic question

2004-11-30 Thread Andrew Robinson
R has compiled with no problem in FreeBSD. It is also included in the ports, thanks to its FreeBSD maintainer. Andrew On Tue, Nov 30, 2004 at 06:58:57AM -0500, Kenneth wrote: > Hi R users: > > I want to know any experience compiling R in other LINUX distributions > besides FEDORA (Red Hat) or

Re: [R] New trellis settings

2004-11-30 Thread David Hinds
On Mon, Nov 29, 2004 at 11:23:35PM -0600, Deepayan Sarkar wrote: > > > > Specifically, I have some code that in the past changed the relative > > proportions of text versus the plot area using: > > > > trellis.par.set('fontsize', list(text=8)) > > > > which caused all text to get smaller, and th

[R] xy_plot

2004-11-30 Thread ebashi
R users; Does anyone have a recommendation for a faster plotting function rather than plot(x,y)? when I use plot(x,y) for a large number of variables, it take couple of moments to plot the graph? Sincerely, Sean __ [EMAIL PROTECTED] mailing list https

Re: [R] impute missing values in correlated variables: transcan?

2004-11-30 Thread Frank E Harrell Jr
Jonathan Baron wrote: I would like to impute missing data in a set of correlated variables (columns of a matrix). It looks like transcan() from Hmisc is roughly what I want. It says, "transcan automatically transforms continuous and categorical variables to have maximum correlation with the best

Re: [R] about the subscript assignment porblem

2004-11-30 Thread Thomas Lumley
On Tue, 30 Nov 2004, kh zu wrote: Hi, Dear, I have a problem for the new version 2.0 of R When I put this code in it give me the error: uis$ivhx3[uis$ivhx>0] <- 1*(uis$ivhx==3) Error: NAs are not allowed in subscripted assignments but this would happen in version 1.91 What would happen? In versio

[R] Course ****R/S System: Advanced Programming, Seattle***December 20-21

2004-11-30 Thread elvis
R/S Advanced Programming course in Seattle on December 20-21 Please check out the full description and Agenda of the 2-day class on the website: www.xlsolutions-corp.com/Radv.htm And let us know if we should hold seats for you. Ask for group discount! Here's the outline: Course outline: - O

Re: [R] A basic question

2004-11-30 Thread asemeria
On gentoo: emerge R and dowload,configuration, compilation of all you need to working with base R is performed. Best! A.S. Alessandro Semeria Models and Simulations Laboratory Montecatini Environmental Research Center (Edison Group), Via Ciro Menotti 48, 4802

[R] RecordPlot

2004-11-30 Thread Jean Coursol
I want to do a zoom with recordPlot(). I have problems with lists. (R-2.0.1 patched 2004-11-30 , various linux). I have problems with RecordPlot class structure. > plot(1:10) > saveP <- recordPlot() > dev.off() > sx <- saveP[[1]][[2]][[2]] > saveP[[1]][[2]][[2]] <- sx Error in "[[<-"(`*tmp*`,

Re: [R] impute missing values in correlated variables: transcan?

2004-11-30 Thread Jonathan Baron
On 11/30/04 11:23, roger koenker wrote: >At the risk of stirring up a hornet's nest , I'd suggest that >means are dangerous in such applications. A nice paper >on combining ratings is: Gilbert Bassett and Joseph Persky, >Rating Skating, JASA, 1994, 1075-1079. Here is the abstract, which seems

[R] about the subscript assignment porblem

2004-11-30 Thread kh zu
Hi, Dear, I have a problem for the new version 2.0 of R When I put this code in it give me the error: uis$ivhx3[uis$ivhx>0] <- 1*(uis$ivhx==3) Error: NAs are not allowed in subscripted assignments but this would happen in version 1.91 there are missing data in the data set and it is represente

Re: [R] impute missing values in correlated variables: transcan?

2004-11-30 Thread roger koenker
At the risk of stirring up a hornet's nest , I'd suggest that means are dangerous in such applications. A nice paper on combining ratings is: Gilbert Bassett and Joseph Persky, Rating Skating, JASA, 1994, 1075-1079. url:www.econ.uiuc.edu/~rogerRoger Koenker email [EMA

[R] (no subject)

2004-11-30 Thread anne berthouly
Dear all, I’ve got problems concerning crossed random effects and nested effects using the function “lme”: I have two crossed random effects and I would like to nest each of them in a fixed factor and also test for the effects of these fixed factors. At the moment my model looks like this: Lme1

Re: [R] Using mimR

2004-11-30 Thread Uwe Ligges
Andreas Cordes wrote: Hi, I am trying to use the mimR Package. According to its help pages it needs the RDCOMClient package to run. When I try to evaluate a model I get > m.sat<-mim("..",data=gm.dat) Error in mim.cmd("clear; clear output") : couldn't find function "COMCreate" > m2.sat<-emfit(m

Re: [R] adding regression curve to xyplot

2004-11-30 Thread Carlisle Thacker
The axes are switched in xyplot. Oceanographers like to have temperature as the vertical axis and salinity horizontal. But temperature is used to predict salinity. Confusing, and makes everything complicated. So panel.abline(thislm) would give the wrong line, as would type=c("p","r"). And neit

[R] Google for scientists: search engine

2004-11-30 Thread Hanke, Alex
A new way to search for documentation on your favourite science topic. http://scholar.google.com/ Alex Hanke Department of Fisheries and Oceans St. Andrews Biological Station 531 Brandy Cove Road St. Andrews, NB Canada E5B 2L9 [[alternative HTML version de

[R] impute missing values in correlated variables: transcan?

2004-11-30 Thread Jonathan Baron
I would like to impute missing data in a set of correlated variables (columns of a matrix). It looks like transcan() from Hmisc is roughly what I want. It says, "transcan automatically transforms continuous and categorical variables to have maximum correlation with the best linear combination of

Re: [R] adding regression curve to xyplot

2004-11-30 Thread Deepayan Sarkar
On Tuesday 30 November 2004 10:18, Austin, Matt wrote: > Inside the panel, you have to only use the parts of x and y for each > unique combination of > factor(lonLabels[whichLon100])*factor(latLabels[whichLat100]) > > > > So use "thislm <- lm(x[subscripts]~y[subscripts])" in the panel > function.

Re: [R] 2k-factorial design with 10 parameters

2004-11-30 Thread Bob Wheeler
On Tue, 30 Nov 2004, Sven wrote: I'd like to apply a 2^k factorial design with k=10 parameters. Obviously this results in a quite long term for the model equation due to the high number of combinations of parameters. How can I specify the equation for the linear model (lm) without writing all

Re: [R] adding regression curve to xyplot

2004-11-30 Thread Deepayan Sarkar
On Tuesday 30 November 2004 10:24, Prof Brian Ripley wrote: > On Tue, 30 Nov 2004, Carlisle Thacker wrote: > > Dear R-listers, > > > > It seems that predict() behaves differently within panel.xyplot. > > Am I doing something stupid? > > You fitted the model lm(x ~ y) and supplied new values for t,

Re: [R] adding regression curve to xyplot

2004-11-30 Thread hadley wickham
If you only want a simple linear regression, you might also want to try xyplot(t~s|factor(lonLabels[whichLon100])*factor(latLabels[whichLat100]), +   data=P100,pch=".", type =c("p","r") - type "r" will automatically fit and plot the regression line for you. Hadley

Re: [R] adding regression curve to xyplot

2004-11-30 Thread Deepayan Sarkar
On Tuesday 30 November 2004 10:09, Carlisle Thacker wrote: > Dear R-listers, > > It seems that predict() behaves differently within panel.xyplot. Am > I doing something stupid? > > Thanks, > > Carlisle > > First, without xyplot(): > > lmtest <- lm(t~s,data=subset(P100,whichLon100==1 & whichLat100=

[R] Using mimR

2004-11-30 Thread Andreas Cordes
Hi, I am trying to use the mimR Package. According to its help pages it needs the RDCOMClient package to run. When I try to evaluate a model I get > m.sat<-mim("..",data=gm.dat) Error in mim.cmd("clear; clear output") : couldn't find function "COMCreate" > m2.sat<-emfit(m.sat) Error in toMIM(mim$

Re: [R] adding regression curve to xyplot

2004-11-30 Thread Prof Brian Ripley
On Tue, 30 Nov 2004, Carlisle Thacker wrote: Dear R-listers, It seems that predict() behaves differently within panel.xyplot. Am I doing something stupid? You fitted the model lm(x ~ y) and supplied new values for t, not y. Using panel.abline would be a bit easier: just call panel.abline(thislm)

[R] combinations min problem

2004-11-30 Thread adrian mincu
Hello, as I was not yet able to find a good solution to the following problem, I'd appreciate some help. I am working on finding the minimum of a list of vectors, each comprised of x parameters to include all combinations of n parameters taken k at a time, x>k thx. and rgds. a/m __

RE: [R] adding regression curve to xyplot

2004-11-30 Thread Austin, Matt
Inside the panel, you have to only use the parts of x and y for each unique combination of factor(lonLabels[whichLon100])*factor(latLabels[whichLat100]) So use "thislm <- lm(x[subscripts]~y[subscripts])" in the panel function. Matt Austin Statistician Amgen One Amgen Center Drive M/S 24-2-

Re: [R] How to know if a bug was recognised

2004-11-30 Thread Heinz Tuechler
Thank you for your answer. BTW for a beginner like me it seems very difficult to ask a question without being criticised. Maybe that for a beginner it is not that easy to find an issue which, as you say, "Indeed, no one reported the problem for 2.0.0 at all, including all the alpha and beta test v

RE: [R] How to know if a bug was recognised

2004-11-30 Thread Heinz Tuechler
Dear Henrik Bengtsson! Thank you for your kind answer. I am sorry that at the time of writing my inital message the last version of R was 2004-11-15. As soon as possible I will try the new version. with many thanks Heinz Tüchler At 14:01 30.11.2004 +0100, you wrote: >> -Original Message---

Re: [R] Attn Heinz Tuechler: Re: problem with sp ecial characters ( ä,ö,ü)

2004-11-30 Thread Heinz Tuechler
Thank you for the information. I did already download the 27-11-2004 version and I have the intention to try it within a day and report on the result, if this is of interest. My address ([EMAIL PROTECTED]) should work but I will check that too. Heinz Tüchler At 10:27 30.11.2004 -0500, you wrote:

Re: [R] Package for multivariate binary logistic regression?

2004-11-30 Thread Prof Brian Ripley
You'll have to help us a bit here. `McFadden' comes up in connection with conditional logistic regression, which package survival supports. Is that what you mean? On Tue, 30 Nov 2004, Lynne Baker wrote: I am trying to find out if someone has implemented a (McFadden-type) multivariate binary l

[R] adding regression curve to xyplot

2004-11-30 Thread Carlisle Thacker
Dear R-listers, It seems that predict() behaves differently within panel.xyplot. Am I doing something stupid? Thanks, Carlisle First, without xyplot(): > lmtest <- lm(t~s,data=subset(P100,whichLon100==1 & whichLat100==1)) > lmtest Call: lm(formula = s ~ t, data = subset(P100, whichLon100 == 1

[R] (no subject)

2004-11-30 Thread jruiz99
Hello, I am trying to estimate a choice model with varying choice set for each individual. I would like to fit different kinds of model (logit ,nested logit, probit...). So far I have found that package *mnp* allows me to estimate a probit model with varying choice set. But for estimation of

[R] Package for multivariate binary logistic regression?

2004-11-30 Thread Lynne Baker
I am trying to find out if someone has implemented a (McFadden-type) multivariate binary logistic regresssion package for R? From what I can tell, this is not available for R. Thank you, Lynne Baker [[alternative HTML version deleted]] __

Re: [R] How to know if a bug was recognised

2004-11-30 Thread Duncan Murdoch
On Tue, 30 Nov 2004 13:24:43 +0100, Heinz Tuechler <[EMAIL PROTECTED]> wrote : >Hello! > >A problem with special characters seemed to me to be a bug. I sent a mail >to [EMAIL PROTECTED] concerning the problem (see below). >How can I find out, if this is considered as a bug or an error of myself? >

Re: [R] plotting data in non-orthogonal coords.

2004-11-30 Thread Deepayan Sarkar
On Tuesday 30 November 2004 04:06, Andreas Franke wrote: [...] > It would be nice if one could just plot data given as F(x,y) where > you supply x and y for every data point seperatly so that you dont > need any specific grid. This is one solution, but it's not exactly what you want: library(l

RE: [R] Info

2004-11-30 Thread Prof Brian Ripley
Also, factanal() does not do `principal component analysis', and it may well be that the data are inappropriate for factor analysis if they are appropriate for PCA. If PCA was really intended, prcomp() and princomp() are appropriate tools. On Tue, 30 Nov 2004, Andy Bunn wrote: It looks like fact

[R] Attn Heinz Tuechler: Re: problem with special characters (ä,ö,ü)

2004-11-30 Thread Duncan Murdoch
[I tried to send this message privately, but the return address bounced.] I think this has been fixed in R-patched, but I doubt if the fix has been tested in Win98. Could you please download a copy from and confirm that it has been fixed?

Re: [R] About ROC curves

2004-11-30 Thread Thomas Lumley
On Mon, 29 Nov 2004, Xin Qi wrote: Hi, Dear all R users: Does someone know whether R can calculate the Receiver Operating Characteristic (ROC) Curves? I didn't find it from the packages. Issue 1 of the R Newsletter this year used ROC curve calculation as an example of simple S3 and S4 classes.

[R] glmmPQL

2004-11-30 Thread Alex
Dear listmembers, I've adjusted a mixed model with glmmPQL: nulo<-glmmPQL(POS~1,random=~1|GRUPO, family=binomial, data=new) and I've reched the following results: Linear mixed-effects model fit by maximum likelihood Data: new AIC BIClogLik 53238.5 53260.74 -26616.25 Random e

Re: [R] Testing for S4 objects

2004-11-30 Thread John Chambers
Let me suggest a different test, because slotNames was written to work differently when given a string or a class definition. With your definition, R> x <- "classRepresentation" R> isS4object(x) [1] TRUE which I assume is not what you wanted. (Given a single string, slotNames() tries to look u

Re: [R] A basic question

2004-11-30 Thread Jari Oksanen
On Tue, 2004-11-30 at 13:58, Kenneth wrote: > Hi R users: > > I want to know any experience compiling R in other LINUX distributions > besides FEDORA (Red Hat) or Mandrake, for example in BSD, Debian, > Gentoo, Slackware, vector LINUX, Knoppix, Yopper or CERN linux? > > Hope this is not a "basic

RE: [R] Info

2004-11-30 Thread Andy Bunn
It looks like factanal is unable to optimize from these starting values (kinda like the error message says). So, factanal.fit.mle isn't converging and you have problems with your analysis. Try putting control = list(trace = T) in your code to see what happenens. E.g., R > R > v1 <- c(1,1,1,1,

Re: [R] A basic question

2004-11-30 Thread Peter Dalgaard
Kenneth <[EMAIL PROTECTED]> writes: > Hi R users: > > I want to know any experience compiling R in other LINUX distributions > besides FEDORA (Red Hat) or Mandrake, for example in BSD, Debian, > Gentoo, Slackware, vector LINUX, Knoppix, Yopper or CERN linux? > > Hope this is not a "basic questio

Re: [R] A basic question

2004-11-30 Thread Douglas Bates
Kenneth wrote: I want to know any experience compiling R in other LINUX distributions besides FEDORA (Red Hat) or Mandrake, for example in BSD, Debian, Gentoo, Slackware, vector LINUX, Knoppix, Yopper or CERN linux? Hope this is not a "basic question" Thank you for your help. There is considerable

[R] Info

2004-11-30 Thread Catherine Dempsey
I am having difficulty obtaining the scores from my principal component analysis. I have used this method before and have had no problems. The data set that I am using this time is similar to what I have used in the past. What do I need to do to my dataset in order for me to obtain these scores?

Re: [R] lme in R-2.0.0: Problem with lmeControl

2004-11-30 Thread Douglas Bates
Pavel Khomski wrote: Hello! One note/question hier about specification of control-parameters in the lme(...,control=list(...)) function call: i tried to specify tne number of iteration needed via lme(,control=list(maxIter=..., niterEM=...,msVerbose=TRUE)) but every time i change the defualt

Re: [R] Opening connection to FTP site

2004-11-30 Thread Prof Brian Ripley
On Tue, 30 Nov 2004, michael watson (IAH-C) wrote: Hi Is it possible to open a connection to an FTP site such that I can read the directory listing? Eg: URL <- url("ftp://ftp.ensembl.org";, open="r") Error in url("ftp://ftp.ensembl.org";, open="r") : unable to open connection URL <- url("f

[R] Opening connection to FTP site

2004-11-30 Thread michael watson \(IAH-C\)
Hi Is it possible to open a connection to an FTP site such that I can read the directory listing? Eg: URL <- url("ftp://ftp.ensembl.org";, open="r") Error in url("ftp://ftp.ensembl.org";, open="r") : unable to open connection URL <- url("ftp://ftp.ensembl.org";) open(URL) Error in open.

RE: [R] How to know if a bug was recognised

2004-11-30 Thread Henrik Bengtsson
> -Original Message- > From: [EMAIL PROTECTED] > [mailto:[EMAIL PROTECTED] On Behalf Of Heinz Tuechler > Sent: Tuesday, November 30, 2004 1:25 PM > To: [EMAIL PROTECTED] > Subject: [R] How to know if a bug was recognised > > > Hello! > > A problem with special characters seemed to me to

Re: [R] How to know if a bug was recognised

2004-11-30 Thread Prof Brian Ripley
You have already been sent an answer. It is a bug in your copy of Windows, and it will be worked around in the current version of R-patched. From the CHANGES file: R 2.0.1 patched === We work around reported bugs in Windows XP as to which characters are printable by attempti

Re: [R] 2k-factorial design with 10 parameters

2004-11-30 Thread Petr Pikal
On 30 Nov 2004 at 12:59, Sven wrote: > Hi, > > I'd like to apply a 2^k factorial design with k=10 parameters. > Obviously this results in a quite long term for the model equation due > to the high number of combinations of parameters. > > How can I specify the equation for the linear model (lm

[R] How to know if a bug was recognised

2004-11-30 Thread Heinz Tuechler
Hello! A problem with special characters seemed to me to be a bug. I sent a mail to [EMAIL PROTECTED] concerning the problem (see below). How can I find out, if this is considered as a bug or an error of myself? Which part of FAQs or documentation did I miss to find the answer? thanks in advance

[R] Newsletter - Live-Jackpot und taegliches Glueckshoroskop auf www.casinos.ch

2004-11-30 Thread www.casinos.ch
Live-Jackpot und tägliches Glückshoroskop auf www.casinos.ch Seit dieser Woche stellt die Nr. 1 Casino-Plattform der Schweiz seinen Besuchern zwei weitere, exklusive Dienstleistungen zur Verfügung. DIE LIVE-JACKPOT-ANBINDUNG Nichts ist ärgerlicher für Kunden, als veraltete oder verspätete Info

Re: [R] 2k-factorial design with 10 parameters

2004-11-30 Thread Prof Brian Ripley
On Tue, 30 Nov 2004, Sven wrote: I'd like to apply a 2^k factorial design with k=10 parameters. Obviously this results in a quite long term for the model equation due to the high number of combinations of parameters. How can I specify the equation for the linear model (lm) without writing all c

Re: [R] 2k-factorial design with 10 parameters

2004-11-30 Thread Dimitris Rizopoulos
Hi Sven, just use: lm(y~(x1+x2+x3+...+x10)^10) e.g., y <- rnorm(5000) x1 <- factor(sample(0:1, 5000, TRUE)) x2 <- factor(sample(0:1, 5000, TRUE)) x3 <- factor(sample(0:1, 5000, TRUE)) x4 <- factor(sample(0:1, 5000, TRUE)) lm1 <- lm(y~(x1+x2+x3+x4)^4) summary(lm1) I hope it helps. Best, Dimitris ---

Re: [R] plotting data in non-orthogonal coords.

2004-11-30 Thread Petr Pikal
On 30 Nov 2004 at 11:06, Andreas Franke wrote: > > Hi ! > Thanks for replying to my post and sorry for not being specific > enough. Maybe I am the one who didnt get the point , but as far as I > understand R plots filled.contour(x,y,z,...) in the following way: > > x,y define a grid in cartes

Re: [R] Building latest version of package

2004-11-30 Thread Prof Brian Ripley
On Mon, 29 Nov 2004, michael watson (IAH-C) wrote: Hi I have a package which was built using R 1.9.1 and everything worked fine. I recently upgraded to R 2.0.1 and tried to re-install my package - and I got: Error in library(mypackage) : 'mypackage' is not a valid package -- installed < 2.0.0? I d

  1   2   >