Re: [R] exporting a table to latex

2006-12-13 Thread Dieter Menne
Rainer M Krug RKrug at sun.ac.za writes: I am using the latex() command from the Hmisc package to export table to latex. For formating of the values, I use the format() function. But I don't manage to get the format of the values right: I don't want scientific notation, but if I say

Re: [R] Problem with sas.get function in Hmisc

2006-12-13 Thread Jean Vidal
The workaround you point seems to be the solution for the moment. I not sure I will have time (and courage) to try to figure out where the bug comes from. Hope a more literate R programmer than me is willing to dive in... Thank you for your answer. 2006/12/12, Agustin Calatroni [EMAIL

[R] Problem in Converting Zoo Objects to Dataframes

2006-12-13 Thread Shubha Vishwanath Karanth
Hi R experts, The below is an RBloomerg command. The object intra here is a zoo object. I need to convert this zoo object into a data frame, called bb. library(zoo) library(chron) library(RDCOMClient) library(RBloomberg)

Re: [R] exporting a table to latex

2006-12-13 Thread Rainer M Krug
Dieter Menne wrote: Rainer M Krug RKrug at sun.ac.za writes: I am using the latex() command from the Hmisc package to export table to latex. For formating of the values, I use the format() function. But I don't manage to get the format of the values right: I don't want scientific

[R] Install R in Linux

2006-12-13 Thread lu kan
Hi, Is it possible to install R in a linux box (Debian) without being a root. I know I can compile the R source code, but there is no F77 compiler on the box. So is it possible to install binary R without being a root? Send instant messages to your online friends http://uk.messenger.yahoo.com

[R] ifelse misusage become more and more frequent...

2006-12-13 Thread Martin Maechler
[EMAIL PROTECTED] == [EMAIL PROTECTED] fr [EMAIL PROTECTED] on Tue, 12 Dec 2006 22:24:33 +0100 writes: [EMAIL PROTECTED] ...ifelse, a function of three **vector** [EMAIL PROTECTED] arguments Yes !! I misunderstood the [EMAIL PROTECTED] functioning of ifelse. Seems to

[R] Questions about saving output files and popup

2006-12-13 Thread u9470002
Hi all: After calculating in R I want to show the answer and some explanations not graphic plots in another new device and then save it as txt.file. However, I couldn't find any package or command to do it yet. I know some commands about generating graphics on different types of display or

[R] persp() problem

2006-12-13 Thread Michael Sun
Dear list, I have a problem on persp() x - u1data #first coloum in attached data y - u2data #second coloum in attached data f - function(x,y){qgev(pnorm(rhoF*qnorm(pnorm((qnorm(y)-rho2*qnorm(x)/sqrt(1-rho2^2 +sqrt(1-rhoF^2)*qnorm(0.95)),-0.3935119, 0.4227890, 0.2701648)}

Re: [R] Sweave, Xfig, pdflatex and \setkeys

2006-12-13 Thread David Lindelöf
On Tue, 2006-12-12 at 13:22 -0500, Kevin E. Thorpe wrote: Trouble is that Sweave defines (with \setkeys) the default width of \includegraphics to be 0.8 times the \textwidth. The result is that the graphic is scaled, but not the text. I was looking for a way to temporarily undefine the

Re: [R] expression()

2006-12-13 Thread javier garcia-pintado
Thanks you, this works right. I just would like to note one thing: although I've found italic() in the help for plotmath and I can see that the italic function: mtext(expression(italic(beta)[max]),side=1,line=2) does not work on the greek beta character. Though not strictly necessary, this

[R] Passing arguments to panels in trellis plots

2006-12-13 Thread Marco Chiarandini
Dear all, I am trying to produce survfit plots in a trellis environment and I would like the plots to be logarithmic. I am trying this: print(Ecdf(~time | size*type, groups=alg,data=B,subscripts=TRUE, panel=function(x,groups,subscripts) { t -

Re: [R] How to sum one column in a data frame keyed on other columns

2006-12-13 Thread Alex Brown
Hi George I'm running ubuntu dapper badger with 2.4.0. add the line deb http://cran.R-project.org/bin/linux/ubuntu dapper/ to your /etc/apt/sources.list -Alex Brown On 13 Dec 2006, at 02:20, George Nachman wrote: Thanks, everyone, for the help! Bill: This looks like a really

Re: [R] Questions about saving output files and popup

2006-12-13 Thread Chuck Cleland
u9470002 wrote: Hi all: After calculating in R I want to show the answer and some explanations not graphic plots in another new device and then save it as txt.file. However, I couldn't find any package or command to do it yet. I know some commands about generating graphics on different

Re: [R] persp() problem

2006-12-13 Thread Duncan Murdoch
On 12/13/2006 6:19 AM, Michael Sun wrote: Dear list, I have a problem on persp() x - u1data #first coloum in attached data y - u2data #second coloum in attached data f - function(x,y){qgev(pnorm(rhoF*qnorm(pnorm((qnorm(y)-rho2*qnorm(x)/sqrt(1-rho2^2

Re: [R] expression()

2006-12-13 Thread Prof Brian Ripley
On Wed, 13 Dec 2006, javier garcia-pintado wrote: Thanks you, this works right. I just would like to note one thing: although I've found italic() in the help for plotmath and I can see that the italic function: mtext(expression(italic(beta)[max]),side=1,line=2) does not work on the

[R] Error to install fMultivar package

2006-12-13 Thread gsmatos1
Hi, I tried to install fMultivar package but an error occurs that I could not understand. (I've been worked with linux / Ubuntu 6.06 LTS) install.packages(fMultivar) Warning in install.packages(fMultivar) : argument 'lib' is missing: using /usr /local/lib/R/site-library --- Please select

[R] (no subject)

2006-12-13 Thread Francisco J Molina
Hi, Let us suppose I have a list x = list () x $ name1 = 1 x $ name2 = 'a' in the work environment. Let us suppose that in the body of a function I want to acces to a component of x by using its name as argument of that function. How can this by done? For instance, I was expecting f =

[R] Install R in Linux

2006-12-13 Thread lu kan
Hi, Is it possible to install R in a linux box (Debian) without being a root. I know I can compile the R source code, but there is no F77 compiler on the box. So is it possible to install binary R without being a root? Send instant messages to your online friends http://uk.messenger.yahoo.com

Re: [R] (no subject)

2006-12-13 Thread Scionforbai
Just read some introductory R tutorial... x = list () x $ name1 = 1 x $ name2 = 'a' x $name1 [1] 1 $name2 [1] a name - name1 x[name] $name1 [1] 1 x[[name]] [1] 1 name - name2 x[name] $name2 [1] a x[[name]] [1] a Bye, Scionforbai __

Re: [R] Install R in Linux

2006-12-13 Thread Peter Dalgaard
lu kan wrote: Hi, Is it possible to install R in a linux box (Debian) without being a root. I know I can compile the R source code, but there is no F77 compiler on the box. So is it possible to install binary R without being a root? (We did see it the first time!!) In a word, no. Either

Re: [R] (no subject)

2006-12-13 Thread Uwe Ligges
Francisco J Molina wrote: Hi, Let us suppose I have a list x = list () x $ name1 = 1 x $ name2 = 'a' in the work environment. Let us suppose that in the body of a function I want to acces to a component of x by using its name as argument of that function. How can this by done?

Re: [R] Error to install fMultivar package

2006-12-13 Thread Scionforbai
Just a few tip, I don't know very much about ubuntu. Relevant is: /usr/bin/ld: cannot find -lblas-3 It means, it doesn't find the library blas-3 (libblas-3). Maybe they are installed on your system but in a non-standard dir. Look around (locate blas, if updatedb is set on your system); if you

Re: [R] ifelse misusage become more and more frequent...

2006-12-13 Thread Michael Dewey
At 10:42 13/12/2006, Martin Maechler wrote: [EMAIL PROTECTED] == [EMAIL PROTECTED] fr [EMAIL PROTECTED] on Tue, 12 Dec 2006 22:24:33 +0100 writes: [EMAIL PROTECTED] ...ifelse, a function of three **vector** [EMAIL PROTECTED] arguments Yes !! I misunderstood the [EMAIL

[R] R_Code

2006-12-13 Thread Ilka Schmidt
Hello! I need the Code of the functions runif and rexp. Where can I get them? Can You help me? I thank you for an answer! bye Ilka -- __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] upside down image/data

2006-12-13 Thread Thomas Harte
rm(list=ls(all=TRUE)) graphics.off() # make a test matrix: nr- 3 nc- 4 # the data: ( m- matrix((1:(nr*nc)), nr, nc) ) [,1] [,2] [,3] [,4] [1,]147 10 [2,]258 11 [3,]36

Re: [R] R_Code

2006-12-13 Thread Duncan Murdoch
On 12/12/2006 11:39 AM, Ilka Schmidt wrote: Hello! I need the Code of the functions runif and rexp. Where can I get them? Can You help me? I thank you for an answer! bye Ilka Just type the names. For example: runif function (n, min = 0, max = 1) .Internal(runif(n, min, max))

Re: [R] upside down image/data

2006-12-13 Thread Jenny Barnes
Thomas, Thank you for this example, makes it easier to see what levelplot does - does this mean that EVERY time I want to plot with levelplot() I have to not only reverse the columns [,ncol(output.temp):1] but also have to transform the matrix as below? I am only suprised as I don't remember

[R] Test for Kurtosis.

2006-12-13 Thread Megh Dal
Dear all R users, I question is on the subject of Statistics rather R related. I am searching some documentation on statistical significance test for Kurtosis especially Glyn test. More specifically I am looking for documentation on: 1: Test statistic for Glyn test 2. Statistical

[R] Skipping Results with Errors in a Loop

2006-12-13 Thread Rick Bilonick
I'm estimating models using lme in a for loop. Sometimes the model doesn't converge or there are other problems. This causes the evaluation to stop prematurely. I can't remember the function name that I need to use to allow the loop to continue until the end. Could someone remind me the name of

Re: [R] Skipping Results with Errors in a Loop

2006-12-13 Thread Antonio, Fabio Di Narzo
?try 2006/12/13, Rick Bilonick [EMAIL PROTECTED]: I'm estimating models using lme in a for loop. Sometimes the model doesn't converge or there are other problems. This causes the evaluation to stop prematurely. I can't remember the function name that I need to use to allow the loop to

Re: [R] upside down image/data

2006-12-13 Thread Thomas Harte
the transform that i provided orientates the data matrix so that when plotted with image or levelplot the result is isomorphic to what you see when you print the matrix at the r prompt. i don't know what your data look like---commented, minimal, self-contained, reproducible code would

[R] counting process and output of survfit

2006-12-13 Thread A. Gunes Koru
Hello, I put my data in conditional counting process format basically by adding a stratum variable to distinguish between observations. I fit a model and then used survfit function. These are the first five lines of output. time n.risk n.event survival std.err lower 95% CI upper 95% CI

Re: [R] upside down image/data

2006-12-13 Thread Jenny Barnes
Thanks Thomas, My data arrays each contain 0.5million data points so I couldn't really reproduce them unfortunately. Next time I will try and offer some exapmle code simplified with comments in order to help you and the others on R-help understand my problem more easily. I appreciate your help

Re: [R] persp() problem

2006-12-13 Thread Duncan Murdoch
On 12/13/2006 9:57 AM, Michael Sun wrote: Dear Duncan, Thank you very much for your quick reply. The comands below really worked testx - seq(min(u1data), max(u1data), len=100) testy - seq(min(u2data), max(u2data), len=100) And you are right, I do need my observed data. can trans3d plot

Re: [R] set up directory for R when I start R

2006-12-13 Thread Earl F. Glynn
Duncan Murdoch [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] On 12/11/2006 2:21 PM, Earl F. Glynn wrote: Dieter Menne [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] Aimin Yan aiminy at iastate.edu writes: I want to set default directory for R when I start R. How to

Re: [R] Error to install fMultivar package

2006-12-13 Thread Dirk Eddelbuettel
On Wed, Dec 13, 2006 at 09:34:04AM -0300, gsmatos1 wrote: Content-Description: Mail message body Hi, I tried to install fMultivar package but an error occurs that I could not understand. (I've been worked with linux / Ubuntu 6.06 LTS) This would be quicker: $ sudo apt-get install

Re: [R] exporting a table to latex

2006-12-13 Thread Rau, Roland
Hi, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Rainer M Krug Sent: Wednesday, December 13, 2006 8:47 AM To: r-help@stat.math.ethz.ch Subject: [R] exporting a table to latex Is there another way of formating the values to three decimal

Re: [R] Obtaining Estimates of the Random Effects Parameters

2006-12-13 Thread Doran, Harold
?VarCorr -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Rick Bilonick Sent: Wednesday, December 13, 2006 11:40 AM To: R Help Subject: [R] Obtaining Estimates of the Random Effects Parameters I'm running simulation using lme and sometimes the

[R] Obtaining Estimates of the Random Effects Parameters

2006-12-13 Thread Rick Bilonick
I'm running simulation using lme and sometimes the estimated variance-covariance matrix is not positive definite so that the intervals function won't work for the random effect coefficients. I've tried varcomp from the ape package but this does not return all the coefficients. How can I extract

Re: [R] exporting a table to latex

2006-12-13 Thread Rainer M Krug
Rau, Roland wrote: Hi, -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Rainer M Krug Sent: Wednesday, December 13, 2006 8:47 AM To: r-help@stat.math.ethz.ch Subject: [R] exporting a table to latex Is there another way of formating the values

[R] why are my multiple box whisker plots so small?

2006-12-13 Thread Ellen Husain
Hi all, I'm sure this is reallly basic, but I just can get it to work. I want to plot six box whisker plots together to make one figure so that they appear one below the next. I can do this using par(mfrow=c(6,1)), but each boxwhisker plots end up vertically compressed to the point where I

Re: [R] why are my multiple box whisker plots so small?

2006-12-13 Thread rolf
I'm sure this is reallly basic, but I just can get it to work. I want to plot six box whisker plots together to make one figure so that they appear one below the next. I can do this using par(mfrow=c(6,1)), but each boxwhisker plots end up vertically compressed to the point where I can't

[R] Slightly off-topic, but probably funny to this group...

2006-12-13 Thread David L. Van Brunt, Ph.D.
Nothing like Grad students expressing themselves through music... *statz rappers*http://video.google.com/videoplay?docid=489221653835413043pr=goog-sl http://video.google.com/videoplay?docid=489221653835413043pr=goog-sl -- --- David L. Van Brunt, Ph.D.

Re: [R] why are my multiple box whisker plots so small?

2006-12-13 Thread Francisco J. Zagmutt
Please remember to add a sample of reproducible code when you post to this list (as recommended in the posting guide). You probably want to use the lattice package i.e. library(lattice) ?bwplot Francisco Ellen Husain wrote: Hi all, I'm sure this is reallly basic, but I just can get it

[R] install.packages

2006-12-13 Thread Aimin Yan
I try to type this in my R-winEdt. but I got these. Do you know? Aimin install.packages('http://rh-mirror.linux.iastate.edu/CRAN/bin/windows/contrib/2.4/plotrix_2.1-6.zip') Warning in download.packages(pkgs, destdir = tmpd, available = available, : no package

Re: [R] expression()

2006-12-13 Thread Paul Murrell
Hi javier garcia-pintado wrote: Thanks you, this works right. I just would like to note one thing: although I've found italic() in the help for plotmath and I can see that the italic function: mtext(expression(italic(beta)[max]),side=1,line=2) does not work on the greek beta

Re: [R] install.packages

2006-12-13 Thread Benilton Carvalho
install.packages(plotrix) b. On Dec 13, 2006, at 1:42 PM, Aimin Yan wrote: I try to type this in my R-winEdt. but I got these. Do you know? Aimin install.packages('http://rh-mirror.linux.iastate.edu/CRAN/bin/ windows/contrib/2.4/plotrix_2.1-6.zip') Warning in download.packages(pkgs,

Re: [R] Passing arguments to panels in trellis plots

2006-12-13 Thread Deepayan Sarkar
On 12/13/06, Marco Chiarandini [EMAIL PROTECTED] wrote: Dear all, I am trying to produce survfit plots in a trellis environment and I would like the plots to be logarithmic. I am trying this: print(Ecdf(~time | size*type, groups=alg,data=B,subscripts=TRUE,

[R] Course***January-February 2007*** R/Splus Fundamentals and Programming Techniques @ 7 cities: New York, San Diego, Washington DC, San Francisco, Chicago, Seattle, Princeton

2006-12-13 Thread elvis
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce our January-February 2007 R/S-plus Fundamentals and Programming Techniques : www.xlsolutions-corp.com/Rfund.htm *** New York / January 8-9, 2007 *** San Diego / January 8-9, 2007 *** Washington DC / January 18-19, 2007 ***

Re: [R] Passing arguments to panels in trellis plots

2006-12-13 Thread Marco Chiarandini
Dear Deepayan, Since you haven't bothered to follow the posting guide at all, it took me a while to figure out that you live in that alternate R universe created by Prof Harrell. Can't help you there, but in the standard universe, things seem fairly simple: I am sorry for not having been

[R] How to avoid test for NAs in foreign function call

2006-12-13 Thread Knut M. Wittkowski
We have packed logical vectors into integers, 32 flags at a time and then want to AND or OR these vectors of integers using other C functions. The problem: occasionally, the packed sequence of 32 logical values resembles NA, causing the error message: Error in bitAND(packed1, packed2, lenx) :

Re: [R] Problem with sas.get function in Hmisc

2006-12-13 Thread Hong Ooi
___ After digging around in the sas.get code, I found a change that seems to be causing the problem: old version: ... status - sys(paste(sasprog, sasin, -log, log.file), output = FALSE) ...

[R] Extracting tolerance in R?

2006-12-13 Thread Keith . Chamberlain
Dear list, How is the tolerance for a model parameter in an lm() call extracted? I did not see a solution in the documentation for lm(), or predict(), nor in the archives using 'tolerance' as the search string. I also checked into the nlme package, though nothing popped out at me. Sincerely,

[R] Obiously obsolete

2006-12-13 Thread J Dougherty
On Tuesday 12 December 2006 07:59, javier garcia-pintado wrote: beta max Sorry, couldn't resist. JD __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] Reverse order of grouping factor in grouppedData

2006-12-13 Thread RMan54
I created the following groupedData object (nlme library): gd - groupedData(Conc ~ Time | Subj, order.groups=T, FUN = myf, data=mydata) The idea of the myf function is to reverse the order of the grouping factor Subj (or better, reorder from largest to smallest). In the mydata data

Re: [R] Reverse order of grouping factor in grouppedData

2006-12-13 Thread Deepayan Sarkar
On 12/13/06, RMan54 [EMAIL PROTECTED] wrote: I created the following groupedData object (nlme library): gd - groupedData(Conc ~ Time | Subj, order.groups=T, FUN = myf, data=mydata) The idea of the myf function is to reverse the order of the grouping factor Subj (or better,

Re: [R] Extracting tolerance in R?

2006-12-13 Thread Chuck Cleland
[EMAIL PROTECTED] wrote: Dear list, How is the tolerance for a model parameter in an lm() call extracted? I did not see a solution in the documentation for lm(), or predict(), nor in the archives using 'tolerance' as the search string. I also checked into the nlme package, though

Re: [R] Passing arguments to panels in trellis plots

2006-12-13 Thread Deepayan Sarkar
On 12/13/06, Marco Chiarandini [EMAIL PROTECTED] wrote: Dear Deepayan, Since you haven't bothered to follow the posting guide at all, it took me a while to figure out that you live in that alternate R universe created by Prof Harrell. Can't help you there, but in the standard universe,

[R] strategy for doing an ANOVA on unbalanced data

2006-12-13 Thread Steven Lacey
Hi, I would like some help deciding if and how to average my data before running an ANOVA. Let me first describe the data and what makes it unique. Hopefully, that will generate some ideas because I am not sure what kind of model I need to use. I don't know how to describe this succinctly, so

[R] Syntax for getting covariance matrix

2006-12-13 Thread Elena
I am new at using R and was wondering if someone can give me the exact syntax for getting the covariance matrix. I used logit to generate the equation. Specifically the syntax was: glm(formula = Status ~ A + B + C + D + E, family = binomial(link = logit), data = logit) Do you know

Re: [R] Extracting tolerance in R?

2006-12-13 Thread KeithC.
Dear Chuck, Thank you. Vif() is exactly what I needed! Sincerely, KeithC. __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

[R] subset question

2006-12-13 Thread Aimin Yan
I have a data set p1982, its structure is the following Then I take 20 observations from this dataset, and assign to pr. in p1982, p has 1982 levels, in dataset pr, p should have 1 levels. But I do str(pr), it shows that p still has 1982 levels. also for these pr$aa [1] ARG THR ASP CYS

Re: [R] subset question

2006-12-13 Thread Richard M. Heiberger
help([.factor) a - factor(letters[1:5]) a a[1:3] a[1:3, drop=T] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented,

Re: [R] Syntax for getting covariance matrix

2006-12-13 Thread Prof Brian Ripley
On Wed, 13 Dec 2006, Elena wrote: I am new at using R and was wondering if someone can give me the exact syntax for getting the covariance matrix. I used logit to generate the equation. Specifically the syntax was: glm(formula = Status ~ A + B + C + D + E, family = binomial(link = logit),

[R] xyplot: discrete points + continuous curve per panel

2006-12-13 Thread RMan54
I have a number of x, y observations (Time, Conc) for a number of Subjects (with subject number Subj) and Doses. I can plot the individual points with xyplot fine: xyplot(Conc ~ Time | Subj, Groups=Dose, data=myData, panel = function(x,y) {

[R] digital length

2006-12-13 Thread XinMeng
Hell sir: How can I control the digital length of data? e.g: 0.1234 is the output of an algorithm. What I want is 0.12 instead. How can I get this ? Thanks a lot! My best __ R-help@stat.math.ethz.ch mailing list