search for: incid

Hello, I have a problem concerning logistic regressions. When I add a quadratic term to my linear model, I cannot draw the line through my scatterplot anymore, which is no problem without the quadratic term. In this example my binary response variable is "incidence", the explanatory variable is "sun": > model0<-glm(incidence~1,binomial) > model1<-glm(incidence~sun,binomial) > anova(model0,model1,test="Chi") Analysis of Deviance Table Model 1: incidence ~ 1 Model 2: incidence ~ sun Resid. Df Resid. Dev Df Devianc...

I am merging two dataframes using a relational key (incident number and incident year), but not all the records match up. I want to be able to review only the records that cannot be merged for each individual dataframe (essentially trying to select records from one dataframe using a multi-value relational key from the other dataframe). The following cod...

Incident Information:- Originator: "G Lalithambika"<lalithambika@mumbai.tcs.co.in> Recipients: samba@samba.org Subject: FOCHAP06 WARNING: The file FOCHAP06.DOC.pif you received was infected with the W32/SirCam@MM virus. The file attachment was not successfully cleaned.

Hi, I'm really struggling with barplot I have a data.frame with 3 columns. The first column represents an "incident" type The second column represents a "month" The third column represents a "time" Code for a sample data.frame incidents <- rep(c('a','b','d','e'), each =25) months <- rep(c(1,2), each =10) times <-rnorm(100) # make m...

What''s the best way to handle a requirement that shows up as a sub-requirement requirement in other features? For example let''s say users can enter dates in various forms throughout my application. There is one set of global rules specifying the formats in which dates may be entered, and how they are interpreted. I put that in one feature. In various other features, like

Dear all, I have a stat question that may not be related to R, but I would like to have your advice. I have just read a medical paper in which the authors report the 1-p (where p is the cumulative survival probability from the Kaplan Meier curve) as incidence of disease. Specifically, the study followed ~12000 women on drug A and ~20000 women on drug B for 12 months. During that period 29 women on drug A and 80 on drug B had the disease. The incidence of disease for A and B was 0.24% and 0.30% respectively. However, instead of reporting the...

Hi, I am wondering if anyone can explain to me if cumulative incidence (CI) is just "1 minus kaplan-Meier survival"? Under what circumstance, you should use cumulative incidence vs KM survival? If the relationship is just CI = 1-survival, then what difference it makes to use one vs. the other? And in R how I can draw a cumulative incidence plot. I kn...

I am designing a package to analyze fire department data. I created an S4 class called CAD (short for Computer Aided Dispatch data) as part of this package that essentially extends the data.frame class. When the user initializes a CAD object, a data.frame object is broken down into seperate incident data and unit response data that are stored in separate slots. This separation is intentional because typically a user would want to analyze either incident data (e.g. plot fires by day of week) or response data (what is the 90th percentile response time for each unit). The two data sets are li...

Im writing a module in the lib directory. If I change the module, the changes are not visible to the application when i releod the web page (i am in the development evironment). I must restart the web server (WEBrick) and so i see the changes... but I lost more time. How can I reload the module in the dirictory lib, without reload WEBrick? Thanks so much --Reis -- Posted via

I would like to compare the incidence rates of three groups. They are supposed to have different risks so I would like to test whether there is a increasing trend in the incidence rates. Does R or any packages provide a trend test for incidence rates? I checked epiR and epitools. It seems they do not have this function. Thank you...

HI Greg, Sorry, I misunderstand your question. I am not sure whether it works with numSummary() from library(Rcmdr). You could use other ways, such as: test<-read.table(text=" ?id year incidents ?100??? 1??????? 0 ?101??? 1??????? 1 ?102??? 1??????? 21 ?100??? 1??????? 27 ?101??? 1??????? 3 ?102??? 1??????? 12 ?100??? 2??????? 5 ?101??? 2??????? 5 ?102??? 2??????? 19 ?100??? 2??????? 10 ?101??? 2??????? 2 ?102??? 2??????? 12 ?100??? 3??????? 0 ?101??? 3??????? 0 ?102??? 3??????? 22 ?10...

...I found a problem with 'lme4'. Basically, once you load the package 'aod' (Analysis of Overdispersed Data), the functions 'lmer' and 'glmer' don't work anymore: library(lme4) (fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)) (gm1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), family = binomial, data = cbpp)) install.packages("aod") library(aod) (fm1 <- lmer(Reaction ~ Days + (Days|Subject), sleepstudy)) (gm1 <- glmer(cbind(incidence, size - incidence) ~ period + (1 | herd), family...

Dear R-helpers, I tried to create a 3D surface showing the interaction between two continuous explanatory variables; the response variable is binary (0/1). The model is: model<-glm(incidence~sun*trees,binomial) then I used "wireframe" to create a 3D plot: xyz<-expand.grid(sun=seq(30,180,1),trees=seq(0,4000,10)) xyz$incidence<-as.vector(predict(model,xyz)) wireframe(incidence~sun*trees,xyz,scales=list(arrows=FALSE)) which gives me a 3D plot, but the scaling of t...

Hi, I have 2 models, how can I save both of them in one transaction? I want to do something like this: begin transaction project.save person.save end transaction The only examples I''ve seen is with saving 2 instances of the same model in a transaction. thanks scott. -------------- next part -------------- An HTML attachment was scrubbed... URL:

Hi I am working on educational assignment to produce an *Incidence matrix *from a *BIB design *using R language software. I found a web page *http://wiki.math.yorku.ca/index.php/R:_Incidence_matrix * about the problem. But it produces Data matrix instead of Incidence matrix. can you please help me out using R software. thanks and regards zaheer *following...

...le, censored Dark,male, censored Fair, female, censored Dark, female, censored Am I correct in thinking that the 2 "(Gray's) Tests" which will be printed by R tell us (i) if there are significant differences between those with fair hair and those with dark hair as regards cumulative incidence of relapse [taking into account sex differences] (ii) if there are significant differences between those with fair hair and those with dark hair as regards cumulative incidence of TRM [taking into account sex differences] ? The 'est'and 'var'values are the same regardless of whe...

Hi, In my analysis of impacts of insecticide-treated bednets on malaria, I look at the relationship between malaria incidence and mosquito behaviors. The condensed data set is copied here. Ordinary regression (lm) shows that Incidence was negatively related to Mortality. This makes sense because the latter reflected the strength of killing mosquitoes by insecticide-treated nets. Since the original data set has a compl...

Dear all, how can I convert an m x n incidence matrix into an m x m adjacency matrix or an n x n adjacency matrix? The current matrix contains binary data, hence the new matrix would contain counts of common occurrences. Thank you for your help. Phil

Hello, I have a problem regarding calculation of Cumulative Incidence Function. The event of interest is failure of bone-marrow transplantation, which may occur due to relapse or death in remission. The data set that I have consists of- lifetime variable, two indicator variables-one for relapse and one for death in remission, and the other variables are donor typ...

Hi In the controller I have the function def convert_sd_ticket_to_incident ServiceDeskTicket.close_ticket(sd_id,created_by_id,modified_by_id) Incident.convert_sd_ticket_to_incident(sd_id,created_by_id,modified_by_id) end Now in the ServiceDeskTicket model def self.close_ticket(sd_id,created_by,modified_by) begin ActiveRecord::Base.transaction...