Re: [R] Convert windows package to Linux source package

2008-10-31 Thread Renaud Gaujoux
Indeed in the case of a lazy-loaded package it looks a bit tricky, but it should still be possible, as I guess one can access the source code of any function defined in a package, once loaded (am I wrong?). When no lazy-load takes place then it is actually quite straight forward as you said Du

Re: [R] Problems installing R-2.8.0 on Suse linux 10.1

2008-10-31 Thread Prof Brian Ripley
On Thu, 30 Oct 2008, William Hudspeth wrote: Hello, I am trying to install R on Suse linux 10.1. To get configure to run correctly, I had to set FPICFLAGS=fPIC. FPICFLAGS=-fPIC would have been correct (or -fpic). It is odd that FFLAGS has been set to -g and not '-g -O2': did you change that?

Re: [R] Suppressing internal grid in filled.contour

2008-10-31 Thread Dieter Menne
Martin Maechler stat.math.ethz.ch> writes: > Yes, indeed, PDFs *are* fine Apologies to the original poster, I was not aware of this. Dieter __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting g

[R] Problem of running R console from Java on linux box

2008-10-31 Thread ZhaoWei
Hi the list: I can not create R process from Java application with my linux box, the process never was created, or just terminated immediately after I called Runtime.getRuntime().exec(...) The command lines I have tried are(the paths are all correct): "/bin/sh /usr/bin/R" "/bin/sh -c /u

Re: [R] package pls

2008-10-31 Thread Bjørn-Helge Mevik
tsn4867 <[EMAIL PROTECTED]> writes: >For the package "pls", I need to understand the algorithm for "simpls.fit" > for Partial Least Squares. I'm not sure if simpls.fit tries to find the > weight vectors (loadings) to maximize which of the two: Cov(Xw, y) or > maximize Cov^2(Xw,y)? Are these o

[R] rimage don't install on mac OS X 10.5.5

2008-10-31 Thread Florencio Gonzalez Blazquez
Hello, i'm trying to install de rimage package form source. I have installed fftw2.1.5 and libjpeg, but during compilation appear this message: * Installing *source* package 'rimage' ... checking for g++... g++ checking for C++ compiler default output... a.out checking whether the C++ compil

[R] loglogistic cumulative distribution used by survreg

2008-10-31 Thread mario ouwens
Dear all, What is the cumulative distribution (with parameterization) used within survreg with respect to the log-logistic distribution? That is, how are the parameters linked to the survivor function? Best regards, Mario [[alternative HTML version deleted]] __

[R] how to retrieve some records from a dataframe

2008-10-31 Thread Kurapati, Ravichandra (Ravichandra)
Hi , Dataframe Df contains more than 30 no of records.but I want only first 30 rows of data. How can I retrieve it. Thanks K.Ravichandra [[alternative HTML version deleted]] __ R-help@r-project.org mailing list htt

Re: [R] FW: Re: linux batch question

2008-10-31 Thread Martin Maechler
> "ML" == Mark Leeds <[EMAIL PROTECTED]> > on Wed, 29 Oct 2008 17:13:04 -0500 (CDT) writes: ML> Hi Phil: That's EXACTLY what it is. Thanks so much. It's nice to know ML> that the R Gods don't hate me. I hope it's okay that I'm going to cc ML> r-help ML> in case this t

Re: [R] rimage don't install on mac OS X 10.5.5

2008-10-31 Thread Prof Brian Ripley
On Fri, 31 Oct 2008, Florencio Gonzalez Blazquez wrote: Hello, i'm trying to install de rimage package form source. I have installed fftw2.1.5 and libjpeg, but during compilation appear this message: Where did you install the fftw headers? If in /usr/local/include, this is apparently not in

Re: [R] how to retrieve some records from a dataframe

2008-10-31 Thread nalluri pratap
Df[1:30,] --- On Fri, 31/10/08, Kurapati, Ravichandra (Ravichandra) <[EMAIL PROTECTED]> wrote: From: Kurapati, Ravichandra (Ravichandra) <[EMAIL PROTECTED]> Subject: [R] how to retrieve some records from a dataframe To: r-help@r-project.org Date: Friday, 31 October, 2008, 2:14 PM Hi ,

