[R] A complicated 'aggregate'

Hi, I have a financial (zoo) time series with prices and volumes (although I can get the coredata as a matrix). Due to the data-source some indices have multiple observations. I want to aggregate these according to a weighted average. 11:00:0134 1000 11:00:0135 500

Re: [R] A simple question about summary.glm

Thanks Uwe for your question, Yesterday I found out what I wanted to know: In the summary table the values refer to the different of each treatment to the first treatment alphabetical (in my case it was the control group and therefore served as Tukey test, which is what I wanted in the first

[R] Goodman Kruskal's tau

Hi I need to know which package in R calculates the Goodman Kruskal's tau statistic for nominal data. Also is there any implementation for multiple classification analysis (Andrews at al 1973) in R? Any information on this would be greatly appreciated. Thank you Upasna --

Re: [R] MTBF Reliability calculations

On Tue, 2007-07-31 at 15:10 -0700, Jennings, Eric wrote: I'm working on a project involving reliability values (known failure rates) for a system with approximately 700 components with a set cconfiguration. I'm looking to compute a parts-count MTBF (mean time between failures) for the

Re: [R] array loop

Hi Dong GUO 郭东 [EMAIL PROTECTED] napsal dne 31.07.2007 15:27:35: Thanks, Petr. I changed the equation mark from = to -, then, it works fine. Dont know what difference it has made between the = and -.. from help page The operators - and = assign into the environment in which they are

Re: [R] Nonlinear optimization with constraints

Make a function and add the constraints to the function eg. if (Ai1 x ) Inf. By making the result Infinite, you have added the constraints manually. Vu Nguyen-4 wrote: Hello R community, I am sorry for the previous accidental posting as I pressed a wrong key. I am using R for

Re: [R] deriv for psigamma

UweL == Uwe Ligges [EMAIL PROTECTED] on Tue, 31 Jul 2007 12:13:45 +0200 writes: UweL francogrex wrote: Hi, 2 questions: [] Question 2: deriv(~gamma(x),x) expression({ .expr1 - gamma(x) .value - .expr1 .grad - array(0,

[R] RWeka cross-validation and Weka_control Parametrization

Hello, I have two questions concerning the RWeka package: 1.) First question: How can one perform a cross validation, -say 10fold- for a given data set and given model ? 2.) Second question What is the correct syntax for the parametrization of e.g.

[R] Custom axis

Dear R users, I would like to draw a plot with a custom scale for the axis. More precisely, instead of plotting y on x, I want to plot y on a monotone function of x (for instance a*x+b). Which command and/or package should I use in order to get this result? Thanks Florent Bresson

[R] need help with pdf-plot

Hello, I'm trying to plot a set of barplots like a matrix (2 rows, 10 columns fromreduced_mat) to a pdf. It works with the following parameters: pdf(test.pdf,width=ncol(reduced_mat)*2, height=nrow(reduced_mat)*2, pointsize = 12) par(mfcol = c(nrow(reduced_mat),ncol(reduced_mat)), oma =

Re: [R] Nonlinear optimization with constraints

If you put in penalties for breaking constraints, it is generally better to make the penalty depend on the size of the violation. This keeps the function continuous (though usually not differentiable at the boundary), and gives the optimizer a hint about which way to go. Patrick Burns [EMAIL

Re: [R] Custom axis

x - 1:10 y - rnorm(10,10,1) x2 - 3*x + 2 plot(y ~ x, xaxt = n) axis(side=1,at = x, labels = x2) Joris Florent Bresson [EMAIL PROTECTED]

[R] Re : Custom axis

Maybe I do not explain well what I would like to do. I do not want to change the labels of the axis, but the scale. What I want is a general procedure for changing the scale. Its like using a logarithmic scale on a plot. Labels are the same, but the increases of x along the x-axis are defined

[R] Re : Custom axis

What about the other way round? x - 1:10 y - rnorm(10,10,1) x2 - 3*x + 2 plot(y ~ x2, xaxt = n) axis(side=1,at = x2, labels = x) Florent Bresson

[R] Cut marks on a plot's y-axis to indica te it is a truncated axis

Hi, I've plotted a Kaplan-Meier curve but the curves only range from 0.7 to 1 on the y-axis. Therefore I have used: - ylim=c(0.7,1) [although I think convention dictates that you plot 0.5 to 1 to show the median? A few papers I've read have done this]?? BUT, I would like a symbol like // (but

[R] Cut marks on a plot's y-axis to indica te it is a truncated axis

Hi, I've plotted a Kaplan-Meier curve but the curves only range from 0.7 to 1 on the y-axis. Therefore I have used: - ylim=c(0.7,1) [although I think convention dictates that you plot 0.5 to 1 to show the median? A few papers I've read have done this]?? BUT, I would like a symbol like // (but

Re: [R] array loop

Thanks again, Petr. Following the reference, that would be true that = only assign values to the top level...So apparently using '-' is the safe all the time to assign values. Dong On 8/1/07, Petr PIKAL [EMAIL PROTECTED] wrote: Hi Dong GUO ¹ù¶« [EMAIL PROTECTED] napsal dne 31.07.2007

[R] clear workspace

Dear all, How can I clear the workspace, as we do in Matlab clear all?? Many thanks in advance. Dong [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] Simple table with frequency variable

Hallo, Im trying to find out how to tabulate frequencies of factors when the data have a frequency variable. e,g: i-rep(1:5,2) j-rep(1:2,5) N-10*i+j table(i,j) gives a table of ones as each combination occurs only once. How does one get a table with the corresponding N's? Thanks! Gerrit. --

Re: [R] Cut marks on a plot's y-axis to indicate it is a truncated axis

David Lloyd wrote: Hi, I've plotted a Kaplan-Meier curve but the curves only range from 0.7 to 1 on the y-axis. Therefore I have used: - ylim=c(0.7,1) [although I think convention dictates that you plot 0.5 to 1 to show the median? A few papers I've read have done this]?? BUT, I

