Re: [R] Registry problem when loading RMySQL

2008-12-10 Thread Prof Brian Ripley
On Wed, 10 Dec 2008, Gerald van der Grijn wrote: Hi all, I am exploring the possibility of using R in my meteorological verification project. I have installed R on my Windows (XP) desktop machine and have also installed the RMySQL and verification packages. However, I run into problems when I t

Re: [R] How can I draw bars

2008-12-10 Thread Jim Lemon
Raymond Balise wrote: I need to make a graphic to show problems on different parts of chromosomes (think of a graphic showing the number of frayed threads as colors along different parts of a worn out rope). I want to draw bars going from left to right across a page and color different parts of

[R] plot Geneland result in a map

2008-12-10 Thread vincenzo landi
Hello, anybody know the procedure to plot the geneland result in a map? thank you ** Vincenzo Landi Post Doctorate student Animal genomic and breeding cell:0039/3395388713 Fax. 075-5857122 [[alternative HTML version del

Re: [R] expected variable name error pos 98349 WInBUGS in R

2008-12-10 Thread Uwe Ligges
Googling for WinBUGS leads you to http://www.mrc-bsu.cam.ac.uk/bugs/winbugs/contents.shtml which contains a link to the BUGS list. Uwe Ligges Anamika Chaudhuri wrote: I did run the same model in WinBUGS, letting WinBUGS generate initials but now it gives me the error 'data value already given

Re: [R] How to add accuracy, sensitivity, specificity to logistic regression output?

2008-12-10 Thread David Hajage
You could calculate this yourself, with the result object ( yourobject$linear.predicators). 2008/12/9 pufftissue pufftissue <[EMAIL PROTECTED]> > Hi, > > Is there a way when doing logistic regression for the output to spit out > accuracy, sensitivity, and specificity? > > I would like to know the

Re: [R] R and Scheme

2008-12-10 Thread Duncan Murdoch
[deletions, including attribution, which I think were Stavros then Luke then Peter:] In R, most data types (including numeric vectors) do not have a standard external representation which can be read back in without evaluation. The default print form is not readable in this sense but dput is a

[R] Registry problem when loading RMySQL

2008-12-10 Thread Gerald van der Grijn
Hi all, I am exploring the possibility of using R in my meteorological verification project. I have installed R on my Windows (XP) desktop machine and have also installed the RMySQL and verification packages. However, I run into problems when I try to load RMySQL. I get the following error:

Re: [R] Logical inconsistency

2008-12-10 Thread emma jane
Thanks Greg, that does make sense.  And I've solved the problem by rounding the variables before taking the difference between them. Thanks to all who replied. Emma Jane  From: Greg Snow <[EMAIL PROTECTED]> .com.br>; Wacek Kusnierczyk <[EMAIL PROTECTED]>;

[R] plot multivariate function with known coefficient.

2008-12-10 Thread YIHSU CHEN
Dear R users -- I have a simple question that I cannot find the answers in the archive email: how to plot a multivariate function if I know the coefficients already? For instance, if I like to plot a 3D fig of Pi against X and Y with the function Pi = a0 + a1X + a2Y, where a0, a1 & a2 are known.

[R] Help with Table headers

2008-12-10 Thread Haroon Malik
Hi, I have a little problem with table header. I am enclosing my TOY example: I have a tab separated file(tsv) with two columns as follow: \\disk(c:) \\processor(1) 322 40 422 60 I read the file using read.table as follow: tab<-read.table("c:/R/file1.

[R] How to Break Axis in Lattice Plot

2008-12-10 Thread Larry Ma
Dear All, I have searched the following problem in the R help achives, and there did not seem to be solutio for it. The problem is how to break the axis in the xyplot (lattice plot). I read the posts using the package plotrix, buy it doesn't work with the xyplot, which is lattice based. A

[R] convert dataframe to matrix for cmdscale

2008-12-10 Thread William Simpson
I have a dataframe like this (toy example): x y z "a" "a" 0 "a" "b" 1 "a" "c" 2 "b" "a" .9 "b" "b" 0 "b" "c" 1.3 "c" "a" 2.2 "c" "b" 1.1 "c" "c" 0 The observations are from a matrix like this: c 2.2 1.1 0.0 b 0.9

Re: [R] plot multivariate function with known coefficient.

2008-12-10 Thread Gavin Simpson
On Tue, 2008-12-09 at 16:07 -0800, YIHSU CHEN wrote: > Dear R users -- > > I have a simple question that I cannot find the answers in the archive > email: how to plot a multivariate function if I know the coefficients > already? For instance, if I like to plot a 3D fig of Pi against X and > Y wi

Re: [R] Corrupt Header<==Help

2008-12-10 Thread Henrique Dallazuanna
Try this: x <- read.table("c:/R/file1.tsv", sep = '\t', header = TRUE, check.names = FALSE) On Wed, Dec 10, 2008 at 1:06 AM, Haroon Malik <[EMAIL PROTECTED]>wrote: > Hi, > > > > I have a little problem with table header. I am enclosing my TOY example: > > > > I have a tab separated file(tsv) w

[R] how to set font size in 'pt'

2008-12-10 Thread Hyunchul Kim
Dear all, I need to make a publication quality figure. Axis labels and xy labels should be written in a given font type and font size in "pt" but not in "px" or "inch". I tried dev2bitmap(... , pointsize= ) and pdf(... , pointsize= ) but I couldn't do what I want. I heard that many people use R t

Re: [R] Corrupt Header<==Help

2008-12-10 Thread Haroon Malik
Thank you it works to some extent, now I get two extra back slashes: > x disk(c:) processor(1) 1 322 40 2 422 60 From: Henrique Dallazuanna [mailto:[EMAIL PROTECTED] Sent: December-10-08 5:18 AM To: Haroon Malik Cc: r-help@r-project.o

[R] Corrupt Header<==Help

2008-12-10 Thread Haroon Malik
Hi, I have a little problem with table header. I am enclosing my TOY example: I have a tab separated file(tsv) with two columns as follow: \\disk(c:) \\processor(1) 322 40 422 60 I read the file using read.table as follow: tab<-read.table("c:/R/file1.

[R] Help with Table headers

2008-12-10 Thread Haroon Malik
Hi, I have a little problem with table header. I am enclosing my TOY example: I have a tab separated file(tsv) with two columns as follow: \\disk(c:) \\processor(1) 322 40 422 60 I read the file using read.table as follow: tab<-read.table("c:/R/file1.

[R] Multpile (45x8) graphs of the same page / device: titles crammed

2008-12-10 Thread tsunhin wong
Dear R users, I'm trying to plot 45x8 graphs on the same pdf / device for the sake of visual comparison. par(mfcol=c(45,8)) par(mai=c(0,0,0,0)) In ?title, I can see there are cex and font settings: I set cex = 0.01 and font = 1: it is still very large, and then I tried setting font < 1, e.g. font

Re: [R] convert dataframe to matrix for cmdscale

2008-12-10 Thread William Simpson
On re-reading ?cmdscale I see that I can also use cmdscale(d) on a full matrix rather than just the lower triangle. d: a distance structure such as that returned by 'dist' or a full symmetric matrix containing the dissimilarities. So how to get a full symmetric matrix from a dataframe l

Re: [R] Stepwise regression

2008-12-10 Thread Frank E Harrell Jr
Maithili Shiva wrote: Hi, I have the response variable 'Y' and four predictors say X1, X2, X3 and X4. Assuming all the assmptions like Y follows normal distribution etc. hold and I want to run linear multiple regression. How do I run the stepwise regression (forward as well as the backward re

Re: [R] Polar coordinates - contour plots

2008-12-10 Thread Todor Kondic
Thanks Bert for the answer! I'll try to find an appropriate package. In the end, i have "fiddling". But, can anyone tell me, were I to create an underlying implementation of cooridnate transforms, perhaps controlled by user via 'par', would I need to change the handling of e.g x,y data by every di

[R] Stepwise regression

2008-12-10 Thread Maithili Shiva
Hi, I have the response variable 'Y' and four predictors say X1, X2, X3 and X4. Assuming all the assmptions like Y follows normal distribution etc. hold and I want to run linear multiple regression. How do I run the stepwise regression (forward as well as the backward regression). >From other

Re: [R] Multpile (45x8) graphs of the same page / device: titles crammed

2008-12-10 Thread baptiste auguie
Hi, "font" should be an integer as described in ?par. I think you want to play with cex.main (possibly cex in combination, depending on what your are plotting) x <- seq(0, 10) pdf(width=8, height=20) par(mfrow=c(45, 8), mai=c(0,0.1,0.1,0)) sapply(1:(45*8), function(ii) { plot(x,

[R] R Graphical parameters

2008-12-10 Thread Metcalfe, John
Dr. Ripley, Thanks much for your reply. I've attached a representative sample; I apologize for not including this with my previous email. Best, John John Z. Metcalfe, M.D., M.P.H. Division of Pulmonary and Critical Care Medicine University of California, San Francisco San Francisco General Hos

Re: [R] convert dataframe to matrix for cmdscale

2008-12-10 Thread Henrique Dallazuanna
Try this: x$zMean <- ave(x$z, apply(x[,1:2], 1, function(x)paste(sort(x), collapse = "")), FUN = mean) xtabs(zMean ~ x + y, data = x) On Wed, Dec 10, 2008 at 11:16 AM, William Simpson < [EMAIL PROTECTED]> wrote: > On re-reading ?cmdscale I

Re: [R] Corrupt Header<==Help

2008-12-10 Thread Henrique Dallazuanna
see this: cat(names(x), '\n') The extra back slashes is a escape character. On Wed, Dec 10, 2008 at 10:15 AM, Haroon Malik <[EMAIL PROTECTED]>wrote: > Thank you it works to some extent, now I get two extra back slashes: > > > x > > disk(c:) processor(1) > > 1 322

Re: [R] Multpile (45x8) graphs of the same page / device: titles crammed

2008-12-10 Thread Prof Brian Ripley
font is an integer number (1,2,3,4,5 on that device): please do read ?par. You are looking of 'pointsize': see ?pdf, but I think rather that you should use a larger plotting surface (see 'height' and 'width') On Wed, 10 Dec 2008, tsunhin wong wrote: Dear R users, I'm trying to plot 45x8 gr

Re: [R] plot Geneland result in a map

2008-12-10 Thread Christoph Heibl
Hi Vincenzo, have you looked at ?PlotTesselation (and ?PostProcessChain)? If you need to export these maps in binary format (e.g., to import it into GRASS GIS) you can use the function 'geneland2ascii' in the haplo package on my website. Best wishes

Re: [R] Polar coordinates - contour plots

2008-12-10 Thread hadley wickham
On Wed, Dec 10, 2008 at 4:33 AM, Todor Kondic <[EMAIL PROTECTED]> wrote: > Thanks Bert for the answer! I'll try to find an appropriate package. > In the end, i have "fiddling". > > But, can anyone tell me, were I to create an underlying implementation > of cooridnate transforms, perhaps controlled

Re: [R] How to get Greenhouse-Geisser epsilons from anova?

2008-12-10 Thread Skotara
Dear John, thank you for the kind offer! Sorry, I just made a mistake anywhere I can not trace back, now it works as you described it. Thank you again! Dear Peter, thank you for the information, I did not know about the quotation marks. It indeed works using "G-G Pr"! The SPSS and R output for

Re: [R] controlling axes in plot.cuminc (cmprsk library)

2008-12-10 Thread Arthur Allignol
Hi, You could try to use the option xaxt = "n" in plot.cuminc, instead of axes. Hope that helps, Arthur Amy Krambrink wrote: Dear R-help list members, I am trying to create my own axes when plotting a cumulative incidence curve using the plot.cuminc function in the CMPRSK library. The defa

Re: [R] Multpile (45x8) graphs of the same page / device: titles crammed

2008-12-10 Thread tsunhin wong
What if cex=0.001 with font=1 will still have the title font too large? (so that the upper halves of title is cut off and cannot be seen?) Will there be a workaround? Can there be a work around so that I can label each small graphs with the info I need in a "smaller" way? Thanks! Regards, J

Re: [R] glm error message when using family Gamma(link=inverse)

2008-12-10 Thread Terry Therneau
John, You have specified a model with E(y) = 1/eta where eta = X beta is the linear predictor and E(y) must be >0, since the family is Gamma and you have a lot of covariates in the model. glm now has to try to find a best linear predictor, but under the

Re: [R] R and Scheme

2008-12-10 Thread Wacek Kusnierczyk
Wacek Kusnierczyk wrote: > Stavros Macrakis wrote: > > >> There is no equivalent to set-car!/rplaca (not even pairlists and >> expressions). For example, r<-pairlist(1,2); r[[1]]<-r does not create a >> circular list. And in general there doesn't seem to be substructure sharing >> at the semant

Re: [R] Multpile (45x8) graphs of the same page / device: titles crammed

2008-12-10 Thread Detlef Steuer
-BEGIN PGP SIGNED MESSAGE- Hash: SHA1 On Wed, 10 Dec 2008 09:22:55 -0500 "tsunhin wong" <[EMAIL PROTECTED]> wrote: > What if cex=0.001 with font=1 will still have the title font too > large? (so that the upper halves of title is cut off and cannot be > seen?) > Will there be a workaround?

[R] xyplot sorted by date

2008-12-10 Thread Ptit_Bleu
Hello I have a dataframe with 3 columns C1, C2 and C3. C1 and C2 are numerical data and C3 is the date of the data (format : %Y-%m-%d %H:%M:%S") With the lattice package, I would like to plot C1 as function of C2 sorted by the dates, that is something like : xyplot(C1 ~ C2 | as.Date(C3)) It plo

Re: [R] R and Scheme

2008-12-10 Thread Stavros Macrakis
On Wed, Dec 10, 2008 at 4:40 AM, Duncan Murdoch <[EMAIL PROTECTED]>wrote: > [deletions, including attribution, which I think were Stavros then Luke > then Peter:] >> >> In R, most data types (including numeric vectors) do not have a standard external representation which can be read back in w

[R] how to merge panel data stored by variable?

2008-12-10 Thread Viktor Nagy
Hi, I have two datasets stored in tab-separated format in the following way file1: country year1year2 Germanyvar1 var1 Hungary var1 var1 file2: country year1year2 Germanyvar2 var2 Hungary var2 va

Re: [R] Multpile (45x8) graphs of the same page / device: titles crammed

2008-12-10 Thread tsunhin wong
I found I have to use font=1, and I found it works when I set: par(mfcol=c(45,8)) par(mar=c(0,0.5,1,0.5)) instead of par(mfcol=c(45,8)) par(mai=c(0,0,0,0)) On Wed, Dec 10, 2008 at 9:22 AM, tsunhin wong <[EMAIL PROTECTED]> wrote: > What if cex=0.001 with font=1 will still have the title font too >

Re: [R] how to merge panel data stored by variable?

2008-12-10 Thread Stefan Grosse
> I have two datasets stored in tab-separated format in the following way > file1: > country year1year2 > Germanyvar1 var1 > Hungary var1 var1 > > file2: > country year1year2 > Germanyvar2 var2 > Hungary v

Re: [R] xyplot sorted by date

2008-12-10 Thread Stefan Grosse
> With the lattice package, I would like to plot C1 as function of C2 sorted > by the dates, that is something like : > xyplot(C1 ~ C2 | as.Date(C3)) > > It plots somethings but not the truth. > I have compared the output with the output of this script > > I don't know whether it is ther perfe

Re: [R] How to get Greenhouse-Geisser epsilons from anova?

2008-12-10 Thread John Fox
Dear Nils, > -Original Message- > From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On > Behalf Of Skotara > Sent: December-10-08 9:22 AM > To: Peter Dalgaard > Cc: r-help@r-project.org; John Fox > Subject: Re: [R] How to get Greenhouse-Geisser epsilons from anova? > > Dear John, > thank

[R] mixed exponential distribution

2008-12-10 Thread Fazekas, Jacob
Good morning, Is there anyway to do Mixed Exponential Distribution in R? I am trying to load some lag-weighted empirical survival distribution into R and run a mixed exponential on that data. Thanks, Jacob Fazekas Jacob Fazekas Assistant Actuary Auto-Owners Insurance Company 517-703-2543 [EMA

[R] 2-Y-axes on same plot

2008-12-10 Thread Joe Trubisz
Hi... Is this possible in R? I have 2-sets of data, that were collected simultaneously using 2- different data acquisition schemes. The x-values are the same for both. The y-values have different ranges (16.4-37.5 using one method, 557-634 using another). In theory, if you plot both plots

[R] RMySQL Keeps crashing

2008-12-10 Thread PeterBvolk
Hi all, I am working on MySQL 5.0 with RMySQL 0.7-2. All commands except for Select statements work. I can create tables, Truncate then etc. When I issue a Select statement with dbSendQuery containing more than one column in the selection then R crashes. I have been looking around but could not f

[R] request information

2008-12-10 Thread Lisbeth Garzon
Dear R-Help team, I'm researcher University Nacional of Colombia, I'll use R for my research, and I would want to know, how i can to estimate Logistic Regression model such as binary logistic model, multinomial logistic model, and logit models with repeated measurement, for i can estimate discrete

[R] .Fortran and memory allocation

2008-12-10 Thread Gunnar Hellmund
Hi I know it is not recommended to use dup=FALSE with .Fortran calls, but ... 1) Is it possible to allocate memory for a large array of type double or integer without initializing the array (and how is this done) 2) If the answer to question 1) is 'Yes': Is it possible to pass the uninitialized

Re: [R] 2-Y-axes on same plot

2008-12-10 Thread Ben Tupper
On Dec 10, 2008, at 9:35 AM, Joe Trubisz wrote: Hi... Is this possible in R? I have 2-sets of data, that were collected simultaneously using 2- different data acquisition schemes. The x-values are the same for both. The y-values have different ranges (16.4-37.5 using one method, 557-634 u

[R] line length in legend

2008-12-10 Thread Giovanni Petris
Hello, Is it possible to set the length of the lines appearing in a legend? For example, with > plot(window(EuStockMarkets, end = 1992), plot.type = "s", lty = c("13", > "6413", "431313", "B4")) > legend(1991.7, 2200, legend = colnames(EuStockMarkets), lty = c("13", "6413", > "431313", "B4"))

Re: [R] request information

2008-12-10 Thread Stefan Grosse
Lisbeth Garzon schrieb: > Dear R-Help team, > > I'm researcher University Nacional of Colombia, I'll use R for my research, > and I would want to know, how i can to estimate Logistic Regression model > such as binary logistic model, multinomial logistic model, and logit models > with repeated measu

[R] tapply within a data.frame: a simpler alternative?

2008-12-10 Thread baptiste auguie
Dear list, I have a data.frame with x, y values and a 3-level factor "group", say. I want to create a new column in this data.frame with the values of y scaled to 1 by group. Perhaps the example below describes it best: x <- seq(0, 10, len=100) my.df <- data.frame(x = rep(x, 3), y=c(3*sin(

Re: [R] Logical inconsistency

2008-12-10 Thread Stephan Kolassa
Hi Emma, unfortunately, rounding variables before taking the difference will not solve your problem, because the *rounded* variables are subject to the same (effectively) random internal representation. Examples: > round(8.3,20)-round(7.3,20)>=1 [1] TRUE > round(2.3,20)-round(1.3,20)>=1 [1] F

Re: [R] tapply within a data.frame: a simpler alternative?

2008-12-10 Thread Chuck Cleland
On 12/10/2008 12:02 PM, baptiste auguie wrote: > Dear list, > > I have a data.frame with x, y values and a 3-level factor "group", say. > I want to create a new column in this data.frame with the values of y > scaled to 1 by group. Perhaps the example below describes it best: > >> x <- seq(0, 10,

Re: [R] tapply within a data.frame: a simpler alternative?

2008-12-10 Thread hadley wickham
On Wed, Dec 10, 2008 at 11:02 AM, baptiste auguie <[EMAIL PROTECTED]> wrote: > Dear list, > > I have a data.frame with x, y values and a 3-level factor "group", say. I > want to create a new column in this data.frame with the values of y scaled > to 1 by group. Perhaps the example below describes i

Re: [R] tapply within a data.frame: a simpler alternative?

2008-12-10 Thread Peter Dalgaard
baptiste auguie wrote: > Dear list, > > I have a data.frame with x, y values and a 3-level factor "group", say. > I want to create a new column in this data.frame with the values of y > scaled to 1 by group. Perhaps the example below describes it best: > >> x <- seq(0, 10, len=100) >> my.df <- da

Re: [R] tapply within a data.frame: a simpler alternative?

2008-12-10 Thread baptiste auguie
Excellent! I completely forgot its name and existence. Perhaps ave should be mentioned on the help page of either by, tapply, split. Many thanks, baptiste On 10 Dec 2008, at 17:20, Chuck Cleland wrote: On 12/10/2008 12:02 PM, baptiste auguie wrote: Dear list, I have a data.frame with x,

Re: [R] exporting rast from R to GRASS

2008-12-10 Thread Roger Bivand
Thybério Luna Freire gmail.com> writes: > > Hi, everybody! > i created a imagem by kriging using geoR package. I imported points > from GRASS("zn", after converted to geodata "zn_geo"), the border > "zn_border" and a raster mask. Then i interpolated the points by > kriging and created a raster

Re: [R] tapply within a data.frame: a simpler alternative?

2008-12-10 Thread baptiste auguie
On 10 Dec 2008, at 17:25, hadley wickham wrote: On Wed, Dec 10, 2008 at 11:02 AM, baptiste auguie <[EMAIL PROTECTED]> wrote: Dear list, I have a data.frame with x, y values and a 3-level factor "group", say. I want to create a new column in this data.frame with the values of y scaled to

Re: [R] line length in legend

2008-12-10 Thread Gabor Grothendieck
Perhaps you could just place the labels right on the graphs: library(zoo) z <- as.zoo(window(EuStockMarkets, end = 1992)) plot(z, screen = 1) text(time(z)[1], z[1], colnames(z), cex = .6) On Wed, Dec 10, 2008 at 11:58 AM, Giovanni Petris <[EMAIL PROTECTED]> wrote: > > Hello, > > Is it possible t

Re: [R] line length in legend

2008-12-10 Thread baptiste auguie
From the code of legend() the length seems to be hard-wired (seg.len = 2). You could copy the code and add this "seg.len" as a free parameter in your own custom function. An alternative is to use the lattice package which has a size argument for this purpose. See ? xyplot in the key section.

[R] First call to constructor fails (R.oo)

2008-12-10 Thread Valentin Ruano
Hi all, I'm trying to do some object-oriented programming in R using the R.oo package. Right from the start I encountered a strange (at least for me) problem. I define a new class/constructor based on the R.oo documentation. However the first attempt to create an object fails: === Code: === li

Re: [R] First call to constructor fails (R.oo)

2008-12-10 Thread Henrik Bengtsson
Hi, Any R.oo Object must accept no arguments, e.g. Object(). Thus, in your case you need to be able to support the call MyClass(). The reason is that the first time you use the class, R.oo will create a static object of the class by calling (). Typically, I would do it this way: setConstructor

Re: [R] Logical inconsistency

2008-12-10 Thread Gustavo Carvalho
Hello, An alternative to round(): isTRUE(all.equal((2.3-1.3),1)) Regards, Gustavo. On Wed, Dec 10, 2008 at 3:12 PM, Stephan Kolassa <[EMAIL PROTECTED]> wrote: > Hi Emma, > > unfortunately, rounding variables before taking the difference will not > solve your problem, because the *rounded* vari

Re: [R] Better way to find distances between points in a set?

2008-12-10 Thread cgw
Hah! I DID do my homework first. From now on I guess I'll have to post the complete list of trails I followed before deciding I was lost. That appears to be the only way to stave off the Harpies. So, before I posted, and before I even wrote my toy script, I found and read thru

[R] read.table with different row lengths

2008-12-10 Thread Chris Poliquin
Hi, I need to read in a series of text files with a time series on each row. The series are of different lengths and I'd like to just use the first row as the length and have R ignore extra values in rows that go over this length. For example: 1 0 3 4 5 1 3 5 6 8 7 7 2 1 1 1 4 7 7 7 So

Re: [R] RCurl::postForm() -- how does one determine what the names are of each form element in an online html form?

2008-12-10 Thread Tony Breyal
Thank you Felix and also to the individual who replied off-list. re: html code -- you are both indeed correct that the form elements are named within the html code for a simple form, and i thank you both for letting me know about this. For simple forms i think i will try and write myself a functio

Re: [R] convert dataframe to matrix for cmdscale

2008-12-10 Thread William Simpson
Thanks Henrique Here's what I came up with temp<-read.table(filename, header=TRUE) tempm<-tapply(temp[,3],temp[,1:2],c) #put dataframe into matrix form tempm<-(tempm + t(tempm))/2 #add matrix to matrix flipped about diagonal and divide by 2 Bill On Wed, Dec 10, 2008 at 1:29 PM, Henrique Dallazu

[R] read.*: How to read from a URL?

2008-12-10 Thread Michael Friendly
The question is how to use a URL in place of a file= argument for read.*.functions that do not support this internally. e.g., utils::read.table() and her family all support a file= argument that can take a URL equally well as a local file. So, if I have a file on the web, I can equally well d

Re: [R] read.table with different row lengths

2008-12-10 Thread Chris Poliquin
For some reason R is reading in the file and creating new rows with fill=TRUE that begin with NA. In other words, the R table I end up with has more rows than my text file. These are very long time series (3500 observations). Each of my text file has about 200 rows. With the example I

Re: [R] read.table with different row lengths

2008-12-10 Thread Marc Schwartz
on 12/10/2008 12:50 PM Chris Poliquin wrote: > Hi, > > I need to read in a series of text files with a time series on each > row. The series are of different lengths and I'd like to just use the > first row as the length and have R ignore extra values in rows that go > over this length. > > For

Re: [R] read.*: How to read from a URL?

2008-12-10 Thread hadley wickham
Hi Michael, In general, I think you should be able to do: gimage <- read.jpeg(url(gimageloc)) or alternatively use the EBImage from bioconductor which will read from a url automatically (it also opens a much wider range of file types) library(EBImage) img <- readImage(gimageloc, TrueColor) Had

Re: [R] Better way to find distances between points in a set?

2008-12-10 Thread Charles C. Berry
On Wed, 10 Dec 2008, [EMAIL PROTECTED] wrote: Hah! I DID do my homework first. From now on I guess I'll have to post the complete list of trails I followed before deciding I was lost. That appears to be the only way to stave off the Harpies. So, before I posted, and before I even wro

Re: [R] how to merge panel data stored by variable?

2008-12-10 Thread Viktor Nagy
2008/12/10 Stefan Grosse <[EMAIL PROTECTED]> > > > I have two datasets stored in tab-separated format in the following way > > file1: > > country year1year2 > > Germanyvar1 var1 > > Hungary var1 var1 > > > > file2: > > country year1

Re: [R] read.*: How to read from a URL?

2008-12-10 Thread Duncan Murdoch
If you use url("http://euclid.psych.yorku.ca/SCS/Gallery/Private/langrens.csv";) then most functions will work, because that creates a connection to read the URL. Most functions that can read files do it through connections. Duncan Murdoch On 12/10/2008 2:17 PM, Michael Friendly wrote: The q

Re: [R] how to merge panel data stored by variable?

2008-12-10 Thread hadley wickham
On Wed, Dec 10, 2008 at 1:54 PM, Viktor Nagy <[EMAIL PROTECTED]> wrote: > 2008/12/10 Stefan Grosse <[EMAIL PROTECTED]> >> >> > I have two datasets stored in tab-separated format in the following way >> > file1: >> > country year1year2 >> > Germanyvar1 var1 >> > Hu

[R] missing argument

2008-12-10 Thread jonas garcia
Dear list, I have a question and I'm going to give an example of my problem f<- function(d1, d2, d3) { d<- d1*d2/d3 return(d) } v1<- 1 v2<- 2 If I try f(v1, v2, v3) Error in f(v1, v2, v3) : object "v3" not found I obviously got the above error message. I would like to add something to my functio

[R] Linking a library with init(argc,argv) function

2008-12-10 Thread Maxime Debon
Hi, I would like to embed the LibMesh library in R. I have already linked C++ code and build libraries as it is explained in the R-extension manual with trivial examples : dyn.load(paste("/home/default/LibFooR",.Platform$dynlib.ext,sep="")) I am facing a difficulty with loading a library d

[R] (no subject)

2008-12-10 Thread Sara Sievert
Hello! I've fitted a time series to the exponential power density (Subbotin) with the following function: library(MASS) dep<- function(x,location,scale,tail) # Exponential power (=Subbotin) { const<- 2*scale*tail^(1/tail) *gamma(1+1/tail) z<- (x-location)/scale exp(-1/tail*abs(z)^tail)/con

[R] Calculating Leverage or studentized residuals with gls

2008-12-10 Thread dt Excellent
Dear colleagues, hi. I am using the function gls in R-package nlme and I want to compute either the leverage values or the studentized residuals. Can any body tell me how can do that? With regards ___ ΧρησιμοποιείÏ[[el

Re: [R] Linking a library with init(argc,argv) function

2008-12-10 Thread Dirk Eddelbuettel
On Wed, Dec 10, 2008 at 09:12:12PM +0100, Maxime Debon wrote: > I am facing a difficulty with loading a library depending on a library > (libMesh) which needs an initialization function (libMesh::init(argc, > argv);). Create 'fake' argc/argv as eg in char *myargv[] = { "somename", "--foo

[R] Confusion with Converting Factors to Dates using as.date

2008-12-10 Thread Josip Dasovic
Dear R-Helpers: I'm having a problem getting dates into the correct format. I have a data frame, which is based on a .csv file that I imported into R via read.table. R has converted my date variables to factors; when I use the as.Date command, most of the values are converted "correctly" (and b

Re: [R] How to Break Axis in Lattice Plot

2008-12-10 Thread Deepayan Sarkar
On Tue, Dec 9, 2008 at 8:04 PM, Larry Ma <[EMAIL PROTECTED]> wrote: > Dear All, > > I have searched the following problem in the R help achives, and there did > not seem to be solutio for it. > > The problem is how to break the axis in the xyplot (lattice plot). I read > the posts using the packag

Re: [R] Better way to find distances between points in a set?

2008-12-10 Thread Stavros Macrakis
On Tue, Dec 9, 2008 at 11:52 PM, Charles C. Berry <[EMAIL PROTECTED]>wrote: > The 'better way' to do almost anything starts with a reading of the > _posting guide_, which reminds you to >Do your homework before posting [Reasons whyfor deleted]]... > > Oh yes, if you are too lazy to look

[R] trouble loading candisc

2008-12-10 Thread Pete Shepard
Hello, I am having trouble loading the package candisc onto my R distribution. I am using 2.7.1-2. I do a "> install.packages("candisc" and get the following output. Warning in install.packages("candisc") : argument 'lib' is missing: using '/usr/local/lib/R/site-library' --- Please select a CRA

Re: [R] Confusion with Converting Factors to Dates using as.date

2008-12-10 Thread Peter Dalgaard
Josip Dasovic wrote: Dear R-Helpers: I'm having a problem getting dates into the correct format. I have a data frame, which is based on a .csv file that I imported into R via read.table. R has converted my date variables to factors; when I use the as.Date command, most of the values are convert

Re: [R] Confusion with Converting Factors to Dates using as.date

2008-12-10 Thread Marc Schwartz
on 12/10/2008 02:41 PM Josip Dasovic wrote: > Dear R-Helpers: > > I'm having a problem getting dates into the correct format. I have a > data frame, which is based on a .csv file that I imported into R via > read.table. > > R has converted my date variables to factors; when I use the as.Date > co

Re: [R] Confusion with Converting Factors to Dates using as.date

2008-12-10 Thread Josip Dasovic
Thank you very much, Peter. As is often the case, R gave me exactly what I asked it to give me, but not what I wanted it to give me. :) Cheers, Josip Research Associate Human Security Report Project School for International Studies Simon Fraser University Suite 7200--515 W. Hastings St. Vancouver

Re: [R] missing argument

2008-12-10 Thread Tony Breyal
?missing never used it myself, but looks like it might help you :-) Tony Breyal. On 10 Dec, 19:09, "jonas garcia" <[EMAIL PROTECTED]> wrote: > Dear list, > I have a question and I'm going to give an example of my problem > > f<- function(d1, d2, d3) > { > d<- d1*d2/d3 > return(d) > > } > > v1<-

Re: [R] missing argument

2008-12-10 Thread Tony Breyal
I think i might of misunderstood the question. I was thinking you meant what if one of the arguments were not used, as in the example below (see ?missing): f<- function(d1, d2, d3) { vec <- c(missing(d1), missing(d2), missing(d3)) if(any(vec)) { return(0) }

Re: [R] R and Scheme

2008-12-10 Thread Wacek Kusnierczyk
Peter Dalgaard wrote: > Johannes Huesing wrote: >> Stavros Macrakis <[EMAIL PROTECTED]> [Wed, Dec 10, 2008 at >> 04:59:25AM CET]: >>> So I conclude that what is really meant by "R semantics are based on >>> Scheme >>> semantics" is "R has functions as first-class citizens and a correct >>> implemen

[R] repeated searching of no-missing values

2008-12-10 Thread Patrizio Frederic
hi all, I have a data frame such as: 1 blue 0.3 1 NA0.4 1 red NA 2 blue NA 2 green NA 2 blue NA 3 red 0.5 3 blue NA 3 NA1.1 I wish to find the last non-missing value in every 3ple: ie I want a 3 by 3 data.frame such as: 1 red 0.4 2 blue NA 3 blue 1.1 I have written a little

Re: [R] missing argument

2008-12-10 Thread Bert Gunter
I believe the usual practice in this case is simply to give default values for arguments: function(x, y, opt.arg1 = 0, opt.arg2 = sin(1),...) If you haven't already done so, perusal of "An Introduction to R" -- especially the "Named Arguments and Defaults" -- would be appropriate. Cheers, Bert G

Re: [R] Better way to find distances between points in a set?

2008-12-10 Thread Carl Witthoft
Ok... then.. I didn't mean to start a flame war, so first off I apologize for firing off my second comment. I think a fair compromise is: I should have at least noted that I got as far as finding dist() (and which.max()), thus showing "some" effort on my part before posting my query. Anyw

Re: [R] Better way to find distances between points in a set?

2008-12-10 Thread Charles C. Berry
On Wed, 10 Dec 2008, Stavros Macrakis wrote: On Tue, Dec 9, 2008 at 11:52 PM, Charles C. Berry <[EMAIL PROTECTED]>wrote: The 'better way' to do almost anything starts with a reading of the _posting guide_, which reminds you to Do your homework before posting [Reasons whyfor deleted]]...

Re: [R] read.*: How to read from a URL?

2008-12-10 Thread Prof Brian Ripley
On Wed, 10 Dec 2008, hadley wickham wrote: Hi Michael, In general, I think you should be able to do: gimage <- read.jpeg(url(gimageloc)) Note that would not be really correct: it would need to be gimage <- read.jpeg(con <- url(gimageloc)) close(con) since it otherwise leaks a connection (w

Re: [R] read.*: How to read from a URL?

2008-12-10 Thread Martin Morgan
Prof Brian Ripley wrote: On Wed, 10 Dec 2008, hadley wickham wrote: Hi Michael, In general, I think you should be able to do: gimage <- read.jpeg(url(gimageloc)) Note that would not be really correct: it would need to be gimage <- read.jpeg(con <- url(gimageloc)) close(con) since it other

[R] How do I multiply labeled vectors of numbers?

2008-12-10 Thread Bill McNeill (UW)
I'm an R newbie trying to do an operation that seems like it should be very simple, but after much playing around with the interface and reading "An Introduction to R" I can figure out how to do it. I have two text files. File1 contains: typen dog2 cat4 horse6 File2 cont

Re: [R] How do I multiply labeled vectors of numbers?

2008-12-10 Thread Jim Holtman
?merge Sent from my iPhone On Dec 10, 2008, at 15:17, "Bill McNeill (UW)" <[EMAIL PROTECTED]> wrote: I'm an R newbie trying to do an operation that seems like it should be very simple, but after much playing around with the interface and reading "An Introduction to R" I can figure out h

[R] bigmemory with reshape?

2008-12-10 Thread Avram Aelony
Hello, I am running into memory boundaries and would like to try to make use of the bigmemory (or any other memory enabling) library. Can anyone help with suggestions as to how this might work? > library(reshape) > s <- melt( d[,1:62], id=c(1) ) Error: cannot allocate vector of size 16.0 Mb >

  1   2   >