Re: [R] Likelihood ratio test using R

For the third one: ?anova.glm test=Chisq will be LRT. For the first two, you can have the answer from ordinary stat book. On Wed, Nov 5, 2008 at 1:11 PM, Maithili Shiva <[EMAIL PROTECTED]> wrote: > Hi! > > I am working on the Logistic Regression using R. My R script is as follows > > > ONS <-

Re: [R] contour plot, failing to interpolate through all data

Folkes, Michael dfo-mpo.gc.ca> writes: > I'm having trouble making contour lines for this attached, sparse > dataset (low data:NA ratio!). Is it the high number of NA's, or funny > layout of the densities, or something else that's causing this? > > start data to be in file='testmatrix.csv

Re: [R] How to extract following data

RON70 yahoo.com> writes: > > - > 2005-01-17T00:00:00+05:30 > 10149 > 1288.40002 > Looks suspiciously like XML, and let's hope the real data are more like this below, without the "-" and with a nice header 2005-01-17T00:00:00+05:30 10149 1288.40002 2005-01-18T0

[R] How to extract following data

Hi everyone, I have this kind of raw dataset : - 2005-01-17T00:00:00+05:30 10149 1288.40002 - 2005-01-18T00:00:00+05:30 10149 1291.69995 - 2005-01-19T00:00:00+05:30 10149 1288.19995 I was looking for some R procedure to extract data from this, that should

[R] Likelihood ratio test using R

Hi! I am working on the Logistic Regression using R. My R script is as follows ONS <- read.csv("ONS.csv",header = TRUE) ONS.glm <- glm(Y ~ Age1+Age2+Sex+Education+Profession+TimeInJob+ TimeInCurrentJob+OwnResidence+Dependents+ApplIncome+FamilyInco+IIR+FOIR+YearsAtBank+SavingsAccount+CurrentAcc

Re: [R] How to create margin area around a graphic made by xyplot

Hi Dr. Alireza Zolfaghari wrote: > Hi list, > Does any body know how I can create a margin area around the graphic created > by xyplot (the same thing we do using par in plot function)? I tried > viewport, but does not work. The trick is to assign the xyplot() result and then print it explicitl

[R] FW: need some help

Hi , > df Session_Setup DCT FwdDataVols_bin countsComp 1 User_Initiated NoRLL 1 5058 User_Initiated+NoRLL+1 2 User_Initiated NoRLL 2584 User_Initiated+NoRLL+2 3 User_Initiated NoRLL 3191 User_Initi

Re: [R] use abline() for regression model in the plot

Hey, it works exactly!thx man. But i have to write comment on what is the command i am using line by line, so i was wondering can you explain the command you use in words for me? many many thanks. -- View this message in context: http://www.nabble.com/use-abline%28%29-for-regression-model-in-th

Re: [R] How to generate a new factor variable by two other factor variables

Thanks. It works :) Shuguang On Wed, Nov 5, 2008 at 12:24 AM, Marc Schwartz <[EMAIL PROTECTED]> wrote: > The easiest way is probably to use interaction(): > >> df > factorA factorB > 1 0 0 > 2 0 0 > 3 1 0 > 4 0 1 > 5 1 1 > > > # Note th

Re: [R] merge 2 txt file

What exactly are you trying to do? You need to provide at least a sample of the input and then what you would expect as output. Exactly what to the columns contain and how you expect to 'merge' them. On Tue, Nov 4, 2008 at 6:52 PM, Alessandro <[EMAIL PROTECTED]> wrote: > Hi all, > > > > I have t

Re: [R] use abline() for regression model in the plot

On Wed, Nov 5, 2008 at 6:47 AM, Swanton0822 <[EMAIL PROTECTED]> wrote: > > after that i have been ask to do the following: > 'The regression equation for EXPEND depending on ECAB is > E(EXPEND)=119.0+1.73xECAB when there is outliners and > E(EXPEND)=97.8+1.94xECAB when the outliners are omitted fro

[R] contour plot, failing to interpolate through all data

I'm having trouble making contour lines for this attached, sparse dataset (low data:NA ratio!). Is it the high number of NA's, or funny layout of the densities, or something else that's causing this? This is a subsample of the dataset, and I get the same problem when using the full data. The only

