[R] package.skeleton from within function: objects not found

2008-01-28 Thread Tineke Casneuf
Hi all, I ran into a strange error: I am trying to create a package skeleton for a new source package from within a function. Objects that are created in this function are to be included in my package, but for some reason, I get an error message saying that these objects cannot be found. Here is

Re: [R] formula for nls

2008-01-28 Thread Jarosław Jasiewicz
Christian Ritz pisze: > Hi Jarek, > > an alternative approach is to provide more precise starting values! > > It pays off to realise that it's possible to find quite good guesses > for some of the parameters in your model function: > > t ~ tr+(ts-tr)/((1+(a*h)^n)^(1-(1/n))) > > The par

Re: [R] package.skeleton from within function: objects not found

2008-01-28 Thread Prof Brian Ripley
You need to set the 'environment' argument (the help file is incomplete). e.g. env <- sys.frames()[[sys.nframe()]] package.skeleton(name = pkgName, list=c("f","e", "myenv"), env=env) On Mon, 28 Jan 2008, Tineke Casneuf wrote: > Hi all, > > I ran into a strange error: I am trying to create a pac

Re: [R] package.skeleton from within function: objects not found

2008-01-28 Thread Tineke Casneuf
Indeed, this works. Thanks! On Jan 28, 2008 9:30 AM, Prof Brian Ripley <[EMAIL PROTECTED]> wrote: > You need to set the 'environment' argument (the help file is incomplete). > e.g. > > env <- sys.frames()[[sys.nframe()]] > package.skeleton(name = pkgName, list=c("f","e", "myenv"), env=env) > > >

Re: [R] How to fill bar plot with textile rather than color

2008-01-28 Thread Domenico Vistocco
If you type > example(barplot) you will find an example. Ciao, domenico CHENYS wrote: > Hi, I'm looking for a tool which can fill bar chart with dash, skewed line, > or grids, rather than pure color. Any one have the idea how to do that in R? > Or maybe in Matlab will also be helpful. > > Than

Re: [R] using facet_grid() from ggplot2 with additional text in labels

2008-01-28 Thread ONKELINX, Thierry
Dear Rainer, I think you'll have to recode the labels of the levels to get the output that you want. levels(ssq$gi) <- c("gi = 1", "gi = 2", "gi = 3") HTH, Thierry ir. Thierry Onkelinx Instituut voor natuur- en bos

Re: [R] Bug in levels() function?

2008-01-28 Thread Peter Dalgaard
Groot, Philip de wrote: > Hello all, > > I am not sure whether it actually is a bug, but it is not the behaviour I > would expect. Please consider this: > > >> Sibships >> > [1] Patient_2400 Patient_2400 Patient_345 Patient_345 Patient_8901 > [6] Patient_8901 Patient_4008 Patient_4

[R] how to skip last lines while reading the data in R

2008-01-28 Thread mrafi
hey all greetings hey all am an engineering student...and am trying to learn R i am trying to automate reading a specific type of file...and perform certain functions...but i want to omit lines in the end of the file.. there is an option for skiping the lines before begining...but how can i ask R

[R] Bug in levels() function?

2008-01-28 Thread Groot, Philip de
Hello all, I am not sure whether it actually is a bug, but it is not the behaviour I would expect. Please consider this: >Sibships [1] Patient_2400 Patient_2400 Patient_345 Patient_345 Patient_8901 [6] Patient_8901 Patient_4008 Patient_4008 Patient_7991 Patient_7991 [11] Patient_8353 Patie

[R] basic spatial query

2008-01-28 Thread jgarcia
Hello; Please coud you advise me of a simple way to select points (x,y coordinates) that fall within a polygon. I've got a set of polygons, each one defined by an arbitrary number of points, and several points inside each polygon. I know this is simple with a GIS, but I'd rather do it inside R. Tha

[R] [R-pkgs] dynlm: new version 0.2-0

2008-01-28 Thread Achim Zeileis
Dear useRs, I've release a new version of the "dynlm" package to CRAN which adds two new features: o instrumental variables regression (two-stage least squares) via formulas like dynlm(y ~ x1 + x2 | z1 + z2 + z3, data = mydata) where z1, z2, z3 are the instruments which can again

Re: [R] how to skip last lines while reading the data in R

2008-01-28 Thread Henrique Dallazuanna
Perhaps: data <- read.table(textConnection(rev(rev(readLines('data.txt'))[-(1:2)]))) will skip the last two lines '-(1:2)' On 28/01/2008, mrafi <[EMAIL PROTECTED]> wrote: > > hey all > greetings > > hey all am an engineering student...and am trying to learn R > i am trying to automate reading a

Re: [R] basic spatial query

2008-01-28 Thread Roger Bivand
ija.csic.es> writes: > > Hello; > Please coud you advise me of a simple way to select points (x,y > coordinates) that fall within a polygon. > I've got a set of polygons, each one defined by an arbitrary number of > points, and several points inside each polygon. > I know this is simple with a

