similar to: exporting high quality graphics from R in Mac OSX

Full_Name: George W. Gilchrist Version: 1.6.2 OS: OS X Submission from: (NULL) (128.239.124.126) Start with a tab-delimited or comma-delimited text file created on the Mac and use read.table("filename.txt", header=T) to read it in. When the first column of the file contains a character vector, and there is a header line, the first letter of the first column of the fifth row is appended

I have a very large dataframe and wish to extract a subset of rows. I have a two column matrix listing the starting and ending indices of one subset on each row. My idea is to create a vector of indices that could be applied to the dataframe and I have a solution using a for loop (below). But surely there is some more elegant way to do this! I looked thorough the archives without

For several years, I have been using Splus to analyze an ongoing series of datasets that have a bimodal distribution. I have used the following functions, in particular the ms() function, to estimate the parameters: two means, two standard deviations, and one proportion. Here is the code I've been using in S: btmp.bi <- function(vec, p, m1, m2, sd1, sd2) {

I have just installed MacPorts 1.7.0 for Leopard using a downloaded .dmg file. Following the download and installation, I looked for the .profile file to change the environment variables, but found that no .profile file was created. When I tried to do so, I got an error message: /Applications/TextEdit.app/Contents/MacOS/TextEdit .profile 2008-12-27 14:19:02.393 TextEdit[47100:10b] ***

Does anyone have an idea of how one might plot male and female symbols on a graph using R? Thanks! .................................................................. George W. Gilchrist Email #1: gwgilc at wm.edu Department of Biology, Box 8795 Email #2: kitesci at cox.net College of William & Mary Phone: (757) 221-7751 Williamsburg, VA

In S-Plus, there is a function ms() [minimum sums] that estimates parameters in a multivariate model and returns, along with the estimates, a loglikelihood score. It's very handy, but I have not found an R version yet. Does anyone know if such a thing exists? Thanks. Cheers, George ================================================================== George W. Gilchrist

I need to extract the degrees of freedom and log likelihoods from a series of mixed models computed using lmer/lme4. If I ask for logLik (lmer.object), I get something like > logLik(lmer.object) 'log Lik.' -177.1000 (df=10) Can I easily get that df from there (or elsewhere) into an object? Thank you for any ideas. George

I am trying to do a series of xyplots plots, using a "for" loop to substitute the appropriate variables for each plot. The basic command works fine by itself and produces a nice plot: > i<-3 > trellis.device(theme="col.whitebg") > xyplot(as.formula(paste(tmp00[2*i], "~ ", tmp00[(2*i)-1], +

Thanks to all who responded. Not all of the advice helped me with that particular problem, but, I have learned a lot just from your answers. 1. R for Mac 10 (Carbon version) uses the Mac convention of using ":" for the usual Unix "/" and thanks to Dr. Lumley for that since I'm new to R, Mac10, Macs in general, and Unix. 2. Interestingly, Dr. Ripley suggested to

Hi guys, I'm new in wine debugging, but i need to get content of textbox to console. I tried Code: WINEDEBUG=+edit ./app.exe 2>&1 It gets content of lineedit, but not textedit. Do you have some idea how read content of textbox from win app in wine? Thanks

In S-Plus, I can obtain polynomial contrasts for an ordered factor with contr.poly(). The function also exists in R, however is limited to factors where the levels are equally spaced. In S-Plus, one can obtain the contrasts for a set of numeric values representing unequally spaced ordered factors. Has anyone implemented this in R? I see that the S-Plus function calls another function (poly.raw())

Dear colleages in R, I have earlier been working with R in Linux, where reading in a table containing Scandinavian letters ("?", "?", and "?") in the header as part of variable names has not caused any problem whatsoever. However, when trying to do the same in R running on new MacOS-X (with an Intel processor) with the same original text table does not seem to

Hello all I'm trying to do a simulation that involves identifying the minimum point between two peaks of a (usually) bimodal distribution. I can do this easily if there are only two peaks: CnBdens<-density(Ys/Xs) #probability density function for ratio of Ys to Xs for(p in 1:512) ifelse(CnBdens$y[p]>CnBdens$y[p-1],peak1<-p,break) #identifies first peak in probability

I'm trying to learn statistics and R at the same time. I have an undergraduate science degree and one year of calculus (30 years ago), but never took a stats course. I hope to take some stats courses in the next year, but thought I would start to see how much I could teach myself. I work for an organization that analyses behavior change communication programs regarding HIV/AIDS and

Hi, I would like to have four plots plotted in the one figure and under each plot the caption (a) ... (d). Using mtext I can get so far: > plot(1:10, exp(1:10), log = "y", xlab = "xlab") > mtext("(a)",side=1,line=4.5) > plot(1:10, exp(1:10), log = "y", xlab = "xlab") > mtext("(b)",side=1,line=4.5) > plot(1:10,

Hello everyone Not strictly an R question but close... hopefully someone will be able to help. I wish to compare the switchpoints in two switchpoint regressions. The switchpoints were estimated using the segmented library running in R, and I have standard errors for the estimates. I initially thought I could just bootstrap confidence intervals for the difference between the switchpoints,

People I apologise for asking a general stats question, but I'm at a bit of a loss as to what to do following some hostile referees' comments. If I have a fully randomised blocked design, with only three blocks, should I treat block as a random or fixed effect? I have read comments about not treating block as a random effect if the number of blocks is less than 6 or 7: is this

Hello everybody, I'm a newbie in R and just went through an introduction class recently. Here's my problem. I have 2 text files (.txt) with plain numbers ('doubles' for those who know c++) ordered into 2 columns as below: coordinate1 value1 coordinate2 value2 coordinate3 value3 ... ... coordinateN valueN I would like to write a small programme in which i would: 1. take

Hi, I want to start using has_many :through. I think that eventually I''ll be removing all uses has_and_belongs_to_many from my project. Here is my situation I''m using SVN I want to use edge rails Reading the wiki about switching to edge rails[1], I can''t figure out which combination of subsections I must perform. Do I want edge rails to reside in the vendor

Hello people I am doing some thinking about how to analyse data on dimorphic animals - where different individuals of the same species have rather different morphology. An example of this is that some male beetles have large horns and small wings, and rely on beating the other guys up to get access to mates, whereas others have smaller horns and larger wings, and rely on mobility to