Re: [R] how to multiply a constant to a matrix?

2006-05-24 Thread Michael
imagine when you have complicated matrix algebra computation using R, you cannot prevent some middle-terms become quadratic and absorbs into one scalar, right? if R cannot intelligently determine this, and you have to manually add drop everywhere, do you think it is reasonable? On 5/23/06,

Re: [R] conditional replacement

2006-05-24 Thread Petr Pikal
Hi if speed is critical you can use (x=30x=60)*x+(x30)*30+(x60)*60 HTH Petr On 23 May 2006 at 16:34, Rogerio Porto wrote: From: Rogerio Porto [EMAIL PROTECTED] To: Sachin J [EMAIL PROTECTED], R-help@stat.math.ethz.ch Date sent: Tue, 23 May

Re: [R] problem with ad.test

2006-05-24 Thread Raymond Wan
Hi, On Tue, 23 May 2006, Bhismadev Chakrabarti wrote: i am a novice and have been trying to use the anderson-darling test on a simple text file with one column of data. ... EQ is the name ( in the top row of the column, imported with read.table( file, header=TRUE) of the column of the data.

Re: [R] normality testing with nortest

2006-05-24 Thread Raymond Wan
Hi all, On Mon, 22 May 2006, Rolf Turner wrote: It is worth noting that if the null hypothesis is true, then the p-value is uniformly distributed on [0,1]. This should be kept in mind when assessing the ``instability'' of p-values. Just wanted

[R] Does R have EGARCH modeling function?

2006-05-24 Thread Michael
I've downloaded fSeries, but looks like it just has an interface to OX(TM) Garch Modeling Software,and that OX(TM) software package is not free. So where can I find an EGARCH function that is truely usable? Thanks a lot! [[alternative HTML version deleted]]

[R] ASC/NZSA 2006 - Detailed Program Available via the Website

2006-05-24 Thread David Scott
Australian Statistical Conference / New Zealand Statistical Association Conference 'Statistical Connections' SKYCITY, Auckland, New Zealand 3 - 6 July 2006 LESS THAN TWO MONTHS TO GO Register Online Now! Register now for ASCNZ06 - an international conference with a full and varied

[R] general Gauss-Newton or support for NSUR: contemporaneously correlated non-linear models

2006-05-24 Thread Thomas Wutzler
Dear r-Help readers, 1) Is there support for NSUR in some R package yet? 2) Is there a general function of applying the Gauss-Newton or Marquard method, in which the function of calculating the partial derivatives can be specified by the user? Contemporaneously correlated non-linear models

[R] Two functions with the same name

2006-05-24 Thread Ana Conesa
__ 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

Re: [R] Two functions with the same name

2006-05-24 Thread Petr Pikal
Function written or loaded last will mask a function written or loaded before. HTH Petr PLEASE do read the posting guide! ^ On 24 May 2006 at 13:47, Ana Conesa wrote: Date sent: Wed, 24 May 2006 13:47:09 +0200 To:

[R] data.frame

