Re: [R] Re: Reasons not to answer very basic questions...

2004-11-30 Thread Uwe Ligges
Jim Lemon wrote: A.J. Rossini wrote: and perhaps the most important reason for the particular socratic form of teaching on this list... Golly, anyone who read Plato's Dialogues would realize that the Socratic method involves patiently leading the questioner stepwise through the solution, not

[R] about cancor.R

2004-11-30 Thread yan zhao
Hello, I'm a beginning user of R, now I have a question about canonical correlation analysis (cca). In R,there is a function cancor.R used for cca; For example X(n*p1) and Y(n*p2)are the two matrix to be analyzed. In the example given by R, when n max(p1,n2), cancor(X,Y) works; but when np1 or

Re: [R] about cancor.R

2004-11-30 Thread Prof Brian Ripley
On Tue, 30 Nov 2004, yan zhao wrote: Hello, I'm a beginning user of R, now I have a question about canonical correlation analysis (cca). In R,there is a function cancor.R used for cca; For example There is no such function in R, and I can't find it in any add-on package. Do you mean the function

Re: Re: [R] systemfit - SUR

2004-11-30 Thread contact
Arne Henningsen [EMAIL PROTECTED] schrieb am 29.11.2004, 17:02:12: On Monday 29 November 2004 16:42, [EMAIL PROTECTED] wrote: Hello to everyone, I have 2 problems and would be very pleased if anyone can help me: 1) When I use the package systemfit for SUR regressions, I get two

[R] Problem with print() and backslashes.

2004-11-30 Thread Kevin Brinkmann
Dear R List I have a small problem concerning the output of print(). My version: R.version _ platform i386-portbld-freebsd5.2 arch i386 os freebsd5.2 system i386, freebsd5.2 status major1 minor9.0 year 2004 month04 day 12 language R Consider this: I

Re: [R] Problem with print() and backslashes.

2004-11-30 Thread Prof Brian Ripley
On Tue, 30 Nov 2004, Kevin Brinkmann wrote: Dear R List I have a small problem concerning the output of print(). My version: R.version _ platform i386-portbld-freebsd5.2 arch i386 os freebsd5.2 system i386, freebsd5.2 status major1 minor9.0 year 2004 month04 day

Re: [R] Problem with print() and backslashes.

2004-11-30 Thread Peter Dalgaard
Kevin Brinkmann [EMAIL PROTECTED] writes: Consider this: I want to print a backslash with an exclamation mark. Here is the output. print( \! ) [1] ! Now I try it differently... print( \\! ) [1] \\! The output contains two backslashes. Why? (Didn't we do that one only

RE: [R] Problem with print() and backslashes.

2004-11-30 Thread Henrik Bengtsson
-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Kevin Brinkmann Sent: Tuesday, November 30, 2004 9:43 AM To: [EMAIL PROTECTED] Subject: [R] Problem with print() and backslashes. Dear R List I have a small problem concerning the output of

[R] Creating a factor from a combination of vectors

2004-11-30 Thread Yves Brostaux
Dear list, Here's a little problem I already solved with my own coding style, but I feel there is a more efficient and cleaner way to write it, but had no success finding the clever solution. I want to produce a factor from a subset of the combination of two vectors. I have the vectors a et b

Re: [R] systemfit - SUR

2004-11-30 Thread Arne Henningsen
On Tuesday 30 November 2004 09:34, [EMAIL PROTECTED] wrote: Arne Henningsen [EMAIL PROTECTED] schrieb am 29.11.2004, 17:02:12: On Monday 29 November 2004 16:42, [EMAIL PROTECTED] wrote: Hello to everyone, I have 2 problems and would be very pleased if anyone can help me: 1)

Re: Re: [R] systemfit - SUR

2004-11-30 Thread contact
Arne Henningsen [EMAIL PROTECTED] schrieb am 30.11.2004, 10:20:29: On Tuesday 30 November 2004 09:34, [EMAIL PROTECTED] wrote: Arne Henningsen schrieb am 29.11.2004, 17:02:12: On Monday 29 November 2004 16:42, [EMAIL PROTECTED] wrote: Hello to everyone, I have 2 problems

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

2004-11-30 Thread Andreas Franke
Am Dienstag, 30. November 2004 10:26 schrieb Petr Pikal: Hallo Andreas There is probably no one who can give you some answer as you did provide almost no facts what you really did and what is wrong. You probably need to transform your coordinates to orthogonal and plot them as you wish. But

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

2004-11-30 Thread Petr Pikal
Hallo Andreas There is probably no one who can give you some answer as you did provide almost no facts what you really did and what is wrong. You probably need to transform your coordinates to orthogonal and plot them as you wish. But you have to do it yourself. I had some data in polar

Re: [R] systemfit - SUR

2004-11-30 Thread Arne Henningsen
On Tuesday 30 November 2004 10:50, [EMAIL PROTECTED] wrote: [ . . .] Please show how you obtained these results. I did not provide the steps, because the data is not public. A general remark: In most cases a questioner can use the example data provided with a package or construct simple

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

2004-11-30 Thread Pavel Khomski
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 values maxIter (e.g.

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

2004-11-30 Thread Eric Lecoutre
Hi Yves, Using your objects, here is a way: cascombo=do.call(paste,c(cas,sep=.)) factor(do.call(paste,c(df,sep=.)),levels=cascombo,labels=rownames(cas)) [1] LowNA NA NA NA NA NA NA NA NA NA NA Medium NA NA [16] NA NA NA NA NA NA NA NA NA High

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

2004-11-30 Thread Gabor Grothendieck
Yves Brostaux brostaux.y at fsagx.ac.be writes: : : Dear list, : : Here's a little problem I already solved with my own coding style, but I : feel there is a more efficient and cleaner way to write it, but had no : success finding the clever solution. : : I want to produce a factor from a

[R] A basic question

2004-11-30 Thread Kenneth
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 question Thank you for your help. Kenneth

[R] 2k-factorial design with 10 parameters

2004-11-30 Thread Sven
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) without writing all combinations explicitly down

[R] (sem assunto)

2004-11-30 Thread cynthia
- gora! [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

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

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

2004-11-30 Thread Petr Pikal
On 30 Nov 2004 at 11:06, Andreas Franke wrote: snip 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

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] 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

[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

[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

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)

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

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 be a bug.

[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

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 -

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

[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

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

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 question BSD is

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 -

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] 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 up

[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

[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 http://cran.r-project.org/bin/windows/base/rpatched.html and confirm that it has been fixed?

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

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:

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? Which

[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]]

[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] 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

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

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

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

[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, xk thx. and rgds. a/m

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] 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$data) :

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==1)) lmtest

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: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, not y.

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: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. That's

[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

[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/ 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

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 neither

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.sat)

[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:

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

[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$ivhx0] - 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 represented

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 to

[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*`, 1, value =

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,

[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: -

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$ivhx0] - 1*(uis$ivhx==3) Error: NAs are not allowed in subscripted assignments but this would happen in version 1.91 What would happen? In version

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

[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

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 the plot area

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] 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 will

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 categorical

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 -

[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

[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=distance19) augPred(fm1, length.out = 2, level = c(0,1))

[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 -

[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 = 6,

[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] 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] 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

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]]

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.

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 described in the

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

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

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 an

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

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 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 expression

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=distance19) augPred(fm1, length.out =

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: D1 -

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(text=8))

[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) { blank - a[i:N] # or whatever else you

[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) *

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

[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 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

  1   2   >