Re: [R] overlaying plots

You can simply do plot(xx) plot(yy,add=T) Ale iberna - Original Message - From: Mike Jones [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Sent: Thursday, June 02, 2005 9:53 PM Subject: [R] overlaying plots surely this not hard, but i can't figure it out. i'd like to overlay the

Re: [R] mplot :how to deal with missing data

If you have to deal with missing values, you might be interested in na.omit(). ?na.omit Regards, Christoph -- Christoph Buser [EMAIL PROTECTED] Seminar fuer Statistik, LEO C13 ETH (Federal Inst. Technology) 8092 Zurich

[R] [R-pkgs] New CRAN package `coin'

Conditional Inference Procedures in a Permutation Test Framework The `coin' package implements a general framework for conditional inference procedures, commonly known as permutation tests, theoretically derived by Strasser Weber (1999). The conditional expectation and covariance for a broad

Re: [R] dot in formula

Adrian Baddeley [EMAIL PROTECTED] writes: gReetings, I want to manipulate a formula object, containing the name . so that . is replaced by a desired (arbitrary) expression. What is a safe way to do this? I'm not vouching for the safety, but something along the lines of mycall - quote(a+b)

Re: [R] mvr function

McGehee, Robert writes: dataSet - data.frame(y = vol[, 12]) dataSet$X - data.matrix(vol[, 1:11]) ans.pcr - pcr(y ~ X, 6, data = dataSet, validation = CV) If there's a more elegant way of doing this without using data frames of matrices, I'd be interested as well. I actually find using

Re: [R] dot in formula

Adrian Baddeley writes: I want to manipulate a formula object, containing the name . so that . is replaced by a desired (arbitrary) expression. How about myf - y ~ . update(myf, . ~ -. + X) -- HTH, Bjørn-Helge Mevik __ R-help@stat.math.ethz.ch

Re: [R] dot in formula

On Fri, 3 Jun 2005, Bjørn-Helge Mevik wrote: Adrian Baddeley writes: I want to manipulate a formula object, containing the name . so that . is replaced by a desired (arbitrary) expression. How about myf - y ~ . update(myf, . ~ -. + X) or just update(myf, . ~ X). Here . is not being

RE: [R] post hoc kruskal wallis

On Thu, 2 Jun 2005, Liaw, Andy wrote: By post hoc I guess you meant pairwise comparisons. You might want to check out the npmc package on CRAN. ?oneway_test in `coin' has an example, essentially ### Length of YOY Gizzard Shad from Kokosing Lake, Ohio, ### sampled in Summer 1984, Hollander

[R] geometric mean regression

Hi, is it possible to perform a geometric mean regression with R ? Thanks. Emmanuel Poizot Cnam/Intechmer B.P. 324 50103 Cherbourg Cedex Phone (Direct) : (00 33)(0)233887342 Fax : (00 33)(0)233887339

[R] Simplify formula for iterative programming

Dear R-ians, I am looking for the simplification of a formula to improve the calculation speed of my program. Therefore I want to simplify the following formula: H = Si (Sj ( sqrt [ (Ai - Aj)² + (Bi - Bj)² ] ) ) where: A, B = two vectors (with numerical data) of length n sqrt = square root

[R] Segmentation fault using gentoo linux and mysql

Dear Listmembers We have a data server with some amounts of data which we try to analyze using R. Now having some MySQL databases with MySQL Ver 12.22 Distrib 4.0.24 and R 2.1.0 compiled from sources on this machine: Fujitsu Siemens CELSIUS V810 4 GB RAM Dual AMD Opteron Processor 246 Harddisk

[R] Need Help - Urgent

Hello, I am a student and some really urgent help.I am using R software and while creating a grid with eight elements its showing the error and while if i do it with seven it says that the memory is not sufficient... here is the line ...with the error

[R] post hoc kruskal wallis

I´m looking for a program with a post hoc kruskall Wallis test like the Tukey- type non parametric test of ZAR(Biostatistical Analysis 2ºEdition). SSPS hasn´t got any non parametric post host test. I don´t know if R has an appropiate post hoc non parametric test. I´m not sure if NDWD test of

Re: [R] Simplify formula for iterative programming

Dear Stef Please check my code below carefully, because it is Friday afternoon and there might be mistakes in it: ## initialize two random vectors (length n = 100) a - rnorm(100) b - rnorm(100) ## computation of your formula (see ?outer which is very ## useful here) (c - sum(sqrt(outer(a,a,-)^2

[R] reading tables into R

I have been using R and Perl. When I read in a text file using the read.table option, and I try to mathematically manipulate the individual elements in the table, I keep getting an object is not subsettable error. If I try to use a different method, it works, but takes too much time(basically,

[R] Creating datasets

Hi folks, Having lost count of the times I have read the relevant section in R-exts.html (and by the way, where has that useful file gone? I had to look at it on CRAN.), searched for creating datasets (this is the first time that Jon Baron's excellent search site has let me down) and read

Re: [R] Creating datasets *IN PACKAGES*

Hi Jim, are you sure you haven't confused *source* packages with *installed* (aka binary) packages ?? Source packages' ./data/ subdirectories have no funny files in them. Martin __ R-help@stat.math.ethz.ch mailing list

Re: [R] Creating datasets

Jim Lemon wrote: Hi folks, Having lost count of the times I have read the relevant section in R-exts.html (and by the way, where has that useful file gone? I had to look at it on CRAN.), searched for creating datasets (this is the first time that Jon Baron's excellent search site has let me

Re: [R] Creating datasets

On 6/3/2005 6:46 PM, Jim Lemon wrote: Hi folks, Having lost count of the times I have read the relevant section in R-exts.html (and by the way, where has that useful file gone? I had to look at it on CRAN.), Should be in R_HOME/doc/manual. searched for creating datasets (this is the first

Re: [R] reading tables into R

On Jun 3, 2005, at 8:40 AM, Wagle, Mugdha wrote: I have been using R and Perl. When I read in a text file using the read.table option, and I try to mathematically manipulate the individual elements in the table, I keep getting an object is not subsettable error. If I try to use a different

RE: [R] Creating datasets

RSiteSearch(creating datasets, restrict=doc) gives me R-exts.html as the first hit... I did pretty much the same thing you did and had no problem. In the data/ directory of the source package I only have the one .rda file and nothing else. Did you try installing the package and test it before

RE: [R] reading tables into R. .

Hi, The file I am reading is a text file, whose contents are a matrix that has 15 rows and 58 columns. The first row has column names, and the first column has row names, so the format is correct as far as using read.table is concerned. The other values in the table are all float values

RE: [R] Need Help - Urgent

You're asking for a matrix with length(x)^2*length(y)^2*length(s)^2*length(t)^2 rows. What are these lengths? Is that what you're expecting? Perhaps you want the distinct factor levels in x,y,s,t rather than x,y,s,t themselves? Reid Huntsinger -Original Message- From: [EMAIL PROTECTED]

[R] RE: GARCH (1 , 1), Hill estimator of alpha, Pareto estimator

Dear R users, Could you please help me out. I am unable to formulate a model in R to graphs a GARCH (1 , 1) model, the Hill estimator (of alpha), and the Pareto estimator. I have a financial time series data from the DAX and CAC40. I just got introduce to R. I am working on a paper which

Re: [R] reading tables into R. .

Hi, a good start is to verify that what you have read meets your expectations. Try str(data1) dim(data1) print(data1) and make sure you read ?read.table. There are plenty of *arguments* that affects the way R reads your data text file. /Henrik Wagle, Mugdha wrote: Hi, The file I am

Re: [R] reading tables into R

You will have to be more specific, if you want specific help. However, when I get the object is not subsettable error, it is because I have accidentally referred to the wrong object. For example, cespl[1,4] Error in cespl[1, 4] : object is not subsettable mode(cespl) [1] function It

[R] plots from batchfile on windows

Hi, On windows I'd like to run a batchfile that leaves me a plot. As a test I have the scriptfile test.r which only contains: x - 1:10; y - sample( 10 ) plot( x, y ) Now I tried the following (but nothing worked): a) R --vanilla test.r in the cmd window, = that doesn't give me a

Re: [R] plots from batchfile on windows

On 6/3/05, Hans-Peter [EMAIL PROTECTED] wrote: Hi, On windows I'd like to run a batchfile that leaves me a plot. As a test I have the scriptfile test.r which only contains: x - 1:10; y - sample( 10 ) plot( x, y ) Now I tried the following (but nothing worked): a) R

[R] How to 'de-cross' a table?

Dear R users, I have received a table in the following format: id a b c1 c2 d1 d2 1 1 1 65 97 78 98 2 1 2 65 97 42 97 3 2 1 65 68 97 98 4 2 2 65 97 97 98 Factors of the design are: a, b, and e, where e has levels c and d. The levels c and d then have 2

[R] usign system() both quietly and quickly

I try to run an external application (A diffusive transport model programmed in Fortran 95) in a loop using the approach for(i in 1:100) { 1. write parameters[i] to inputfile for application 2. Run application using system(application.exe) 3. Read output written to a file by application to a

[R] Reading biplot function's source code

Hi everybody, Excuse me for this silly question, but after searching the help archives I'm still unable to find my way to read the source code of the biplot function in R. I have installed the mvbutils package, and tried: fixr(biplot) which only gives me: function (x, ...)

Re: [R] plots from batchfile on windows

On 6/3/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: This was discussed previously this week. Look right at the end of: https://www.stat.math.ethz.ch/pipermail/r-help/2005-June/071147.html I had to add a win.graph... command, now it works: x - 1:10; y - sample( 10 )

RE: [R] plots from batchfile on windows

Just open a device before you plot: pdf( plotfile.pdf ) plot( x, y ) dev.off() also have a look at: ?Devices Best Bendix -- Bendix Carstensen Senior Statistician Steno Diabetes Center Niels Steensens Vej 2 DK-2820 Gentofte Denmark tel: +45 44 43 87 38 mob: +45 30 75 87 38

Re: [R] Simplify formula for iterative programming

Dear Christoph, Thanks for your help! I checked it in R and it works if we extend a and b with one element for each run. Unfortunately, I actually want to merge two vectors and then calculate the H for the merge. It is consequently no addition of 1 element but an addition of x elements. I did

Re: [R] Reading biplot function's source code

Hector Villalobos Ortiz wrote: Hi everybody, Excuse me for this silly question, but after searching the help archives I'm still unable to find my way to read the source code of the biplot function in R. I have installed the mvbutils package, and tried: fixr(biplot) which only gives me:

Re: [R] Need Help - Urgent

Those sizes are smaller than I thought they would be, but you repeat the variables more than once in expand.grid, so its still large overall. Maybe you can try a simpler example first and then build up to what you really want once you get the simpler one working. The underlying problem is that

Re: [R] How to 'de-cross' a table?

On Fri, Jun 03, 2005 at 04:33:44PM +0200, Sander Oom wrote : Dear R users, I have received a table in the following format: id a b c1 c2 d1 d2 1 1 1 65 97 78 98 2 1 2 65 97 42 97 3 2 1 65 68 97 98 4 2 2 65 97 97 98 Factors of the design are: a,

[R] Load package using R code rather than GUI Package menu

Is there any way to load a package using R code rather than the Load Package command in the Package menu? I would like to load the boot package. R.2.1.0 Patched Win 2k. Thanks, John John Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics Baltimore VA Medical Center GRECC and University of

[R] Problems of performance between linux and freebsd

Dear All, I posted the following message to the Freebsd-questions mailing = snip On my laptop I've being using linux for some years now and landed at last to the gentoo distribution which I tuned for working with the statistical

Re: [R] Load package using R code rather than GUI Package menu

library(boot) On Friday 03 June 2005 18:00, John Sorkin wrote: Is there any way to load a package using R code rather than the Load Package command in the Package menu? I would like to load the boot package. R.2.1.0 Patched Win 2k. Thanks, John John Sorkin M.D., Ph.D. Chief,

[R] Histogram of multiple series on one histogram

Hello, I have three sample, for example a-c(10,20,10,20,30) b-c(10,20,20,30,30) c-c(20,20,10) I would like to have only one histogram with these series, I try the following code: hist(a) hist(b,add=TRUE,col=red) hist(c,add=TRUE,col=green) but it's

[R] help with guiDlgOpen

Hi, I am trying to open a file with guidlgOpen. on windows. When the file is in a directory whose name contains a space, the name is splitted into 2, i.e. Afile - guiDlgOpen(title= Open file,defaultFile=,defaultDir=, multi=FALSE, filters = c(All files (*.*), *.*)) file [1] C:/Program

[R] Hclust question

Hey, I am running hclust on several different distance matrices and I have a question thats more about labeling. I've been looking for a way to label the edge values on the graph with their distances between them. I've been looking through the documentation and I haven't found anything yet.

Re: [R] How to 'de-cross' a table?

On Fri, 3 Jun 2005, Sander Oom wrote: Dear R users, I have received a table in the following format: id a b c1 c2 d1 d2 1 1 1 65 97 78 98 2 1 2 65 97 42 97 3 2 1 65 68 97 98 4 2 2 65 97 97 98 Factors of the design are: a, b, and e, where e has levels c

Re: [R] Hclust question

On Fri, 03 Jun 2005 12:22:52 -0400, Mike Schuler (MS) wrote: Hey, I am running hclust on several different distance matrices and I have a question thats more about labeling. I've been looking for a way to label the edge values on the graph with their distances between them. I've

Re: [R] RE: GARCH (1 , 1), Hill estimator of alpha, Pareto estimator

First, I doubt if you need library(stats), because it is normally automatically attached when R starts. To confirm, request search(). When I just did this with R 2.1.0 patched, package:stats was the third item. Second, if you've already produced qq-plots and histograms, then I doubt

Re: [R] help with guiDlgOpen

On 6/3/2005 12:17 PM, Marco Zucchelli wrote: Hi, I am trying to open a file with guidlgOpen. on windows. When the file is in a directory whose name contains a space, the name is splitted into 2, i.e. Afile - guiDlgOpen(title= Open file,defaultFile=,defaultDir=, multi=FALSE, filters =

[R] Dirty Rotten Hack. (reversing tickmarks on axes?)

I feel dirty. I have some graphs I'm building to communicate chargeback rates and service usage for our backup system here at the University of Florida. These come down to daily data points on a graph of number-of-bytes transferred and stored. Since we chargeback on the same basis (price per

RE: [R] Dirty Rotten Hack. (reversing tickmarks on axes?)

?par tck Make tck a positive fraction .5 (e.g. .02) in your call. -- Bert Gunter Genentech Non-Clinical Statistics South San Francisco, CA The business of the statistician is to catalyze the scientific learning process. - George E. P. Box -Original Message- From: [EMAIL

Re: [R] Dirty Rotten Hack. (reversing tickmarks on axes?)

Allen S. Rout [EMAIL PROTECTED] writes: I feel dirty. I have some graphs I'm building to communicate chargeback rates and service usage for our backup system here at the University of Florida. These come down to daily data points on a graph of number-of-bytes transferred and stored.

RE: [R] Problems of performance between linux and freebsd

From: [EMAIL PROTECTED] Dear All, I posted the following message to the Freebsd-questions mailing = snip On my laptop I've being using linux for some years now and landed at last to the gentoo distribution which I tuned

Re: [R] geometric mean regression

Poizot Emmanuel wrote: Hi, is it possible to perform a geometric mean regression with R ? Thanks. As has been said on this list before, This is R, there is no if, only how, but if you actually wanted to ask how it is possible, it would help if you explained what is geometric mean

[R] R IRC channel?

Is there an IRC channel for R? I have a bunch of relatively simple questions that I can't find answers to. So it'd be nice to send them to an IRC channel rather than clutter up the mailing list. I'm sure after a few examples and pointers to the right docs I'll be able to find my way.

Re: [R] R IRC channel?

On 6/3/2005 1:20 PM, Robert Citek wrote: Is there an IRC channel for R? I have a bunch of relatively simple questions that I can't find answers to. So it'd be nice to send them to an IRC channel rather than clutter up the mailing list. I'm sure after a few examples and pointers to the

RE: [R] Histogram of multiple series on one histogram

I think you can use barplot for what you want. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Depire Alexandre Sent: June 3, 2005 12:14 PM To: r-help@stat.math.ethz.ch Subject: [R] Histogram of multiple series on one histogram Hello, I have three

[R] Errors on bar plots

Hi R Gurus, I am a new R user. I did my home work by reading the manuals and searching the internet for examples but failed to find an example before writing this message for help. I would like to create bar plots with the standard error of the mean (or 95% confidence intervals) on the bar.

[R] noise poser spectral density

Dear Signal Processing Expert, I would like to generate a random stationary signal of gaussian probability density function to simulate narrow band noise at the output of an IF amplifier. I know the receiver's system temperature (Ts) and IF bandwidth (B) therefore I assume that my narrow band

Re: [R] Errors on bar plots

On 6/3/05, Fikret Isik [EMAIL PROTECTED] wrote: Hi R Gurus, I am a new R user. I did my home work by reading the manuals and searching the internet for examples but failed to find an example before writing this message for help. I would like to create bar plots with the standard error of

RE: [R] Histogram of multiple series on one histogram

Something like this!?!? Rnew a [1,] 10 1 [2,] 20 1 [3,] 10 1 [4,] 20 1 [5,] 30 1 [6,] 10 2 [7,] 20 2 [8,] 20 2 [9,] 30 2 [10,] 30 2 [11,] 20 3 [12,] 20 3 [13,] 10 3 Rbarplot(table(new[,2],new[,1]),beside=T,legend.text=c(a,b,c)) -Original Message- From: [EMAIL PROTECTED]

Re: [R] Histogram of multiple series on one histogram

Hello alexandre, what you are trying to do is *not* an histogram (as a density estimator), if you divide each bar in 3, the surfaces of a won't sum to 1. However a barplot or a barplot2 (in package gplots, bundle gregmisc) would do the trick. See graph 54 on the graph gallery :

[R] factor vector manipulation

Hi, I have one question on factor vector. I have 3 factor vectors: a-factor(c(1, 2, 3)) b-factor(c(a, b, c)) c-factor(c(b, a, c)) what I want is like: c x 1 b 2 2 a 1 3 c 3 which means, I use b as keys and vector a as values and I find values for c. I used the following codes: x-c()

[R] p-value 1 in fisher.test()

The following contingency table generates p-value 1 from fisher.test() ff = c(0,10,250,5000); dim(ff) = c(2,2); fhisher.test(ff)$p.value Sean [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list

Re: [R] p-value 1 in fisher.test()

Yi-Xiong Zhou [EMAIL PROTECTED] writes: The following contingency table generates p-value 1 from fisher.test() ff = c(0,10,250,5000); dim(ff) = c(2,2); fhisher.test(ff)$p.value On this system (PIII, Fedora Core 3): ff = c(0,10,250,5000); dim(ff) = c(2,2); fisher.test(ff)$p.value -1

Re: [R] p-value 1 in fisher.test()

At 10:01 PM +0200 6/3/05, Peter Dalgaard wrote: Yi-Xiong Zhou [EMAIL PROTECTED] writes: The following contingency table generates p-value 1 from fisher.test() ff = c(0,10,250,5000); dim(ff) = c(2,2); fhisher.test(ff)$p.value On this system (PIII, Fedora Core 3): ff =

[R] ts.intersect a multivariate and univariate ts

This seems like a FAQ, but I can't figure it out. I have a mv ts object: R tsp(pg) [1] 1982 20031 R dim(pg) [1] 22 12 and a univariate ts: R tsp(rw) [1] 1690 19961 Yet, when I try to intersect them: R tsp(ts.intersect(rw, pg)) [1] 1982 21761 the process goes awry. How to I

[R] Lattice xyplot -- footnote font size / mtext

Is there a way of controlling the font size and alignment of a footnote in an xyplot, or alternatively of using 'mtext' to place a footnote at the bottom of a graph? --Rich Richard Kittler Advanced Micro Devices, Inc. Sunnyvale, CA __

Re: [R] geometric mean regression

I presume the reference is to the 'geometric mean functional regression' or the 'line of organic correlation' or 'reduced major axis regression'. If so, this is relatively easy alsmost trivial to implement in R. Maybe it's in a package, but I never looked. I worked from Helsel's description in

FW: [R] p-value 1 in fisher.test()

The problem is from the most recent version. Here is the platform: R 2.1.0 2005-04-18, on winXP SP2, Sean -Original Message- From: Douglas Grove [mailto:[EMAIL PROTECTED] Sent: Friday, June 03, 2005 12:49 PM To: Yi-Xiong Zhou Subject: Re: [R] p-value 1 in fisher.test() You need to

Re: [R] Lattice xyplot -- footnote font size / mtext

On Friday 03 June 2005 15:47, Kittler, Richard wrote: Is there a way of controlling the font size and alignment of a footnote in an xyplot, or alternatively of using 'mtext' to place a footnote at the bottom of a graph? Your best bet is using 'page', e.g. library(grid) library(lattice)

Re: [R] R IRC channel?

/server irc.freenode.net /join #R :) On Fri, 3 Jun 2005, Robert Citek wrote: Is there an IRC channel for R? I have a bunch of relatively simple questions that I can't find answers to. So it'd be nice to send them to an IRC channel rather than clutter up the mailing list. I'm sure

RE: [R] Lattice xyplot -- footnote font size / mtext

Great! That should do it. Thanks. --Rich -Original Message- From: Deepayan Sarkar [mailto:[EMAIL PROTECTED] Sent: Friday, June 03, 2005 2:16 PM To: r-help@stat.math.ethz.ch Cc: Kittler, Richard Subject: Re: [R] Lattice xyplot -- footnote font size / mtext On Friday 03 June 2005

Re: [R] p-value 1 in fisher.test()

On 03-Jun-05 Peter Dalgaard wrote: Yi-Xiong Zhou [EMAIL PROTECTED] writes: The following contingency table generates p-value 1 from fisher.test() ff = c(0,10,250,5000); dim(ff) = c(2,2); fhisher.test(ff)$p.value On this system (PIII, Fedora Core 3): ff = c(0,10,250,5000);

[R] rearrange data

Dear all: I have this: A1 B1 C1 D1 E1 A2 B2 C2 D2 E2 A3 B3 C3 D3 E3 And I want this A1 E1 B1 E1 C1 E1 D1 E1 A2 E2 B2 E2 C2 E2 D2 E2 A3 E3 B3 E3 C3 E3 D3 E3 Example: m- matrix(1:15,nrow=3,byrow=T) m v- unlist(list(t(m[,1:4]))) u- rep(c(5,10,15),c(4,4,4)) data.frame(v,u)

[R] lattice: change point color in densityplot()?

When using the densityplot function from lattice library, is there a way to pick a different color *for the points only*? (I'm using densityplot(..., plot.points=TRUE, ...)) Using col parameter changes the color for both, the points and the curve. __

Re: [R] RE: GARCH (1 , 1), Hill estimator of alpha, Pareto estimator]

Ukech U. Kidi wrote: dax- diff(log(DAX_CAC$DAX[1:1865])) m1- garch(dax) Error: couldn't find function garch m1- garch(dax[1:1865]) Error: couldn't find function garch m1- garch(dax[1:1865]) I am sorry, but I forgot to change the addres to r-help in the reply. Well, I am not sure,

RE: [R] rearrange data

reshape() should do it: d V1 V2 V3 V4 V5 1 A1 B1 C1 D1 E1 2 A2 B2 C2 D2 E2 3 A3 B3 C3 D3 E3 d2 - reshape(d, varying=list(names(d)[1:4]), direction=long)[c(3,1)] d2[order(d2$V5),] V1 V5 1.1 A1 E1 1.2 B1 E1 1.3 C1 E1 1.4 D1 E1 2.1 A2 E2 2.2 B2 E2 2.3 C2 E2 2.4 D2 E2 3.1 A3 E3 3.2 B3 E3 3.3

Re: [R] lattice: change point color in densityplot()?

On Friday 03 June 2005 17:49, M. K. wrote: When using the densityplot function from lattice library, is there a way to pick a different color *for the points only*? (I'm using densityplot(..., plot.points=TRUE, ...)) Using col parameter changes the color for both, the points and the curve.

Re: [R] p-value 1 in fisher.test()

On 03-Jun-05 Ted Harding wrote: And on mine (A: PII, Red Had 9, R-1.8.0): ff - c(0,10,250,5000); dim(ff) - c(2,2); 1-fisher.test(ff)$p.value [1] 1.268219e-11 (B: PIII, SuSE 7.2, R-2.1.0beta): ff - c(0,10,250,5000); dim(ff) - c(2,2); 1-fisher.test(ff)$p.value [1]

Re: RE: [R] rearrange data

perfect. thank you Andy! reshape() should do it: d V1 V2 V3 V4 V5 1 A1 B1 C1 D1 E1 2 A2 B2 C2 D2 E2 3 A3 B3 C3 D3 E3 d2 - reshape(d, varying=list(names(d)[1:4]), direction=long)[c(3,1)] d2[order(d2$V5),] V1 V5 1.1 A1 E1 1.2 B1 E1 1.3 C1 E1 1.4 D1 E1 2.1 A2 E2 2.2 B2 E2 2.3 C2

[R] How to change the value of a class slot

I defined an S4 class with a slot i. Then I wrote a regular function that attempted to increment i. This didn't work, apparently because of the general rule that a function can't change the values of its arguments outside the function. I gather there are ways around it, but the Green book

Re: [R] factor vector manipulation

On 6/3/05, Weiwei Shi [EMAIL PROTECTED] wrote: Hi, I have one question on factor vector. I have 3 factor vectors: a-factor(c(1, 2, 3)) b-factor(c(a, b, c)) c-factor(c(b, a, c)) what I want is like: c x 1 b 2 2 a 1 3 c 3 which means, I use b as keys and vector a as values and I

Re: [R] ts.intersect a multivariate and univariate ts

Andy: I have a mv ts object: R tsp(pg) [1] 1982 20031 R dim(pg) [1] 22 12 and a univariate ts: R tsp(rw) [1] 1690 19961 Yet, when I try to intersect them: R tsp(ts.intersect(rw, pg)) [1] 1982 21761 the process goes awry. How to I get rw and pg to be one ts that

[R] locator() via tcltk

Hello, I'm trying to write a function using tcltk to interactively modify a plot and gather locator() data. I've read Peter's articles in Rnews, the help pages in tcltk, http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/, plus a post in R-help sometime ago, but haven't found a solution. The

Re: [R] lattice: change point color in densityplot()?

--- Deepayan Sarkar [EMAIL PROTECTED] wrote: Indirectly. ?panel.densityplot says that there's a parameter for line color, so you could do densityplot(rnorm(100), col = 'red', col.line = 'cyan') Ah, great, thanks. __ R-help@stat.math.ethz.ch

Re: [R] locator() via tcltk

Sebastian Luque wrote: Hello, I'm trying to write a function using tcltk to interactively modify a plot and gather locator() data. I've read Peter's articles in Rnews, the help pages in tcltk, http://bioinf.wehi.edu.au/~wettenhall/RTclTkExamples/, plus a post in R-help sometime ago, but haven't

[R] the test result is quite different,why?

data:http://fmwww.bc.edu/ec-p/data/wooldridge/CRIME4.dta a$call lm(formula = clcrmrte ~ factor(year) + clprbarr + clprbcon + clprbpri + clavgsen + clpolpc, data = cri) bptest(a,st=F) Breusch-Pagan test data: a BP = 34.4936, df = 10, p-value = 0.0001523 bptest(a,st=T)

[R] locator()

Dear Sebastian: Here is a snippet of code that you may find interesting for labeling points on the plane: plane1 - function(x,main1=) { plot(x,main=main1,ylab= ) zz - locator() text(zz$x,zz$y, labels=paste((,round(zz$x,1), ,,round(zz$y,1),),sep=), pos=2,cex=.55)

[R] Re: locator() via tcltk

That's very interesting, thanks a lot Duncan. I should re-read the R Language manual in more detail. Sebastian On 6/3/05, Duncan Murdoch [EMAIL PROTECTED] wrote: Sebastian Luque wrote: Hello, I'm trying to write a function using tcltk to interactively modify a plot and gather locator()