Re: [R] Using dll with Visual Studio Compiler

2007-06-13 Thread Prof Brian Ripley
On Tue, 12 Jun 2007, Ian McCarthy wrote: Hi, I have created a dll in Fortran and used the Visual Studio 2005 Compiler. I've read that certain problems can arise based on the Since Visual Studio does not create Fortran, that's a bit vague. And let us assume you are on Win32, but you failed

[R] passing (or obtaining) index or element name of list to FUN in lapply()

2007-06-13 Thread Stephen Tucker
Hello everyone, I wonder if there is a way to pass the index or name of a list to a user-specified function in lapply(). For instance, my desired effect is something like the output of L - list(jack=4098,sape=4139) lapply(seq(along=L),function(i,x) if(i==1) jack else sape,x=L) [[1]] [1] jack

Re: [R] PCA for Binary data

2007-06-13 Thread Prof Brian Ripley
On Tue, 12 Jun 2007, Spencer Graves wrote: The problem with applying prcomp to binary data is that it's not clear what problem you are solving. The standard principal components and factor analysis models assume that the observations are linear combinations of unobserved common

Re: [R] Viewing a data object

2007-06-13 Thread Stephen Tucker
Hi Horace, I have also thought that it may be useful but I don't know of any Object Explorer available for R. However, (you may alread know this but) (1) you can view your list of objects in R with objects(), (2) view objects in a spreadsheet-like table (if they are matrices or data frames)

[R] if statement

2007-06-13 Thread Jiong Zhang, PhD
Hi all, I have a rather naive question. I have the height of 100 individuals in a table and I want to assign the tallest 30% as Case=1 and the bottom 30% as Case=0. How do I do that? thanks. jiong The email message (and any attachments) is for the sole use of the intended recipient(s) and

[R] Panel data

2007-06-13 Thread Megh Dal
Dear all R users, I have a small doubt about panel data analysis. My basic understanding on Panel data is a type of data that is collected over time and subjects. Vector Autoregressive Model (VAR) model used on this type of data. Therefore can I say that, one of statistical tools used for

[R] Trouble making JRI.jar with Ubuntu and Java6

2007-06-13 Thread taivo
Hi, Forum newb here, looking for some help. Have been trying to install an R-Java interface to make R calls from Java. JRI's configure script runs fine, but when it comes to make, I get the error: too few arguments to function 'R_ParseVector' Java runs fine. R runs fine. But I can't get this

Re: [R] PCA for Binary data

2007-06-13 Thread ssls sddd
Dear Prof Brian Ripley, Would you also recommend some packages for non-binary data to do variable and feature selection? Thanks a lot! Alex On 6/12/07, Prof Brian Ripley [EMAIL PROTECTED] wrote: On Tue, 12 Jun 2007, Spencer Graves wrote: The problem with applying prcomp to binary

Re: [R] if statement

2007-06-13 Thread Simon Blomberg
My solutions are usually too baroque, but does this do what you want? x - rnorm(100) quants - quantile(x, c(.3, .7)) Case - rep(2, length(x)) # 2 lies in the middle of the distribution Case[x = quants[1]] - 0 Case[x = quants[2]] - 1 Case Cheers, Simon. On Mon, 2007-06-11 at 15:14 -0700, Jiong

[R] Confusion with sapply

2007-06-13 Thread Patnaik, Tirthankar
Hi, I have some confusion in applying a function over a column. Here's my function. I just need to shift non-March month-ends to March month-ends. Initially I tried seq.dates, but one cannot give a negative increment (decrement) here.

Re: [R] if statement

2007-06-13 Thread Ken Knoblauch
Simon Blomberg s.blomberg1 at uq.edu.au writes: My solutions are usually too baroque, but does this do what you want? x - rnorm(100) quants - quantile(x, c(.3, .7)) Case - rep(2, length(x)) # 2 lies in the middle of the distribution Case[x = quants[1]] - 0 Case[x = quants[2]] - 1 Case

Re: [R] passing (or obtaining) index or element name of list to FUN in lapply()

2007-06-13 Thread Prof Brian Ripley
On Tue, 12 Jun 2007, Stephen Tucker wrote: Hello everyone, I wonder if there is a way to pass the index or name of a list to a user-specified function in lapply(). For instance, my desired effect is something like the output of L - list(jack=4098,sape=4139)

[R] extractor rows from a matrix

