similar to: Multiple lon lat points in the map with ggplot2

If I wanted to fit a logit model and account for clustering of observations, I would do something like: library(Design) f <- lrm(Y1 ~ X1 + X2, x=TRUE, y=TRUE, data=d) g <- robcov(f, d$st.year) What would I do if I wanted to do the same thing with a probit model? ?robcov says the input model must come from the Design package, but the Design package appears not to do probit? Thanks very

Dear all, I am a R beginner, and I am looking for a way to do the same thing for all levels of a column in a table. Basically, I have a bunch of protein sequences composed of different amino acid residues, and each residue is represented by an uppercase letter. I want to calculate the ratio of different amino acid residues at each position of the proteins. Here is an example table: Proteins

suddenly a few days ago it takes about a minute for local help to come up . i type : ?hist and get starting httpd help server ... done but only after a minute does http://127.0.0.1:15686/library/graphics/html/hist.html appear . Lee De Cola, PhD, MCP. DATA to Insight LDECOLA@COMCAST.NET Reston, Virginia 703 709 6972 571 315 0577 mobile [[alternative HTML version deleted]]

Hi, I want to make a plot similar to sm1 (attached). The code I tried is: dcols <- densCols(x,y) smoothScatter(x,y, col = dcols, pch=20,xlab="A",ylab="B") abline(h=0, col="red") But it turned out to be s1 (attached) with big dots. I was wondering if anything wrong with my code. Thanks,Zhengyu -------------- next part -------------- A non-text

Hi people, I'm having problems of connection with a Linksys SPA IP PHONE 942 when I use the WAN port, most of the time when I try to connect to the network or restart the IP Phone I can't get internet connectivity . I tried using both static IP and DHCP, but the problem is the same. Some have had similar problem with this brand of IP Phone? Thanks for the help and attention. Hugs

Hi all, I have two vectors (columns) called "efinal" and "efinal 2". I want to plot them on the same plot and "draw" a shaded area beween the two lines using function polygon I have tried all but I don ?t understand the polygon area, can you help me with examples? plot(efinal,type="l",ylim=range(min(efinal2),

I am lousy at simple regex and I have not found a solution to a simple problem. I have a vector with some character values that I want to split. Sample data dd1 <- c( "XXY (mat harry)","XXY (jim bob)", "CAMP (joe blow)", "ALP (max jack)") Desired result dd2 <- data.frame( xx = c("XXY", "XXY", "CAMP",

hi i try to do a graph of a time series which shows in red the values > -0.05, in blue the values >0.05 and in white the values between -0.05 and 0.05 for un exemple :http://www.appinsys.com/globalwarming/enso.htm thanks !!!!! denisse -- View this message in context: http://r.789695.n4.nabble.com/how-to-do-a-graph-with-tree-different-colors-tp4635206.html Sent from the R help mailing

Hi, Kbytes RSS rNo 2689632 1450876 1433788 2689632 1450876 1433788 2689632 1450876 1433788 2689632 1450876 1433788 2689632 1450876 1433788 2689632 1450876 1433788 2689632 1450876 1433788 2689632 1460168 1443084 plot(data$rNo,data$RSS,pch=0,type="b",col=36,axes=FALSE, ylab="", xlab="",las=2,lwd=2.5) title(" PID -

Dear All, I have a data sets as shown below A (Patient ID ), B and C are the Concentration of drug in blood on day 1 and day 4, D is the difference in conc. To do this in R I have written a code as follows, identified the number of patients who have more concentration on day 4 . Here I want to identify specifically the patient ID (is he patient 1 or 2 or 5 and 7), whose concentration is more. How

I would like to plot two data sets (frequency (y-axis) of mean values for 0-1(x=axis)) on a single histogram for comparison. The hist() only allow the overlay of two histograms, and although barplot() allows beside=TRUE, it does not show frequency values (like hist) but rather all of the values. Is there any way that I can use the hist() to plot two data sets similar to the barplot(). Any help or

I am following a document teaching how to use regression and right at the onset I get an R error. I understand that variables "conc" and "signal" should have the same length but I am using the what R manual suggests to drop the error and it is not cooperating. What gives? The manual says that the default is na.ombit by the way, not that it's addition does anything How

A stupid question but does anyone know how to express the actual colours used by the main Canadian political parties? I want to do a couple of ggplot2 plots and have lines or rectangles that accurately reflect the party colours. I can probably play around with RColorBrewer or something to figure it out but if some some already has got them it would save me some time especially with the NDP

Hi all, Sorry for the VP8 question on the Theora list, but I think all the relevant people are here :) I'm using Icecast to distribute some Ogg Theora streams at the moment with a view to adding VP8 along side these in future. Ideally I'd just use Ogg VP8+Vorbis, I know there is a mapping for this already but how much support from the browsers can I expect for this configuration?

Hi all, I am using quantile regression with svy design. I want to retrieve summary regression statistics (std error, p-value), since I don't have any in my output: Commands: clus1_d<- svydesign(id=~cd002_co, weights=~wtper, strata=~str, data=data) bclus1<-as.svrepdesign(clus1_d,type="bootstrap",replicates=100) fit1<-

Hi: I have a dataframe containing comma seperated group of words such as milk,bread bread,butter beer,diaper beer,diaper milk,bread beer,diaper I want to output the frequency of occurrence of comma separated words for each row and collapse duplicate rows, to make the output as shown in the following dataframe: milk,bread 2 bread,butter 1 beer,diaper 3 milk,bread 2 Thanks for help! deb

Dear Group, I want to do a t-test calculation on a large data set. I am pasting some part of it structure(list(Name = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("CTK100", "CTK103", "CTK121"), class = "factor"), Score = c(236.9726, 207.0055, 237.3464, 224.4774, 236.5034,

Gracias Carlos por la clara explicación, si alguien más quiere aportar, bienvenido El 28 de enero de 2014, 12:07, Carlos Ortega <cof@qualityexcellence.es>escribió: > Hola Marta, > > Lo que quiere decir (lo que yo entiendo que quiere decir) es que > directamente el código es la misma función de Google a la que se la recubre > ("wrapper") con código "R"

Dear all, I am trying to plot a matrix I have? as an image str(matrixToPlot) ?num [1:21, 1:66] 0 0 0 0 0 0 0 0 0 0 . ?that contains only 0s and 1s, where the xlabel will be Labeled as str(xLabel) ?num [1:66] 1e+09 1e+09 1e+09 1e+09 1e+09 ... and the yLabels will be labeled as str(yLabel) ?num [1:21] -88 -87 -86 -85 -84 -83 -82 -81 -80 -79 ... I have found on the internet that I can do

Hi, I have developed one application using ggmap package.It is based on google map. I am a bit confused regarding its licensing terms.I want to know that can i use it with my other applications. Is it legal? -- With Regards, Narendra Pratap Principal Research Developer Hitech Robotic Systemz Ltd. [[alternative HTML version deleted]]