similar to: Contrasts with an interaction. How does one specify the dummy variables for the interaction

Dear r-community, Today I have completed my PhD.  I would like to take this opportunity to thank the r-community for helping me with the r-coding.  I use r to do data manipulation during my PhD and I benefit a lot through the discussion in the r-forum. I will continue using R and help the others too. Thank you so much. Regards, Roslina UniSA. [[alternative HTML version deleted]]

I am running my code in a loop and it does not work but when I run it outside the loop I get the values I want. n <- 1000; # Sample size fitglm <- function(sigma,tau){ x <- rnorm(n,0,sigma) intercept <- 0 beta <- 0 ystar <- intercept+beta*x z <- rbinom(n,1,plogis(ystar)) xerr <- x + rnorm(n,0,tau) model<-glm(z ~ xerr, family=binomial(logit))

I am new to R and I am trying to do a monte carlo simulation where I generate data and interject error then test various cut points; however, my output was garbage (at x equal zero, I did not get .50) I am basically testing the performance of classifiers. Here is the code: n <- 1000; # Sample size fitglm <- function(sigma,tau){ x <- rnorm(n,0,sigma) intercept <- 0 beta

R 2.5 windows XP I am getting an error from glm() that I don't understand. Any help or suggestions would be appreciated. N.B. 1<=AAMTCAREJ<=327900 > summary(data$AAMTCAREJ) Min. 1st Qu. Median Mean 3rd Qu. Max. 1.0 404.3 1430.0 6567.0 5457.0 327900.0 > fitglm<-glm(AAMTCAREJ~sexcat+H_AGE+SmokeCat+InsuranceCat+MedicadeCat+ +

library(ROCR) n <- 1000 fitglm <- function(iteration,intercept,sigma,tau,beta){ x <- rnorm(n,0,sigma) ystar <- intercept+beta*x z <- rbinom(n,1,plogis(ystar)) xerr <- x + rnorm(n,0,tau) model<-glm(z ~ xerr, family=binomial(logit)) *int*<-coef(model)[1] *slope*<-coef(model)[2] # when add error you are suppose to get slightly bias slope. However when I change

Hi. Excuse my English. I've installed Samba+OpenLDAP as a PDC. Everything works fine but Samba ignores completely group information. Linux is ok. Any clue? I'm going crazy here! Here's the sittuation: user: fish1 home dir: /home/reaml/swim/fish1 primary group: swimmers other groups: smokers Directory of smoker's group: /home/realm/smokers Here's an 'ls -l' on

Hi, We are running into a problem with a Samba setup and would like to know if a current fix or workaround is at all possible. Our setup is a NetApp filer serving NFS v4 that is mounted by Solaris and Linux servers. On those servers we are using Samba to create shares of those NFSv4 mounted filesystems. We are migrating to this NFSv4 setup from an existing Solaris NFSv3+Posix ACL setup that also

Hi all, I tried last week to do traffic shaping on a production system, object of the exercise was simply to throttle the outgoing traffic. tc qdisc add dev eth4 root tbf rate 300mbit burst 300kb latency 50ms But the server became rather instable, crashing repeatedly without anything in the logs. Can anybody spot glaring mistakes in the tc command above, or tell me what I should have done

Hello there, I have two groups (men and women) and I know per group how many of them smoke or don't smoke (women 40 of 200; men 100 of 300). I would like to know how I can compare in R if men and women differ significantly in their smoking. However, because there are more men in the sample than women I cannot just compare the number of smokers and non-smokers in both groups, right?! (I would

Dear UseRs, Suppose I have data regarding smoking habits of a prospective cohort and wish to determine the risk ratio of colorectal cancer in the smokers compared to the non-smokers. What do I do at the end of the study with people who die of heart disease? Can I just censor them exactly the same as people who become uncontactable or who die in a plane crash? If not, why not? I'm thinking

Hi all! Just a word of warning: after updating a few of our x86_64 based web frontend boxes to the new kernel, we began to get weird MySQL timeouts. The problem went away again when we downgraded to the previous kernel-2.6.18-53.1.6.el5.x86_64.rpm regards, Bent Terp

Hi. Can anyone suggest a simple way to obtain in R a list of vector coordinates of the following form? The code below is Mathematica. In[5]:= Flatten[Table[{i,j,k},{i,3},{j,4},{k,5}], 2] Out[5]= {{1,1,1},{1,1,2},{1,1,3},{1,1,4},{1,1,5},{1,2,1},{1,2,2},{1,2,3},{1 ,2,4},{1,2, 5},{1,3,1},{1,3,2},{1,3,3},{1,3,4},{1,3,5},{1,4,1},{1,4,2},{1,4,3}, {1,4,

I was downloading offices from sourceforge.net to test on Wine and to see if I liked one better that the ones that I was using now. Liberty Office launches then says that there is no parent window. Terp office prints this in the console log. Modules: Module Address Debug info Name (11 modules) PE 400000- 44e000 Export terp PE 60b50000-60b54000 Deferred advapi32 PE

Hi, Does anyone know if: with R can you take a set of numbers and aggregate them like you can in SPSS? For example, could you calculate the percentage of people who smoke based on a dataset like the following: smoke = 1 non-smoke = 2 variable 1 1 1 2 2 1 1 1 2 2 2 2 2 2 When aggregated, SPSS can tell you what percentage of persons are smokers based on the frequency of 1's and 2's. Can

Hi list! While regular backup solutions like amanda or bacula are very good at their job, ie keeping point2point copies of the files currently on disk, I find them less suited for archiving - having unused files move to tape in duplo and stay there until requested. I've even read of multi-tier solutions - move to slower disks after a week and further on to tape after a month. Does anyone

On Mon, Jun 2, 2008 at 2:03 PM, Johnny Hughes <johnny at centos.org> wrote: > I would also not use XFS in production ... but that is just me. Interesting, I thought that XFS was fairly safe for use. What would you recommend for filesystems in the 50-500 terabyte range? (And yes, we do actually run a 70 TB at the moment, so I'm not asking just to annoy you; I'm genuinely

Hi, i am trying to read a .txt file, do a couple of select if statements on my data, and then finally use the ?table function to get frequency counts on the data. Specifically, i am looking at answering the following question: What is the frequency of Grade 7 students in the province of Alberta who are smokers? I am having some problems: 1)i cannot get the column names to show up when print

Dear All, I'am looking for examples showing that correlation does not imply causality, the targeted audience consists of undergraduate students (their first year at the university but in the BioMathStat track). All practicals are under R. I was able to extract this from R datasets: ### begin data(sunspots) data(lynx) spots <- window(sunspots, freq = 1, start = 1880, end = 1900) lnx <-

On Jan 30, 2008 8:25 AM, Johnny Hughes <johnny at centos.org> wrote: > Joshua Baker-LePain wrote: > > On Wed, 30 Jan 2008 at 10:18am, Johnny Hughes wrote > > > >> Bent Terp wrote: > > > >>> Has something changed with regard to the mount options? We use > >>>

Here is an even more awkward property of tapply with factors. What we want to do is to determine the modal value of the first factor given the level of the second factor. We do that by applying table in a tapply call then finding the maximum count then ... Turns out that table gives unexpected results in tapply. > table(Machines$Machine) A B C 18 18 18 > tapply( Machines$Machine,