similar to: Looking for R analogue of S-plus function ms()

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

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

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

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

[This was filtered (looking like spam) and manually approved by your list maintainer, MM ] EvolDir wrote: > Colleagues, > > A few days ago I posted a message inquiring what folks were using for > teaching biostats. I received a stunning 38 replies offering various > opinions, some quite surprising, along with links to web pages, pointers for > coverage, commentaries on

Can Samba be used to generated WinPopUp RPC calls on TCP port 135? I am familiar with using the "smbclient -M" command to send messages using the Win Popup protocol, otherwise known as Windows Messenger Service, or what a Microsoft user does from a DOS prompt as a "net send" command. WinPopup can be generated either via Microsoft file sharing protocols, using the

Hi, My branch office as got a 256Kbps b/w from their service provider at a very very high rate per Mb. They don''t require 256Kbps at all but the ISP does not offer anything low. Can we restrict the bandwith to say 64Kbps nothing fancy? How do I go about it? With warm regards, -Payal

Hi, I have 2 interfaces - one for adsl and other for LAN on my Linux gateway machine. The IP addresses are 10.10.10.3 & 192.168.10.101 respectively. Now my routing tables show this particular entry. What exactly is this? 169.254.0.0/16 dev eth0 scope link Or by traditional route -n, 169.254.0.0 0.0.0.0 255.255.0.0 U 0 0 0 eth0 What exactly is this 169.254.0.0/16

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

Greetings all, While I haven''t given up on getting a dead gateway detection script a la nano.txt to work on my T1/DSL router box, I have taken a break and played with the mpath2.sh script that Robert Kurjata wrote back in 2001. This shows more promise in that I can actually get it to work - is it foolish to try and use this script for this purpose? Here''s my setup: My router

Greetings all I'm running Samba 3.0 on an Ubuntu box as a PDC and I'm having trouble changing passwords with XP clients - here's my smb.conf http://pastebin.com/m1bb6d4a6 I've played with a variety of "passwd chat" settings but no joy. I am trying to use pam_cracklib.so - here's my /etc/pam.d/common-password file: http://pastebin.com/m1a1d5f89 I've tried the

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