similar to: Categorical Variables and glm()

In the following code, the only difference between the two plots is the order the variables are plotted. In this case, the plot of "cdata.den" in plot #1 is different from its plot in #2. Specifically, "cdata.den" spans the x-axis from -5 to 30 in plot #1 and from 0 to 20 in plot #2. Does anyone understand why these two plots do not yield the same result? #Make density

I am simulating sickness among a group of families. Part of the task is to randomly draw who in the family will be sick, randomly drawing from family ID's where Dad =1, Mom = 2, Kid1 = 3, Kid2 = 4., etc. My census of Dads is of the form shown below. Dad_ID Spouse (Y=1;N=0) #Kids #People_Becoming_Sick 1 1 0 1 2

All - I am evaluating an arima(2,1,3) and arima(3,1,3) and notice the log-likelihood of the restricted model is higher than the log-likelihood of the unrestricted. Since these are nested models, I thought the unrestricted model would have a log-likelihood at least as large as that of the restricted model. Am I interpreting the "loglik" output incorrectly? Regards, Stephen

I've seen with other software the capability for the optimizer to switch algorithms if it is not making progress between iterations. Is this capability available in optim()? Thanks, Stephen Collins, MPP | Analyst Health & Benefits | Aon Consulting [[alternative HTML version deleted]]

I currently run simulations on an IBM x31 Thinkpad (laptop), with an Intel Pentium M, 1.3GHz processor. I am planning to buy a desktop to help speed up my work, but I am wondering if R is compatible with some of the newer technology that has come to market (i.e., Intel Core 2 Duo, Quad Core processors, 64-bit operating systems). Does anyone know the threshold or limit on R's

Hi all, simple Q: how do I extract the upper and lower CI for predicted probabilities directly for a glm - I'm sure there's a one line to do it but I can't find it. the predicted values I get with the predict (.. "response") Thanks Gerard ********************************************************************************** The information transmitted is intended only for

Dear R-experts, I am running R version 2.7.1 on Windows Vista. I have a small dataset which consists of ?chick ID?, ?year (0, 1)?, ?hatching order [HO, defined as first, second and third-hatched chick]?, and the binary outcome of interest ?death (0, 1)?. So a subset of my dataset looks like this on a txt file: y ID Yr HO 1 1 1 First 0 2 1 First 0 3 1 Second 0 4 1 First 1 5 1 First 0 6 1 Third

Hello, I'm struggling with an elementary problem with R. I have a simple data frame such as this one giving the number of accidents subdivided by sex, age and region. sex age region no_of_accidents F young north 10 F young south 12 F old north 5 F old south 7 M young north 24 M

Hi How can i enable this feature because: Running /opt/xensource/libexec/mpp-rdac --enable it says: No module mppUpper found for kernel 2.6.32.12-0.7.1.xs1.0.0.298.170582xen, aborting. There is no MPP modules... Terveisin/Regards, Pekka Panula, Sofor Oy - Jatkuvat palvelut _______________________________________________ Xen-users mailing list Xen-users@lists.xensource.com

Hi all, I have a multiple queing situation I'd like to simulate to get some idea of the distributions - waiting times and allocations etc. Does R has a package available for this - many years ago there used to be a language called "simscript" for discrete event simulation and I was wondering if R has an equivalent (or hopefully with graphics, something better!). Apologies if there

Dear all, sorry to bother you all with this but I've been trying to use the loglm in MASS package (v2.8.0) and cannot get any sensible output. I'm wondering am I doing something very foolish or missing something obvious. For example, I tried the documentation help(loglm) example - here's the code # Case 1: frequencies specified as an array. sapply(minn38,

Dear all, I've replicated the cheese tasting example on p175 of GLM's by McCullagh and Nelder. This is a 4 treatment (rows) by 9 ordinal response (cols) table. Here's my simple code: #### cheese library(MASS) options(contrasts = c("contr.treatment", "contr.poly")) y = c(0,0, 1, 7, 8,8,19, 8,1, 6,9,12,11, 7,6, 1, 0,0, 1,1, 6, 8,23,7,

Hi all! I need to compile de IBM Mpp-rdac driver in XCP 1.0 to use with the IBM DS4700 storage. Is there any way to do this? -- Rogério da Costa _______________________________________________ xen-api mailing list xen-api@lists.xensource.com http://lists.xensource.com/mailman/listinfo/xen-api

Hi all, apologies if this is obvious - but I can't see it and would appreciate some quick help! the matrix mhouse is 26x3 and I'm computing odds ratios. The simple code below "should" compute the odds vector for every pair (325) i.e. 26C2 in cols 1 and 2. On the first i=1 outer loop the inner j loop runs from 2 to 26 ok and then I get the error (Error: subscript out of bounds)

Has anyone had any experience with the Dell PowerConnect 2724 Gigabit Switch and VoIP traffic? It seems to support Class of Service (CoS), but not full Quality of Service (QoS). Will CoS help prioritize VoIP traffic at all? I am setting up an asterisk PBX for 22 users with very light traffic, but I don't want to have the call start skipping just because someone is pulling a large

I am using mle as a wrapper from optim( ). How would I extract the convergence code, to know that optim( ) converged properly? Thanks, Stephen Collins, MPP | Analyst Global Strategy | Aon Benfield [[alternative HTML version deleted]]

Say I have a formula Y ~ 1 + X, where X is a categorical variable. A previous thread showed how to evaluate this model using the mle package from "stats4" (see below). But, the user had to create the data matrix, X, including the column of one's for the regression constant. Is there a way to nest the linear formula in the code below, so the data matrix doesn't explicitly

Hi all, I read the setup at http://www.blackh0le.net/articles/vpn-dun-howto.html to setup my VPN. However, I'm having a problem which I think is proxy-ARP not working. I like to ask you to see if you know what's going on. When I ping 10.77.1.1 from windows XP machine the packets get to the 10.77.1.1 machine, but they don't have a return path to get back. When I do ping the windows

http://bugs.freedesktop.org/show_bug.cgi?id=25795 Summary: windows not redrawn after resize Product: xorg Version: 7.4 Platform: x86-64 (AMD64) OS/Version: Linux (All) Status: NEW Severity: normal Priority: medium Component: Driver/nouveau AssignedTo: nouveau at lists.freedesktop.org

Hello, I have IBM Storwize V3700 storage, connected to 2 IBM x3550 M4 servers via fiber channel. The servers are with QLogic ISP2532-based 8Gb Fibre Channel to PCI Express HBA cards and run Centos 5.10 When I export a volume to the servers, each of them sees the volume twice, i.e /dev/sdb and /dev/sdc, with the same size. Previously I have installed many systems with IBM DS3500 series of