[R] Re: locator() via tcltk

2005-06-03 Thread Sebastian Luque
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

[R] locator()

2005-06-03 Thread Erin Hodgess
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

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

2005-06-03 Thread ronggui
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)

Re: [R] locator() via tcltk

2005-06-03 Thread Duncan Murdoch
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

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

2005-06-03 Thread M. K.
--- 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

[R] locator() via tcltk

2005-06-03 Thread Sebastian Luque
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 ide

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

2005-06-03 Thread Achim Zeileis
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 a

Re: [R] factor vector manipulation

2005-06-03 Thread Gabor Grothendieck
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 a

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

2005-06-03 Thread Ross Boylan
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 admoni

Re: RE: [R] rearrange data

2005-06-03 Thread jose silva
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.

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

2005-06-03 Thread Ted Harding
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.val

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

2005-06-03 Thread Deepayan Sarkar
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

RE: [R] rearrange data

2005-06-03 Thread Liaw, 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 E2 2.4 D2 E2 3.1 A3 E3 3.2 B3 E

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

2005-06-03 Thread Wojtek Slusarski
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

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

2005-06-03 Thread M. K.
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. __

[R] rearrange data

2005-06-03 Thread jose silva
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,

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

2005-06-03 Thread Ted Harding
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 =

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

2005-06-03 Thread Kittler, Richard
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 15:

Re: [R] R IRC channel?

2005-06-03 Thread Dan Bolser
/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 s

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

2005-06-03 Thread Deepayan Sarkar
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) xypl

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

2005-06-03 Thread Yi-Xiong Zhou
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] geometric mean regression

2005-06-03 Thread Michael Grant
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 hi

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

2005-06-03 Thread Kittler, Richard
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 __ R-help@stat.m

[R] ts.intersect a multivariate and univariate ts

2005-06-03 Thread Andy Bunn
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

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

2005-06-03 Thread Don MacQueen
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 = c(0,10,250

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

2005-06-03 Thread Peter Dalgaard
"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.

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

2005-06-03 Thread Yi-Xiong Zhou
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 https://sta

[R] factor vector manipulation

2005-06-03 Thread Weiwei Shi
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 c

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

2005-06-03 Thread Romain Francois
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 : http://addicted

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

2005-06-03 Thread Jim Brennan
Something like this!?!? R>new 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 R>barplot(table(new[,2],new[,1]),beside=T,legend.text=c("a","b","c")) -Original Message- From: [EMAIL PROT

Re: [R] Errors on bar plots

2005-06-03 Thread Gabor Grothendieck
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 e

[R] noise poser spectral density

2005-06-03 Thread DIM PATRIDGE
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 no

[R] Errors on bar plots

2005-06-03 Thread Fikret Isik
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.

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

2005-06-03 Thread Jim Brennan
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 sample,

Re: [R] R IRC channel?

2005-06-03 Thread Duncan Murdoch
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 ri

[R] R IRC channel?

2005-06-03 Thread Robert Citek
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. R

Re: [R] geometric mean regression

2005-06-03 Thread Kjetil Brinchmann Halvorsen
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 regress

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

2005-06-03 Thread Liaw, Andy
> From: [EMAIL PROTECTED] > > Dear All, > I posted the following message to the Freebsd-questions mailing > = > > On my laptop I've being using linux for some years now and > "landed" at last > to > the gentoo distribution which I

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

2005-06-03 Thread Peter Dalgaard
"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 > s

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

2005-06-03 Thread Berton Gunter
?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 PRO

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

2005-06-03 Thread Allen S. Rout
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] help with guiDlgOpen

2005-06-03 Thread Duncan Murdoch
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 =

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

2005-06-03 Thread Spencer Graves
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

Re: [R] Hclust question

2005-06-03 Thread Friedrich . Leisch
> 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

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

2005-06-03 Thread Thomas Lumley
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 a

[R] Hclust question

2005-06-03 Thread Mike Schuler
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. A

[R] help with guiDlgOpen

2005-06-03 Thread Marco Zucchelli
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

[R] Histogram of multiple series on one histogram

2005-06-03 Thread Depire Alexandre
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

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

2005-06-03 Thread Miguel A. Arranz
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] Problems of performance between linux and freebsd

