Re: [R] Background color in a grid plot seems to interfere with plot()

2012-09-23 Thread Marius Hofert
Dear Jeff, Thanks for helping. I thought that gridBase would take care of using simple 'graphics' and grid plots (?) I tried grid.grill, but it's not so trivial to set it up such that the grid is drawn at the axis ticks. Cheers, Marius Jeff Newmiller writes: > Same reason. > > grid() is pa

[R] GGally diagonal axis: changing font size

2012-09-23 Thread Andrew Digby
ggpairs() from GGally is fantastic, but I'm stuck trying to change the diagonal labels and increase the size of the font in the diagonal axis plots. I thought using ggally_diagAxis(... labelSize) would work, but it doesn't appear to. I've tried the following three approaches, all of which don't

Re: [R] Background color in a grid plot seems to interfere with plot()

2012-09-23 Thread Marius Hofert
Dear Paul, Thank you for helping. This works great. I then tried to put in a grid (via grid()). Why does that fail? Cheers, Marius require(grid) require(gridBase) pdf(file="Rplot.pdf", width=8, height=8, onefile=FALSE) ## set up the grid layout plot.new() # start (empty) new page with 'gr

[R] Regression Problem

2012-09-23 Thread Vignesh Prajapati
Hello all, I am new to R, I am learning regression and logistic modeling with categorical predictor variables, when there is only one predictor categorical variable i can use as.numeric() but when more than two variable then what is solution? can anyone suggest me? Thanks vignesh

Re: [R] Vegan vegdist Jaccard matrix problem