2007-06-13 Thread billycorg
hi! i have a little problem: my data's matrix has 1093 rows and 3 columns. i'd like to extract each rows.. something like this: ht= my matrix Dt=(???)=a vector with t=1,2...1093 what can i do? thank you! Vincenzo -- View this message in context:

[R] export to a dat file that SAS can read

2007-06-13 Thread Rina Miehs
Hello i have a data frame in R that some SAS users need to use in their programs, and they want it in a dat file, is that possible? and which functions to use for that? my data frame is like this: out13[1:100,] faridniveau1 niveau3 p1

Re: [R] export to a dat file that SAS can read

2007-06-13 Thread Uwe Ligges
Rina Miehs wrote: Hello i have a data frame in R that some SAS users need to use in their programs, and they want it in a dat file, is that possible? What is a dat file? and which functions to use for that? I *guess* write.table() will do the trick, given dat is what I guess it is...

Re: [R] export to a dat file that SAS can read

2007-06-13 Thread Christophe Pallier
One solution: Use 'write.csv' (or write.csv2) in R to output a csv file. Then import it in SAS (by default, it imports into the work library; then you need to copy into another libray to have the SAS file). Christophe On 6/13/07, Rina Miehs [EMAIL PROTECTED] wrote: Hello i have a data

Re: [R] extractor rows from a matrix

2007-06-13 Thread Uwe Ligges
billycorg wrote: hi! i have a little problem: my data's matrix has 1093 rows and 3 columns. i'd like to extract each rows.. something like this: ht= my matrix Dt=(???)=a vector with t=1,2...1093 A vector with t=1,2, ..., 1093 and 3 columns is the matrix itself, isn't it? If you want to

Re: [R] Design library installation problem

2007-06-13 Thread Uwe Ligges
Ian Watson wrote: Dear Listers I have tried to install Frank Harrell's two libaries: Hmisc and Design. I found that Hmisc was listed in the list of packages from the Install Packages command on the Packages menu, but Design was not. I installed Hmisc from this list, and when I

Re: [R] Cause of error message in cov function?

2007-06-13 Thread Uwe Ligges
Matthew Keller wrote: Hi all, I have written a script in R that simulates genetically informative data - it is posted on my website and available to the public. This is my first time to write a script for use by others and am learning that it isn't as easy as it seems. To the issue. My

Re: [R] Read Windows-like .INI files into R data structure?

2007-06-13 Thread Vladimir Eremeev
One more question, inspired by this one, just to increase my R skill level. Earl F. Glynn wrote: I need to process some datasets where the configuration information was stored in .INI-like files, i.e., text files with sections like this: [Section1] var1=value1 var2=value2 [Section2]

Re: [R] Trouble making JRI.jar with Ubuntu and Java6

2007-06-13 Thread Prof Brian Ripley
You have a version mismatch. You haven't told us your version of JRI, but it is not recent enough for R 2.5.0. Also, you told us you are using java-6-sun-1.6.0.00, but the commands you show indicate otherwise. On Tue, 12 Jun 2007, taivo wrote: Hi, Forum newb here, looking for some help.

[R] Fitted value

2007-06-13 Thread livia
I would like to fit a Pareto Distribution and I am using the following codes. I thought the fitted (fit1) should be the fitted value for the data, is it correct? As the result of the fitted turns out to be a single value for all. fit=vglm(ycf1 ~ 1, pareto1(location=alpha), trace=TRUE, crit=c)

Re: [R] Confusion with sapply

2007-06-13 Thread Uwe Ligges
Patnaik, Tirthankar wrote: Hi, I have some confusion in applying a function over a column. Here's my function. I just need to shift non-March month-ends to March month-ends. Initially I tried seq.dates, but one cannot give a negative increment (decrement) here.

Re: [R] Awk and Vilno

2007-06-13 Thread Ted Harding
On 13-Jun-07 01:24:41, Robert Wilkins wrote: In clinical trial data preparation and many other data situations, the statistical programmer needs to merge and re-merge multiple input files countless times. A syntax for merging files that is clear and concise is very important for the

Re: [R] Linux equivalent to windows menus and script window

2007-06-13 Thread Stefan Grosse
JGR is pretty close to the windows R-GUI (with the menus) but even better since it is system independent and e.g. has syntax highlightning in the editor. You could have a look at the screenshots at http://rosuda.org/JGR/index.htm What you need is a Java 5 SDK installed at the system and

Re: [R] Confusion with sapply