2005-06-03 Thread v . demartino2
Dear All, I posted the following message to the Freebsd-questions mailing = 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 software

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

2005-06-03 Thread John Sorkin
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 M

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

2005-06-03 Thread Charles Plessy
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 desig

Re: [R] Need Help - Urgent

2005-06-03 Thread roger bos
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 th

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

2005-06-03 Thread Gavin Simpson
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] Simplify formula for iterative programming

2005-06-03 Thread Stefaan Lhermitte
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] plots from batchfile on windows

2005-06-03 Thread BXC (Bendix Carstensen)
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] plots from batchfile on windows

2005-06-03 Thread Hans-Peter
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 ) win.g

[R] Reading biplot function's source code

2005-06-03 Thread Hector Villalobos Ortiz
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, ...) UseMethod("bi

[R] usign system() both quietly and quickly

2005-06-03 Thread Henrik Andersson
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 l

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

2005-06-03 Thread Sander Oom
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 rep

Re: [R] plots from batchfile on windows

2005-06-03 Thread Gabor Grothendieck
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):

[R] plots from batchfile on windows

2005-06-03 Thread Hans-Peter
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

Re: [R] reading tables into R

2005-06-03 Thread Don MacQueen
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" I

Re: [R] reading tables into R. .

2005-06-03 Thread Henrik Bengtsson
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] RE: GARCH (1 , 1), Hill estimator of alpha, Pareto estimator

2005-06-03 Thread Ukech U. Kidi
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 mu

RE: [R] Need Help - Urgent

2005-06-03 Thread Huntsinger, Reid
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] [m

RE: [R] reading tables into R. .

2005-06-03 Thread Wagle, Mugdha
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 (numer

RE: [R] Creating datasets

2005-06-03 Thread Liaw, Andy
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 befor

Re: [R] reading tables into R

2005-06-03 Thread Sean Davis
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

2005-06-03 Thread Duncan Murdoch
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 firs

Re: [R] Creating datasets

2005-06-03 Thread Douglas Bates
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

Re: [R] Creating datasets *IN PACKAGES*

2005-06-03 Thread Martin Maechler
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 https://stat.ethz.ch/

[R] Creating datasets

2005-06-03 Thread Jim Lemon
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 t

[R] reading tables into R

2005-06-03 Thread Wagle, Mugdha
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

Re: [R] Simplify formula for iterative programming

2005-06-03 Thread Christoph Buser
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,"-

[R] post hoc kruskal wallis

2005-06-03 Thread amberfer
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 th

[R] Need Help - Urgent

2005-06-03 Thread Nivesh Pawar
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 > Z<-as.matrix(expand.grid(x,y,x,y,t,s,t,

[R] Segmentation fault using gentoo linux and mysql

2005-06-03 Thread Dubravko Dolic
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] Simplify formula for iterative programming

2005-06-03 Thread Stefaan Lhermitte
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] geometric mean regression

2005-06-03 Thread Poizot Emmanuel
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 --

RE: [R] post hoc kruskal wallis

2005-06-03 Thread Torsten Hothorn
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, Hollande

Re: [R] dot in formula

2005-06-03 Thread Prof Brian Ripley
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 bein

Re: [R] dot in formula

2005-06-03 Thread Bjørn-Helge Mevik
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

Re: [R] "mvr" function

2005-06-03 Thread Bjørn-Helge Mevik
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 fin

Re: [R] dot in formula

2005-06-03 Thread Peter Dalgaard
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 <-

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

2005-06-03 Thread Torsten Hothorn
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 cl