Re: [R] "Goto" loop

2008-10-31 Thread Peter Dalgaard
megh wrote: > Is there anything like "goto" loop, which exists in most computer programs? > e.g. I am looking for this kind of stuff : > > if(i < 6) goto "step-02" > > Any idea? > > Regards, It doesn't exist, but it can always be replaced by if() {} else {} constructs. (You don't usually see go

[R] Row and Column positions

2008-10-31 Thread Shubha Vishwanath Karanth
Hi R, m=data.frame(a=c(1,NA,5,5),b=c(4,5,6,7),c=c(NA,NA,NA,5)) I want to know the methods of getting row and column positions of NA in the above dataframe. How do I do this? Thanks, Shubha This e-mail may contain confidential and/or privileged i...{{dropped:13}}

Re: [R] Row and Column positions

2008-10-31 Thread Claudia Beleites
Am Freitag 31 Oktober 2008 12:17:30 schrieb Shubha Vishwanath Karanth: > m=data.frame(a=c(1,NA,5,5),b=c(4,5,6,7),c=c(NA,NA,NA,5)) ? which HTH Claudia -- Claudia Beleites Dipartimento dei Materiali e delle Risorse Naturali Università degli Studi di Trieste Via Alfonso Valerio 6/a I-34127 Trieste

[R] *** buffer overflow detected ***: /usr/lib64/R/bin/exec/R terminated on R 2.6.2 to 2.8.0: logging a bug?

2008-10-31 Thread Corrado
Dear friends, there is a serious problem with RODBC on R from 2.6.2 to 2.8.0 on a postgresql 8.3 database. Let's say we define an empty dsn called "test", with user "postgres" and password "none". We test the connection with another application, and it works properly. When you call the databa

[R] "Goto" loop

2008-10-31 Thread megh
Is there anything like "goto" loop, which exists in most computer programs? e.g. I am looking for this kind of stuff : if(i < 6) goto "step-02" Any idea? Regards, -- View this message in context: http://www.nabble.com/%22Goto%22-loop-tp20263445p20263445.html Sent from the R help mailing list

Re: [R] "Goto" loop

2008-10-31 Thread Barry Rowlingson
2008/10/31 megh <[EMAIL PROTECTED]>: > > Is there anything like "goto" loop, which exists in most computer programs? Really? Not since 1968, I think: http://www.cs.utexas.edu/users/EWD/ewd02xx/EWD215.PDF > e.g. I am looking for this kind of stuff : > > if(i < 6) goto "step-02" > > Any idea? Th

Re: [R] how to retrieve some records from a dataframe

2008-10-31 Thread Henrique Dallazuanna
See head: head(Df, n = 30) On Fri, Oct 31, 2008 at 6:44 AM, Kurapati, Ravichandra (Ravichandra) < [EMAIL PROTECTED]> wrote: > Hi , > > > > Dataframe Df contains more than 30 no of records.but I want only > first 30 rows of data. > > > >How can I retrieve it. > > > > Thanks > > K.Ravicha

Re: [R] Row and Column positions

2008-10-31 Thread Henrique Dallazuanna
Try: which(is.na(m), arr = T) On Fri, Oct 31, 2008 at 9:17 AM, Shubha Vishwanath Karanth < [EMAIL PROTECTED]> wrote: > Hi R, > > > > m=data.frame(a=c(1,NA,5,5),b=c(4,5,6,7),c=c(NA,NA,NA,5)) > > > > I want to know the methods of getting row and column positions of NA in > the above dataframe. How

Re: [R] rimage don't install on mac OS X 10.5.5

2008-10-31 Thread Florencio Gonzalez Blazquez
Thanks for your rapid response. I tried your sugerences and set C_INCLUDE_PATH=/usr/local/include. I also look in /usr/local/include and here it appears fftw.h and also fftw3.f, fftw3.h, fftw_threads.h and rfftw.h rfftw_threads.h but it don't work also. I will try to update to R 2.8.0 Thanks

[R] dotplot with a log scale

