Re: [R] with for objects

Hi Gavin, Thanks for your suggestion - it wasn't quite what I wanted, but at least it got me thinking in the right direction. I realised I didn't really need automatic access to all the slots/methods of a function, just to some I used commonly. So I wrote a function that generates a list

[R] Table to pdf

Hi all, I am new in R world and I haven't been able to find the answer to my question in the documentation I looked up so far. I hope someone can help. In the R function I am writing I have set the graphical output to be saved into a pdf file. I would also like to include in this

AW: [R] Table to pdf

Hi, See ?Sweave in library(tools). HTH Thomas -Ursprüngliche Nachricht- Von: Ana Conesa [mailto:[EMAIL PROTECTED] Gesendet: Donnerstag, 4. Dezember 2003 09:34 An: [EMAIL PROTECTED] Betreff: [R] Table to pdf Hi all, I am new in R world and I haven't been able to find the

Re: [R] Table to pdf

On Thu, 4 Dec 2003, Ana Conesa wrote: I am new in R world and I haven't been able to find the answer to my question in the documentation I looked up so far. I hope someone can help. In the R function I am writing I have set the graphical output to be saved into a pdf file. I would

[R] Matrix Decomposition

Hi, How can one decompose the matrix M into M= W * t(W) I know of the chol() command but I am looking more for hypersphere decomposition and spectral decomposition. Thanks. Coomaren P. Vencatasawmy - [[alternative HTML version deleted]]

Re: [R] Matrix Decomposition

I think perhaps you are looking for something like symsqrt which is given on page 285 of S Poetry. Patrick Burns Burns Statistics [EMAIL PROTECTED] +44 (0)20 8525 0696 http://www.burns-stat.com (home of S Poetry and A Guide for the Unwilling S User) Coomaren Vencatasawmy wrote: Hi, How can one

[R] extracting p value from GEE

Dear R users, If anyone can tell me how to extract the p values from the output of gee? Many thanks in advance. Yu-Kang _ MSN http://www.msn.com.tw/english/ __ [EMAIL PROTECTED]

Re: [R] Matrix Decomposition

I should add that there is a shar file of functions from S Poetry, so you need not type them in yourself. Pat I think perhaps you are looking for something like symsqrt which is given on page 285 of S Poetry. Patrick Burns Burns Statistics [EMAIL PROTECTED] +44 (0)20 8525 0696

[R] Selecting subsamples

Hi all, I?m working with a dataset with 9 columns and 2000 rows. Each row represents an individual and one of the columns represents the volume of that individual (measured in cubic meters). I?d like to select a sample from this dataset (without considering any probability of the rows) in which

RE: [R] Selecting subsamples

CM, maybe s - which(data.frame$attribute = 100) is a starting point!? regards,christian -Ursprüngliche Nachricht- Von: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Auftrag von [EMAIL PROTECTED] Gesendet: Donnerstag, 4. Dezember 2003 13:18 An: [EMAIL PROTECTED] Betreff: [R] Selecting

Re: [R] Selecting subsamples

Hallo I assume you want equal size samples x-runif(1000) this construction gives you y to be set if sum(y) 5 while(sum(y-sample(x,10))5) y-sample(x,10) Cheers On 4 Dec 2003 at 8:18, [EMAIL PROTECTED] wrote: Hi all, I?m working with a dataset with 9 columns and 2000 rows. Each row

RE: [R] get mean of several rows

1. Using rowMeans() is more efficient for computing row means. 2. Mean of five rows is the same as mean of the five row means (exception: if you have NAs, and set na.rm=TRUE). So you can just do something like: k - ceiling(nrow(x) / 5) ktimes - if (rem - nrow(x) %% 5) c(rep(5, k-1), rem) else

[R] Power-function of the F-Test

Hello, does anybody know the procedure how to write a function to generate the power of the two sample F-Test ? Would be glad if somebody can help me. S.Zank [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list

Re: [R] HMisc describe -- error with dates

On Wed, 3 Dec 2003 11:53:44 -0500 Tanya Murphy [EMAIL PROTECTED] wrote: Thank you Frank and Gabor for the fixes and checking and rechecking! Everything seems to work well with the Hmisc functions tried--upData, describe and summary. To summarize: 1. Add the testDateTime and formatDateTime

[R] regression with limited range response

Dear R experts How can you perform a regression analysis in R when the dependent variable is countiuous but bounded, say between 0 and 100? I would be grateful for pointers to R-functions but also for hints to relavant litterature since I have never worked with this problem before. Thanks

RE: [R] Power-function of the F-Test

If you really only have two samples, power.t.test will do. (An F with ndf=1 and ddf=p is the same as the square of a t with df=p.) Andy From: Sebastian Zank Hello, does anybody know the procedure how to write a function to generate the power of the two sample F-Test ? Would be

Re: [R] Comparing Negative Binomial Regression in Stata and R. Constants differ?

On Thu, 4 Dec 2003, Paul E. Johnson wrote: I looked for examples of count data that might interest the students and found this project about dropout rates in Los Angeles High Schools. It is discussed in the UCLA stats help pages for the Stata users:

Re: [R] regression with limited range response

What are your assumptions about the distribution of the dependent variable and its relationship to independent variables? Have you considered glm with family = quasibinomial? hope this helps. spencer graves Kim Mouridsen wrote: Dear R experts How can you perform a regression

Re: [R] Comparing Negative Binomial Regression in Stata and R. Constants differ?

Without going into details, did you remember that the intercept depends very much on the **contrasts** you use (for your factors)? Regards, Martin __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] get mean of several rows

On Thu, 4 Dec 2003, Jan Wantia wrote: Dear all! After hours of trying around, I gave up: I have a 2-dimensional array, and I know how to split it into its rows and how to get the mean for every row using 'sapply'. But what I want is to calculate the mean over the first n rows, and then

RE: [R] get mean of several rows

The following: c( tapply( x, (row(x)-1)%/%5, mean ) ) gives a vector whose first element is the mean of every element in the rows 1 through 5 inclusive, whose second element is the mean of every element rows 6 through 10 inclusive, etc. --- Date: 4 Dec 2003 14:03:57 +0100 From: Jan

Re: [R] R performance--referred from Bioconductor listserv

Please see below. Michael Benjamin wrote: Hi, all-- I wanted to start a (new) thread on R speed/benchmarking. There is a nice R benchmarking overview at http://www.sciviews.org/other/benchmark.htm, along with a free script so you can see how your machine stacks up. Looks like R is

Re: [R] regression with limited range response

Kim This would depend a lot on how the DV is distributed within the 0-100 range. It's even conceivable that the usual linear model could be used without much problem (e.g. if the mean is 50 and the sd is small, then the bounds of 0 and 100 will exist more in theory than in practice; of course,

RE: [R] Irregular Time series

You could take a look at the irregular time-series (its) package on CRAN. Your series is certainly irregular, and possibly a little er... sparse. Anyway, the following might get you started: require(its) mydates - c(28.8.1962,27.6.1977,19.7.1989,26.6.1995,26.7.1999) data -

Re: [R] Comparing Negative Binomial Regression in Stata and R. Constants differ?

On Thu, 4 Dec 2003, Peter Dalgaard wrote: Prof Brian Ripley [EMAIL PROTECTED] writes: report. In the Stata output, they report alpha, same as 1/theta from the R glm.nb output. Except for minor differences in standard errors, only the intercept estimates markedly differ. What are

[R] Comparing Negative Binomial Regression in Stata and R. Constants differ?

I looked for examples of count data that might interest the students and found this project about dropout rates in Los Angeles High Schools. It is discussed in the UCLA stats help pages for the Stata users: http://www.ats.ucla.edu/stat/stata/library/count.htm and See:

RE: [R] regression with limited range response

... the dependent variable is countiuous but bounded, say between 0 and 100. Do you know why the dependent variable is 'bounded'? ... hints to relavant litterature If the dependent var is a percent John Cornell's book on mixture data might be of interest. Is it truncated in the sense that

[R] Font-style

Hi, how to change the font-style (e.g. Arial or Times) in Plots?? Thank Lars Peters - Lars Peters University of Konstanz Limnological Institute D-78457 Konstanz Germany phone: +49 (0)7531 88-2930 fax: +49 (0)7531 88-3533 e-mail: [EMAIL PROTECTED]

Re: [R] add a point to regression line and cook's distance

One way of implementing some Bayesian techniques is to add data points based on prior knowledge. E.g., see Gelman, Carlin, Stern Rubin, in Bayesian Data Analysis (1997) for how a prior on a regression parameter can be interpreted as an additional data point. (Section 8.9 in my 2000

Re: [R] convert data

Muhammad Subianto [EMAIL PROTECTED] writes: Dear R-helper, I have a data set like: OLDa ALL OLDc OLDa OLDb NEW OLDb OLDa ALL . . . ALL OLDc NEW I want to convert that data as OLDa=1, OLDb=2, OLDc=3, NEW=4 and ALL=5 or the result like: 1 5 3 1 2 4 2 1 5 . . . 5 3

[R] bug in as.POSIXct ?

I think that there is a bug in the as.POSIXct function on Windows. Here is what I get on Win2000, Pentium III machine in R 1.8.1. dd1 - ISOdatetime(2003, 10, 26, 0, 59, 59) dd2 - ISOdatetime(2003, 10, 26, 1, 0, 0) dd2 - dd1 Time difference of 1.000278 hours Now, the 26th of October was the

RE: [R] convert data

Do you just want to out the data vectors together? Try ?scan, ?c, ?data.frame my.data.frame - data.frame( names = c(OLDa, ALL, OLDc, OLDa, OLDb, NEW..), values = c(1,5,3,1,2,4.) HTH, -Andy __ [EMAIL PROTECTED] mailing list

Re: [R] Font-style

That depends on the plotting device: I answered this _yesterday_, so please consult the list archives. Date: Wed, 3 Dec 2003 11:05:20 + (GMT) From: Prof Brian Ripley [EMAIL PROTECTED] To: Savano [EMAIL PROTECTED] Subject: Re: [R] FONT TYPE On Thu, 4 Dec 2003, Lars Peters wrote: how to

RE: [R] performance gap between R 1.7.1 and 1.8.0

From: Martin Maechler [mailto:[EMAIL PROTECTED] [...] A very first step of diagnosis might be to activate trace(read.dcf) trace(library) options(verbose = TRUE) A step further might be to patch read.dcf such that it prints info if(getOption(verbose)) { print what I am

Re: [R] performance gap between R 1.7.1 and 1.8.0

Liaw, Andy [EMAIL PROTECTED] writes: I'm really out of ideas at this point. Can anyone see the problem given the above, or at least tell me what to try next? Hmm. Looking into the actual Rprof.out file is sometimes revealing. Could you provide a link to it (I assume it might be a bit large

Re: [R] bug in as.POSIXct ?

[EMAIL PROTECTED] wrote: I think that there is a bug in the as.POSIXct function on Windows. Here is what I get on Win2000, Pentium III machine in R 1.8.1. dd1 - ISOdatetime(2003, 10, 26, 0, 59, 59) dd2 - ISOdatetime(2003, 10, 26, 1, 0, 0) dd2 - dd1 Time difference of 1.000278 hours Now, the

[R] assigning colors to barplot when beside=TRUE

dear list, i am having trouble coloring the bars in a barplot. my data have two groups, which i would like to plot side by side. within each group i want to sort the observations in decreasing order, like a pareto chart. the bar colors would relfect the value of a third variable. below i have

RE: [R] bug in as.POSIXct ?

Not sure but perhaps this is related to bug PR#3646 reported in: http://maths.newcastle.edu.au/~rking/R/devel/03b/0195.html --- Date: Thu, 4 Dec 2003 11:20:15 -0600 From: [EMAIL PROTECTED] To: [EMAIL PROTECTED] Subject: [R] bug in as.POSIXct ? I think that there is a bug in the

Re: [R] Interrupt handling

The short answer is `no', nothing has changed since 1.8,0. There was a change at 1.8.0, and that was deliberate and described in the NEWS file. You seem rather confused about this! A slighter longer answer is that this depended on your unnamed OS having such interrupts, so for many R users

[R] Code for Hodrick-Prescott Filter

Hi, I saw ur reply for Hofdrick-Prescott filter in a web site. I thought maybe u can get some help me. Can I find a patch for excel to run Hodrick-Prescott filter? I want to automate it for forecast experiments. Thanx for your time and attenton Ayhan Tosuner [[alternative HTML

Re: [R] assigning colors to barplot when beside=TRUE

Hi Alejandro Munoz Del Rio wrote: dear list, i am having trouble coloring the bars in a barplot. my data have two groups, which i would like to plot side by side. within each group i want to sort the observations in decreasing order, like a pareto chart. the bar colors would relfect the value of

RE: [R] volume of an irregular grid

On 03-Dec-03 Karim Elsawy wrote: I have a 3d irregular grid of a surface (closed surface) I would like to calculate the volume enclosed inside this surface can this be done in R any help is very much appreciated best regards karim Hi Karim, You should be able to create your own function

Re: [R] assigning colors to barplot when beside=TRUE

Quoting Paul Murrell [EMAIL PROTECTED]: I haven't figured out exactly what your code is supposed to produce, but I suspect that you are calculating colv incorrectly. In your example, you get a colv with 5 elements. There are a total of 10 bars plotted (some are zero height because the

[R] Processing calendar dates with R

I am a beginner in R with a background in SAS. Are there built-in R methods of reading dates for calculating elapsed days between two calendar dates? If so, are there any examples I can browse? Thanks in anticipation. John Byrne. Lecturer in Information Systems. Australian Catholic

[R] passing par() options to pdf under Sweave

I'm using Sweave to generate LaTeX output from R. Thank-you very much for this capability. It is very useful. Could someone please tell me how to pass non-default graphics parameters from par() to the pdf driver when it is called by Sweave? (I need to change some of the margin defaults. )

[R] simple sphericity test using bartlett and cov?

Hey, Is this right for a sphericity test? x - mymatrix mycov - cov(x) bartlett.text (array (mycov), myconditions) __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

[R] Re: passing par() options to pdf under Sweave

Sorry for wasting the list for this. Had a blind spot. It is very easy. Just include the par() commands inside the R chunks that generate the plots. Anne On Thu, 4 Dec 2003, Anne York wrote: I'm using Sweave to generate LaTeX output from R. Thank-you very much for this capability. It is

Re: [R] Selecting subsamples

[EMAIL PROTECTED] wrote [that he has a data set with 9 variables (columns) measured on 2000 individuals (rows) and wants a sample] in which the sum of the volume of the individuals in that sample = 100 cubic m. Let's suppose that this information is held in d, a data frame, and that

[R] grid packages since R-1.7.1

I'm having a spot of bother using code that worked with R-1.7.1 but will not work with 1.8.1. The beginning of the saga is with grid.polygon ostensibly not findable. One does exist in /R-1.8.1/library/grid/R, and when I specifically load the grid package (which probably isn't a good idea),

Re: [R] get mean of several rows

Jan Wantia [EMAIL PROTECTED] asked: I have a 2-dimensional array, and I know how to split it into its rows and how to get the mean for every row using 'sapply'. But what I want is to calculate the mean over the first n rows, and then the second n rows, etc., so

Re: [R] grid packages since R-1.7.1

On Thursday 04 December 2003 21:19, Patrick Connolly wrote: I'm having a spot of bother using code that worked with R-1.7.1 but will not work with 1.8.1. The beginning of the saga is with grid.polygon ostensibly not findable. One does exist in /R-1.8.1/library/grid/R, and when I