2012-09-23 Thread Jari.Olsanen
kbannarm hotmail.com> writes: > > This is my script: > > Madr<-read.csv("MadRegSppSite.csv", sep=",", row.names=1) > rownames(Madr) <- Madr[,1] > Madr <- Madr[,-1] > MadrJaccard<-vegdist(Madr, method="jaccard",binary=FALSE) > MadrJaccard<-as.matrix(dist(MadrJaccard)) > write.table(MadrJaccard,

[R] Rows not common in dataframes

2012-09-23 Thread Chintanu
Hi, I have two dataframes (Dataframe_A, Dataframe_B) with the same no. of columns. The first column of both the dataframes contains unique names. I wish to have Dataframe_A with the rows that are NOT common to Dataframe_B. With merge (), it is possible to get the common rows or to merge rows, but

Re: [R] Confused by code?

2012-09-23 Thread Peter Alspach
Tena koe I think you probably meant: x[as.logical(z)] <- y[as.logical(z)] i.e., choosing those elements of × and y where z is 1 (TRUE as logical). Whereas what you have written: ×[z] <- y[z] references the 0th (by default indexing starts at 1 so this is empty (see ×[0]) and the first element

Re: [R] Return level plots

2012-09-23 Thread Roy Mendelssohn
Hi Michelle: On Sep 23, 2012, at 3:38 PM, MichelleNCSU wrote: > On Fri, Sep 21, 2012 at 3:45 PM, David Winsemius [via R] > wrote: >> >> On Sep 21, 2012, at 7:17 AM, MichelleNCSU wrote: >> >>> Hello, >>> >>> First of all, let me apologize that my statistics background is modest at >>> best. >

[R] Vegan vegdist Jaccard matrix problem

2012-09-23 Thread kbannarm
This is my script: Madr<-read.csv("MadRegSppSite.csv", sep=",", row.names=1) rownames(Madr) <- Madr[,1] Madr <- Madr[,-1] MadrJaccard<-vegdist(Madr, method="jaccard",binary=FALSE) MadrJaccard<-as.matrix(dist(MadrJaccard)) write.table(MadrJaccard,file="MadrJaccard.csv",sep=",",row.names=T) The in

[R] Confused by code?

2012-09-23 Thread Bazman76
x<-matrix(c(1,0,0,0,1,0,0,0,1),nrow=3) > y<-matrix(c(0,0,0,1,0,0,1,1,0),nrow=3) > z<-matrix(c(0,1,0,0,1,0,1,0,0),nrow=3) > x[z]<-y[z] The resultant matrix x is all zeros except for the last two diagonal cells which are 1's. While y is lower triangualr 0's with the remaining cells all ones. I real

Re: [R] Background color in a grid plot seems to interfere with plot()

2012-09-23 Thread Paul Murrell
Hi On 24/09/12 09:36, Marius Hofert wrote: Hi, Why does the upper left panel (in the plot below) not have a gray background? That is a symptom of the conflict that occurs when the 'graphics' package and the 'grid' package both try to initialise a new page. A good rule of thumb is to start a

[R] problems with function geese() in geepack

2012-09-23 Thread André Gabriel
Good evening, In my research I am studying the marginal models, where the main goal is on the structure of the association. My practical example has cluster with up to 600 observations and with this database, the function geese() return me the following message: This application has reque

[R] How to Write a Model in R that has the Log taken of the Intercept

2012-09-23 Thread Cheryl Johnson
Hi, I know that +1 is used to specify an intercept in a R model. An example of this would be: y~x+1 If I want to have a model where the log of the intercept is taken, the equation y~x+log(1) will not take the log of the intercept. Any suggestions on how to take the log of the intercept will be a

Re: [R] Return level plots

2012-09-23 Thread MichelleNCSU
On Fri, Sep 21, 2012 at 3:45 PM, David Winsemius [via R] wrote: > > On Sep 21, 2012, at 7:17 AM, MichelleNCSU wrote: > >> Hello, >> >> First of all, let me apologize that my statistics background is modest at >> best. >> >> I am doing some extreme value analysis on model output (WRF) which have >>

[R] Background color in a grid plot seems to interfere with plot()

2012-09-23 Thread Marius Hofert
Hi, Why does the upper left panel (in the plot below) not have a gray background? Cheers, Marius require(grid) require(gridBase) pdf(file="Rplot.pdf", width=8, height=8, onefile=FALSE) ## set up the grid layout gl <- grid.layout(5, 5, widths=unit(c(1.8, 8, 0.8, 8, 0.8), "cm"),

Re: [R] Exactly Replicating Stata's Survey Data Confidence Intervals in R

2012-09-23 Thread Thomas Lumley
On Sat, Sep 22, 2012 at 2:51 AM, Anthony Damico wrote: > Survey: Mean estimation > > Number of strata = 1 Number of obs= 183 > Number of PSUs = 15 Population size = 9235.4 > Design df= 14 > > ---

Re: [R] prune in rpart: choose number terminal nodes

2012-09-23 Thread Prof Brian Ripley
On 23/09/2012 17:32, Uwe Ligges wrote: On 21.09.2012 11:49, agent dunham wrote: Dear community, I've an rpart object, and I know the CP I want. I'd like to know if it's possible also to fix the number of terminal nodes I want. I don't know of any existing function, but you can write it you

Re: [R] aggregate help

2012-09-23 Thread Sam Steingold
Thanks. Why does > aggregate(z, list(id=z$id),FUN=list) id id a1 a2 1 10 10, 10, 10 a, a, b x, x, z 2 20 20, 20b, by, y 3 30 30 c z work, but aggregate(z, list(id=z$id),FUN=function(l) { t <- sort(table(l),decreasing=TRUE) list(length(t),t[

Re: [R] prune in rpart: choose number terminal nodes

2012-09-23 Thread Uwe Ligges
On 21.09.2012 11:49, agent dunham wrote: Dear community, I've an rpart object, and I know the CP I want. I'd like to know if it's possible also to fix the number of terminal nodes I want. I don't know of any existing function, but you can write it yourself: All the information you need is

Re: [R] grid: How to merge cells in grid.layout?

2012-09-23 Thread Marius Hofert
Ahh, now I see what you mean... Thanks, that indeed works. Cheers, Marius Marius Hofert writes: > Bert Gunter writes: > >> Inline below. >> >> On Sun, Sep 23, 2012 at 1:41 AM, Marius Hofert >> wrote: >>> Dear grid expeRts, >>> >>> I would like to create a layout with grid that looks like th

Re: [R] BRugs has a bug to use "OpenBUGS_PATH"

2012-09-23 Thread Uwe Ligges
On 21.09.2012 20:49, R. Michael Weylandt wrote: On Fri, Sep 21, 2012 at 7:00 PM, Shuguang Sun wrote: Hi, When I used BRugs (Version 0.8.0), I found a bug about the findOpenBUGS way using "OpenBUGS_PATH". I don't know how to contact the developer, so I think someone here may help me. The ea

Re: [R] grid: How to merge cells in grid.layout?

2012-09-23 Thread Marius Hofert
Bert Gunter writes: > Inline below. > > On Sun, Sep 23, 2012 at 1:41 AM, Marius Hofert > wrote: >> Dear grid expeRts, >> >> I would like to create a layout with grid that looks like the following, but >> with cells (1,1), (1,4), (4,1), and (4,4) removed and cells (2,1) and (3,1) >> (and (4,2)

Re: [R] grid: How to merge cells in grid.layout?

2012-09-23 Thread Bert Gunter
Inline below. On Sun, Sep 23, 2012 at 1:41 AM, Marius Hofert wrote: > Dear grid expeRts, > > I would like to create a layout with grid that looks like the following, but > with cells (1,1), (1,4), (4,1), and (4,4) removed and cells (2,1) and (3,1) > (and (4,2) and (4,3)) combined to one cell (s

[R] grid: How to merge cells in grid.layout?

2012-09-23 Thread Marius Hofert
Dear grid expeRts, I would like to create a layout with grid that looks like the following, but with cells (1,1), (1,4), (4,1), and (4,4) removed and cells (2,1) and (3,1) (and (4,2) and (4,3)) combined to one cell (so that contents can easily be centered. How can this be achieved? requir

Re: [R] scraping with session cookies

2012-09-23 Thread Heramb Gadgil
This may be because connection to the site via R is taking a lot of time. I too faced this problem for the site "Social-Mention". I tried very primitive approach. I put the 'if' condition in the loop. if(length(output)==0){getURL(site) }else{continue with the code} It might help you. Best, Hera

Re: [R] If Command in Plot

2012-09-23 Thread arun
HI, I guess you want to plot with a subset of your data z1. z1[with(z1,e>0),]  plot(y~x,data=z1[with(z1,e>0),]) A.K. - Original Message - From: Bhupendrasinh Thakre To: r-help Cc: Sent: Saturday, September 22, 2012 11:48 PM Subject: [R] If Command in Plot Hi Team, I am trying to