Re: [R] A complicated 'aggregate'

Use wtd.mean from Hmisc and by: Lines - 11:00:0134 1000 11:00:0135 500 11:00:0135 1000 11:00:0234 500 11:00:0235 500 library(Hmisc) # for wtd.mean library(zoo) library(chron) # replace with DF - read.table(mytable.dat) DF -

Re: [R] clear workspace

On Wed, 2007-08-01 at 14:06 +0200, Dong GUO 郭东 wrote: Dear all, How can I clear the workspace, as we do in Matlab clear all?? Many thanks in advance. Dong ?rm E.g.: rm(list = ls()) will remove everything shown by ls(). Look at ?ls to see possible arguments to that function to fine

Re: [R] Cut marks on a plot's y-axis to indicate it is a truncated axis

--- David Lloyd [EMAIL PROTECTED] wrote: Hi, I've plotted a Kaplan-Meier curve but the curves only range from 0.7 to 1 on the y-axis. Therefore I have used: - ylim=c(0.7,1) [although I think convention dictates that you plot 0.5 to 1 to show the median? A few papers I've read have

Re: [R] plot matrix data- lattice?

Try this: mat - matrix(1:24, 6, dimnames = list(year = 2001:2006, region = letters[1:4])) library(lattice) xyplot(Freq ~ year | region, as.data.frame.table(mat)) On 8/1/07, Dong Guo [EMAIL PROTECTED] wrote: Dear all, I have a matrix, dim = (years, regions) I would like to plot the data

[R] plot matrix data- lattice?

Dear all, I have a matrix, dim = (years, regions) I would like to plot the data in a lattice so that each panel is region's plot with y-axis based on values, x-axis based on year. how can I do that? Many thanks in advance. Dong [[alternative HTML version deleted]]

Re: [R] plot matrix data- lattice?

?matrix On 8/1/07, Dong Guo [EMAIL PROTECTED] wrote: Thanks, Gabor. My matrix is from a big array(year, regions,variables). so, matrix does not have row names or col names, how could i add the col names or row names?? Thanks again. Dong On 8/1/07, Gabor Grothendieck [EMAIL PROTECTED]

Re: [R] plot matrix data- lattice?