[R] parcoord

Hi All, Could you please advise me how to put xlab for var1,var2,var3,.. in parcoord()? With regards, Abu _ Win £1000 John Lewis shopping sprees with BigSnapSearch.com [[alternative HTML version deleted]] _

[R] merge 2 txt file

Hi all, I have two txt file with X,Y,Z column and I need to merge together I tried file_all <- merge("file1.txt","file2.txt") but I don't sure about the result. Is It this code correct? Thanks Ale [[alternative HTML version deleted]] _

[R] perform Kruskal-Wallis test without using the built-in command in R

Hi, again i am stuck in my presentation, and i have never learn R before in my life but need this to be done, so please help me out for a favour: http://www.nabble.com/file/p20333155/kew.dat kew.dat run this in R and these comes up: Month Year Rain 1 Jan 1900 74.40 2

Re: [R] parcoord

Dear Abu, One option would be use title() as follows: library(MASS) parcoord(state.x77[, c(7, 4, 6, 2, 5, 3)]) title(xlab='Your title here') HTH, Jorge On Tue, Nov 4, 2008 at 6:07 PM, Abu Naser <[EMAIL PROTECTED]> wrote: > > Hi All, > > Could you please advise me how to put xlab for var1,var

[R] use abline() for regression model in the plot

http://www.nabble.com/file/p20332968/spe.dat spe.dat Hi, i have a problem in sorting out some command in R and i am really hoping some expert can help me out please! i have the spe.dat file which i upload here, and when u read into R u got something like this: A5 <- read.table('spe.dat' ,header

[R] neural package

Hello, I am very happy to see that the Rprop neural network training algorithm is now available in the neuralnet package. I have a few questions about the package and a few feature requests also. According to the documentation the startweights argument of the neuralnet functio

Re: [R] Problem: Sweave and math environments

David Scott wrote: On Tue, 4 Nov 2008, Duncan Murdoch wrote: On 11/4/2008 9:15 AM, Sarah Goslee wrote: I'm sure there are a variety of ways around this. My chosen solution is to replace mydata$somevar with mydata[, "somevar"] Sarah You should almost never need to modify valid

Re: [R] [OT] factorial design

On 4/11/2008, at 6:29 PM, Edna Bell wrote: Dear R Gurus: I vaguely remember reading that if interaction was present in a factorial design, then the main effect results were suspect. However, I was reading a text which now uses the tests for main effects even if interaction is present. Which

Re: [R] AIC in time series

Vasileios Ismyrlis gmail.com> writes: > > library("fGarch") > fit = garchFit(formula=~ar(1)+~garch(1,1), data=garat) > > Now I want to count AIC for this model. How can I do it? I cannot do > it with the AIC function of stats package, because R tells me: > "Error in UseMethod("logLik") : no ap

Re: [R] Problem: Sweave and math environments

On Wed, 5 Nov 2008, David Scott wrote: On Tue, 4 Nov 2008, Duncan Murdoch wrote: On 11/4/2008 9:15 AM, Sarah Goslee wrote: > I'm sure there are a variety of ways around this. My chosen solution > is to replace > > mydata$somevar > with > mydata[, "somevar"] > > Sarah > You should

Re: [R] Problems with scan

David C. Howell wrote: I have been having problems with using scan(). I searched the archives and found someone with the same problem several years ago, but did not find a solution. I want to prompt the user for input, scan in that input, and then go on to use it as a variable. The simplest

Re: [R] Problem: Sweave and math environments

On Tue, 4 Nov 2008, Duncan Murdoch wrote: On 11/4/2008 9:15 AM, Sarah Goslee wrote: I'm sure there are a variety of ways around this. My chosen solution is to replace mydata$somevar with mydata[, "somevar"] Sarah You should almost never need to modify valid R code to work with Sweave. (Th

Re: [R] response surfaces

thank you for all of your help! On Tue, Nov 4, 2008 at 3:13 PM, Bert Gunter <[EMAIL PROTECTED]> wrote: > Take a look at the AlgDesign package. > > Of course, lm(), or whatever be your favorite version of multiple > regression, will fit the results. > > Obviously (I think) you'll have to write you

