similar to: Odd graphic device behavior

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",

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

Hi again all, Several replied ASAP that I also needed reshape loaded and not just reshape2. Hmmm tried that and I had some output but not the correct format. What I need is to run simulations of time overlap between species as per the simulation program data input constraints: The basis for the simulations is a species by _time-use matrix in which species are arranged in rows, and time

Hi, does anyone know how to make decimal points midline when plotting? The journal to which we are going to submit a manuscript require this particular formatting, and it gives direction how to do this on PC: hold down ALT key and type 0183 on the number pad Thanks John [[alternative HTML version deleted]]

Dear R software faculty:       There is a question bothering me. That is why we can not set the property of row.names=false when we output csv data. But there exists that parameter in R help. Hoping for your reply! Best regards!                                                                         Li Nan [[alternative HTML version deleted]]

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

Hi, I am actually working on some auto-routine to import XML file, run some analysis on them and create graph as jpeg. The files are in different language french/english/danish even chinese. At the moment I'm focusing on the European language. I import them using the XML package and specify encoding="UTF-8" which seems to work pretty well when I write the text in the console, the

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

Given that we occasionally run into problems with comparing Excel results to R results, and other spreadsheet-induced errors, I thought this might be of interest. http://www.nextnewdeal.net/rortybomb/researchers-finally-replicated-reinhart-rogoff-and-there-are-serious-problems The punchline: "If this error turns out to be an actual mistake Reinhart-Rogoff made, well, all I can hope is that

As I've mentioned a number of times. I find it very useful to have lowess() become a generic function so that a lowess.formula() can be defined. Below is a patch that makes both changes, as well as updating the corresponding help documentation. Gregory R. Warnes Manager, Non-Clinical Statistics Pfizer Global Research and Development Tel: 860-715-3536 ? DESCRIPTION ? Makefile ?

Hi there, I'm plotting some glass RI values just by plotting plot(x) then I put on my lowess smoother lines(lowess(x)) now I want to put on some 95% Confidence Interval bands of the lowess smoother, but don't know how?? Thanks -- Gareth Campbell PhD Candidate The University of Auckland P +649 815 3670 M +6421 256 3511 E gareth.campbell@esr.cri.nz gcam032@gmail.com [[alternative

Dear All, I have a question, suppose X is a dataframe, with column names as "x1", "x2", "x3", ..... And I would like to use the i-th column by X[,'xi']. But it seems the single quote and double quote are different. So if I run X[, names(X)[i]], it has some error. Please use the below example code X = matrix(rnorm(50),ncol = 5) X = data.frame(X)

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

I'm using this code to plot a smoothed line. These two columns of data really represent 4 groups and I'd like to plot a separate line for each group but have them all in the same plot. The R-Docs for lowess do not seem to indicate some type of "GROUPS=var_name" option. What would be the syntax for this? plot(AWGT ~ lipid ) lines(lowess(lipid , AWGT, f=.8)) -- Dean

Hi all, I am trying to do something simple, but which is deviling me. I want to add an up-arrow to the x-axis, pointing to a specific location on the axis. I want the arrow to be as large as possible. Here is where I'm at with my code: #this part if from the ?plot example: require(stats) plot(cars) lines(lowess(cars)) #I am trying to use mtext to add the arrow... but I don't know how

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 R users I used the function line(x,y) and line(lowess(x,y)) to see the correlation between 2 variables (x,y). Here is my question: is there a way to ask R to tell me the equation of -this line : line(x,y) -this curve: line(lowess(x,y)) Best regards -- View this message in context: http://www.nabble.com/equation-of-a-line-or-curve-tp20628845p20628845.html Sent from the R help mailing list

Hi, I was trying as well as looking for an answer without success (a bit strange since it should be an easy problem) and therefore I will appreciate you help: My simple script is: # Loadings data of 5 columns and 100 rows of data data1<-read.csv("C:/?/MyPCA.csv") pairs(data1[,1:4]) pca1 <- princomp(data1[,1:4], score=TRUE, cor=TRUE) biplot(pca1) The biplot present the data

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]]