Re: [R] maptools no such file

2008-01-28 Thread Roger Bivand
gmail.com> writes: > > No, I get the same error message there... Please do not needlesly delete the thread content. Your original question was: "I'm having problems reading a shapefile with read.shape (maptools). I'm absolutely sure my file is there, but I get "no such file". The wd is ok, si

Re: [R] package.skeleton from within function: objects not found

2008-01-28 Thread Tineke Casneuf
Professor Ripley, do you have an idea why it works for me to save different types of objects, but not for the environment object I construct in the function: ### Example of a function that works without error: fun2 <- function(myname){ f <- function(x,y) x+y g <- function(x,y) x-y d <- da

Re: [R] maptools no such file

2008-01-28 Thread Roger Bivand
Roger Bivand nhh.no> writes: > > gmail.com> writes: > > > > > No, I get the same error message there... > > Your original question was: > > "I'm having problems reading a shapefile with read.shape (maptools). I'm > absolutely sure my file is there, but I get "no such file". The wd is ok, >

Re: [R] double-click in RData file versus load( file )

2008-01-28 Thread Cleber Nogueira Borges
hello Duncan, Gabor Many thanks for your help! I think that the line: if( chartr("/","\\",getwd() )==R.home() ) setwd("C:\\") can solve my problem in this moment!! :-) Thanks again Cleber > On 26/01/2008 9:03 PM, Cleber Nogueira Borges wrote: >> hi Gabor and Duncan, >> >> I make a test a

Re: [R] KM estimation for interval censoring?

2008-01-28 Thread David Winsemius
"gallon li" <[EMAIL PROTECTED]> wrote in news:[EMAIL PROTECTED]: > Does anybody know if there is such a function to estimate the > distribution for interval censored data? > > survfit doesn't work for this type of data as I tried various > references. My reading of R-help posts by Therneau and

[R] x-axis

2008-01-28 Thread mohamed nur anisah
Hi, I want to plot a graph and here is my code: ec<-rep(0,length(e)) fc<-rep(0,length(f)) plot(e,ec,type="p",col=1,pch=19) points(f,fc,col=2,pch=20) legend(1.0e+08,1.0,c("dog", "human"),text.col="green4",pch=c(19,20),col=c(1,2)) my major problem here is the x-axis is too lar

Re: [R] using facet_grid() from ggplot2 with additional text in labels

