similar to: aggregate binary response data

I have a large dataset organised in choice groups see below. Each choice represents a separate occasion with 1 product chosen out of the 6 offered.            +-------------------------------------------------------------------------------------------------+      | dn   obs   choice      acid   br                 date       cdate   situat~n   mth   year   set |     

Dear List members I have a large dataset organised in choice groups see sample below +-------------------------------------------------------------------------------------------------+ | dn obs choice acid br date cdate situat~n mth year set |

Dear List I am trying to fit a multinomial model using the mlogit package. Attempting to load the data into mlogit presents the following error. MLOG<-mlogit.data(Mult3,shape="long",choice="CHOICE",alt.var="mode.ids",indivs = "set3",chid.var = "obs") Error in `row.names<-.data.frame`(`*tmp*`, value = c("1.1", "1.2",

Dear List Has anyone used R to distnguish between alternative forecasting models? In particular is the Diebold Mariano test available for use within R. Any assistance would be greatly appreciated. Graham Leask Lecturer in Strategy Economics & Strategy Group Aston University Aston Triangle Birmingham B4 7ET Tel: 0121 204 3150 E Mail: g.leask at aston.ac.uk

Has anyone implemented Tibrishani's gap statistic in R or S plus? If so I would greatly appreciate the relevant script file. Any help will be much appreciated Kind regards Dr Graham Leask Economics and Strategy Group Aston Business School Aston University Aston Triangle Birmingham B4 7ET Tel: Direct line 0121 204 3150 email g.leask@aston.ac.uk [[alternative HTML version deleted]]

On Saturday 15 September 2007, Tanya Lattner wrote: > 2) Download llvm-2.1, llvm-test-2.1, and the llvm-gcc4.0 source. > Compile everything. Run "make check" and the full llvm-test suite > (make TEST=nightly report). I tried to do this, but ran into trouble. LLVM itself compiled fine. GCC compilation is broken when Java support is enabled. I can file a detailed bug report

Dear List I am looking for a script that will calculate the Box M test to test the homogeneity of the variance/covariance matrix between two matrices. If anyone could send me the script I would appreciate it. I am aware of the scepticism about this test, where due to extreme sensitivity a p value of 0.01 is recommended. Despite this however Box's M test is the established

Hi, I have constructed a dataframe as follows: Oil <- rep(c("Oil1","Oil2","Oil3"),8) Comp <- rep(rep(c("C1","C2"),c(4,4)),3) Mth <- rep(c("M1","M1","M2","M2"),6) Meas <- rep(c(1,2),12) Result <- rnorm(24,mean=100, sd=5) df <- data.frame(Oil, Comp, Mth, Meas, Result) The same compound

Is there a function available in R that implements Mojena's Upper Tail Rule or that draws a Mojena plot? I would also like to find a function that implements Duda and Hart's stopping rule. Finally with function cophenet how can I achieve a straightforward Cophenet correlation coefficient (i.e. one number for example 0.876) that gives the agreement between the structure of the

Summary ======= After an unclean shutdown the file /var/run/dovecot/master.pid remained behind. This prevented dovecot from starting, and gave a misleading error message. To be more resilient and fault-tolerant, I recommend that dovecot also check the validity of the PID in /var/run/dovecot/master.pid in order to determine whether or not another dovecot process is running. Detail ====== In

Hi there, I am trying to run simulations using R with linear mixed model (lme). There are two factors in my fixed effect model, educ (treatment and control) and mth (visit 1, 2, and 3). What I want to obtain is the estimated treatment difference (treatment - control) at visit 3, plus the standard error and p-value. This can be easily obtained in SAS using lsmeans or estimate statements, but I

Hi, I am traying to implement an STAR model, but I have some problems. I am following the instruction of the model, that they are in: http://bm2.genes.nig.ac.jp/RGM2/R_current/library/tsDyn/man/star.html that they are from: http://bm2.genes.nig.ac.jp/RGM2/pkg.php?p=tsDyn The model is: star(x, m=2, noRegimes, d = 1, steps = d, series, rob = FALSE, mTh, thDelay, thVar, sig=0.05, trace=TRUE,

Dear List I am having a bit of a problem getting a program to work. For each of i=1 to n persons I have a matrix (different for each person) with m rows. What I want to do, is create m new data sets such that the first is made up of the first row for each person from the original matrices, the second contains the second row for each person from the original matrices etc etc up to the mth

For a list say; list1<-{1,2,3,4,5,2,1} How do I remove the duplicates please? My real list is 20,000 obs long of dates with many duplicates Regards Glenn [[alternative HTML version deleted]]

Warnings about /metno/internweb/fouweb/htdocs/teknisk/verktoy/R/faq.html BAD LINK: http://stat.ethz.ch/R/manual/ BAD LINK: http://nlme.stat.wisc.edu/MEMSS/ BAD LINK: http://stats.mth.uea.ac.uk/Rcgi/ BAD LINK: http://software.biostat.washington.edu/statsoft/ess/ ------------------------------------------------------------ Warnings about

I have a large dataset (~1 million rows) of three variables: ID (patient's name), DATE (of appointment) and DIAGNOSIS (given on that date). Patients may have been assigned more than one diagnosis at any one appointment - leading to two rows, same ID and DATE but different DIAGNOSIS. The diagnoses may change between appointments. I want to subset the data in two ways: - define groups

Sensei H. Peter Anvin, I have successfully modified/enhanced memdisk so that one can chainload a local operating system after running diskless FreeDOS. I would like your advice regarding the appropriate interface for invoking this functionality. In addition, I would like to know what steps to take so that you feel comfortable considering including this in the standard syslinux distribution.

Removing rows with NAs, using na.omit(), doesn't seem to be working for me. Dataset: > str ( ex10s ) 'data.frame': 2189576 obs. of 5 variables: $ LOPNR : int 58 58 58 58 64 64 64 64 64 64 ... $ DIAGNOS: Factor w/ 173 levels "F20","F200","F2000",..: 128 128 128 128 105 105 105 160 105 105 ... $ X_DATE : int 20060821 20061207 20080102 20090904

Hi, I am using Rattle 2.6.4 with R 2.12.2 on win64, is this a bug ? Following is the content after execute the associate analysis process: Summary of the Apriori Association Rules: Number of Rules: 23351 Summary of the Measures of Interestingness: support confidence lift Min. :0.1250 Min. :1 Min. :2.667 1st Qu.:0.1250 1st Qu.:1 1st Qu.:2.667 Median :0.1250

Thanks. I was "believing" an Excel import of the file - as you say, a re-import into R still has 333 columns. I shall look elsewhere for the cause of the mysterious truncation. Stuart > ----- Original Message ----- > From: "Thomas Lumley" <tlumley at u.washington.edu> > To: "Stuart Leask" <stuart.leask at nottingham.ac.uk> > Cc: