similar to: Lattice plot within a "for" loop does not happen?

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

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

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

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

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/> <%=

Hello, Excuse me for posting two questions in one day, but I figured it would be better to ask my questions in separate emails. I will again give the caveat that I'm not a statistician by training, but have a fairly decent understanding of probability and likelihood. As before, I'm trying to fit a nonlinear model to a dataset which has two main factors using nlme. Within the dataset

On UNIX, I am executing an rsync command, from within a script. The command goes something like this: /usr/bin/rsync --verbose --progress --stats --compress --recursive --times --perms --links --safe-links source_dir/ user at target_machine:/parent_path/source_dir In other words, I am replicating source_dir on a remote machine. It ends up next to a lot of sibling, directories, like this: On

Hello, I'm looking to upgrade my Samba from 3.0.4 to 3.0.7. Samba is my PDC currently. Are there any issues that I need to be aware of before I upgrade? Will I need to re-join the domain after the upgrade? Rohan Gilchrist rohan@e-mailme.org http://www.e-mailme.org/~rohan/ 0412 648 909 *************************************************************************** This e-mail and any files