2008-01-28 Thread Rainer M Krug
Hi Hadley On 27/01/2008, hadley wickham <[EMAIL PROTECTED]> wrote: > > > I am using ggplot2 at the moment and I must say it is definitely better > > then ggplot - good work. > > > > My problem is that I am using facet_grid() in the following way: > > > > > p <- ggplot(ssq, aes(x=year, y=-log(ssq)

Re: [R] how to skip last lines while reading the data in R

2008-01-28 Thread Barry Rowlingson
Henrique Dallazuanna wrote: > Perhaps: > > data <- read.table(textConnection(rev(rev(readLines('data.txt'))[-(1:2)]))) > Euurgh! Am I the only one whose sense of aesthetics is enraged by this? To get rid of the last two items you reverse the vector, remove the first two items, then reverse t

Re: [R] Bug in levels() function?

2008-01-28 Thread hadley wickham
> And subsetting a factor retains the original factor levels. To drop > unused levels, just use factor(f[index]) or f[index, drop=TRUE]. The > opposite behaviour can be even more annoying/dangerous because it leads > to empty cells dropping out of tables and bars disappearing from barplots. Of cou

Re: [R] using facet_grid() from ggplot2 with additional text in labels

2008-01-28 Thread hadley wickham
> > ggopt(strip.text = function(variable, value) paste(variable, value, sep="= > ")) > > That's exactly what I was looking for - thanks. One thing that I should mention is that this is likely to change at some point in the future. Eventually it will become: + facet_grid(strip.text = function(var

[R] How to extract Hopach MSS value

2008-01-28 Thread Prateek Singh
HI, I was running clustering through Hopach method. I used hopach library and got the output. But the output has lots other things which we doesn't need. I want only MSS(Median Split Silhouette) value. for my further work. please suggest me way how do I extract only MSS value. Or you can say ho

[R] help with checking out-of-range values in each column in data frame

2008-01-28 Thread Tom Cohen
Dear list, I have following data, where I want to check if any value in each column is out of range. For example, column "f1" can only take values 1-5, so if any values less than 1 or > 5 will be defined as missing value (i.e. NA), column "f4" can only take values of 1-3 and any values th

Re: [R] how to skip last lines while reading the data in R

2008-01-28 Thread Prof Brian Ripley
On Mon, 28 Jan 2008, Barry Rowlingson wrote: > Henrique Dallazuanna wrote: >> Perhaps: >> >> data <- read.table(textConnection(rev(rev(readLines('data.txt'))[-(1:2)]))) >> > > Euurgh! Am I the only one whose sense of aesthetics is enraged by > this? To get rid of the last two items you reverse th

[R] Grouping data via an index

2008-01-28 Thread Tobin, Jared
Hello r-help, I have a lengthy vector of data (with values anywhere from 1-200), and another index vector of 'groups' representing values 0-2, 3-5, 6-8, ... of length 67. The index vector has the structure (1, 4, 7, ... , 196, 199), where each value is the midpoint of each respective group. I'm

[R] R on an eeePC

2008-01-28 Thread Dr. Walter H. Schreiber
Dear list, I wonder if somebody has succeeded in installing R on an eeePC (Xandros desktop). Searching via Rseek (term eeePC) and in eeePC forums (term Cran) left me without proper hits. Best wishes, Walter. __ R-help@r-project.org mailing list

[R] R loops

2008-01-28 Thread cvandy
I'm a new user and am having trouble with loops. In the following, I'm trying to add the results of "test" and the loops are not working. I've simplified the loop. What am I doing wrong? Thanks! > test<-numeric(20) > tot<-numeric(20) > for(i in 1:20){test[i]<-1} > for (i in 1:20){tot[i]<-(test[i]

Re: [R] help with checking out-of-range values in each column in data frame

2008-01-28 Thread Gabor Grothendieck
Try: transform(data, f1 = factor(f1, levels = 1:5), f4 = factor(f4, 1:3)) On Jan 28, 2008 8:12 AM, Tom Cohen <[EMAIL PROTECTED]> wrote: > Dear list, > > I have following data, where I want to check if any value in each > column is out of range. For example, column "f1" can only take values >

Re: [R] how to skip last lines while reading the data in R

2008-01-28 Thread Gabor Grothendieck
If you don't mind reading it in twice its just: DF <- read.table("xy.dat", header = TRUE, nrow = length(readLines("xy.dat")) - 3) tail(DF) # or DF <- read.table("xy.dat", header = TRUE, nrow = length(count.fields("xy.dat")) - 3) tail(DF) # or DF <- read.table("xy.dat", header = TRUE, nrow = nrow(r

[R] [R-pkgs] New package: gRain - gRAphical Independence Networks

2008-01-28 Thread Søren Højsgaard
Dear useRs I have uploaded a new package, gRain, for propability propagation in graphical independence networks; sometimes also called probabilistic expertsystems and Bayesian networks. Regards Søren Højsgaard ___ R-packages mailing list [EMAIL PRO

[R] Odp: x-axis

2008-01-28 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 28.01.2008 13:06:16: > Hi, > > I want to plot a graph and here is my code: > >ec<-rep(0,length(e)) > fc<-rep(0,length(f)) Strange. If I understand correctly, you repeat zero length(e) or length(f) times and then you plot those zeroes at x axis = e value

Re: [R] how to skip last lines while reading the data in R

2008-01-28 Thread Prof Brian Ripley
On Mon, 28 Jan 2008, mrafi wrote: > but then the number of levels would reamain the same...!! Please explain: the levels of factors are taken from the data which is actually read. > > > Prof Brian Ripley wrote: >> >> On Mon, 28 Jan 2008, Barry Rowlingson wrote: >> >>> Henrique Dallazuanna wrote

Re: [R] How to fill bar plot with textile rather than color

2008-01-28 Thread yaosheng CHEN
Thank you, and also thanks to John Kane. I did some tests last night, if you plot something like height <- t(t(c(1,-1,1))) bardensity <- t(t(c(10,10,0))) barangle <- t(t(c(45,135,0))) barplot(height, density = bardensity, angle = barangle) You can get grids. Otherwise, only slanted lines or colo

Re: [R] How to fill bar plot with textile rather than color

2008-01-28 Thread John Kane
density ?barplot aa <- c(4,5,6) barplot(aa, density=2, col='red', border="blue") --- CHENYS <[EMAIL PROTECTED]> wrote: > Hi, I'm looking for a tool which can fill bar chart > with dash, skewed line, > or grids, rather than pure color. Any one have the > idea how to do that in R? > Or maybe in M

Re: [R] how to skip last lines while reading the data in R

2008-01-28 Thread mrafi
but then the number of levels would reamain the same...!! Prof Brian Ripley wrote: > > On Mon, 28 Jan 2008, Barry Rowlingson wrote: > >> Henrique Dallazuanna wrote: >>> Perhaps: >>> >>> data <- >>> read.table(textConnection(rev(rev(readLines('data.txt'))[-(1:2)]))) >>> >> >> Euurgh! Am I the

Re: [R] KM estimation for interval censoring?

2008-01-28 Thread Robert Gentleman
The Icens package implements a number of different estimation paradigms. Briefly, for interval censored data, the KM does not apply, but Turnbull's method does. Relevant citations for Icens are in it. For a Cox model approach, there was at one time a package, not at CRAN, by Commenges (and a Bi

Re: [R] how to skip last lines while reading the data in R

2008-01-28 Thread Henrique Dallazuanna
Ok, I agree, but I assume that her don't know the number of rows of file. In your example: rev2 <- function(x,n=100,m=5){ for(i in 1:n){ y=x[1:(length(x)-m)] } return(y) } is needed that open other textConnection - if use the example posted by me. Is there other option? On 28/01

Re: [R] Grouping data via an index

2008-01-28 Thread Peter Dalgaard
Tobin, Jared wrote: > Hello r-help, > > I have a lengthy vector of data (with values anywhere from 1-200), and > another index vector of 'groups' representing values 0-2, 3-5, 6-8, ... > of length 67. The index vector has the structure (1, 4, 7, ... , 196, > 199), where each value is the midpoint

Re: [R] R on an eeePC

2008-01-28 Thread Zembower, Kevin
Doing 'RSiteSearch("eee")' yields some hits. I knew that the ASUS eeePC had come up on r-help. -kevin -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Dr. Walter H. Schreiber Sent: Monday, January 28, 2008 9:32 AM To: r-help@r-project.org Subject: [R] R on

Re: [R] R loops

2008-01-28 Thread Duncan Murdoch
On 1/28/2008 8:48 AM, cvandy wrote: > I'm a new user and am having trouble with loops. > In the following, I'm trying to add the results of "test" and the loops are > not working. > I've simplified the loop. What am I doing wrong? > Thanks! >> test<-numeric(20) >> tot<-numeric(20) >> for(i in 1:20

[R] How to get out the t-test value matrix for a linear regression ?

2008-01-28 Thread leo_aries
Hi, all I've written some R script to calculate the linear regression of a matrix. Here below is my script: >x<-matrix(scan("h:/data/xxx.dat",0),nrow=46,ncol=561,byrow=TRUE) >year <- NULL >year <- cbind(year,as.matrix(x[,1])) >lm.sol<-lm(x~year) >xtrend<-coef(lm.sol)[2,]# get the matri

Re: [R] double-click in RData file versus load( file )

2008-01-28 Thread Cleber Nogueira Borges
hi Duncan, I can't access file greater than 3 mb in my work! :-( :-( ((( unfortunately, is not a joke! I need a new job :-D ))) When I return home, I will do download the patch and report the results to you! :-( Thanks for your attention! Cleber > It is now fixed in the version availabl

Re: [R] Poisson Maximum Likelihood Estimation

2008-01-28 Thread Paul Sweeting
Many thanks for your suggestions on the forum (and to those that have emailed me directly) - I'll make a start on the poisson mle approach shortly! Paul -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Gavin Simpson Sent: 25 January 2008 15:14 To: David Wi

Re: [R] KM estimation for interval censoring?

2008-01-28 Thread Martin Maechler
> "RobG" == Robert Gentleman <[EMAIL PROTECTED]> > on Mon, 28 Jan 2008 08:18:06 -0800 writes: RobG> The Icens package implements a number of different estimation paradigms. RobG> Briefly, for interval censored data, the KM does not apply, but RobG> Turnbull's method does

Re: [R] matrix creation

2008-01-28 Thread Gavin Simpson
hits=-2.6 tests=BAYES_00 X-USF-Spam-Flag: NO Hi Michelle, You don't show your read.csv or read.table call, nor the output of str(obj) where obj is the name of the object you read the data into. I notice that you have explicit 0 and NA. Is there a chance that you have entered NA into the cells th

[R] matrix creation

2008-01-28 Thread Michelle DePrenger-Levin
Hello, I am trying to create multiple matrices (to run a PVA) but can't import all of them from a .csv without the numbers treated as labels and not factors. I can enter the matrix slowly: Site05_96 <- matrix(c(0.07,0,0.03,0.00,NA,0.00, 0.09,0.16667,0.31,0.42,NA,0.00, 0.00,0,0.00,0.0

Re: [R] sub-plot in a plot

2008-01-28 Thread Henrique Dallazuanna
See 'fig' argument in ?par example: x <- rnorm(100) plot(x) par(fig=c(2/3,1,2/3,1), new=T) hist(x, main="") On 28/01/2008, lamack lamack <[EMAIL PROTECTED]> wrote: > > > > > From: [EMAIL PROTECTED]: [EMAIL PROTECTED]: sub-plot in a plotDate: Mon, 28 > Jan 2008 15:43:40 + > > > Dear all, how

[R] Odp: How to get out the t-test value matrix for a linear regression ?

2008-01-28 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 28.01.2008 17:16:37: > Hi, all > > I've written some R script to calculate the linear regression of a matrix. > Here below is my script: > > >x<-matrix(scan("h:/data/xxx.dat",0),nrow=46,ncol=561,byrow=TRUE) > > >year <- NULL > >year <- cbind(year,as.matrix

Re: [R] how to skip last lines while reading the data in R

2008-01-28 Thread Gabor Grothendieck
or even: head(read.table("xy.dat", header = TRUE), -2) On Jan 28, 2008 10:52 AM, Gabor Grothendieck <[EMAIL PROTECTED]> wrote: > If you don't mind reading it in twice its just: > > DF <- read.table("xy.dat", header = TRUE, nrow = > length(readLines("xy.dat")) - 3) > tail(DF) > # or > DF <- read.

Re: [R] using facet_grid() from ggplot2 with additional text in labels

2008-01-28 Thread Rainer M Krug
Dear Thierry On 28/01/2008, ONKELINX, Thierry <[EMAIL PROTECTED]> wrote: > > Dear Rainer, > > I think you'll have to recode the labels of the levels to get the output > that you want. I haven't thought about using levels - that was probably the reason why it threw the sorting out. Thanks Rainer

[R] Come join me on MaLaYsiaN BesT…

2008-01-28 Thread khai
Come join me on MaLaYsiaN BesT! Mari kita berkongsi cerita mengenai apa sahaja mengenai malaysia. cth: gambar lucu, gambar aneh, gambar hantu, artis malaysia, info bola sepak malaysia, sejarah bola sepak, berita sukan bola sepak terkini dan arena bola sepak . Click here to join: http://malaysia

Re: [R] sub-plot in a plot

2008-01-28 Thread Greg Snow
See the subplot function in the TeachingDemos package. Hope this helps, From: [EMAIL PROTECTED] on behalf of lamack lamack Sent: Mon 1/28/2008 8:51 AM To: r-help@r-project.org Subject: [R] sub-plot in a plot From: [EMAIL PROTECTED]: [EMAIL PROTECTED]: sub

Re: [R] Function for translation of a list into a matrix as used by ordination?

2008-01-28 Thread Billy_Schweiger
Thanks for all who commented on this question. It turns out that there is a nice set of functions in library BiodiversityR that translate a stack or list input data file to a community matrix: import.from.Excel(file = file.choose(), sheet = "community", sitenames = "sites", column = "species"

Re: [R] Bug in levels() function?

2008-01-28 Thread Thomas Lumley
This is not a bug; it is deliberately designed this way. There are circumstances when you want to drop levels on subsetting and other circumstances where you don't, so the default behaviour can't make everyone happy. However, there is an option to get the behaviour you want > x<-as.factor(LETT

Re: [R] [OT] vernacular names for circular diagrams

2008-01-28 Thread Rolf Turner
Thanks very much for the enlightenment. Very interesting indeed, and I am glad to find Nightingale exonerated of her purported crime. cheers, Rolf On 29/01/2008, at 8:25 AM, Greg Snow wrote: > I had heard the same thing about Florence Nightingale, but it seems > t

Re: [R] [OT] vernacular names for circular diagrams

2008-01-28 Thread Rolf Turner
Thank you! cheers, Rolf On 29/01/2008, at 8:38 AM, roger koenker wrote: > Howard Wainer (Graphical Discovery, PUP, 2005, p 20) gives > this dubious honor to Playfair (1759- 1823). Nightingale (1820- > 1910) was far too enlightened for this sort of thing, see for examp

Re: [R] Integer vs numeric

2008-01-28 Thread Roland Rau
Christophe Genolini wrote: > Hi the list. > > I do not understand the philosophy behind numeric and integer. > - 1 is numeric (which I find surprising) > - 2 is numeric. > - 1:2 is integer. > Why is that ? > I hope I can answer your question at least partly: Numeric means "double", i.e. intern

Re: [R] [OT] vernacular names for circular diagrams

2008-01-28 Thread Peter Jepsen
If anyone is interested in seeing "Nightingale's Coxcomb" a.k.a. "Nightingale's Rose", it can be seen at http://www.economist.com/world/europe/displaystory.cfm?story_id=10278643. Best regards, Peter. -Oprindelig meddelelse- Fra: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] På vegne af R

Re: [R] survey: estimating a covariance matrix

2008-01-28 Thread Daniel Oberski
Thomas and David, Thanks for your answers. The svyvar() function does the trick. It does not provide a vcov() method for the variance-covariance matrix of these estimates themselves. But I guess I can bootstrap them using library(boot)... Thanks again, daniel On Jan 25, 2008 6:13 PM, Thomas Lu

r-help@r-project.org

2008-01-28 Thread Matthieu Stigler
Hello I'm trying to implement a grid search for a threshold autoregressive model, it is a model in which the regression coefficients are different according to the regimes (under the lower threshold, between lower and upper, over the upper threshold). Estimation of the threshold is made with C

Re: [R] R on an eeePC

2008-01-28 Thread Peter Dalgaard
Dr. Walter H. Schreiber wrote: > Dear list, > > I wonder if somebody has succeeded in installing R on an eeePC > (Xandros desktop). Searching via Rseek (term eeePC) and in eeePC > forums (term Cran) left me without proper hits. > Try just looking for 'eee' at http://tolstoy.newcastle.edu.au/~rki

Re: [R] How to fill bar plot with textile rather than color

2008-01-28 Thread Greg Snow
Filling bars with lines/grids/points is legacy back to the days when the only way to get high quality plots was to use a pen plotter (on the screen you would see bars made of '*' or similar). The pen plotter would use a mechanical arm to draw the lines using a pen/marker. it was easy to have t

Re: [R] basic spatial query

2008-01-28 Thread Rolf Turner
On 28/01/2008, at 11:14 PM, [EMAIL PROTECTED] wrote: > Hello; > Please coud you advise me of a simple way to select points (x,y > coordinates) that fall within a polygon. > I've got a set of polygons, each one defined by an arbitrary number of > points, and several points inside each polygon. > I

Re: [R] How to get out the t-test value matrix for a linear regression ?

2008-01-28 Thread Peter Dalgaard
leo_aries wrote: > Hi, all > > I've written some R script to calculate the linear regression of a > matrix. > Here below is my script: > > >> x<-matrix(scan("h:/data/xxx.dat",0),nrow=46,ncol=561,byrow=TRUE) >> > > >> year <- NULL >> year <- cbind(year,as.matrix(x[,1])) >> > >

Re: [R] matrix creation

2008-01-28 Thread Michelle DePrenger-Levin
I was asked for the following information and hope it might help those who could answer my question... To import the table I used: AsMi05test=read.csv("C:/AsMi_Site05_1998.csv") > str(AsMi05test) `data.frame': 12 obs. of 8 variables: $ X : Factor w/ 6 levels

Re: [R] [OT] vernacular names for circular diagrams

2008-01-28 Thread Rolf Turner
On 28/01/2008, at 12:07 AM, Peter Dalgaard wrote: > Jean lobry wrote: >>> >>> about an hour North of Paris. Her father inquired - >>> coincidentally during the cheese course - what work I was >>> doing in Paris; I replied that I was researching the >>> activities of a Scot, William

Re: [R] Integer vs numeric

2008-01-28 Thread Ted Harding
On 28-Jan-08 21:23:12, Roland Rau wrote: > Christophe Genolini wrote: >> Hi the list. >> >> I do not understand the philosophy behind numeric and integer. >> - 1 is numeric (which I find surprising) >> - 2 is numeric. >> - 1:2 is integer. >> Why is that ? >> > I hope I can answer your question

Re: [R] [OT] vernacular names for circular diagrams

2008-01-28 Thread Greg Snow
I had heard the same thing about Florence Nightingale, but it seems that this is a confusion of different graphs. Nightingale developed a graph based on a circle, but all the angles were equal and the different values were encoded by using different radii of the slices (and she did the right

Re: [R] Package simex

2008-01-28 Thread Liaw, Andy
The original paper is: Cook, J. R. & Stefanski L. A. (1994) Simulation--extrapolation estimation in parametric measurement error models. Journal of the American Statistical Association, 89, 1314-1328. Andy From: Michael Kubovy > > Dear R-helpers, > > It is not clear to me how you get measurem

[R] LOOPS - GIS

2008-01-28 Thread Johan Lassen
Dear R-user community, I need to solve a problem related to river management. The task is to calculate the distance from specified points in a river to the outlet of the river. The river is defined by a database (here called "M"): LinkID |NodeID1 | NodeID2 |Distout |Length

[R] Integer vs numeric

2008-01-28 Thread Christophe Genolini
Hi the list. I do not understand the philosophy behind numeric and integer. - 1 is numeric (which I find surprising) - 2 is numeric. - 1:2 is integer. Why is that ? Christophe __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinf

Re: [R] [OT] vernacular names for circular diagrams

2008-01-28 Thread roger koenker
Howard Wainer (Graphical Discovery, PUP, 2005, p 20) gives this dubious honor to Playfair (1759- 1823). Nightingale (1820- 1910) was far too enlightened for this sort of thing, see for example her letter to Galton about endowing an Oxford professorship in social statistics (reprinted in Karl Pea

Re: [R] maptools no such file

2008-01-28 Thread Rainer Hurling
Roger, I tried your suggestions and ran into same problems as Pieter did before. On 28.01.2008 12:35 (UTC+1), Roger Bivand wrote: > gmail.com> writes: > >> No, I get the same error message there... > > Please do not needlesly delete the thread content. Your original question was: > > "I'm ha

Re: [R] R for Mac 10.3.9

2008-01-28 Thread Benilton Carvalho
http://cran.fhcrc.org/bin/macosx/old/powerpc/R-2.2.0.dmg But that is sooo old. Maybe he could compile from the source (he'd need to install XCode and a fortran compiler). I assume that once these requirements are installed, the compilation would be as simple as: ./configure make May

[R] Very Important Message

2008-01-28 Thread David Ellis
Hello, I am Mr. Ellis David by name and I?m the CEO/OWNER of Davidellis Textiles company Inc. in West Africa. We supply textiles to over hundred of customers spread all over the world.We need someone that will be representing my company as a Manager/Representative, who will attend to customers in

Re: [R] R on an eeePC

2008-01-28 Thread Gabor Csardi
RSiteSearch("R asus") G. On Mon, Jan 28, 2008 at 03:31:45PM +0100, Dr. Walter H. Schreiber wrote: > Dear list, > > I wonder if somebody has succeeded in installing R on an eeePC (Xandros > desktop). Searching via Rseek (term eeePC) and in eeePC forums (term > Cran) left me without proper hit

Re: [R] Fourier Analysis and Curve Fitting in R

2008-01-28 Thread Carson Farmer
Rolf Turner wrote: >> >> On 26/01/2008, at 10:54 AM, Carson Farmer wrote: >> >>> Dear List, >>> >>> I am attempting to perform a harmonic analysis on a time series of snow >>> depth, in which the annual curve is essentially asymmetric (i.e. snow >>> accumulates slowly over time, and the subsequent

Re: [R] how to skip last lines while reading the data in R

2008-01-28 Thread Romain Francois
Hi, What about this : data <- read.table( pipe('head -n-2 data.txt') ) However, the aesthetics will only appeal to people having head (maybe it needs to be a recent version as well ...). Cheers, Romain Barry Rowlingson wrote: > Henrique Dallazuanna wrote: > >> Perhaps: >> >> data <- read

[R] (no subject)

2008-01-28 Thread Loai Mahmoud Awad Alzoubi
Hi all I am trying to generate a normal unbalanced data to estimate the coefficients of LM, LMM, GLM, and GLMM and their standard errors. Also, I am trying to estimate the variance components and their standard errors. Further, I am trying to use the likelihood ratio test to test H0: sigma^2_b

Re: [R] Integer vs numeric

2008-01-28 Thread Ted Harding
On 28-Jan-08 22:40:02, Peter Dalgaard wrote: > [...] > AFAIR, space is/was more of an issue. If you do something like > > for i in 1:1e7 > some.silly.simulation() > > then you have 40 MB sitting there doing nothing, and 80 MB if > it had been floating point. Hmmm ... there's something to be

Re: [R] maptools no such file

2008-01-28 Thread Roger Bivand
On Mon, 28 Jan 2008, Rainer Hurling wrote: > Roger, > > I tried your suggestions and ran into same problems as Pieter did before. > Rainer: I don't think that we know what platform Pieter was using, apart from the working directory that looked like Windows, and where the problem could not be r

Re: [R] Integer vs numeric

2008-01-28 Thread Peter Dalgaard
(Ted Harding) wrote: > Further to the above: The help > > ?":" > > says: > > Value: > For numeric arguments [as opposed to factors], > a numeric vector. This will be of type 'integer' > if 'from' and 'to' are both integers and > representable in the integer type, otherwise of > type 'nu

[R] trouble installing quantreg packages

2008-01-28 Thread Arjun Adhikari
I am trying to install the "quantreg Package" on MacOSX. When I choose CRAN mirror, it gives me the error: Error in read.dcf(file = tmpf) : Line longer than buffer size In addition: There were 44 warnings (use warnings() to see them) Is there any known current problem with this package? Sincerely

Re: [R] Integer vs numeric

2008-01-28 Thread Christophe Genolini
> Further to the above: The help > > ?":" > > says: > > Value: > For numeric arguments [as opposed to factors], > a numeric vector. This will be of type 'integer' > if 'from' and 'to' are both integers and > representable in the integer type, otherwise of > type 'numeric' ??? This is

[R] R for Mac 10.3.9

2008-01-28 Thread Paul Johnson
I know nothing of Macintosh, so please be patient. My student has a Macintosh with OSX 10.3.9. The R for Mac says it is for 10.4.4 or higher. Aside from saying "get a new Mac," what can be said to my student? Can you point me at the newest version of R that did work on 10.3 ? pj -- Paul E. Jo

[R] Problems compiling shogun 0.4.4

2008-01-28 Thread Brian McNally
Hello, I'm trying to compile Shogun 0.4.4 from source on Redhat Enterprise Linux 4 (update 5). I've compiled and installed R-2.6.1 already and that seems to work fine. I had some initial problems finding R development files, which were fixed by setting CFLAGS='-I/tmp/R-2.6.1/ include' and C

[R] [R-pkgs] eRm: new version 0.9-6

2008-01-28 Thread Patrick Mair
Dear useRs, new and updated features in eRm 0.9-6 for extended Rasch modeling: - infit and outfit mean-square statistics added to functions personfit() and itemfit(). - new method: plotPImap() plots a person-item map (cf. Bond & Fox, 2007), i.e., a map of locations of item (and threshold) para

[R] [OT] - standard errors for parameter estimates under ridge regression and lasso?

2008-01-28 Thread Andrew Robinson
Dear R community, I'm curious to know how people go about estimating standard errors for parameter estimates after model selection by ridge regression and the lasso. Do you have any practical or theoretical advice? Warmly, Andrew -- Andrew Robinson Department of Mathematics and Statistics

Re: [R] matrix creation

2008-01-28 Thread Gavin Simpson
hits=-2.6 tests=BAYES_00 X-USF-Spam-Flag: NO On Mon, 2008-01-28 at 12:17 -0700, Michelle DePrenger-Levin wrote: > I was asked for the following information and hope it might help those who > could answer my question... That looks fine to me Michelle. You will have problems with as.matrix on thi

Re: [R] [OT] vernacular names for circular diagrams

2008-01-28 Thread hadley wickham
On Jan 28, 2008 1:25 PM, Greg Snow <[EMAIL PROTECTED]> wrote: > I had heard the same thing about Florence Nightingale, but it seems that this > is a confusion of different graphs. > > Nightingale developed a graph based on a circle, but all the angles were > equal and the different values were en

[R] question about order

2008-01-28 Thread Waverley
I have a data vector as following: > z [1] 183.1370 201.9610 113.7250 140.7840 156.2750 42.1569 42.1569 42.1569 [9] 240.1960 308.4310 42.1569 42.1569 42.1569 42.1569 42.1569 42.1569 [17] 42.1569 42.1569 42.1569 42.1569 279.8040 42.1569 42.1569 when I sort, it gave me the right ord

Re: [R] question about order

2008-01-28 Thread Benilton Carvalho
that seems right order() gives you the indexes "idx" such that x[idx] == sort(x) > set.seed(123) > x <- rnorm(10) > idx <- order(x) > identical(x[idx], sort(x)) [1] TRUE best b On Jan 28, 2008, at 8:19 PM, Waverley wrote: I have a data vector as following: z [1] 183.1370 201.9610 113.72

Re: [R] question about order

2008-01-28 Thread Achim Zeileis
On Mon, 28 Jan 2008, Waverley wrote: > I have a data vector as following: > > z > [1] 183.1370 201.9610 113.7250 140.7840 156.2750 42.1569 42.1569 42.1569 > [9] 240.1960 308.4310 42.1569 42.1569 42.1569 42.1569 42.1569 42.1569 > [17] 42.1569 42.1569 42.1569 42.1569 279.8040 42.1569

Re: [R] Fourier Analysis and Curve Fitting in R

2008-01-28 Thread stephen sefick
well if you want to find the spectral density aka what frequencies explain most of the variance then I would suggest the spectral density. This can be implemented with spec.pgram(). This is conducted with the fast fourier transform algorithm. > a<-ts(data, frequency = 1) #make the time series w

Re: [R] Package simex

2008-01-28 Thread Michael Kubovy
On Jan 28, 2008, at 11:30 AM, Liaw, Andy wrote: > The original paper is: > > Cook, J. R. & Stefanski L. A. (1994) Simulation--extrapolation > estimation in parametric measurement error models. Journal of the > American Statistical Association, 89, 1314-1328. > > From: Michael Kubovy >> >> Dea

[R] add/subtract matrices, ignoring NA or missing values

2008-01-28 Thread Ng Stanley
Hi, For example, given two 2x2 matrices m1 and m2. I would like to add/subtract element by element > m1 [,1] [,2] [1,] NA NA [2,]12 > m2 [,1] [,2] [1,]1 NA [2,] NA2 > m1 + m2 [,1] [,2] [1,] NA NA [2,] NA4 How can I ignore the NA, and g

[R] how to simulate seasonal cointegrated series

2008-01-28 Thread tom soyer
Hi, does anyone know how to simulate two seasonal data series that are cointegrated? Thanks! -- Tom [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the pos

Re: [R] add/subtract matrices, ignoring NA or missing values

2008-01-28 Thread Rolf Turner
On 29/01/2008, at 3:34 PM, Ng Stanley wrote: > Hi, > > For example, given two 2x2 matrices m1 and m2. I would like to add/ > subtract > element by element > >> m1 > [,1] [,2] > [1,] NA NA > [2,]12 > >> m2 > [,1] [,2] > [1,]1 NA > [2,] NA2 > >> m1 + m2 >

  1   2   >