Re: [R] response surfaces

Take a look at the AlgDesign package. Of course, lm(), or whatever be your favorite version of multiple regression, will fit the results. Obviously (I think) you'll have to write your own plotting code or perhaps make use of anything available in the other packages Ben mentioned. -- Bert Gunter

[R] Problems with scan

I have been having problems with using scan(). I searched the archives and found someone with the same problem several years ago, but did not find a solution. I want to prompt the user for input, scan in that input, and then go on to use it as a variable. The simplest version of my program i

[R] AIC in time series

Hi everybody, I have fitted an ar(1),Garch(1,1) model to some observations with the help of the garchFit function which is in the fGarch package. Here what I've done: library("fGarch") fit = garchFit(formula=~ar(1)+~garch(1,1), data=garat) Now I want to count AIC for this model. How can I do it?

Re: [R] ordered logistic regression of survey data with missing variables

You can analyse multiple imputations with the survey and mitools package, and there is a toy example including ordinal logistic regression at http://faculty.washington.edu/tlumley/survey/svymi.html If am I reading their documentation correctly, 'mice' creates what Rubin calls 'proper imputa

Re: [R] Building with MKL on Ubuntu

I can see a couple of problems: 1) No "-lpthread" in --with-blas 2) No "-L" prefix in the library path in --with-lapack In addition, I don't think you need to add -lmkl to --with-lapack, although that is probably harmless. Martyn Quoting Prof Brian Ripley <[EMAIL PROTECTED]>: > Look in config.

Re: [R] response surfaces

Erin Hodgess gmail.com> writes: > > Dear R People: > > Are there functions or packages for response surfaces, please? > RSiteSearch("response surface") does bring up a lot of hits, although most of them are from people asking questions about plotting 3D response surfaces in lattice etc..

Re: [R] How to suppress errors from htmlTreeParse() function in XML package?

Martin Morgan wrote: Hi Tony -- Tony Breyal <[EMAIL PROTECTED]> writes: Dear R-help, The following code downloads an html document into variable 'doc' and then stores an internal representation into variable 'html.tree'. Even if the html code is malformed, this still works which is fantasti

[R] Extract Residual Variance from lmer using VarCorr

I am trying to extract the residual variance from an lmer object. I have tried using VarCorr() which allows me to see the estimate but can not save it by itself. How can I get this from VarCorr() or is there another function available for this? Thanks, ben __

[R] Courses*** R/Splus Fundamentals & R Advanced Programming: November - December 2008 at 3 locations by XLSolutions Corp

XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce our*** R/Splus Fundamentals and Programming Techniques and R Advanced Programming***courses at 3 USA locations for November - December 2008. (1) R/Splus Fundamentals and Programming Techniques http://www.xlsolutions-corp.co

Re: [R] IWLS vs direct ML estimation

sandsky <[EMAIL PROTECTED]> wrote: > I think that IWLS provides the computational flexibility. Because when there > exists no analytical solution, IWLS makes it possible to solve via the > numerical solution. > > Do you have any idea for this? My impression is the same as yours. M.H.P. __

Re: [R] Zoo seems to be running slow in R 2.8.0 windows

2.3-24 chron and 1.5-4 zoo in both 2.8 and 2.7.2 On Tue, Nov 4, 2008 at 1:43 PM, Gabor Grothendieck <[EMAIL PROTECTED]> wrote: > One other question. Are you running the same version of the chron > package on both 2.7.2 and 2.8.0? > > packageDescription("chron")$Version > > will tell you. If the

Re: [R] Zoo seems to be running slow in R 2.8.0 windows