2006-05-24 Thread Sara Mouro
Dear all, Does any one knows why should I get the following error message, when trying to do a simple data.frame?? DataF-data.frame(Subject,BiomR,Spp,Capas,Litter,Herbs,LitterD,MaxCanH,DDifS p,DSSp,Slope, CanDens,NearestSp) Erro em data.frame(Subject, BiomR, Spp, Capas, Litter, Herbs,

Re: [R] data.frame

2006-05-24 Thread Petr Pikal
Hi Your arguments has different length and therefore the error message data.frame is an object, which resembles a table from Excel, it has the same number of rows in each column From help page: A data frame, a matrix-like structure whose columns may be of differing types (numeric, logical,

Re: [R] Two functions with the same name

2006-05-24 Thread Petr Pikal
Hi you probably created the second function in .Global Environment e.g. directly in console. see ls() and look for the name of the function. If it is there the only way I know about is to remove the function from environment by rm(name.of.the.function) then the masked function shall by

[R] median of a survfit object

2006-05-24 Thread David Hajage
Hello R users ! Here a survfit object : library(survival) essai - aml[aml$x == Maintained,] calc - survfit(Surv(essai$time, 1 - essai$status)) calc Call: survfit(formula = Surv(essai$time, 1 - essai$status)) n events median 0.95LCL 0.95UCL 11 4 103 28 Inf I

Re: [R] Cross correlation/ bivariate/ mantel

2006-05-24 Thread Marcelino de la Cruz
McClatchie, Sam (PIRSA-SARDI mcclatchie.sam at saugov.sa.gov.au writes: Colleagues I have two spatial datasets (latitude, longitude, fish eggs) and (latitude, longitude, fish larvae) at the same 280 stations (i.e. 280 cases). I want to determine if the 2 datasets are spatially

[R] data transformation

2006-05-24 Thread gabriela escati peñaloza
Hi, I have great problems with my work in R. I look for to model the growth of fish. I have Longitudinal data, a serie of repeated measures for each individual. Using the corresponding packages nlme in R. I treat to fit to the data different growth functions, wich were entered by me. Unfortunately

Re: [R] median of a survfit object

2006-05-24 Thread Heinz Tuechler
see: [R] How to access results of survival analysis Xiaochun Li (06 May 2006) At 15:03 24.05.2006 +0200, David Hajage wrote: Hello R users ! Here a survfit object : library(survival) essai - aml[aml$x == Maintained,] calc - survfit(Surv(essai$time, 1 - essai$status)) calc Call:

Re: [R] median of a survfit object

2006-05-24 Thread David Hajage
It works ! thank you ! 2006/5/24, Heinz Tuechler [EMAIL PROTECTED]: see: [R] How to access results of survival analysis Xiaochun Li (06 May 2006) At 15:03 24.05.2006 +0200, David Hajage wrote: Hello R users ! Here a survfit object : library(survival) essai - aml[aml$x ==

Re: [R] dendrogram plotting problem

2006-05-24 Thread Marcelino de la Cruz
j.joshua thomas researchjj at gmail.com writes: Dear List RGui Version : 2.3.0 User : 1 month I am having the *dendrogram plotting problem * The code i tried: library(cluster) DD-DataSetS01022 # 575 x 2 matrix VC-hclust(dist(DD),ave) *Warning message: NAs introduced by

Re: [R] data.frame

2006-05-24 Thread P Ehlers
Sara, You didn't read your data into R correctly. If your data are really in the form you posted (one long column of mixed data types and lots of blank lines), then I would: 1. remove the blank lines with my text editor and save (say, mydata.txt) 2. scan() the variable names into a vector

Re: [R] glmmADMB and the GPL -- formerly-- How to buy R.

2006-05-24 Thread Marc Schwartz
On Tue May 23 18:01:03 CEST 2006 Dave Fourier wrote: Dear List, Some of you have been following the discussion of the GPL and its inclusion in the glmmADMB package we created for R users. I would like to provide a bit of background and include an email we received from Prof. Ripley so

Re: [R] multiple plots with par mfg

2006-05-24 Thread Yan Wong
On 23 May 2006, at 21:47, Romain Francois wrote: Hi, An other possibility might be to use two devices and use dev.set to go from one to another : Thanks. Actually I did try that, but there are quite a lot of points to plot, and the switching between plots slowed the whole simulation

[R] converting hexadecimal

2006-05-24 Thread Romain Lorrilliere
Hi, do you know, a method to convert an hexadecimal value to the corresponding integer value (decimal) ? thinks for help. Romain -- Lorrillière Romain UMR 8079 Laboratoire Ecologie, Systématique et Evolution Bât. 362 Université Paris-Sud 91405 Orsay cedex France tel : 01 69 15

[R] How to make attributes persist after indexing?

2006-05-24 Thread Heinz Tuechler
Dear All! For descriptive purposes I would like to add attributes to objects. These attributes should be kept, even if by indexing only part of the object is used. I noted that some attributes like levels and class of a factor exist also after indexing, while others, like comment or label vanish.

Re: [R] How to make attributes persist after indexing?

2006-05-24 Thread Gabor Grothendieck
You could create your own child class with its own [ method. [.myfactor - function(x, ...) { attr - attributes(x) x - NextMethod([) attributes(x) - attr x } gx - structure(fx, class = c(myfactor, class(fx))) attributes(gx[1]) On 5/24/06, Heinz Tuechler [EMAIL

Re: [R] converting hexadecimal

2006-05-24 Thread Duncan Murdoch
On 5/24/2006 11:01 AM, Romain Lorrilliere wrote: Hi, do you know, a method to convert an hexadecimal value to the corresponding integer value (decimal) ? as.double (which is called by as.numeric) automatically converts hex to numeric values. You need to make sure the hex strings have a 0x

Re: [R] How to make attributes persist after indexing?

2006-05-24 Thread Marc Schwartz (via MN)
On Wed, 2006-05-24 at 17:20 +0100, Heinz Tuechler wrote: Dear All! For descriptive purposes I would like to add attributes to objects. These attributes should be kept, even if by indexing only part of the object is used. I noted that some attributes like levels and class of a factor exist

[R] multiple destinations in duration (survival) analysis

2006-05-24 Thread Dimitri Szerman
Hi, I'm trying to estimate a (parametric) competing risks model in the context of duration (or survival, if you wish) analysis. That is, instead of studying the transition of subjects to death, I wish to study the transitions to multiple *destinations* (which is different from studying multiple

[R] R Reference Card (especially useful for Newbies)

2006-05-24 Thread Berton Gunter
Newbies (and others!) may find useful the R Reference Card made available by Tom Short and Rpad at http://www.rpad.org/Rpad/Rpad-refcard.pdf or through the Contributed link on CRAN (where some other reference cards are also linked). It categorizes and organizes a bunch of R's basic, most used

[R] glmmADMB and the GPL -- formerly-- How to buy R.

2006-05-24 Thread dave fournier
I promise the list that this will be my last posting regarding this matter. The package glmmADMB is not and never was the property or a product of Otter Research Ltd. It was simply a probram I wrote using ADMB-RE and gave away. I do this all the time to help various people. I don't monitor or

Re: [R] glmmADMB and the GPL -- formerly-- How to buy R.

2006-05-24 Thread Jason Barnhart
- Original Message - From: dave fournier [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Sent: Tuesday, May 23, 2006 9:01 AM Subject: [R] glmmADMB and the GPL -- formerly-- How to buy R. Dear List, Some of you have been following the discussion of the GPL and its inclusion in the

Re: [R] normality testing with nortest

2006-05-24 Thread Spencer Graves
see inline Raymond Wan wrote: Hi all, On Mon, 22 May 2006, Rolf Turner wrote: It is worth noting that if the null hypothesis is true, then the p-value is uniformly distributed on [0,1]. SG: This is true in many but not all cases, and I think is a reasonable

Re: [R] Statistical Power

2006-05-24 Thread Charles C. Berry
On Tue, 23 May 2006, Christopher Brown wrote: How can I compute a power analysis on a multi-factor within-subjects design? If you are capable of installing source packages and if you know what a general linear hypothesis test is, you can: download 'hpower' to your computer

Re: [R] UBTA - Wall Street Journal Agree - Ref. wj277963

2006-05-24 Thread Alma Lyles
Kelsey, P2P Sports Betting Software May Change Casino Sports Betting VANCOUVER, BRITISH COLUMBIA - UBA Technology Inc. (UBTA), has entered into initial negotiations to install its proprietary betting exchange software in traditional land-based casinos. Read the whole story:

[R] Regression line limited by the rage of values

2006-05-24 Thread Andreas Svensson
Hi In R, using plot(x,y) followed by abline(lm(y~x)) produces a graph with a regression line spanning the whole plot . This means that the line extends beyond the swarm of data points to the defined of default plot region. With par(xpd=T) it will span the entire figure region. But how can

Re: [R] Regression line limited by the rage of values

2006-05-24 Thread Marc Schwartz (via MN)
On Wed, 2006-05-24 at 18:51 +0200, Andreas Svensson wrote: Hi In R, using plot(x,y) followed by abline(lm(y~x)) produces a graph with a regression line spanning the whole plot . This means that the line extends beyond the swarm of data points to the defined of default plot region.

Re: [R] Statistical Power

2006-05-24 Thread Rick Bilonick
On Wed, 2006-05-24 at 09:22 -0700, Charles C. Berry wrote: On Tue, 23 May 2006, Christopher Brown wrote: How can I compute a power analysis on a multi-factor within-subjects design? If you are capable of installing source packages and if you know what a general linear hypothesis test

[R] changing font size in plot(effect())

2006-05-24 Thread Emilie Berthiaume
I can't seem to be able to change the font size in an effect display. I've tried the following: par(cex.lab=4) plot(effect (alti,reg8), ylab=detection probability) and plot(effect (alti,reg8), ylab=detection probability, cex=4) but nothing changes. Can anyone help me? thanks. Emilie

[R] FastICA and matching variance of inputs

2006-05-24 Thread Andrew Moroz
I have been playing around with the FastICA algorithm in R, and I have come across the following issue. I'm new to this topic, so I'm hoping someone will be able to shed some insight. We know that the variance of the Independent Components is constrained to be 1 and they are uncorrelated, so the

Re: [R] changing font size in plot(effect())

2006-05-24 Thread Chuck Cleland
Try this: trellis.par.set(par.ylab.text = list(cex = 4)) plot(effect (alti,reg8), ylab=detection probability) ?trellis.par.set Emilie Berthiaume wrote: I can't seem to be able to change the font size in an effect display. I've tried the following: par(cex.lab=4) plot(effect

Re: [R] Statistical Power

2006-05-24 Thread Ioannis Dimakos
On Wed, May 24, 2006 20:22, Rick Bilonick wrote: Downloading, unzipping the package on to FC4 running R 2.3.0 gives an error message about a bad description file when trying to install hpower. R CMD INSTALL hpower Indeed, it does so on both my Ubuntu and Win XP boxes. The problem is that

[R] problem-nlme

2006-05-24 Thread gabriela escati peñaloza
Hi, I have great problems with my work in R. I look for to model the growth of fish. I have Longitudinal data, a serie of repeated measures for each individual. Using the corresponding packages nlme in R. I treat to fit to the data different growth functions, wich were entered by me. Unfortunately

[R] Joining variables

2006-05-24 Thread Guenther, Cameron
Hello, If I have two variables that are factors or characters and I want to create a new variable that is the combination of both what function can I use to accomplish this? Ex. Var1Var2 SA100055113 19851113 And I want NewVar SA10005511319851113 Thanks in

Re: [R] Joining variables

2006-05-24 Thread Marc Schwartz (via MN)
On Wed, 2006-05-24 at 14:45 -0400, Guenther, Cameron wrote: Hello, If I have two variables that are factors or characters and I want to create a new variable that is the combination of both what function can I use to accomplish this? Ex. Var1 Var2 SA100055113

[R] LARS error

2006-05-24 Thread Cal Stats
Hi All, when i run cv.lars(x3,y3) it runs fine. but when i run cv.lars(x3,y3,fraction=seq(0,0.1,100)) I get the following error. Error in apply((y[omit] - fit)^2, 2, mean) : dim(X) must have a positive length Any help/suggestions will be appreciated.

Re: [R] Joining variables

2006-05-24 Thread Guenther, Cameron
Thanks to all who responded. The paste function is what I was looking for. Thanks again. Cameron Guenther, Ph.D. Associate Research Scientist FWC/FWRI, Marine Fisheries Research 100 8th Avenue S.E. St. Petersburg, FL 33701 (727)896-8626 Ext. 4305 [EMAIL PROTECTED] -Original Message-

Re: [R] multiple destinations in duration (survival) analysis

2006-05-24 Thread Greg Snow
I have been working on a similar project and here is how I approached it (though I would be very happy to hear other ideas): Our situation was that we wanted to predict length of stay in the emergency room of the hospital, there are multiple competing places to go from the ER: released to go

Re: [R] Joining variables

2006-05-24 Thread Sundar Dorai-Raj
Guenther, Cameron wrote: Hello, If I have two variables that are factors or characters and I want to create a new variable that is the combination of both what function can I use to accomplish this? Ex. Var1 Var2 SA100055113 19851113 And I want

Re: [R] Joining variables

2006-05-24 Thread Berton Gunter
What does combination of both mean exactly. I can think of two interpretations that have two different answers. If you give a small example (as the posting guide suggests) it would certainly help me provide an answer. -- Bert Gunter Genentech Non-Clinical Statistics South San Francisco, CA

[R] Problem with Statistics::R Save plots

2006-05-24 Thread Filipe Garrett
__ 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

Re: [R] Joining variables

2006-05-24 Thread Dimitrios Rizopoulos
look at ?paste(). Best, Dimitris Dimitris Rizopoulos Ph.D. Student Biostatistical Centre School of Public Health Catholic University of Leuven Address: Kapucijnenvoer 35, Leuven, Belgium Tel: +32/(0)16/336899 Fax: +32/(0)16/337015 Web: http://med.kuleuven.be/biostat/

[R] Decimal places

2006-05-24 Thread Robert Mcfadden
Hello, I would like to force R to print 2 decimal places. I use options(digits=2) but as far as I know it prints significant digits so in certain cases is more then 2. Is there other way to do it? Best, Robert [[alternative HTML version deleted]]

Re: [R] Regression line limited by the range of values

2006-05-24 Thread Andreas Svensson
Thankyou very much Marc for that nifty little script. When I use it on my real dataset though, the lines are fat in the middle and thinner towards the ends. I guess it's because lines draw one fitted line for each x, and if you have hundreds of x, this turns into a line that is thicker that it

Re: [R] Joining variables

2006-05-24 Thread mark.lopresti
Hello, Although I cannot speak for the original individual posting this question but I would like to see an example, if possible, of how to construct a factor of several factors. I think this would be more efficient than paste if such a procedure were possible and, of course, would also

Re: [R] Regression line limited by the range of values

2006-05-24 Thread Marc Schwartz (via MN)
On Wed, 2006-05-24 at 21:53 +0200, Andreas Svensson wrote: Thankyou very much Marc for that nifty little script. When I use it on my real dataset though, the lines are fat in the middle and thinner towards the ends. I guess it's because lines draw one fitted line for each x, and if you

Re: [R] Survey proportions... Can I use population as denominator?

2006-05-24 Thread Thomas Lumley
On Tue, 23 May 2006, David L. Van Brunt, Ph.D. wrote: Just giving the survey package a spin... I'm accustomed to stata, and it seems very similar in many respects. One thing is throwing me, however. I've gotten my data in, and specified the design. Looks like the weighting is right (based

[R] biplot with grouped observations

2006-05-24 Thread Kim Milferstedt
Hi, I'd like to plot the results of PCA in biplots, grouping the observations by using different symbols. With xyplot in lattice I can easily differentiate the groups (see for example DAAG, p 285) but it seems difficult (or impossible?) to plot the vectors for the variables as well as xyplot

[R] legend title in effects plot

2006-05-24 Thread Viktor Tron
Could someone tell how I can change/remove the legend title in a (multiline) effect plot? Thanks V In general, how is one supposed to find out answers to such questions? It is virtually impossible to keep track of what graphical parameters are passed on to which low level function and which

Re: [R] Decimal places

2006-05-24 Thread Prof Brian Ripley
On Wed, 24 May 2006, Robert Mcfadden wrote: Hello, I would like to force R to print 2 decimal places. I use options(digits=2) but as far as I know it prints significant digits so in certain cases is more then 2. Is there other way to do it? ?round, possibly format(rounf(x, 2), nsmall =

Re: [R] problem-nlme

2006-05-24 Thread Andrew Robinson
Gabriela, it may be that the model that you are trying to fit is vey complex. Try something simpler. Note that the default setting for the random effects is the same as the fixed effects. In my experience this is a very difficult model to fit, especially for the Richards and von Bertalanffy

Re: [R] Decimal places

2006-05-24 Thread Robert Mcfadden
Thank you very much. Round(x, digids=2) works perfect. Rob -Original Message- From: Prof Brian Ripley [mailto:[EMAIL PROTECTED] Sent: Wednesday, May 24, 2006 10:47 PM To: Robert Mcfadden Cc: r-help@stat.math.ethz.ch Subject: Re: [R] Decimal places On Wed, 24 May 2006, Robert Mcfadden

Re: [R] Regression line limited by the rage of values

2006-05-24 Thread Greg Snow
For single lines the usual method is to just use the lines or segments function, however I have been thinking that there may be some other uses for clipping within a plot region when adding info (abline, but others as well). So here is a first stab at a function to clip within the region, it

[R] Logistic Regression - Results?

2006-05-24 Thread Wojciech Gryc
Hi, I use SPSS at work and have R installed both at work and on my home machine. I've been trying to do some logistic regressions in R and SPSS, but the results I'm getting are different. I've followed a few R tutorials, and with most of them, I get the following instructions: result - glm(z ~ x

Re: [R] How to make attributes persist after indexing?

2006-05-24 Thread Heinz Tuechler
Thank you for your answer, Gabor. I will see, if I understood it. Heinz At 11:31 24.05.2006 -0400, Gabor Grothendieck wrote: You could create your own child class with its own [ method. [.myfactor - function(x, ...) { attr - attributes(x) x - NextMethod([) attributes(x) -

[R] plot xy data with error bars

2006-05-24 Thread Tibi Codilean
Dear All, I have x,y data with errors associated with both x and y. What would be the easiest way of creating plots with both x and y ? Thanks __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] plot xy data with error bars

2006-05-24 Thread Tibi Codilean
Dear All, I have x,y data with errors associated with both x and y. What would be the easiest way of creating plots with both x and y error bars? Thanks __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] Does R have EGARCH modeling function?

2006-05-24 Thread Spencer Graves
I just got 39 hits for RSiteSearch(egarch). The one that looks to me like it comes the closest to answering your question is http://finzi.psych.upenn.edu/R/Rhelp02a/archive/49245.html;. I don't know the definition of egarch, but I would expect that it would not be too hard

Re: [R] general Gauss-Newton or support for NSUR: contemporaneously correlated non-linear models

2006-05-24 Thread Spencer Graves
Have you looked at nlme and the book Pinheiro and Bates (2000) Mixed-Effects Models in S and S-Plus (Springer)? I'm not familiar with NSUR, but it sound like nlme might handle it very well. hope this helps, Spencer Graves Thomas Wutzler wrote: Dear r-Help

Re: [R] How to make attributes persist after indexing?

2006-05-24 Thread Marc Schwartz
On Thu, 2006-05-25 at 00:43 +0100, Heinz Tuechler wrote: Thank you for your answer, Gabor. I will see, if I understood it. Heinz At 11:31 24.05.2006 -0400, Gabor Grothendieck wrote: You could create your own child class with its own [ method. [.myfactor - function(x, ...) { attr

Re: [R] general Gauss-Newton or support for NSUR: contemporaneously correlated non-linear models

2006-05-24 Thread David Scott
On Wed, 24 May 2006, Spencer Graves wrote: Have you looked at nlme and the book Pinheiro and Bates (2000) Mixed-Effects Models in S and S-Plus (Springer)? I'm not familiar with NSUR, but it sound like nlme might handle it very well. hope this helps, Spencer Graves

[R] Ploting two datasets

2006-05-24 Thread Anderson de Rezende Rocha
Dears, I need your help. I have two sets A and B each one containing n lines and being composed by the atributes x and y. Three questions: 1) How can I plot A and B both on the same graphic? I need to show that the classes A and B are separable. In this way, I need to plot the points

[R] boosting

2006-05-24 Thread stephenc
Hi I am using boosting for a classification and prediction problem. For some reason it is giving me an outcome that doesn't fall between 0 and 1 for the predictions. I have tried type=response but it made no difference. Can anyone see what I am doing wrong? Screen output shown below:

Re: [R] exporting long character vectors to dbf

2006-05-24 Thread Eduardo Leoni
I am trying to use RODBC, but still run into some problems. If I understood correctly, I have to change the default type for char vectors to text (or varchar(3000)). But when I try something like library(RODBC) channel - odbcConnect(supremo, uid=R, case=tolower) setSqlTypeInfo(MySQL,