similar to: polynomial contrasts in R

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) {

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

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

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

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

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've been experimenting with drop1 for my biostatistics class, to obtain the so-called Type III sums of squares. I am fully aware of the deficiencies of this method, however I feel that the students should be familiar with it. What I find baffling is that when applied to a fully balanced design, you obtain different sums of squares. I've used this for several years in Splus and R and never

My apologies to the list for sending this without adequate research. I have found my answer; please ignore! Thanks. I've been experimenting with drop1 for my biostatistics class, to obtain the so-called Type III sums of squares. I am fully aware of the deficiencies of this method, however I feel that the students should be familiar with it. What I find baffling is that when applied to a fully

Hi there The default option for saving graphics from R (1.9.1) on my Mac is as a pdf file. If I open the file in Acrobat reader it looks really good and crisp, and is obviously saved as vector graphics, since I can zoom in as much as I like and it continues to look really nice. If I import it into MS Word (from office 2000), or Textedit, however, it imports it as a bitmap and unless I save

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

Hello [R]-help I am trying to find > a package where you can do ANOVA based trend analysis on grouped data > using orthogonal polynomial contrasts coefficients, for unequally > spaced factor levels. The closest hit I've had is from this web site: >(http://webcache.googleusercontent.com/search?q=cache:xN4K_KGuYGcJ:www.datavis.ca/sasmac/orpoly.html+Orthogonal+polynomial >l but I

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

For block effects with small variance, lmer will sometimes estimate the variance as being very close to zero and issue a warning. I don't have a problem with this -- I've explored things a bit with some simulations (see below) and conclude that this is probably inevitable when trying to incorporate random effects with not very much data (the means and medians of estimates are plausibly

Dear All I would like to produce interaction boxplots and this seems to work: par(mfrow=c(2,2)) A=sample(rnorm(50,50,10)) B=sample(rnorm(50,100,10)) Test=merge(A,B,by=0)#by=0 where 0 is the row.names TreatA=(gl(2,50,100,labels=c("High","Low"))) TreatB=rep(gl(2,25,50,labels=c("High","Low")),2) Newdata=data.frame(TreatA,TreatB,Test)

Hi, I would like to figure out how to determine what the id of the last element in a table is. Is there an equivalent to the Table.find() method which will tell me the number of elements in a table? also, maybe more importantly, where are the available methods on tables documented? What i am trying to do is to increment the id and work through a table and check whether the id i am currently

Hello, I am currently using R to run an external program and then read the results the external program sends to the stdout which are tsv data. When R reads the results in it converts it to to a list of strings which I then have to maniuplate with a whole slew of commands (which, figuring out how to do was a reall challenge for a newbie like myself)--see below. Here's the code I'm

Hello, I was wondering if anyone else has experienced this. >From an XP machine printing to a SAMBA printer, printing is painfully slow. Selecting the printer, and even sending the job through windows makes the app appear non-responsive. The print job is completed fine, but just takes a long time to print. XP boxes with SP1 print like a dream. I've managed to replicate this on 4 boxes

I have a view with a form and i want to be able to process the button selected in a controller - not sure how to do this. The form does not map directly to a model - all I want to do is be able to pass back the selected button and know which one it is. So, i have a form and several radio buttons - the submit (not shown) takes me to the process_answer action where I want to determine which

Still the same problems I'm afraid. I read the release notes and followed them to the letter. I'm now seeing: Finished Dependency Resolution lvm2-2.02.46-8.el5.i386 from base has depsolving problems --> Missing Dependency: device-mapper >= 1.02.32-1 is needed by package lvm2-2.02.46-8.el5.i386 (base) mkinitrd-5.1.19.6-54.i386 from base has depsolving problems --> Missing

Hi, I am trying to create a basic form in which I want one of the fields initialized before displaying it. My code (in a view) is as follows: <% if ( session[:user_id] != nil ) then logged_in_user = User.find(session[:user_id]) end %> <% form_for :suggestion do |form| %> <label for="suggestion_title">Topic:</label><br/> <%=