One other question. Are you running the same version of the chron package on both 2.7.2 and 2.8.0? packageDescription("chron")$Version will tell you. If there is a version difference does changing the chron version change the performance? On Tue, Nov 4, 2008 at 1:37 PM, Gabor Grothendieck <[E

Re: [R] Prevent read.table from converting "+" and "-" to 0

Daren Tan hotmail.com> writes: > I am using read.table("data.txt", sep="\t") to read in a tab-limited text file. However, two columns of > data were read wrongly. read.table converts "+" and "-" in the two columns to 0. I have tried setting other > parameters but to no avail. Looks like http:

Re: [R] Zoo seems to be running slow in R 2.8.0 windows

On Tue, Nov 4, 2008 at 1:06 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > this should be right- I have looked at it with summaryRprof() and it > is spending most of its time on things I would expect... > read.production and read.zoo. Is this the best way to communicate the It indicates that chr

Re: [R] Building with MKL on Ubuntu

Look in config.log to see what's wrong. (E.g. is /opt/intel/mkl/10.0.2.018/lib/em64t in the ld.so cache?) And note the warnings in the manual about using --with-lapack: it is most definitely not recommended. R-devel would be a better place to ask questions about this. On Tue, 4 Nov 2008, An

Re: [R] how to cite R data sets

The data sets should be documented... library(MASS) ?cats and then the source should tell you how to site it and specifically ?airquality On Tue, Nov 4, 2008 at 11:03 AM, dxc13 <[EMAIL PROTECTED]> wrote: > > Hi, > > As part of a project I am doing, I am using the "airquality" data set that > c

Re: [R] Zoo seems to be running slow in R 2.8.0 windows

this should be right- I have looked at it with summaryRprof() and it is spending most of its time on things I would expect... read.production and read.zoo. Is this the best way to communicate the results? Or do you want it copy and pasted into the email? thank you very much On Tue, Nov 4, 2008 a

[R] how to cite R data sets

Hi, As part of a project I am doing, I am using the "airquality" data set that comes built into R as a means of showing an example of how my analysis will be carried out. I know that citation("package") will produce a citation for any R package, but is there a proper way to cite this data set?

Re: [R] IWLS vs direct ML estimation

Mike, Thank you for your valuable reply. I have read Charnes et al. (1976) and Bradley (1973). However, they just showed the equivalence of IWLS and ML solutions, but didn't mentioned the advantage (or disadvantage) of IWLS comparing with ML estimation. Charnes, A., Frome, E. L., and Yu, P. L.,

Re: [R] Zoo seems to be running slow in R 2.8.0 windows

Here is the output: > summaryRprof() $by.self self.time self.pct total.time total.pct index.search 0.30 55.6 0.30 55.6 .C 0.18 33.3 0.18 33.3 inDL 0.06 11.1

[R] Building with MKL on Ubuntu

Hi all, I'm trying to build R from subversion with MKL 10.0.2 on Ubuntu. I tried: ./configure --with-blas='-L/opt/intel/mkl/10.0.2.018/lib/em64t -lmkl -lguide' --with-lapack='/opt/intel/mkl/10.0.2.018/lib/em64t -lmkl -lmkl_lapack' --enable-R-shlib and got: R is now configured for x86_64-unknown-

Re: [R] Zoo seems to be running slow in R 2.8.0 windows

Are you sure? summaryRprof() says that based on your Rprof.out file that 55% of the time is being spent in index.search which is for searching help files. On Tue, Nov 4, 2008 at 12:36 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > this is from the read.production command > > On Tue, Nov 4, 200

[R] response surfaces

Dear R People: Are there functions or packages for response surfaces, please? Thanks in advance, Sincerely, Erin -- Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailto: [EMAIL PROTECTED]

Re: [R] ESS Toolbar missing after Ubuntu Update

Andreas, I am late with this but in case it might help, I installed Ubuntu 8.10 this weekend and used "apt-get install emacs" and "apt-get install ess" to place Emacs version 22.2.1 and ess 5.3.9. Simple tests seemed to indicate that menus come up, R files load, and code could be executed, as expe

[R] Clustering In R. (rookie)

Hi all. I have alrge microarray dat set that i would like to analyze using hierarchical clustering. The problem is when i use the command below, > hc<- hclust(dist(array), "ave") i get get this feedback... Error in as.vector(x, mode) : cannot coerce type 'closure' to vector of type 'any' Can som

Re: [R] How to suppress errors from htmlTreeParse() function in XML package?

Hi Tony -- Tony Breyal <[EMAIL PROTECTED]> writes: > Dear R-help, > > The following code downloads an html document into variable 'doc' and > then stores an internal representation into variable 'html.tree'. Even > if the html code is malformed, this still works which is fantastic. > However, as

Re: [R] How to generate a new factor variable by two other factor variables

... and also note the "drop" argument. Adding "drop = TRUE" to the call is usually what is wanted, e.g. df$factorC <- with(df, interaction(factorA, factorB), drop = TRUE ) -- Bert Gunter -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Marc Schwartz Se

Re: [R] Zoo seems to be running slow in R 2.8.0 windows

See ?Rprof On Tue, Nov 4, 2008 at 12:01 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > R version 2.8.0 (2008-10-20) > i386-pc-mingw32 > > locale: > LC_COLLATE=English_United States.1252;LC_CTYPE=English_United > States.1252;LC_MONETARY=English_United > States.1252;LC_NUMERIC=C;LC_TIME=English_Uni

Re: [R] Zoo seems to be running slow in R 2.8.0 windows

as a follow up I am using the same things as above in 2.7.2 and it works just fine... On Tue, Nov 4, 2008 at 12:01 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > R version 2.8.0 (2008-10-20) > i386-pc-mingw32 > > locale: > LC_COLLATE=English_United States.1252;LC_CTYPE=English_United > States.125

[R] Zoo seems to be running slow in R 2.8.0 windows

R version 2.8.0 (2008-10-20) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 attached base packages: [1] stats graphics grDevices utils datasets metho

Re: [R] How to generate a new factor variable by two other factor variables

The easiest way is probably to use interaction(): > df factorA factorB 1 0 0 2 0 0 3 1 0 4 0 1 5 1 1 # Note the default separator of '.' df$factorC <- with(df, interaction(factorA, factorB)) > df factorA factorB factorC 1 0

Re: [R] getting "small" graphs with Sweave using Rnews.sty

On 11/4/2008 11:15 AM, Felipe Carrillo wrote: \setkeys{Gin}{width=8} gets the desired width but shifts the graph to the right. Is there another latex command to place the graph so that the margins on left and right of the page are the same? \centering or \begin{center} doesn't seem to get it right

Re: [R] getting "small" graphs with Sweave using Rnews.sty

\setkeys{Gin}{width=8} gets the desired width but shifts the graph to the right. Is there another latex command to place the graph so that the margins on left and right of the page are the same? \centering or \begin{center} doesn't seem to get it right. I have about one inch margin to my right

Re: [R] How to generate a new factor variable by two other factor variables

Dear Shuguang, Here are two ways. Perhaps they are not efficient enough, but the work: # Data mydata=read.table(textConnection(" factorA factorB 0 0 0 0 1 0 0 1 1 1"),header=TRUE) closeAllConnections() # Option 1 mydata$factorC=as.factor( apply(mydata

Re: [R] IWLS vs direct ML estimation

sandsky <[EMAIL PROTECTED]> wrote: > I am thinking about IWLS vs ML estimation. When I use glm() for a > 2-parameter distribution (e.g., Weibull), I can otain the MLE of scale > parameter given shape parameter through IWLS. Because this scale parameter > usually converges to the MLE. > > In this

Re: [R] ggplot & annotating charts

Note that even without substantial changes the code could be shortened by making use of read.zoo and other items illustrated below: library(zoo) vix.url <- "http://www.cboe.com/publish/ScheduledTask/MktData/datahouse/vixcurrent.csv"; v <- read.zoo(vix.url, header = TRUE, sep = ",", format = "%m/

Re: [R] writing R extensions

If you are going to contribute it then all you need to do is upload it to the ftp site found in writing R extensions and very nice person will oversee the building of windows and mac binaries. that is as far as my knowledge goes. good luck stephen On Tue, Nov 4, 2008 at 10:14 AM, micha_ <[EMAIL

Re: [R] writing R extensions

It sounds like you simply uncompressed your .tar.gz file and then zipped it up. If so, it should not work correctly. You need to compile it for windows. Try something like Rcmd build --binary myRpackageDir and you may need to include "--force" option in the command above. Also

Re: [R] writing R extensions

On 11/4/2008 10:14 AM, micha_ wrote: Hi, I'm working on a package and got some problems. After I've done R CMD check and build I get the package.tar.gz which I can install under Linux without any problems. Now I wanted to have a Windows version. I heard that I only have to zip the package folder

[R] writing R extensions

Hi, I'm working on a package and got some problems. After I've done R CMD check and build I get the package.tar.gz which I can install under Linux without any problems. Now I wanted to have a Windows version. I heard that I only have to zip the package folder. That worked once, but now the packag

[R] sample selection mixed-effects model with nlme

Dear all, I want to estimate this equation with repeated data. log(y)=b0+b1*log(x1)+b2*log(x1)+B'Z+u if y<=0 then log(y) is undefined.  It seem to be a tobit model with panel data, can I do this in R with nlme, lme or plm ? If not is it valid to use a two step procedure ? Justin BEM BP

[R] How to generate a new factor variable by two other factor variables

How to generate a new factor variable by two other factor variables? For example, if I have two factor variables, factorA and factorB, factorA factorB 0 0 0 0 1 0 0 1 1 1 Is there a simple way to generate a new 4-levels factor variable as factorC

Re: [R] TIme Series AR to MA and (viceversa)

On Tue, 4 Nov 2008, [EMAIL PROTECTED] wrote: Hi, I am new to using R for Time series analysis. I was wondering if there are any functions that can convert ARMA or ARIMA time series into their corresponding AR or MA time series representations (by calculating the corresponding AR or MA coeffi

Re: [R] TIme Series AR to MA and (viceversa)

Sounds like a homework problem... In general, the AR or MA representation of an ARMA process will be of infinite order. GP > Date: Tue, 04 Nov 2008 07:10:03 -0800 (PST) > From: "[EMAIL PROTECTED]" <[EMAIL PROTECTED]> > Sender: [EMAIL PROTECTED] > Precedence: list > DomainKey-Signature: a=rsa-s

Re: [R] Prevent read.table from converting "+" and "-" to 0

On Tue, 4 Nov 2008, Duncan Murdoch wrote: On 11/4/2008 9:08 AM, Daren Tan wrote: I am using read.table("data.txt", sep="\t") to read in a tab-limited text file. However, two columns of data were read wrongly. read.table converts "+" and "-" in the two columns to 0. I have tried setting other p

[R] TIme Series AR to MA and (viceversa)

Hi, I am new to using R for Time series analysis. I was wondering if there are any functions that can convert ARMA or ARIMA time series into their corresponding AR or MA time series representations (by calculating the corresponding AR or MA coefficients). Thanks a lot Kris.

Re: [R] ggplot & annotating charts

Hi Simeon, Could you publish vixarchive.csv too please? Do you have a start on the version in ggplot2? It would be useful to see that code too. Hadley On Tue, Nov 4, 2008 at 8:49 AM, simeon duckworth <[EMAIL PROTECTED]> wrote: > Dear "R-listers" > > I've been trying to figure out how to annotat

Re: [R] ggplot2 scatterplot matrix

#this worked (I changed the dataframe in this code so it is not tested...but should work) plotmatrix(d[,3:6], mapping=aes(shape=d$Site))+scale_shape_manual(values=c(1:12)) On Tue, Nov 4, 2008 at 9:51 AM, hadley wickham <[EMAIL PROTECTED]> wrote: >> plotmatrix(d[,3:6]) >> >> I would like to either

Re: [R] reshape bug?

I've figured out the problem and just pushed a new version up to CRAN. It will (hopefully) be available from mirrors in a couple of days time. Hadley On Tue, Nov 4, 2008 at 8:29 AM, stephen sefick <[EMAIL PROTECTED]> wrote: > no problem- keep up the good work. > > On Tue, Nov 4, 2008 at 9:25 AM,

Re: [R] ggplot2 scatterplot matrix

> plotmatrix(d[,3:6]) > > I would like to either have a colour or more preferably a shape to > correspond to Site > this is what I have tried (along with many other combinations) > plotmatrix(nmds.bug[,3:6])+scale_shape(nmds.bug, shape=Site, pch=1:12) There's not really any easy way to do this in

[R] ggplot & annotating charts

Dear "R-listers" I've been trying to figure out how to annotate charts in ggplot (ie add text to line charts, highlighted boxes etc). By and large, I can get close to what i want with base graphics, but would ideally like to use ggplot whenever possible (additionally, i would like to add text lab

Re: [R] reshape bug?

Hi Stephen, Thanks for the bug report (and apologies for the delay getting back to you). I'm working on a fix today and will hopefully release a new version in the very near future. Hadley On Fri, Oct 31, 2008 at 12:51 PM, stephen sefick <[EMAIL PROTECTED]> wrote: > Hadley et al., > I was using

Re: [R] Problem: Sweave and math environments

On 11/4/2008 9:15 AM, Sarah Goslee wrote: I'm sure there are a variety of ways around this. My chosen solution is to replace mydata$somevar with mydata[, "somevar"] Sarah You should almost never need to modify valid R code to work with Sweave. (There are a few exceptions: code that illus

[R] ggplot2 scatterplot matrix

#question is below data library(ggplot2) d <- (structure(list(Site = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 4L, 4L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 7L, 7L, 8L, 8L, 8L, 8L, 8L, 8L, 9L, 9L, 9L, 9L, 9L, 9L, 1

Re: [R] reshape bug?

no problem- keep up the good work. On Tue, Nov 4, 2008 at 9:25 AM, hadley wickham <[EMAIL PROTECTED]> wrote: > Hi Stephen, > > Thanks for the bug report (and apologies for the delay getting back to > you). I'm working on a fix today and will hopefully release a new > version in the very near futu

Re: [R] Prevent read.table from converting "+" and "-" to 0

On 11/4/2008 9:08 AM, Daren Tan wrote: I am using read.table("data.txt", sep="\t") to read in a tab-limited text file. However, two columns of data were read wrongly. read.table converts "+" and "-" in the two columns to 0. I have tried setting other parameters but to no avail. Set the column c

Re: [R] fine grain tick marks for zoo plots

Many thanks all. The following does the trick for me, taken out of the vignette: plotmonths<-function(z,...){ plot(z,xaxt="n",...) tt <- time(z) m <- unique(as.Date(as.yearmon(tt))) jan <- format(m, "%m") == "01" mlab <- substr(months(m[!jan]), 1, 1) Axis(side = 1, at = m[!jan], labels = ml

Re: [R] Problem: Sweave and math environments

I'm sure there are a variety of ways around this. My chosen solution is to replace mydata$somevar with mydata[, "somevar"] Sarah On Tue, Nov 4, 2008 at 6:16 AM, Christian Reisinger <[EMAIL PROTECTED]> wrote: > Dear list, ... > I want Sweave to "echo" my code, so that the R source code of the cod

Re: [R] Problem: Sweave and math environments

On 11/4/2008 6:16 AM, Christian Reisinger wrote: Dear list, I using R as the statistical package for my dissertation in political sciences and I am using Sweave for documenting the work done for the creation of my dataset. However, I have a problem with Sweave (which is great, by the way) whi

[R] HDF5 and R 2.8.0patched

Hi Everyone, I have updated my R version from 2.7 to 2.8 (patch) and since then I can not open HDF5 files that I saved previously with the command hdf5save from the hdf5 library. As anyone had that problem yet? Any help would be appreciated. Thanks Emmanuel

[R] Prevent read.table from converting "+" and "-" to 0

I am using read.table("data.txt", sep="\t") to read in a tab-limited text file. However, two columns of data were read wrongly. read.table converts "+" and "-" in the two columns to 0. I have tried setting other parameters but to no avail. TIA __

Re: [R] R News, volume 8, issue 2 is now available

Dear r-help list members, A couple of people have requested the direct link to the new issue of R News. It is . John -- John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario, Canada web:

Re: [R] fine grain tick marks for zoo plots

And there are additional examples in vignette("zoo-faq") and example(xyplot.zoo) On Tue, Nov 4, 2008 at 8:26 AM, Gabor Grothendieck <[EMAIL PROTECTED]> wrote: > example(plot.zoo) has an example. > > On Tue, Nov 4, 2008 at 8:10 AM, <[EMAIL PROTECTED]> wrote: >> Dear R Users, >> >> I am trying to g

Re: [R] fine grain tick marks for zoo plots

I think it is one of the fancy x-axis examples. If you have some trouble I may be able to dig up some code, but no promises. Good Luck On Tue, Nov 4, 2008 at 8:26 AM, Gabor Grothendieck <[EMAIL PROTECTED]> wrote: > example(plot.zoo) has an example. > > On Tue, Nov 4, 2008 at 8:10 AM, <[EMAIL PRO

[R] Problem: Sweave and math environments

Dear list, I using R as the statistical package for my dissertation in political sciences and I am using Sweave for documenting the work done for the creation of my dataset. However, I have a problem with Sweave (which is great, by the way) which is (hopefully) easy to solve: I want Sweave t

[R] How to suppress errors from htmlTreeParse() function in XML package?

Dear R-help, The following code downloads an html document into variable 'doc' and then stores an internal representation into variable 'html.tree'. Even if the html code is malformed, this still works which is fantastic. However, as in the example below, i do get some ouput from R in the console

[R] R and real time applications

Is R currently used in real time applications (large scale)? What is the typical architecture, say for real time transaction scoring? Where is the best place to ask this type of questions? I've also asked the question on AnalyticBridge (3000 members), check the R-Project group at www.analyticb

Re: [R] fine grain tick marks for zoo plots

example(plot.zoo) has an example. On Tue, Nov 4, 2008 at 8:10 AM, <[EMAIL PROTECTED]> wrote: > Dear R Users, > > I am trying to get plot.zoo to place monthy tickmarks/labels for a time > series which spans daily data going back a bit over a year. Right now, I > am getting only one tick mark on th

[R] fine grain tick marks for zoo plots

Dear R Users, I am trying to get plot.zoo to place monthy tickmarks/labels for a time series which spans daily data going back a bit over a year. Right now, I am getting only one tick mark on the x-axis for the beginning of 2008. How can I force plot.zoo to place more regular x-axis tick marks

Re: [R] Help needed using 3rd party C library/functions from within R (Nvidia CUDA)

Please see the R posting guide (the footer of this and every R-help message). - This is not an R-help question but an R-devel one. - You have not told us your OS, nor any of the 'at a minimum' information requested in the posting guide. But look at nm -g myFunc.so to see what symbols are expo

[R] Fill vector filed NAs data

Good afternoon, I am dealing with a 2d vector field on a regular grid. The field has sparse NA values (t near the domain borders and not isolated), which I am trying to substitute with interpolated values of the nearest point. What I am looking for is to replace the val

Re: [R] descretizing xy data

Rubén Roa-Ureta wrote: Jon A wrote: Hello, I have a dataset with a continuous independent variable (fish length, range: 30-150 mm) and a binary response (foraging success, 0 or 1). I want to discretize fish length into 5 mm bins and give the proportion of individuals who successfully foraged

Re: [R] a "date" x-axis in a xyplot

Is your variable actually a date? R doesn't recognise dates automatically - you need to create them with as.Date or similar. Hadley On Tue, Nov 4, 2008 at 3:22 AM, <[EMAIL PROTECTED]> wrote: > > Dear all, > > In the x-axis of a xyplot I define dates ,for example 01-10-2007 instead of > numbers

[R] Help needed using 3rd party C library/functions from within R (Nvidia CUDA)

Hello, I'm trying to combine the parallel computing power available through NVIDIA CUDA (www.nvidia.com/cuda) from within R. CUDA is an extension to the C language, so I thought it would be possible to do this. If I have a C file with an empty function which includes a needed CUDA library (cutil

[R] R News, volume 8, issue 2 is now available

Dear R users, The October 2008 issue of 'R News' is now available on CRAN under the Documentation/Newsletter link. John (on behalf of the R News Editorial Board) -- John Fox, Professor Department of Sociology McMaster University Hamilton, Ontario, Canada web: socserv.

  1   2   >