2007-06-13 Thread hadley wickham
On 6/13/07, Patnaik, Tirthankar [EMAIL PROTECTED] wrote: Hi, I have some confusion in applying a function over a column. Here's my function. I just need to shift non-March month-ends to March month-ends. Initially I tried seq.dates, but one cannot give a negative increment (decrement) here.

[R] Subscription

2007-06-13 Thread Lucy Namu
I would like to subscribe and get free software for statistical analysis. Lucy ___ now. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] Read Windows-like .INI files into R data structure?

2007-06-13 Thread Christophe Pallier
var1=value1, A=value3 is almost pure R code. Is it possible to use this feature to solve the problem? Along the same lines: you may write a short script that converts the ini file into R code that can be sourced. From your example, you can generate the following R code: Section1 - list()

Re: [R] Read Windows-like .INI files into R data structure?

2007-06-13 Thread Vladimir Eremeev
Christophe Pallier wrote: var1=value1, A=value3 is almost pure R code. Is it possible to use this feature to solve the problem? Along the same lines: you may write a short script that converts the ini file into R code that can be sourced. From your example, you can generate the

[R] equivalent of windialog on unix??

2007-06-13 Thread BaKaLeGuM
I have on a script something like this toto = winDialog(yesno, Do you want to install the package) if (toto==YES){ but it dont work on unix because of the winDialog i think.. how can i do to change this for unix please? best regards vincent [[alternative HTML version deleted]]

Re: [R] Read Windows-like .INI files into R data structure?

2007-06-13 Thread Uwe Ligges
Vladimir Eremeev wrote: Christophe Pallier wrote: var1=value1, A=value3 is almost pure R code. Is it possible to use this feature to solve the problem? Along the same lines: you may write a short script that converts the ini file into R code that can be sourced. From your example, you

Re: [R] update packages with R on Vista: error

2007-06-13 Thread Stefan Grosse
Dear all, I figured out that the best solution for me is to install R in the /User/me/Documents/R/ directory instead of the program files directory since I am the only user on my notebook and R has then full write access and I do not need to change any environment variables... Thanks for the

Re: [R] Problems with Vista, R 2.5.0 and function save

2007-06-13 Thread Stefan Grosse
Hi Maja, if you are working on your own computer an alternative could be to install R in your user directory: c:\Users\your_loginname_\R or c:\Users\your_loginname_\Documents\R I also had a problem with the writing permissions when I tried updating R packages. Btw. when I work with files I set

Re: [R] Linux equivalent to windows menus and script window

2007-06-13 Thread Peter McMahan
Note: The screenshots at http://rosuda.org/JGR/ don't show the menu bar because they're taken on a mac. But the menu items do exist. On Jun 13, 2007, at 11:02 AM, Stefan Grosse wrote: JGR is pretty close to the windows R-GUI (with the menus) but even better since it is system independent

Re: [R] equivalent of windialog on unix??

2007-06-13 Thread Prof Brian Ripley
?readline On Wed, 13 Jun 2007, BaKaLeGuM wrote: I have on a script something like this toto = winDialog(yesno, Do you want to install the package) if (toto==YES){ but it dont work on unix because of the winDialog i think.. how can i do to change this for unix please? best regards

Re: [R] export to a dat file that SAS can read

2007-06-13 Thread Michael Dewey
At 09:05 13/06/2007, Rina Miehs wrote: Hello i have a data frame in R that some SAS users need to use in their programs, and they want it in a dat file, is that possible? and which functions to use for that? Does library(foreign) ?write.foreign get you any further forward? my data frame is

[R] installing GRASS-R packages

2007-06-13 Thread Philippe LAGACHERIE
Hi, I tried to install R packages required for the GRASS-R interface by using the following command (copied from http://www.geog.uni-hannover.de/grass/statsgrass/grass6_r_install.html): /install.packages (c(sp, spgrass6,rgdal,maptools), dependencies =TRUE)/ rgdal package was installed

[R] installing Rgraphviz under fedora 5

2007-06-13 Thread marco.R.help marco.R.help
Dear list, I have a lot of troubles installing Rgraphviz. I installed graphviz 2.13 from graphviz-2.13.20061222.0540.tar I installed the library Rgraphviz getBioC(Rgraphviz) Running biocinstall version 2.0.8 with R version 2.5.0 Your version of R requires version 2.0 of Bioconductor. trying

[R] Fitted Value Pareto Distribution

2007-06-13 Thread livia
I would like to fit a Pareto Distribution and I am using the following codes. I thought the fitted (fit1) should be the fitted value for the data, is it correct? As the result of the fitted turns out to be a single value for all. fit=vglm(ycf1 ~ 1, pareto1(location=alpha), trace=TRUE, crit=c)

Re: [R] Viewing a data object

2007-06-13 Thread Liaw, Andy
I believe JGR has an object browser. See the screenshots at the bottom of http://rosuda.org/JGR/. Andy From: Stephen Tucker Hi Horace, I have also thought that it may be useful but I don't know of any Object Explorer available for R. However, (you may alread know this but) (1)

Re: [R] specify constraints in maximum likelihood

2007-06-13 Thread Ben Bolker
adschai at optonline.net writes: Hi,I know only mle function but it seems that in mle one can only specify the bound of the unknowns forming the likelihood function. But I would like to specify something like, a = 2b or a = 2b where 'a' and 'b' could be my parameters in the likelihood

Re: [R] Design library installation problem

2007-06-13 Thread John Kane
I cannot find the posting but I believe Brian Ripley posted something here a day or so ago that said that Design had not passed some of the 2.5.0 tests. --- Ian Watson [EMAIL PROTECTED] wrote: Dear Listers I have tried to install Frank Harrell's two libaries: Hmisc and Design. I

Re: [R] Design library installation problem

2007-06-13 Thread Frank E Harrell Jr
Uwe Ligges wrote: Ian Watson wrote: Dear Listers I have tried to install Frank Harrell's two libaries: Hmisc and Design. I found that Hmisc was listed in the list of packages from the Install Packages command on the Packages menu, but Design was not. I installed Hmisc from this list,

Re: [R] installing Rgraphviz under fedora 5

2007-06-13 Thread Prof Brian Ripley
The problem is that libgvc.so.3 is not in your dynamic library search path. That's true by default for FC5, it is in /usr/lib64/graphviz on a x86_64 FC5 box. The proper solution is to fix the FC5 graphviz installation to tell ld.so about the location. A simpler way is to add

Re: [R] installing Rgraphviz under fedora 5

2007-06-13 Thread Michael Lawrence
On 6/13/07, marco.R.help marco.R.help [EMAIL PROTECTED] wrote: Dear list, I have a lot of troubles installing Rgraphviz. I installed graphviz 2.13 from graphviz-2.13.20061222.0540.tar I installed the library Rgraphviz getBioC(Rgraphviz) Running biocinstall version 2.0.8 with R version

Re: [R] if statement

2007-06-13 Thread Gabor Grothendieck
See ?quantcut in the gtools package. On 6/11/07, Jiong Zhang, PhD [EMAIL PROTECTED] wrote: Hi all, I have a rather naive question. I have the height of 100 individuals in a table and I want to assign the tallest 30% as Case=1 and the bottom 30% as Case=0. How do I do that? thanks. jiong

Re: [R] Confusion with sapply

2007-06-13 Thread Gabor Grothendieck
Try this. It takes a Date class date and in the first line month.day.year converts unclass(x) to chron. In the last line of the function we convert back to Date class. Its already vectorized so sapply is unneeded: library(chron) f - function(x) with(month.day.year(unclass(x)), { month -

Re: [R] equivalent of windialog on unix??

2007-06-13 Thread BaKaLeGuM
this is the answer :) toto - tkmessageBox(title = Rcmdr, message = Do you want to install Rcmdr?, icon = question, type = yesno) if (as.character(toto)==yes){ it work on windows and unix (and i think mac) 2007/6/13, Prof Brian Ripley [EMAIL PROTECTED]: ?readline On Wed, 13 Jun 2007,

[R] Differences of correlations testing in R

2007-06-13 Thread rafael
Good morning everybody, I need to compare intercorrelations between variables. My data were collected from 4 samples (with various number of subject) containing 4 variables scores. Some of my hypothesis are about the strength of relations (the sign doesn't matter) between variables across

Re: [R] Confusion with sapply

2007-06-13 Thread Patnaik, Tirthankar
Hi, Many thanks for this Hadley, and Uwe, the packages I used were chron, and zoo. Later I'm also using Hadley's reshape. I was able to run the code for a vector thus: dc - seq(as.Date(2006-01-01),len=10,by=month) dc [1] 2006-01-01 2006-02-01 2006-03-01 2006-04-01 2006-05-01

Re: [R] installing GRASS-R packages

2007-06-13 Thread Roger Bivand
On Wed, 13 Jun 2007, Philippe LAGACHERIE wrote: Hi, I tried to install R packages required for the GRASS-R interface by using the following command (copied from http://www.geog.uni-hannover.de/grass/statsgrass/grass6_r_install.html): /install.packages (c(sp, spgrass6,rgdal,maptools),

[R] How to install RMySQL package in R 2.5 in Windows OS?

2007-06-13 Thread Ruixin ZHU
Dear R-users, It seems that install.packages( ) doesn't work to RMySQL package. Would anybody have the experience of that? Thanks _ Dr.Ruixin ZHU Shanghai Center for Bioinformation Technology [EMAIL PROTECTED] [EMAIL PROTECTED] 86-21-13040647832

Re: [R] export to a dat file that SAS can read

2007-06-13 Thread Thomas Lumley
On Wed, 13 Jun 2007, Uwe Ligges wrote: Rina Miehs wrote: Hello i have a data frame in R that some SAS users need to use in their programs, and they want it in a dat file, is that possible? What is a dat file? and which functions to use for that? I *guess* write.table() will do the

[R] Obtaining the cross validation coefficient of determination

2007-06-13 Thread Pietrzykowski, Matthew (GE, Research)
Hello - I am currently working through some MLR models and wanted to perform cross validation to assess the model's prediction error. As a Minitab user prior to discovering R, I was able to view the PRESS statistic and the resultant cross validation R^2. Is there a canned function already

Re: [R] R Book Advice Needed

2007-06-13 Thread ngottlieb
Roland: Thanks for your reply. I have sort of pay my dues with statistics and doing the hard math reading of Proofs. Years ago reading lots of books on Multi-variate Methods such As Principal Components, Cluster Analysis, Discriminant Analysis, Multi Dimensional Scaling(MDS), Optimization both

Re: [R] R is not a validated software package..

2007-06-13 Thread Terry Therneau
I've been on vacation and so come late to this interesting discussion. Let me add two minor points. 1. I have run into a lot of statements that x is required when dealing with pharma, and in particular wrt NDAs (new drug application). Almost all were false. But I also understand a bit of

Re: [R] ievent.wait

2007-06-13 Thread Greg Snow
Does locator(type='l') (or type ='b') Work for you? -Original Message- From: ryestone [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch r-help@stat.math.ethz.ch Sent: 6/8/07 10:19 AM Subject: [R] ievent.wait I am working on a plot and would be like to click on a few points and then

Re: [R] R Book Advice Needed

2007-06-13 Thread Roland Rau
Hi, [EMAIL PROTECTED] wrote: I am new to using R and would appreciate some advice on which books to start with to get up to speed on using R. My Background: 1-C# programmer. 2-Programmed directly using IMSL (Now Visual Numerics). 3- Used in past SPSS and Statistica. I put together a

Re: [R] barplot and map overlay

2007-06-13 Thread Greg Snow
Look at the subplot function in the TeachingDemos package. -Original Message- From: Héctor Villalobos [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch r-help@stat.math.ethz.ch Sent: 6/11/07 5:48 PM Subject: [R] barplot and map overlay Hi, I wonder if it is possible with the graphics

Re: [R] installing Rgraphviz under fedora 5

2007-06-13 Thread Seth Falcon
marco.R.help marco.R.help [EMAIL PROTECTED] writes: Dear list, I have a lot of troubles installing Rgraphviz. I installed graphviz 2.13 from graphviz-2.13.20061222.0540.tar I installed the library Rgraphviz I'm pretty sure that you will have problems with graphviz 2.13 and that you will

Re: [R] R Book Advice Needed

2007-06-13 Thread ngottlieb
Thanks Alain. Guess bite the bullet with limited budget buy bunch From Amazon and see what reads best and return the rest!. One ends up collecting so many books (most of bought 5 books on Bayesian analysis years ago), still like browsing shelfs! Regards, Neil -Original Message-

Re: [R] Fitted Value Pareto Distribution

2007-06-13 Thread Markus Jäntti
livia wrote: I would like to fit a Pareto Distribution and I am using the following codes. I thought the fitted (fit1) should be the fitted value for the data, is it correct? As the result of the fitted turns out to be a single value for all. fit=vglm(ycf1 ~ 1, pareto1(location=alpha),

Re: [R] R Book Advice Needed

2007-06-13 Thread ngottlieb
Cody: Think you might have asked the question for me Neil. I do time series analysis of return data in finance. I will be creating a factor model based on PCA Or Single Value Decomposition to get Eigenvectors Of the correlation matrix (tends to work better for finance data Than covariance).

Re: [R] R Book Advice Needed

2007-06-13 Thread Roland Rau
Hi Neil, [EMAIL PROTECTED] wrote: At this point, want to jump in avoiding all the Mathematical proofs and just apply R and the packages for what I want to do. I'd still recommend Venables/Ripley: Modern Applied Statistics with S (or often abbrev. MASS, which is also name of the package

Re: [R] equivalent of windialog on unix??

2007-06-13 Thread Prof Brian Ripley
On Wed, 13 Jun 2007, BaKaLeGuM wrote: this is the answer :) toto - tkmessageBox(title = Rcmdr, message = Do you want to install Rcmdr?, icon = question, type = yesno) if (as.character(toto)==yes){ it work on windows and unix (and i think mac) Only if you have tcl/tk installed and an

Re: [R] R Book Advice Needed

2007-06-13 Thread Peter Dalgaard
Roland Rau wrote: Hi, [EMAIL PROTECTED] wrote: I am new to using R and would appreciate some advice on which books to start with to get up to speed on using R. My Background: 1-C# programmer. 2-Programmed directly using IMSL (Now Visual Numerics). 3- Used in past SPSS and Statistica.

Re: [R] Awk and Vilno

2007-06-13 Thread Rogerio Porto
Hey, What we should really compare is the four situations: R alone R + awk R + vilno R + awk + vilno and maybe R + SAS Data step and see what scripts are more elegant (read 'short and understandable') what do you guys think of creating a R-wiki page for syntax comparisons among the

Re: [R] R Book Advice Needed

2007-06-13 Thread ngottlieb
Thanks Roland, fortunately I dug up MASS by Venables/Ripley buried under all my econometric and statistic books. Will be reading it today and order a few of the R books for additional support. Thanks for your suggestions... Regards, Neil -Original Message- From: Roland Rau

Re: [R] ievent.wait

2007-06-13 Thread Sundar Dorai-Raj
Hi, Greg, type = 'b' won't work according to ?locator. Try type = 'o'. HTH,x --sundar Greg Snow said the following on 6/13/2007 7:27 AM: Does locator(type='l') (or type ='b') Work for you? -Original Message- From: ryestone [EMAIL PROTECTED] To:

Re: [R] How to install RMySQL package in R 2.5 in Windows OS?

2007-06-13 Thread Prof Brian Ripley
On Wed, 13 Jun 2007, Ruixin ZHU wrote: Dear R-users, It seems that install.packages( ) doesn't work to RMySQL package. Under Windows, yes. You need the MySQL client libraries for your version of MySQL (or something very close to the same version), so the only safe way is to install from

Re: [R] pretty report

2007-06-13 Thread Weiwei Shi
I think my initial intention is to write multiple worksheets for multiple data frames. write.csv or write.table cannot do that. On 6/12/07, Robert A LaBudde [EMAIL PROTECTED] wrote: At 09:13 PM 6/12/2007, Don wrote: At 5:01 PM -0400 6/12/07, Weiwei Shi wrote: Dear Listers: I have a

Re: [R] Fitted Value Pareto Distribution

2007-06-13 Thread J. R. M. Hosking
livia wrote: I would like to fit a Pareto Distribution and I am using the following codes. I thought the fitted (fit1) should be the fitted value for the data, is it correct? As the result of the fitted turns out to be a single value for all. fit=vglm(ycf1 ~ 1, pareto1(location=alpha),

Re: [R] R Book Advice Needed

2007-06-13 Thread Gabor Grothendieck
There are some online sources that you might find useful. You could get started on those while you decide what books to get: - CRAN contributed documentation http://cran.r-project.org/other-docs.html - S Poetry http://www.burns-stat.com/pages/spoetry.html - Zoonekynd book

Re: [R] VGAM Pareto

2007-06-13 Thread Ben Bolker
livia yn19832 at msn.com writes: I would like to fit a Pareto Distribution and I am using the following codes fit=vglm(ycf1 ~ 1, pareto1(location=alpha), trace=TRUE, crit=c) fitted(fit) But the fitted values turn out to be the same for each observation. I guess the problem is with

Re: [R] R Book Advice Needed

2007-06-13 Thread ngottlieb
Pat: I have done PCA to extract eigenvectors on return series for equities. Rotation does help and does make factors more understandable, have had success doing this. You are right, when doing pure statistical factors, one tends to find first factor which explains most of the variance is the

Re: [R] Viewing a data object

2007-06-13 Thread Horace Tso
Stephen and Christophe, I'm aware of fix and edit and the few issues with fix. Thus my reluctance to use them. Emacs may be the way to go, but from what I heard here it has a steep learning curve. The autocompletion feature in 2.5.1 is great. Andy Liaw points me to JGR which I'm just about to

Re: [R] passing (or obtaining) index or element name of list to FUN in lapply()

2007-06-13 Thread Stephen Tucker
Hi Professor Ripley, Thanks for the response. I apologize, my examples were not too real (though your solutions are indeed clever)... I was trying to ask more generally whether the element name or index of 'listObj' could be obtained by the user-function 'myfunction' when used in

Re: [R] barplot and map overlay

2007-06-13 Thread Héctor Villalobos
Thank you Greg, It works! On 13 Jun 2007 at 8:27, Greg Snow wrote: Look at the subplot function in the TeachingDemos package. -Original Message- From: Héctor Villalobos [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch r-help@stat.math.ethz.ch Sent: 6/11/07 5:48 PM Subject: [R]

Re: [R] ievent.wait

2007-06-13 Thread ryestone
With locator( ) does it only work in a regular R plot or can I use it with iPlot? I am having difficulty getting it to be used with Iplots, it just calls up a new screen when the function is called. Stone. Sundar Dorai-Raj wrote: Hi, Greg, type = 'b' won't work according to ?locator. Try

Re: [R] Subscription

2007-06-13 Thread John Kane
--- Lucy Namu [EMAIL PROTECTED] wrote: I would like to subscribe and get free software for statistical analysis. Lucy It looks like you are subscribed. To download R http://www.r-project.org/. Some other sources of free software http://www.epidata.dk/

Re: [R] extractor rows from a matrix

2007-06-13 Thread John Kane
ex1 - ht[,1]and so on? Have a look at Chapter 5 in the Introduction to R --- billycorg [EMAIL PROTECTED] wrote: hi! i have a little problem: my data's matrix has 1093 rows and 3 columns. i'd like to extract each rows.. something like this: ht= my matrix Dt=(???)=a vector with

[R] Removing Inf and Inf values from a fata frame

2007-06-13 Thread Judith Flores
Hi, I have a csv file with empty values, when I apply the different functions (mean, std, etc.) I create a new data frame, the empty values generate Inf and -Inf values. How can I remove those Inf and -Inf values from the new data frame? I already specified na.rm in the mean and std

Re: [R] Removing Inf and Inf values from a fata frame

2007-06-13 Thread Henrique Dallazuanna
Hi, try with df[!is.infinite(your_column_in_d.f.),] -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22 Ohttp://maps.google.com/maps?f=qhl=enq=Curitiba,+Brazillayer=ie=UTF8z=18ll=-25.448315,-49.276916spn=0.002054,0.005407t=kom=1 On 13/06/07, Judith Flores [EMAIL PROTECTED]

Re: [R] Read Windows-like .INI files into R data structure?

2007-06-13 Thread Earl F. Glynn
[EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED]... .Ini files are, for lack of a better description, ancient. In this case a device is creating the INI files as part of an experiment, so the file format cannot be changed (at least easily). I've looked at XML files from time to time

Re: [R] Removing Inf and Inf values from a fata frame

2007-06-13 Thread Petr Klasterecky
?is.finite for removing the Infs *from the dataframe*, this does not aviod creating them by mean() or std(). Petr Judith Flores napsal(a): Hi, I have a csv file with empty values, when I apply the different functions (mean, std, etc.) I create a new data frame, the empty values

[R] lme() doesn't converge on IGF example

2007-06-13 Thread David Daniel
Running the Chapter 4 examples in Pinheiro Bates' Mixed-Effects Models in S and S-PLUS (2000), I get a message that the default optimizer doesn't converge, but using optim for the optimizer results in convergence: library(nlme) fm1IGF.lis - lmList(IGF) fm1IGF.lme - lme(fm1IGF.lis)

[R] Normal and Poisson tail area expectations in R

2007-06-13 Thread kavindra malik
I am interested in R functions for the following integrals / sums (expressed best I can in text) - Normal: G_u(k) = Integration_{Lower limit=k}^{Upper limit=infinity} [(u -k) f(u) d(u)], where where u is N(0,1), and f(u) is the density function. Poisson: G(lambda,k) = Sum_{Lower

Re: [R] lme() doesn't converge on IGF example

2007-06-13 Thread Douglas Bates
On 6/13/07, David Daniel [EMAIL PROTECTED] wrote: Running the Chapter 4 examples in Pinheiro Bates' Mixed-Effects Models in S and S-PLUS (2000), I get a message that the default optimizer doesn't converge, but using optim for the optimizer results in convergence: library(nlme)

Re: [R] Read Windows-like .INI files into R data structure?

2007-06-13 Thread ngottlieb
Earl: Really depends on the need. XML yes can get crazy (having had to deal with some ugly XML). One can do a correctly formatted XML, that parses via the DOM which does not mean well formatted XML. It's all a matter of design and data structures. XML advantages: one can define own data types

Re: [R] Awk and Vilno

2007-06-13 Thread Tim Churches
Rogerio Porto wrote: Hey, What we should really compare is the four situations: R alone R + awk R + vilno R + awk + vilno and maybe R + SAS Data step and see what scripts are more elegant (read 'short and understandable') I don't think that short and understandable necessarily go

[R] how to optionally include variables in a data.frame at assignment

2007-06-13 Thread Jenny Bryan
I am creating a data.frame inside a function and the set of variables to include depends on the current value of other variables. Is there a way to accomplish this in the original assignment? Or must I first create the core data.frame with the variables I always want and then use if

Re: [R] passing (or obtaining) index or element name of list to FUN in lapply()

2007-06-13 Thread Charles C. Berry
This sounds like a job for mapply, viz: a.list - list(a=cbind(1:4,rnorm(4)),b=cbind(4:1,rnorm(4))) plot.x.main.y - function(x,y,...) plot( x, main=y, ... ) mapply( plot.x.main.y , a.list, names(a.list) ) Chuck On Wed, 13 Jun 2007, Stephen Tucker wrote: Hi Professor Ripley, Thanks for

Re: [R] Normal and Poisson tail area expectations in R

2007-06-13 Thread Charles C. Berry
On Wed, 13 Jun 2007, kavindra malik wrote: I am interested in R functions for the following integrals / sums (expressed best I can in text) - Normal: G_u(k) = Integration_{Lower limit=k}^{Upper limit=infinity} [(u -k) f(u) d(u)], where where u is N(0,1), and f(u) is the density function.

Re: [R] R vs. Splus in Pharma/Devices Industry

2007-06-13 Thread Cody_Hamilton
I should have also noted that Sweave is available for use with R. This is offset, however, by the fact that I will probably never be able to convince anyone to use Latex. This is a pity as I often find myself admiring reports done in Latex as opposed to the ones I have worked on in MS Word.

Re: [R] Normal and Poisson tail area expectations in R

2007-06-13 Thread kavindra malik
Thank you very much. This solves the problem I was trying to solve. I am new to R and am learning. A great lesson in the power of R... Charles C. Berry [EMAIL PROTECTED] wrote: On Wed, 13 Jun 2007, kavindra malik wrote: I am interested in R functions for the following integrals / sums

Re: [R] Formatted Data File Question for Clustering -Quickie Project

2007-06-13 Thread AA
if you look at the data USArrests by doing data(USArrets) USArrets you will see that it is a data.frame. so by analogy you could do the following: Probably you have this data in Excel (I guess from the format in your mail). Have the data in a sheet as: convertsshortBais

[R] Difficulties With Posting To Ongoing Threads on the R Mailing List

2007-06-13 Thread Robert Wilkins
A number of people are having the same problem as me, when you post as a response to an ongoing thread, in place of your message, the following message appears: An embedded charset-unspecified text was scrubbed ... and a link is given that leads to the desired message. It's better than nothing

[R] Offline ? Searching for James Wettenhall's TclTk Examples

2007-06-13 Thread Peter Ruckdeschel
Hi, as a starting point for using Tcl/Tk in R, I used to refer to James Wettenhall's nicely presented TclTk Examples formerly hosted at http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/ These days I have been trying to reach these pages but without success. Does anyone know (James,

[R] Where to Find Data Transformation Software

2007-06-13 Thread Robert Wilkins
Hello All, Here is the requested information. Most of it was on the original post for the Tools For Preparing Data For Analysis thread from last week, but it got overlooked. They are all given under an open source license. Check 'em out! ***

Re: [R] Normal and Poisson tail area expectations in R

2007-06-13 Thread Ravi Varadhan
More interesting is the Poisson convolution. I don't know if there is an analytic solution to this. I looked at Jolley's Summation of Series and Abramowitz and Stegun, but no help there. It seems that discrete FFT technique should work. Does anyone know the answer? Ravi.

  1   2   >