2008-10-31 Thread glaporta
Dear all, I would like to create a dotplot (lattice) with a log scale as the function http://www.nabble.com/dotchart-with-log-scale--td5411625.html#a5411625 mydotchart of Johannes Huesing. Is it possible? Thanx Gianandrea -- View this message in context: http://www.nabble.com/dotplot-with-a-lo

Re: [R] Fitting weibull and exponential distributions to left censoring data

2008-10-31 Thread Terry Therneau
Use the survreg function. There are many different ways to parameterize a Weibull. The survreg function imbeds it a general location-scale familiy, which is a different parameterization than the rweibull function. > y <- rweibull(1000, shape=2, scale=5) > survreg(Surv(y)~1, dist="weibull"

[R] locfit and latest R version

2008-10-31 Thread Stuart Barker (Prof)
I have been using some old R scripts that were prepared for me using R Version 1.9.0. In these, there is a call to library(locfit) and with this version of R, I have no problem. But needing to upgrade to the latest version of R, I find these scripts no longer work, and with help(locfit) I ge

Re: [R] PAM: how to get the best number of clusters

2008-10-31 Thread Maura E Monville
My problem is that I have already a distance (similarity) matrix generated outside R through a C++ code because the criteria to calculate the "distance" between pairs of objects are none of the standard criteria implemented in R. If I got it right, but I might be mistaken, stepFlexclust() performs

[R] chi-square area analysis

2008-10-31 Thread Arctica
Hi I was wondering if anyone have done an area analysis using chi-aquare? I have three areas (area1, area2 and area3) I need to compare them to each other by looking at a binomial distribution factor (0=negative, 1=positive). I have already done a logistiv regression analysis comparing these are

[R] how to compute a roc curve

2008-10-31 Thread Pau Marc Munoz Torres
Hi, I'm trying to set up a prediction software, now i testing the performance of my method, so i need to calculate a ROC curve, specially auc, cut-off, sens and spec, i just looking at ROCH package, but it's a mass for me, i'm not a math guy and I'm getting lost Could any of you recommend me an

Re: [R] stratified kappa (measure agreement or interrater reliability)?

2008-10-31 Thread David Winsemius
On Oct 30, 2008, at 11:38 PM, FD wrote: Hi All: Could anyone point me to a package that can calculate stratified kappa? My design is like this, 4 raters, 30 types of diagnosis scores, 20 patients. Each rater will rate each patient for each type of diagnosis score. The rater's value is

Re: [R] locfit and latest R version

2008-10-31 Thread David Winsemius
My guess is that you have not configured you new installation to see (and update) your collection of packages. As a quick fix, what happens when you go to your favorite CRAN mirror and install a new copy of locfit? My version is 1.5-4 and runs with recent versions of R. -- David Winsemius,

Re: [R] loglogistic cumulative distribution used by survreg

2008-10-31 Thread Terry Therneau
> What is the cumulative distribution (with parameterization) used within > survreg with respect to the log-logistic distribution? > That is, how are the parameters linked to the survivor function? Reference: Kalbfleisch and Prentice, The statistical analysis of failure time data, chapter 2. Th

[R] subset with variable number of subset arguments

2008-10-31 Thread Brar Piening
I've got a big data.frame from which I need to extract data based on a variable number of Id's (InfCode). Until now I've been using the following dull solution as I never needed to search for more than 5 codes. Now my needs have increased faster than my R skills did and I need to call my functi

Re: [R] subset with variable number of subset arguments

2008-10-31 Thread Erik Iverson
I think subset(set, InfCode %in% x) does it, where 'x' is your vector of ids (e.g., c(1,3,9) in your example at the bottom) Brar Piening wrote: I've got a big data.frame from which I need to extract data based on a variable number of Id's (InfCode). Until now I've been using the following

[R] set the language

2008-10-31 Thread Gustave Lefou
Hello, I would like to change my R language from German to French or English. I have read messages in the archives and also R Installation and Administration paragraph 7.2 but I did not get it all. I have to set LANGUAGE=en (what would it be for French, "fr" ?), but I don't know where. I am unde

[R] getting the p-value from lm as a list object

2008-10-31 Thread eric lee
Hi, I'm trying to get the p-value from the 'lm' regression function as a list object. For example, I can get r^2 from the following code by entering summary(fm)$r.squared. Is there a way to get the p-value? If not, is there a function where I can enter the f-value and degrees of freedom to get

Re: [R] how to compute a roc curve

2008-10-31 Thread Saeed Abu Nimeh
Try library("ROCR") Pau Marc Munoz Torres wrote: > Hi, > > I'm trying to set up a prediction software, now i testing the performance > of my method, so i need to calculate a ROC curve, specially auc, cut-off, > sens and spec, i just looking at ROCH package, but it's a mass for me, i'm > not a m

Re: [R] getting the p-value from lm as a list object

2008-10-31 Thread Thomas Petzoldt
eric lee wrote: Hi, I'm trying to get the p-value from the 'lm' regression function as a list object. For example, I can get r^2 from the following code by entering summary(fm)$r.squared. Is there a way to get the p-value? If not, is there a function where I can enter the f-value and degrees

Re: [R] getting the p-value from lm as a list object

2008-10-31 Thread Marc Schwartz
on 10/31/2008 10:59 AM eric lee wrote: > Hi, > > I'm trying to get the p-value from the 'lm' regression function as a list > object. For example, I can get r^2 from the following code by entering > summary(fm)$r.squared. Is there a way to get the p-value? If not, is there > a function where I c

[R] Help needed with Waterfall plot

2008-10-31 Thread Philip Twumasi-Ankrah
Hi friends, I need suggestions/directions on how to producing a waterfall plot for present extend of change in tumour size for a set of respondents in a study.  Example of use of waterfall plot is in the following slides presented at ASCO 2007 by Axel Grothey. Link is http://media.asco.org/pla

Re: [R] set the language

2008-10-31 Thread Prof Brian Ripley
On Fri, 31 Oct 2008, Gustave Lefou wrote: Hello, I would like to change my R language from German to French or English. On Windows, we presume. I have read messages in the archives and also R Installation and Administration paragraph 7.2 but I did not get it all. I have to set LANGUAGE=en (

Re: [R] Help needed with Waterfall plot

2008-10-31 Thread Thomas Petzoldt
Philip Twumasi-Ankrah schrieb: Hi friends, I need suggestions/directions on how to producing a waterfall plot for present extend of change in tumour size for a set of respondents in a study. Example of use of waterfall plot is in the following slides presented at ASCO 2007 by Axel Grothey. Link

[R] Question about betareg using R

2008-10-31 Thread Jordi Garcia
Good morning, I am using R to try to model the proportion of burned area in Portugal. The dependent variable is the proportion (P). The family used is binomial and the epsilon would be binary. They have recommended me to use the Betareg package as it allows to work with proportions. When I use

[R] [ifelse] how to maintain a value from original matrix without probs?

2008-10-31 Thread Diogo André Alagador
Dear all, I have a matrix with positive and negative values. >From this I would like to produce 2 matrices: 1st - retaining positives and putting NA in other positions 2nd - retaining negatives and putting NA in other positions and then apply rowMeans for both. I am trying to use the function

[R] Quantile Regression for Longitudinal Data:error message

2008-10-31 Thread Helen Chen
Quantile Regression for Longitudinal Data. Hi, I am trying to estimate a quantile regression using panel data. I am trying to use the model that is described in Dr. Koenker's article. So I use the code the that is posted in the following link: http://www.econ.uiuc.edu/~roger/research/panel/rq

[R] Kalman Filter

2008-10-31 Thread Sandrine LUNVEN
Hi, I am studying Kalman Filter and it seems to be difficult for me to apply the filter on a simple ARMA. It is easy to construct the state-space model, for instance: dlmModARMA(ar=c(0.4,-0.2),ma=c(0.2,-0.1, sigma2=1) but applying the dlmFilter on it, it doesn't work... I don't know if my problem

Re: [R] "Goto" loop

2008-10-31 Thread Barry Rowlingson
2008/10/31 Peter Dalgaard <[EMAIL PROTECTED]>: > It doesn't exist, but it can always be replaced by if() {} else {} > constructs. (You don't usually see goto in the class of "functional > programming languages" to which R belongs. See also > http://en.wikipedia.org/wiki/Edsger_W._Dijkstra .) Al

Re: [R] [ifelse] how to maintain a value from original matrix without probs?

2008-10-31 Thread Henrik Bengtsson
> A <- matrix(-4:4, ncol=3) > A [,1] [,2] [,3] [1,] -4 -12 [2,] -303 [3,] -214 > Apos <- A; Apos[A <= 0] <- NA; > Apos [,1] [,2] [,3] [1,] NA NA2 [2,] NA NA3 [3,] NA14 > Aneg <- A; Aneg[A >= 0] <- NA > Aneg [,1] [,2] [,3] [1,] -

Re: [R] [ifelse] how to maintain a value from original matrix without probs?

2008-10-31 Thread Erik Iverson
Hello - Diogo André Alagador wrote: Dear all, I have a matrix with positive and negative values. From this I would like to produce 2 matrices: 1st - retaining positives and putting NA in other positions 2nd - retaining negatives and putting NA in other positions and then apply rowMeans for

Re: [R] [ifelse] how to maintain a value from original matrix without probs?

2008-10-31 Thread Greg Snow
The ifelse function sometimes has unwanted side effects. Your problem can easily be done with simple subsetting. Try: > pos <- neg <- A > pos[ pos <= 0 ] <- NA > neg[ neg >= 0 ] <- NA # remove the ='s if you don't want strict positive/negative Hope this helps, -- Gregory (Greg) L. Snow Ph.D.

[R] R-WinEdt unexcpected Asian characters

2008-10-31 Thread Saskia Sandring
Dear list! I have encountered an odd problem with R-WinEdt. The program is mostly running fine and I had no problems with the installation. But several times a day I'm typing in R-WinEdt there appears a menu with Chinese characters where the curser is and I'm not able to type normally. I don’t k

Re: [R] smooth function advice

2008-10-31 Thread tolga . i . uzuner
Thanks to Mark Leeds for an implementation of exponentially weighted moving averages as follows which solves this problem. This also resmooths recursively. ewma<-function(x,lambda = .5, init = (1-lambda)*.raw[good.ind][1],order=1) { .raw <- unclass(coredata(x)) good.ind <- !is.na(.raw) # dete

Re: [R] dotplot with a log scale

2008-10-31 Thread Deepayan Sarkar
On 10/31/08, glaporta <[EMAIL PROTECTED]> wrote: > > Dear all, > I would like to create a dotplot (lattice) with a log scale as the function > http://www.nabble.com/dotchart-with-log-scale--td5411625.html#a5411625 > mydotchart of Johannes Huesing. Is it possible? The usual lattice syntax for

Re: [R] [ifelse] how to maintain a value from original matrix without probs?

2008-10-31 Thread Marc Schwartz
on 10/31/2008 09:59 AM Diogo André Alagador wrote: > Dear all, > > I have a matrix with positive and negative values. >>From this I would like to produce 2 matrices: > 1st - retaining positives and putting NA in other positions > 2nd - retaining negatives and putting NA in other positions > > a

[R] dot colour in dotchart2 function (Hmisc-package)

2008-10-31 Thread Rainer Tischler
Hi, I'm using the dotchart2-function from the Hmisc-package to plot multiple dots on horizontal lines using the add=TRUE option. My question is very simple: How can I change the color of the dots? I couldn't find the right option in the help-pages of the dotchart2-function or the plot.default-fu

Re: [R] smooth function advice

2008-10-31 Thread Bert Gunter
RSiteSearch("ewma",rest="func") would have told you there already exists an ewma implementation in the qcc package. -- Bert Gunter -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of [EMAIL PROTECTED] Sent: Friday, October 31, 2008 10:04 AM To: [EMAIL PR

[R] Digest setting won't "take"

2008-10-31 Thread Farley, Robert
I've been trying to set my email setting to Digest Mode for a while, but it doesn't seem to "take". I log through to the "Unsubscribe or edit options" button on the Web interface to the R-help mailing list page. I scroll down to the "Set Digest Mode" box and change the radio button to "On"

Re: [R] how to compute a roc curve

2008-10-31 Thread Jorge Ivan Velez
Dear Pau, Take a look at lroc in the epicalc package. HTH, Jorge On Fri, Oct 31, 2008 at 10:25 AM, Pau Marc Munoz Torres <[EMAIL PROTECTED]>wrote: > Hi, > > I'm trying to set up a prediction software, now i testing the performance > of my method, so i need to calculate a ROC curve, specially

[R] AIC for quasipoisson link

2008-10-31 Thread Antonio.Gasparrini
Dear fellows, I'm trying to extract the AIC statistic from a GLM model with quasipoisson link. The formula I'm referring to is AIC = -2(maximum loglik) + 2df * phi with phi the overdispersion parameter, as reported in: Peng et al., Model choice in time series studies os air pollution and

[R] help with contrasts for a binomial 3-way GLM

2008-10-31 Thread Susana Zuloaga
Hi I am a new user the R and I am very grateful for all your help but... I have a problem and I can't resolve yet. I am trying to get the contrasts for a binomial 3-way GLM (T= 4 temperature, t= 2 time and c= 2 substrate levels, plus treatment control) in total they are 17 treatments.

Re: [R] Quantile Regression for Longitudinal Data:error message

2008-10-31 Thread roger koenker
If you are going to insist on doing such things you will have to learn to read the documentation. In this case if you do a traceback() you will see that the error is occurring in rq.fit.slm and when you do ?rq.fit.slm you will see that there are several storage sizes that can

[R] reshape bug?

2008-10-31 Thread stephen sefick
Hadley et al., I was using the cast function to reshape some data (aggregate a melted data frame) and I did not put in the fill and for the most part the values that came out were fine, but there were value great than an order of magnitude from the actual value. When I put in the fill argument eve

Re: [R] Compiling R Packages

2008-10-31 Thread Cahill, Joshua
I'm currently using R version 2.7.2. Not the newest, but not even six months old either. Thanks for the idea. -Original Message- From: Erin Hodgess [mailto:[EMAIL PROTECTED] Sent: Thursday, October 30, 2008 4:06 PM To: Cahill, Joshua Cc: r-help@r-project.org Subject: Re: [R] Compiling

Re: [R] R-WinEdt unexcpected Asian characters

2008-10-31 Thread Uwe Ligges
Saskia Sandring wrote: Dear list! I have encountered an odd problem with R-WinEdt. The program is mostly running fine and I had no problems with the installation. But several times a day I'm typing in R-WinEdt there appears a menu with Chinese characters where the curser is and I'm not able

Re: [R] Digest setting won't "take"

2008-10-31 Thread Marianne Promberger
On Friday, 31 October 2008, 10:46 (UTC-0700), Farley, Robert wrote: > I've been trying to set my email setting to Digest Mode for a while, but > it doesn't seem to "take". I think I had such a problem at some point, and I think it was because I didn't have cookies enabled in my browser. m. --

Re: [R] subset with variable number of subset arguments

2008-10-31 Thread Brar Piening
Thank you very much! Brar Erik Iverson schrieb: I think subset(set, InfCode %in% x) does it, where 'x' is your vector of ids (e.g., c(1,3,9) in your example at the bottom) Brar Piening wrote: I've got a big data.frame from which I need to extract data based on a variable number of Id's

[R] replace() error: new columns would leave holes after existing columns

2008-10-31 Thread Kirsten M Barrett
Hello, I have a problem with using replace() to convert a vector of dates from -mm-dd to julian date. For example, I type replace(x,2004-05-14,134) and I receive an error: Error in `[<-.data.frame`(`*tmp*`, list, value = 134) : new columns would leave holes after existi

Re: [R] AIC for quasipoisson link

2008-10-31 Thread Marc Schwartz
on 10/31/2008 01:07 PM [EMAIL PROTECTED] wrote: > Dear fellows, > > I'm trying to extract the AIC statistic from a GLM model with quasipoisson > link. > The formula I'm referring to is > > AIC = -2(maximum loglik) + 2df * phi > > with phi the overdispersion parameter, as reported in: > >

Re: [R] how to compute a roc curve

2008-10-31 Thread Greg Snow
There are various packages that do ROC curves (others have mentioned some, you can search for the others), but ROC curves are not that hard to do by hand in R. Using the classic iris dataset, we can use Sepal.Length as the test/score and species=='virginica' as the true pos/neg that we are tryi

Re: [R] dotplot with a log scale

2008-10-31 Thread glaporta
Perfectly... thanx Deepayan Deepayan Sarkar wrote: > > On 10/31/08, glaporta <[EMAIL PROTECTED]> wrote: >> >> Dear all, >> I would like to create a dotplot (lattice) with a log scale as the >> function >> http://www.nabble.com/dotchart-with-log-scale--td5411625.html#a5411625 >> mydotchart

Re: [R] Formula for Xi

2008-10-31 Thread Johannes Hüsing
Am 30.10.2008 um 16:53 schrieb [EMAIL PROTECTED]: Hallo! "urn:schemas-microsoft-com:office:office" /> [...] [[alternative HTML version deleted]] This is a rare case where I would have liked to see the HTML version in order to see how these tags are interpreted (paragraphs?

[R] Is there a way to vectorize this?

2008-10-31 Thread Nutter, Benjamin
I apologize up front for this being a little long. I hope it's understandable. Please let me know if I need to clarify anything. Several months ago I wrote a series of functions to help me take my R analyses and build custom reports in html files. Each function either builds or modifies a string

Re: [R] AIC for quasipoisson link

2008-10-31 Thread Ben Bolker
Marc Schwartz comcast.net> writes: > > on 10/31/2008 01:07 PM Antonio.Gasparrini lshtm.ac.uk wrote: > > I'm trying to extract the AIC statistic from a GLM model > >with quasipoisson link. > > The formula I'm referring to is > > > > AIC = -2(maximum loglik) + 2df * phi > > > > with phi t

Re: [R] why does sample(x, n) give the same n items in every separate runs?

2008-10-31 Thread Shuhua Zhan
Hello, I'd like to think all respondents for your advises and suggestions. I am sorry for not including the ordered gene expression data frame (largedf) since it is large. The problem is solved now by Greg Down by including the command: rm(.Random.seed) in the code, perhaps I used some R packges

Re: [R] AIC for quasipoisson link

2008-10-31 Thread Bonenfant Christophe
Ben, I think the reference you're searching for is the one below @ARTICLE{Lebreton1992, author = {Lebreton, J.-D. and Burnham, K. P. and Clobert, J. and Anderson, D. R.}, title = {Modeling survival and testing biological hypotheses using marked animals: a unified approach with

Re: [R] Is there a way to vectorize this? [with correction]

2008-10-31 Thread Nutter, Benjamin
** Sorry to repost. I forgot to include a function necessary to make the example work ** I apologize up front for this being a little long. I hope it's understandable. Please let me know if I need to clarify anything. Several months ago I wrote a series of functions to help me take my R analyse

[R] ctree to rules conversion function

2008-10-31 Thread Sharma, Dhruv
Hi, Similar to the R function of converting rpart into rules (see example: http://books.huihoo.org/data-mining-desktop-survival-guide/Convert_Tree. html_ is there any open code to convert ctree (from party) into rule format? This would be very useful. thanks Dhruv Dhruv Sharma

Re: [R] how to compute a roc curve

2008-10-31 Thread Frank E Harrell Jr
Pau Marc Munoz Torres wrote: Hi, I'm trying to set up a prediction software, now i testing the performance of my method, so i need to calculate a ROC curve, specially auc, cut-off, sens and spec, i just looking at ROCH package, but it's a mass for me, i'm not a math guy and I'm getting lost

Re: [R] PAM: how to get the best number of clusters

2008-10-31 Thread Dave Roberts
Maura, No, in general it is not possible. Depending on your goodness-of-clustering metric (and there are many besides silhouette width), the results may demonstrate multiple peaks and be severely non-monotonic. Perhaps more problematic, silhouette width is especially difficult. It is lo

Re: [R] Digest setting won't "take"

2008-10-31 Thread Farley, Robert
No Joy. :-( I've tried with both IE (7.0) and Firefox (3.0.3) and the same thing happens. Firefox has cookies enabled, and I set IE to explicitly allow http://www.r-project.org/ cookies. Robert Farley Metro www.Metro.net -Original Message- From: [EMAIL PROTECTED] [mailto:[EM

[R] Why are these different?

2008-10-31 Thread Farley, Robert
I want to import some "target values" to normalize elements of a dataframe. I'm summarizing the existing totals using as.matrix(xtabs(~tripid_nu+lineon, data=SurveyData)) and reading a .csv table of target values. I can't even get them to list out as similar data items. My idea is to do

[R] Estimating R2 value for a unit-slope regression

2008-10-31 Thread J. Sebastian Tello
Dear all, I am in need to estimate the amount of variation explained in a variable after simulations that produce a predictor which is in the same units as the dependent variable (numbers of species). Since the dependent and predictor variables are the same, I would think the most appropriate ana

Re: [R] replace() error: new columns would leave holes after existing columns

2008-10-31 Thread Kellie Wills
Hi Kirsten, The second argument of replace() is actually an index, so it's trying to replace element number 2004-5-14 = 1985 with 134! See ?replace. I would need more of your code to understand your error exactly. I think you would be better off converting your dates to date format (?Date) and

[R] trellis: conditional statements within a panel function

2008-10-31 Thread erwann rogard
hello, what i'm trying to achieve is shown in the comment (# ...) below, which is surely very wrong, but probably intuitively more meaningful that i can provide in words. xyplot(y~x2|x1+which,data=make.groups( "a"=data.a, "b"=data.b, "c"=data.c ), panel=function(...){ panel.abline(a=0,b

Re: [R] Error Running TinnR with R

2008-10-31 Thread Yihui Xie
In the 'etc' directory of your installation path, e.g. C:\Program Files\R\etc Regards, Yihui -- Yihui Xie <[EMAIL PROTECTED]> Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086 Mobile: +86-15810805877 Homepage: http://www.yihui.name School of Statistics, Room 1037, Mingde Main Building, Renmin Univ

Re: [R] trellis: conditional statements within a panel function

2008-10-31 Thread Deepayan Sarkar
On 10/31/08, erwann rogard <[EMAIL PROTECTED]> wrote: > hello, > > what i'm trying to achieve is shown in the comment (# ...) below, which is > surely very wrong, but probably intuitively more meaningful that i can > provide in words. > > xyplot(y~x2|x1+which,data=make.groups( > "a"=data.a, >

[R] Help about the plot funtion

2008-10-31 Thread German Urrus
  Hello, my name is German, I am from Argentina, and we are working at university with R. This work is about sounds and how they affect out ears, we are using R to analice different recordings of sounds and we need to get the intensity of this recordings in dB.  In order to get the inte

Re: [R] Help about the plot funtion

2008-10-31 Thread Gabor Grothendieck
Try: ?locator On Fri, Oct 31, 2008 at 10:06 PM, German Urrus <[EMAIL PROTECTED]> wrote: >  > Hello, my name is German, I am from Argentina, and we are working at > university with R. This work is about sounds and how they affect out ears, > we are using R to analice different recordings o

[R] count data with some conditions

2008-10-31 Thread sandsky
Hi there, I have a data set: a=cbind(5,2,4,7,8,3,4,11,1,20) I want to count # of data, satistfying a[1]http://www.nabble.com/count-data-with-some-conditions-tp20275722p20275722.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r

Re: [R] count data with some conditions

2008-10-31 Thread David Winsemius
Do you want the count of remaining elements which are strictly greater than the first element? > length(which(a[1] < a[2:10])) [1] 4 or perhaps a bit more deviously: > sum( a[1] Hi there, I have a data set: a=cbind(5,2,4,7,8,3,4,11,1,20) I want to count # of data, satistfying a[1]http://

Re: [R] getting the p-value from lm as a list object

2008-10-31 Thread David Winsemius
In your example, you could have noted that summary had an element named coefficients which was a 2 x 4 object, with names "(Intercept)" and "x". Although you could have asked for the values with coefficients[2,4], the use of names makes the intent more clear. snipped from the str(summary(