Thanks, Gabor. My matrix is from a big array(year, regions,variables). so, matrix does not have row names or col names, how could i add the col names or row names?? Thanks again. Dong On 8/1/07, Gabor Grothendieck [EMAIL PROTECTED] wrote: Try this: mat - matrix(1:24, 6, dimnames = list(year

[R] Reading Matrices

Dear all, I have been successful so far in plotting matrices and getting the regression line. But, as the matrices contains values like 1 and 0 (diagonal line) which is actually not needed as they are for the same gene, is creating bias in my regression line. I wish to neglect that part of

[R] Two-way ANOVA

I've got this dataframe X12 X14 X17 X19 vitigni years 11 4 78 54 rie 2005 21 4 7 4 rie 2005 31 4 75 5 rie 2005 42 5 66 5croa 2005 51 4 4 46croa 2005 62 5 7 6croa 2005 73 2 56 5 rie 2006 83

[R] passing args to R CMD BATCH in win 2000

Hello and sorry to bother. Please help. I searched the archives but could not find out why --args is being ignored on Windows 2000. I try R CMD BATCH --slave 11.R 11.Rout --args 12 and 11.R has x=commandArgs(trail=T) print(x) a=x[length(x)] write.csv(a,file=13.out) q(no) the argument is not

[R] lattice: densityplot: improper length of lim when more than one conditioning variable, scales free, and empty panels

I am using R 2.5.0 and lattice 0.15-5. plot1, plot2, and plot3 all work fine. plot4 gives the error. x-1:12 f-gl(3,4) g-gl(4,3) plot1-xyplot(y~x|f*g) plot2-xyplot(y~x|f*g,scales=list(relation=free)) plot3-densityplot(~x|f*g) plot4-densityplot(~x|f*g,scales=list(relation=free)) Thanks Ben

Re: [R] [R-pkgs] New R package sqldf

On 8/1/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote: Extremely cool and useful I immediately saw uses for it in some of the work I do. So went to look at it. I especially liked the examples with joins! I observed a few small bugs in the documentation: - The help references The sqldf

Re: [R] Two-way ANOVA

I forget to specify that the data in the table are random. and I want to know how it is possible to carry out all the two way ANOVA test in a single run for all the columns. thank you -- Initial Header --- From : [EMAIL PROTECTED] To : Cc : r-help

[R] cross-correlations

Dear R-users, I would like a suggestion. I have two time series covering the period 1000-2000 a.D. and I would like to understand if there are some time significant correlations between them. Samples in the first time series are quite compact because about a sample per year is available.

[R] shadow between two lines in plot()

Dear R users, I used the following code to draw a scatter plot. plot(x,y,type=n) points(x,y,pch=1) And then I used the abline functions to draw two lines. I want to add the shadow between those two lines. abline(h=200) abline(h=300) Any suggestions? Thanks Rebecca

Re: [R] how to sort dataframe levels

Emilio- One possible way to do this is to create a new factor and put the levels in the order you want them in. For example: related.differences$header2 - factor(related.differences$header, levels=c(spontaneous recovery, negative reinforcer, ... etc. ) where you want spontaneous recovery to

Re: [R] [R-pkgs] New R package sqldf

It works after rm(list=ls()) I had a function called 'fn'; if I had paid close attention when I loaded gsubfn, I would have seen the warning. My fault. Thanks again for a most useful package! David L. Reiner Rho Trading Securities, LLC -Original Message- From: Gabor Grothendieck

[R] Warning generated by Panda GateDefender Integra.

08/01/2007 18:46:14 [GMT+0100] For security reasons certain items found in an email with your address as the sender have not been accepted. File name: TRANSCRIPT.SCR Filtered by: Malformed messages Sender: r-help@stat.math.ethz.ch Recipients: [EMAIL PROTECTED] CC:

[R] Splom custom superpanels

I thought one nice addition to a splom figure would be to have the scatterplots in the upper triangle and a color-coordinated correlation matrix on the bottom. So I tried my hand at customizing panel.pairs(), and was rebuffed. Many times. Four hours of fruitless debugging later, I turn to

Re: [R] shadow between two lines in plot()

see ?rect, or, for more general shapes, ?polygon ## EXAMPLES plot(c(0,500),c(0,500),type=n,las=1) rect(par(usr)[1],200,par(usr)[2],300,col=grey90) points(seq(0,500,length=3),seq(0,500,length=3)) plot(c(0,500),c(0,500),type=n,las=1) polygon((par(usr)[1:2])[c(1,1,2,2)],

[R] t-distribution

If I have a calculated t can I get the probability associated with it using an R function by giving it the df and t? I know I can do the whole calculation using t.test() or get the t-distribution using qt(). If t=1.11 and df =9 can I get the probability? Thanks../Murli

Re: [R] t-distribution

for the upper tail: 1-pt(1.11, 9) [1] 0.1478873 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Nair, Murlidharan T Sent: Thursday, 2 August 2007 4:43 AM To: r-help@stat.math.ethz.ch Subject: [R] t-distribution If I have a calculated t can I get the

[R] constrOptim

Hi, I'm having trouble using the constrOptim function to generate the 9-component vector argmin of the function ELfsds: ELfsds - function(pvechat){ LG=0 for(i in 1:9){ LG=LG+log(pvechat[i]) } return(-LG) } with accompanying gradient function:

Re: [R] t-distribution

if you mean the area to the left of the 1.11 point on the x axis of a t dist with 9 degrees of freedom, Then you need to use pt(1.11,9). See ?pt for more info. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Nair, Murlidharan T Sent: Wednesday, August

Re: [R] t-distribution

Bill.Venables at csiro.au writes: for the upper tail: 1-pt(1.11, 9) [1] 0.1478873 wouldn't pt(1.11, 9, lower.tail=FALSE) be more accurate? __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] t-distribution

Well, is t = 1.11 all that accurate in the first place? :-) In fact, reading beween the lines of the original enquiry, what the person probably wanted was something like ta - pt(-1.11, 9) + pt(1.11, 9, lower.tail = FALSE) which is the two-sided t-test tail area. The teller of the parable will

Re: [R] t-distribution

?pt is what you want. Hope this is helpful, Dan Daniel Nordlund Bothell, WA -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Nair, Murlidharan T Sent: Wednesday, August 01, 2007 11:43 AM To: r-help@stat.math.ethz.ch Subject: [R] t-distribution

Re: [R] t-distribution

see ?polygon function. -- Henrique Dallazuanna Curitiba-Paraná-Brasil 25° 25' 40 S 49° 16' 22 O On 01/08/07, Nair, Murlidharan T [EMAIL PROTECTED] wrote: Indeed, this is what I wanted, I figured it from the function you and Mark pointed me. Thank you both. I am trying to plot it to

Re: [R] Reading Matrices

If want you want to do is to delete the diagonal values in your matrix, here is one way of doing it assuming that your matrix is 'x': # ignore the diagonal value in each column z - lapply(1:ncol(x), function(a) x[-a, a]) # put back into a matrix z - do.call('cbind', z) On 8/1/07, Urmi Trivedi

Re: [R] t-distribution

On 01-Aug-07 19:18:05, [EMAIL PROTECTED] wrote: Well, is t = 1.11 all that accurate in the first place? :-) In fact, reading beween the lines of the original enquiry, what the person probably wanted was something like ta - pt(-1.11, 9) + pt(1.11, 9, lower.tail = FALSE) which is the

[R] a question about d F(x, y; rho) / d rho = f(x, y; rho)

Hello, My name is IKJIN LEE and I am studying in University of Iowa. Recently, I have been trying to prove the fact d F(x, y; rho) / d rho = f(x, y; rho), but I could not. In the mean time, I found that Professor S. Le Cessie from Netherlands kindly gave you an elegant proof of d F(x, y; rho)

Re: [R] t-distribution

Indeed, this is what I wanted, I figured it from the function you and Mark pointed me. Thank you both. I am trying to plot it to illustrate the point and I tried this plot(function(x) dt(x, df = 9), -5, 5, ylim = c(0, 0.5), main=t - Density, yaxs=i) Is there an easy way to shade the area under

Re: [R] lattice grayscale theme

On 7/30/07, Patrick Drechsler [EMAIL PROTECTED] wrote: Deepayan Sarkar [EMAIL PROTECTED] writes: On 7/30/07, Patrick Drechsler [EMAIL PROTECTED] wrote: The Gmane interface seems to have some lag at the moment... Deepayan Sarkar [EMAIL PROTECTED] writes: On 7/28/07, Patrick Drechsler

Re: [R] add custom strip to lattice plot

On 7/30/07, Patrick Drechsler [EMAIL PROTECTED] wrote: Hi, what is the recommended way of adding a strip to a lattice plot? See ?strip.default. In the example below I would like to add the value of mean(y) to a new strip.: --8---cut here---start-8---

Re: [R] Simple table with frequency variable

I am not exactly sure what you are asking for. I am assuming that you want a vector that represent the combinations that are given combinations that are present: N [1] 11 22 31 42 51 12 21 32 41 52 table(i,j) j i 1 2 1 1 1 2 1 1 3 1 1 4 1 1 5 1 1 z - table(i,j) which(z==1)

[R] new user question on dataframe comparisons and plots

I'm coming from the scipy community and have been using R on and for the past week or so. I'm still feeling out the language structure, but so far so good. I apologize in advance if I pose any obvious questions, due to my current lack of diction when searching for my issue, or recognizing it if

[R] Extracting a website text content using R

Dear useR, Just wandering whether it is possible that there is any function in R could let me get the text contents for a certain website. Thanks a lot! Best, Leon [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing

Re: [R] Extracting a website text content using R

Yes, there are. (Please see and follow the posting guide if you wish to obtain something more specific) Bert Gunter Genetech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Am Stat Sent: Wednesday, August 01, 2007 2:19 PM To:

[R] Problem to remove loops in a routine

Dear R-users, I have written the following code to generate some trellis plots. It works perfectly fine except that it is quite slow when it is apply to my typical datasets (over several thousands of lines). I believe the problem comes from the loops I am using to subset my data.frame. I read

[R] Moving data from one workspace to another

Hi, Suppose I have a dataframe in one workspace (a .RData file) dedicated to one project. I then create a new workspace (a new .RData file) for another project but I want to move or copy the dataframe to the new workspace. How can I do this efficiently? Just do a dump and

[R] Predict using SparseM.slm

Hi, I am trying out the SparseM package and had the a question. The following piece of code works fine: ... fit = slm(model, data = trainData, weights = weight) ... But how do I use the fit object to predict the values on say a reserved testDataSet? In the regular lm function I would do

Re: [R] Splom custom superpanels

On 8/1/07, Jonathan Williams [EMAIL PROTECTED] wrote: I thought one nice addition to a splom figure would be to have the scatterplots in the upper triangle and a color-coordinated correlation matrix on the bottom. So I tried my hand at customizing panel.pairs(), and was rebuffed. Many times.

Re: [R] Extracting a website text content using R

All right, my question is, if there is(are) such function(s), what is(are) it(they) ? Best, Leon 2007/8/1, Bert Gunter [EMAIL PROTECTED]: Yes, there are. (Please see and follow the posting guide if you wish to obtain something more specific) Bert Gunter Genetech Nonclinical

Re: [R] Problem to remove loops in a routine

First thing to do is to use Rprof to determine where the time is being spent and then you can pinpoint what section of code is taking the time. A quick look says to do all your subsetting at once. You might look into using 'split' to create the subsets and then access the subsets with [[...]].

Re: [R] Predict using SparseM.slm

If you are feeling altruistic you could write a predict method for slm objects, it wouldn't be much work to adapt what is already available and follow the predict.lm prototype. On the other hand if you are looking for something quick and dirty you can always resort to newX %*%

Re: [R] Extracting a website text content using R

work with it as text. for text mining use: 1- http://wwwpeople.unil.ch/jean-pierre.mueller/ 2- tm by Ingo F. Am Stat wrote: Dear useR, Just wandering whether it is possible that there is any function in R could let me get the text contents for a certain website. Thanks a lot! Best,

Re: [R] Moving data from one workspace to another

Walter, Here is what I do in similar situations. I am on WinXP but this should be similar on other systems (I hope). 1. I start R with the new .RData workspace (usually by double-clicking on it). 2. I go to File Change Dir menu item. 3. I change the directory to where the old .RData is. 4.

Re: [R] Extracting a website text content using R

-Original Message- From: [EMAIL PROTECTED] on behalf of Am Stat Sent: Wed 8/1/2007 2:19 PM To: r-help@stat.math.ethz.ch Subject: [R] Extracting a website text content using R Dear useR, Just wandering whether it is possible that there is any function in R could let me get the text

Re: [R] Extracting a website text content using R

Perhaps more fun is library(XML) res = htmlTreeParse(http://www.omegahat.org/RSXML/;, useInternalNodes=TRUE) xpathApply(res, //h1, xmlValue) [[1]] [1] An XML package for the S language Martin Quoting Steven McKinney [EMAIL PROTECTED]: -Original Message- From: [EMAIL PROTECTED]

[R] how to plot a differential equation?

Hi, I would like to plot the following equation: dF(x)/dx=(k1+k2F(x))(1-F(x)) where k1 and k2 are parameters that I have estimated already. How can I plot the curve in R? Thanks! Montserrat Rue Universitat de Lleida (Spain) [[alternative HTML version deleted]]

Re: [R] how to plot a differential equation?

Hi Montserrat, What exactly would you like to plot? Your differential equation can be easily integrated so that you can get an implicit expression for F(x), i.e. expression like G(c,x,F(x)) = 0 where G is a known function and c is an arbitrary constant. For every value of c and each value of x

Re: [R] new user question on dataframe comparisons and plots

Hi Conor, I hope I interpreted your question correctly. I think for the first one you are looking for a conditioning plot? I am going to create and use some nonsensical data - 'iris' comes with R so this should be reproducible on your machine: library(lattice) data(iris) x - iris # make some

[R] Using 'diff' on zoo vs zooreg classes (possible bug?)

Hello, Can anyone explain the following behaviour? To me it seems a bug, but maybe it is intentional. It seems that a diff on a zooreg class that is not _strictly_ regular only considers those entries that are 'deltat' apart. In the following, diff on the zooreg class only returns values where