Re: [R] ANOVA Table for Full Linear Model?

2007-02-06 Thread talepanda
summary.lm tells you F, df, and p value. try: summary(regression9) On 2/6/07, Jason R. Finley [EMAIL PROTECTED] wrote: Hello, I have spent a good deal of time searching for an answer to this but have come up empty-handed; I apologize if I missed something that is common knowledge. I am

Re: [R] Rconsole - setting the size and location of Windows help files (Rgui)

2007-02-06 Thread john.gavin
Hi, # Dimensions(in characters) of the internal pager. pgrows = 48 pgcolumns = 128 Thanks, that worked. Regards, John. -Original Message- From: talepanda [mailto:[EMAIL PROTECTED] Sent: 06 February 2007 06:37 To: Gavin, John Cc: R-help@stat.math.ethz.ch Subject: Re: [R]

Re: [R] ANOVA Table for Full Linear Model?

2007-02-06 Thread Peter Dalgaard
Jason R. Finley wrote: Hello, I have spent a good deal of time searching for an answer to this but have come up empty-handed; I apologize if I missed something that is common knowledge. I am trying to figure out how to get an ANOVA table that shows the sum of squares. degrees of

Re: [R] Metapost device driver

2007-02-06 Thread Prof Brian Ripley
Please discuss this on R-devel: it is way off-topic for R-help. '11' is the number of arguments, it is not 'magic'. On Tue, 6 Feb 2007, Geoff Russell wrote: Hi All, I've started work on a MetaPost device driver (please don't hold your breath). I've copied the XFig driver and renamed

Re: [R] Rconsole - setting the size and location of Windows help files (Rgui)

2007-02-06 Thread Prof Brian Ripley
On Tue, 6 Feb 2007, talepanda wrote: For size, maybe: # Dimensions(in characters) of the internal pager. #pgrows = 25 #pgcolumns = 80 pgrows = 48 pgcolumns = 128 in Rconsole, but location cannot be handled. For the very good reason that you can have multiple pagers and I at least do

Re: [R] ar function in stats

2007-02-06 Thread Prof Brian Ripley
On Mon, 5 Feb 2007, Leeds, Mark (IED) wrote: I had a couple of questions about the ar function that i was hoping someone could answer. [...] I call the ar function with aic=TRUE as below so that it picks the order of the ar model based on BIC. It actually does as it says it does (BIC !=

Re: [R] Rconsole - setting the size and location of Windows help files (Rgui)

2007-02-06 Thread talepanda
For the very good reason that you can have multiple pagers and I at least do not want them exactly on top of one another. Yes. and I knew for the first time now that location of new pager is randomly determined. I think it is cool way. On 2/6/07, Prof Brian Ripley [EMAIL PROTECTED] wrote: On

[R] How-To construct a cov list to use a covariance matrix in factanal?

2007-02-06 Thread Alistair Campbell
Hi, I have a set of covariance matrices but not the original data. I want to carry out some exploratory factor analysis. So, I am trying to construct a covariance matrix list as the input for factanal. I can construct a list which includes the cov, the centers, and the n.obs. But it doesn't

[R] How to do moran's I test?

2007-02-06 Thread Jian Zhang
I want to do moran's I test in R language. I try to use gearymoran in Package ade4,moran in Package spdep, and Moran.I in Package ape. But I do not know how to do it because data format is different. My data: xy dbh 111.0310.7 7 118.110.28 1.2

Re: [R] any implementations for adaptive modeling of time series?

2007-02-06 Thread Peter Nimda
Hi Ansel, thank you for the response generally speaking, wavelets are known to be good at extracting signal from noisy data and are adaptive but I am not familiar with any R implementation of wavelets. I used wavelets before. 1. wavelets is just a particular case of orthogonal function

[R] R equivalent to Matlab spline(x,y,xx)

2007-02-06 Thread Maria Vatapitakapha
Hello I am currently translating some Matlab code into R and would like to know what people would recommend as the equivalent in R for spline(x,y,xx) in matlab. The latter performs a cubic spline interpolation where x,y are the given data vectors and xx is the new abscissa vector. It allows for

Re: [R] LyX 1.4.3-5 problem generating PDF documents

2007-02-06 Thread John Kane
--- Stefan Grosse [EMAIL PROTECTED] wrote: You are posting in the wrong mailing list, mail your question to the lyx users list: http://www.lyx.org/internet/mailing.php OOPs I wondered why I wasn't getting a reply. I use the same account for both lists. My appologies to the list John

Re: [R] ANOVA Table for Full Linear Model?

2007-02-06 Thread Frank E Harrell Jr
Peter Dalgaard wrote: Jason R. Finley wrote: Hello, I have spent a good deal of time searching for an answer to this but have come up empty-handed; I apologize if I missed something that is common knowledge. I am trying to figure out how to get an ANOVA table that shows the sum of

Re: [R] How to do moran's I test?

2007-02-06 Thread Lyndsey Stanfill
First, I use lattice data rather than point pattern data, so this may not be applicable. The following is an example of the Moran's I test using the south dataset, analyzing the variable HR60 for dependence. I use GeoDa to create my spatial weight (https://www.geoda.uiuc.edu/), and I call the

[R] heatmap from xyz data

2007-02-06 Thread Trevor Graham
Hi, I've got some data in a data frame arranged like this: x y z othervariables 0.1 0.2 1.7 0.01 0.2 0.2 1.3 0.23 0.2 0.3 1.1 0.43 etc I'd like to plot a heatmap of this data, with x and y as the (x,y) co- ords

[R] Questions on counts by case

2007-02-06 Thread Serguei Kaniovski
Hi all, for the data below I would like to 1. generate a dummy variable for each group gr of the same composition by people, then save each portion in a separate file, 2. compute the frequency of 1's in x for each person by group gr. So, mike will have freq=2/3, as he has two 1 and one 0 in

Re: [R] lme in R and Splus-7

2007-02-06 Thread Dieter Menne
yyan liu zhliur at yahoo.com writes: I used the function lme in R and Splus-7. With the same dataset and same argument for the function, I got quite different estimation results from these two software. Anyone has this experience before? From the FAQ: R by default uses treatment contrasts

Re: [R] heatmap from xyz data

2007-02-06 Thread Vladimir Eremeev
You can use plot(y~x,col=color.index.in.palette.defined.from(z),pch=20,type=p) where color.index.in.palette.defined.from(z) is a function or an expression, returning either a color index in a predefined palette or any other color representation, suitable for R. This is described in ?par. I

Re: [R] How-To construct a cov list to use a covariance matrix in factanal?

2007-02-06 Thread Prof Brian Ripley
The help page says covmat: A covariance matrix, or a covariance list as returned by 'cov.wt'. Of course, correlation matrices are covariance matrices. and there is an example of a covariance list (ability.cov). factanal(factors = 2, covmat = ability.cov)

[R] Query about merging two tables

2007-02-06 Thread lalitha viswanath
Hi I have table1 which has the foll. columns id age rate and table2 which has the foll. columns id count I wish to get data from table1 for all the ids which are persent in table2 and where the rate is not equal to 999. The ids in table2 are a subset of those in table1 and every id in table2 has

[R] installing packages and windows vista

2007-02-06 Thread Daniel O'Shea
I installed R (R-2.4.1-win32.exe) on a new computer with Windows Vista and a 64 bit operating system (hp dv9000 with intel core t7200). The base R runs fine, but I can not get any of the packages to load. From within R I choose install packages choose a site then a package. I tried installing

[R] Query about merging two tables

2007-02-06 Thread lalitha viswanath
Hi I have table1 which has the foll. columns id age rate and table2 which has the foll. columns id count I wish to get data from table1 for all the ids which are persent in table2 and where the rate is not equal to 999. The ids in table2 are a subset of those in table1 and every id in table2 has

Re: [R] Query about merging two tables

2007-02-06 Thread Wensui Liu
subset(table1, rate != 999id == table2$id) On 2/6/07, lalitha viswanath [EMAIL PROTECTED] wrote: Hi I have table1 which has the foll. columns id age rate and table2 which has the foll. columns id count I wish to get data from table1 for all the ids which are persent in table2 and where

Re: [R] Query about merging two tables

2007-02-06 Thread Vladimir Eremeev
Wensui Liu wrote: subset(table1, rate != 999id == table2$id) On 2/6/07, lalitha viswanath [EMAIL PROTECTED] wrote: Hi I have table1 which has the foll. columns id age rate and table2 which has the foll. columns id count I wish to get data from table1 for all the ids which are

[R] Conditional Sum

2007-02-06 Thread Bert Jacobs
Hi, I would like to make a conditional sum for certain columns in a dataframe. This is how my dataframe looks like: Clinic Rep ColM1 ColM2 ColM3 ... ColM40 A 1 10 01 B 1 0-1 00 C 1 0-1 1-1 A 2 11 -1

Re: [R] installing packages and windows vista

2007-02-06 Thread Duncan Murdoch
On 2/6/2007 10:33 AM, Daniel O'Shea wrote: I installed R (R-2.4.1-win32.exe) on a new computer with Windows Vista and a 64 bit operating system (hp dv9000 with intel core t7200). The base R runs fine, but I can not get any of the packages to load. From within R I choose install packages

Re: [R] installing packages and windows vista

2007-02-06 Thread Peter Dalgaard
Daniel O'Shea wrote: I installed R (R-2.4.1-win32.exe) on a new computer with Windows Vista and a 64 bit operating system (hp dv9000 with intel core t7200). The base R runs fine, but I can not get any of the packages to load. From within R I choose install packages choose a site then a

Re: [R] Conditional Sum

2007-02-06 Thread Kaskelma, Heikki
Yes, consider: df=data.frame(Clinic=rep(c(A, B, C), 2), Rep=c(1, 1, 1, 2, 2, 2), colM1=c(1, 0, 0, 1, -1, 1), ColM2=c(0, -1, -1, 1, 0, 0), ColM3=c(0, 0, 1, -1, 0, 1), ColM40=c(1, 0, -1, 0, 1, -1) ) Clinic=1; Rep=2

[R] RES: RdbiPgSQL in R 2.4.1

2007-02-06 Thread Eduardo Dutra de Armas
Hi Seth I'm running R on WinXP. On the first time I've installed Rdbi and RdbiPgSQL from zip files not downloaded from Bioconductor. No problems during installation. Now I tried to install from biocLite and a message was showed for RdbiPgSQL. biocLite(c(Rdbi, RdbiPgSQL)) Running getBioC version

[R] R in Industry

2007-02-06 Thread Doran, Harold
The other day, CNN had a story on working at Google. Out of curiosity, I went to the Google employment web site (I'm not looking, but just curious). In perusing their job posts for statisticians, preference is given to those who use R and python. Other languages, S-Plus and something called SAS

Re: [R] How to do moran's I test?

2007-02-06 Thread Roger Bivand
On Tue, 6 Feb 2007, Jian Zhang wrote: I want to do moran's I test in R language. I try to use gearymoran in Package ade4,moran in Package spdep, and Moran.I in Package ape. But I do not know how to do it because data format is different. Please consider posting on the R-sig-geo mailing list.

Re: [R] RES: RdbiPgSQL in R 2.4.1

2007-02-06 Thread Seth Falcon
Eduardo Dutra de Armas [EMAIL PROTECTED] writes: Hi Seth I'm running R on WinXP. On the first time I've installed Rdbi and RdbiPgSQL from zip files not downloaded from Bioconductor. No problems during installation. Now I tried to install from biocLite and a message was showed for RdbiPgSQL.

[R] Naming bring to top graphical Device

2007-02-06 Thread Sebastien Durand
Dear all, Here is my questions: 1- Under a WINDOWS installation of R-2.4.1, can we change the naming of a new ploting device open by the command windows()?. Instead of the default name e.g.: Device 2 I would like to use something like Density plot or whatever! 2- Under a MAC OS X

[R] glm gamma scale parameter

2007-02-06 Thread WILLIE, JILL
I would like the option to specify alternative scale parameters when using the gamma family, log link glm. In particular I would like the option to specify any of the following: 1. maximum likelihood estimate 2. moment estimator/Pearson's 3. total deviance estimator Is this easy? Possible?

[R] Confirmation Of Results For Electronic Email Winners

2007-02-06 Thread Angela Frenz
Confirmation Of Results For Electronic Email Winners MICROWORLD INTERNATIONAL PROMOTIONS EUROPEAN CENTRAL OFFICE 162 ZUIDOOST 1109 AMSTERDARM THE NETHERLANDS We happily announce to you the draw of the Microworld International promotions organized for email users all over the World as part of

Re: [R] glm gamma scale parameter

2007-02-06 Thread Prof Brian Ripley
I think you mean 'shape parameter'. If so, see the MASS package and ?gamma.shape. glm() _is_ providing you with the MLE of the scale parameter, but really no estimate of the shape (although summary.glm makes use of one). On Tue, 6 Feb 2007, WILLIE, JILL wrote: I would like the option to

[R] A question regarding cex and pch=. in lattice

2007-02-06 Thread Saptarshi Guha
Hello, I'm using lattice and opened an X11 device with the following call X11(width=5,height=5,pointsize=1) I then ran the following code library(lattice) x-rnorm(30,sd=2) y-runif(30) xyplot(y~x,pch=.,col=black,cex=1) If i

Re: [R] glm gamma scale parameter

2007-02-06 Thread Prof Brian Ripley
On Tue, 6 Feb 2007, Prof Brian Ripley wrote: I think you mean 'shape parameter'. If so, see the MASS package and ?gamma.shape. Also http://www.stats.ox.ac.uk/pub/MASS4/#Complements leads to several pages of discussion. glm() _is_ providing you with the MLE of the scale parameter, but

Re: [R] A question regarding cex and pch=. in lattice

2007-02-06 Thread Prof Brian Ripley
How do you know the points are not being plotted? All you can tell is that your graphics card / monitor / X11 driver is not displaying them, not at all the same thing. My guess is that the rectangle is being displayed only if the differences in its discretized x and y coordinates are both

Re: [R] glm gamma scale parameter

2007-02-06 Thread WILLIE, JILL
Thank you. You are correct, the shape parameter is what I need to change I think I see how to use the MASS package to do it...or if not, at least I have enough now to figure it out. A question to reconcile terminology which will speed me up, if you have time to help me a bit more: phi = 'scale

Re: [R] A question regarding cex and pch=. in lattice

2007-02-06 Thread Saptarshi Guha
On Feb 6, 2007, at 4:05 PM, Prof Brian Ripley wrote: How do you know the points are not being plotted? All you can tell is that your graphics card / monitor / X11 driver is not displaying them, not at all the same thing. Yes, that is what I meant but I suppose I conveyed it wrongly.

Re: [R] glm gamma scale parameter

2007-02-06 Thread Prof Brian Ripley
On Tue, 6 Feb 2007, WILLIE, JILL wrote: Thank you. You are correct, the shape parameter is what I need to change I think I see how to use the MASS package to do it...or if not, at least I have enough now to figure it out. A question to reconcile terminology which will speed me up, if you

Re: [R] R in Industry

2007-02-06 Thread Kuhn, Max
As someone who has (reluctantly) sent job postings to R Help, I think that a SIG would be a good idea. Max -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Doran, Harold Sent: Tuesday, February 06, 2007 2:08 PM To: R-help@stat.math.ethz.ch Subject: [R] R

Re: [R] R in Industry

2007-02-06 Thread ¨Tariq Khan
Our (financial) firm has struggled immensely with finding appropriate candidates and I believe this would be *tremendously* useful. I really hope that your suggestion is implemented. Importantly, I believe it would help create a commercial presence for R, and secondly reduce use of products like

Re: [R] R in Industry

2007-02-06 Thread Bert Gunter
... two main drawbacks of R at our firm (as viewed by our IT dept) are lack of guaranteed support as well as the difficulty in finding candidates. -- Just an aside: lack of guaranteed support -- absolutely true in theory, absolutely false in practice. I doubt that the voluntary support found on

[R] formula for test statistics of arima coefficients

2007-02-06 Thread Leeds, Mark \(IED\)
I have an object, ARestopt, that is the result of a call to the arima function. I want to calculate test statistics ( null is coeffs are zero ) for each of the estimated coefficients. My coefficents are in a vector called ARestopt$coef and my covariance matrix is ARestopt$var.coef. I thought

Re: [R] Questions on counts by case

2007-02-06 Thread jim holtman
On 2/6/07, Serguei Kaniovski [EMAIL PROTECTED] wrote: Hi all, for the data below I would like to 1. generate a dummy variable for each group gr of the same composition by people, then save each portion in a separate file, This gives back a list that you can go through and write out each

[R] abbreviate dataframe for Sweave output

2007-02-06 Thread stubben
I wanted to print the first and last rows of some dataframes in Sweave using dots in columns to separate the two parts. Head and tail almost work, but I have problems with factors and row names. z-data.frame(id=letters[1:26], x=sample(1:26,26)) rbind(head(z,3), ., tail(z,1)) id x 1

Re: [R] R in Industry

2007-02-06 Thread Muenchen, Robert A (Bob)
That sounds like a good idea. The name R makes it especially hard to find job postings, resumes or do any other type of search. Googling resume+sas or job opening+sas is quick and fairly effective (less a few airline jobs). Doing that with R is of course futile. At the risk of getting flamed, it's

[R] 2 Courses: (1) R/Splus Advanced Programming - San Francisco - (2) Data Mining: Practical Tools and Techniques in R/Splus - Salt Lake City

2007-02-06 Thread Sue Turner
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce our*** R/S Advanced Programming ***course in San Francisco on March 15-16 and *** Data Mining: Practical Tools and Techniques in R/Splus *** course in Salt Lake City on March 26-27 Ask for group discount and reserve your

Re: [R] How-To construct a cov list to use a covariance matrix in factanal?

2007-02-06 Thread Alistair Campbell
Thanks for that Brian, I have worked through the examples. They work because the covmat were produced by the cov.wt which provides output as a list object. I am trying to construct my own list object to use as the covmat. There are no obvious instructions on how to do this. So, here is what I

Re: [R] abbreviate dataframe for Sweave output

2007-02-06 Thread Charles C. Berry
Chris, Not precisely what you want (does both rows and columns), but can be specialized to your case easily. Use with xtable in a chunk with 'results=tex' as in xtable( dot.matrix( diag(10) ) ) Chuck dot.matrix - function(x , rws=NULL, cls=NULL, in.math=FALSE ) { ## Purpose:

[R] How to reorder rows in dataframe by text flag

2007-02-06 Thread Dale Steele
Given two columns of type character in a dataframe of the form: col1col2 31* 66 0 0* 102*66 71* 80 31 2* 66 31* 47 38* How do I generate the following dataframe? Ie. col1 contains row item with * and col2 contains row member without * col1col2

Re: [R] R in Industry

2007-02-06 Thread Wensui Liu
I've been looking for job that allows me to use R/S+ since I got out of graduate school 2 years ago but with no success. I am wondering if there is something that can be done to promote the use of R in industry. It's been very frustrating to see people doing statistics using excel/spss and even

Re: [R] How-To construct a cov list to use a covariance matrix in factanal?

2007-02-06 Thread William Revelle
Alistair, I have worked through the examples. They work because the covmat were produced by the cov.wt which provides output as a list object. I am trying to construct my own list object to use as the covmat. There are no obvious instructions on how to do this. So, here is what I have done

Re: [R] How to reorder rows in dataframe by text flag

2007-02-06 Thread talepanda
I think several ways can do that. my code is: data.frame(t(apply(dat,1,function(x)as.numeric(sub(\\*,,x[sort(1:2,grep(\\*,x)==2)]) HTH On 2/7/07, Dale Steele [EMAIL PROTECTED] wrote: Given two columns of type character in a dataframe of the form: col1col2 31* 66 0 0*

Re: [R] How to reorder rows in dataframe by text flag

2007-02-06 Thread talepanda
I'm sorry, I forgot my definition of dat. dat in code of previous post is: dat-data.frame(col1=c(31*,0,102*,71*,31,66,47),col2=c(66,0*,66,80,2*,31*,38*)) dat col1 col2 1 31* 66 20 0* 3 102* 66 4 71* 80 5 31 2* 6 66 31* 7 47 38* On 2/7/07, talepanda [EMAIL PROTECTED]

Re: [R] multinomial logistic regression with equality constraints?

2007-02-06 Thread Roger Levy
Hi again Jasjeet, Walter, I have a further question about an error message I get when running multinomRob. I am simulating a dataset where I look at the effect of making a previous categorical choice on the probability of making the same choice later on. Given the following code: n - 20

Re: [R] How to reorder rows in dataframe by text flag

2007-02-06 Thread Gabor Grothendieck
Try this: DF - data.frame( col1 = factor(c(31*, 0, 102*, 71*, 31, 66, 47)), col2 = factor(c(66, 0*, 66, 80, 2*, 31*, 38)) ) replace(DF, TRUE, as.numeric(sub(*, , as.matrix(DF), fixed = TRUE))) On 2/6/07, Dale Steele [EMAIL PROTECTED] wrote: Given two columns of type character

Re: [R] R in Industry

2007-02-06 Thread Matthew Keller
Bob, Far from flaming you, I think you made a good point - one that I imagine most people who use R have come across. The name R is a big impediment to effective online searches. As a check, I entered R software, SAS software, SPSS software, and S+ software into google. The R 'hit rate' was only

Re: [R] multinomial logistic regression with equality constraints?

2007-02-06 Thread Walter Mebane
Roger, Error in if (logliklambda loglik) bvec - blambda : missing value where TRUE/FALSE needed In addition: Warning message: NaNs produced in: sqrt(sigma2GN) That message comes from the Newton algorithm (defined in source file multinomMLE.R). It would be better if we

[R] hi

2007-02-06 Thread Havel Masterson
Hi, Viragra 1,80 Ciralis 3,00 Lervitra 3,35 http://officecop.printeryml-com Replace - with . in the above link -- would bring him here, beyond the reach of Dumbledores help and protection, and into my waiting arms. And here he is ...the boy you all believed had been my downfall...

[R] Histogram of ones.

2007-02-06 Thread Dong H. Oh
Dear expeRts, I'd like to picture histogram of ones. For example, hist(rep(1, 100), col=lightblue, border=black) A bin is pictured ranging on [0,1] and range of x axis is also [0,1]. How can I adjust the bin range from [0,1] to [0.9, 1] ? I am looking forward to hearing from expeRt.

Re: [R] Histogram of ones.

2007-02-06 Thread talepanda
If I correctly understand, hist(rep(1, 100), col=lightblue, border=black,breaks=0:10*0.1) see: breaks arg in ?hist On 2/7/07, Dong H. Oh [EMAIL PROTECTED] wrote: Dear expeRts, I'd like to picture histogram of ones. For example, hist(rep(1, 100), col=lightblue, border=black) A bin is

[R] boxplot statistics in ggplot

2007-02-06 Thread Vikas Rawal
I need to make weighted boxplots. I found that ggplot makes them. I would however like to label them with the boxplot statistics (the median, q1 and q3). In the boxplot function in r-base, I could output boxplot statistics and then write a text on the plot to place the labels. How would one do it

[R] R Search

2007-02-06 Thread İbrahim Mutlay
I know that two of the search engine for R is available: http://www.rseek.org/ http://www.dangoldstein.com/search_r.html -- Ýbrahim Mutlay [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list