search for: infert

Hi R-users, I have the following problem with CrossTable function within ?gmodels? package: the output of the function (format ?spss? and asresid=T) can not be stored within another object. For example: >library(gmodels) >data(infert, package = "datasets") > CrossTable(infert$education, infert$induced)->aa # the function prints everything ok on the screen > aa # works just fine $t y x 0 1 2 0-5yrs 4 2 6 6-11yrs 78 27 15 12+ yrs 61 39 16 ?????. But when I call >CrossTable(infe...

...seem to get my brain around it. If I want to produce a set of tables for all the variables in my data, how can I do that without having to type them into the table command one by one. So, I would like to use (t? s? r?)apply to use one command instead of the following set of table commands: data(infert, package = "datasets") attach(infert) table.education<-table(education) table.age<-table(age) table.parity<-table(parity) etc. To make matters worse, what I subsequently need is the chi-square for each and all of the pairs of variables. Such as: chi.education.age<-chisq.t...

Dear R-helpers, I am very new to R and trying to run the conditional logit model using "clogit " command. I have more than 4000 observations in my dataset and try to predict the dependent variable from 14 independent variables. My command is as follows clmtest1 <- clogit(Pin~Income+Bus+Pop+Urbpro+Health+Student+Grad+NE+NW+NCC+SCC+CH+SE+MRD+strata(IDD),data=clmdata) However, it

Hi, I am playing with clogit(case~spontaneous+induced+strata(stratum),data=infert) from clogit help file. This line works. 1. But, why strata(stratum) doesn't have a coefficient like spontaneous and induced? 2. When I remove strata(stratum) from the command, this function seems to keep running forever. Why? 3. I think the equation for clogit looks like P=1/(1+ exp(-1*(...

--begin inclusion -- I have a matched-case control dataset that I'm using conditional logistic regression (clogit in survival) to analyze. I'm trying to conduct k-folds cross validation on my top models but all of the packages I can find (CVbinary in DAAG, KVX) won't work with clogit models. Is there any easy way to do this in R? -end inclusion -- The clogit funciton is simply a

Hi all: My data has 3 variables: age(3levels : <30y=1 30-50y=2, >50y=3) gender(Male=0, Female=1) CD4 cell count(raw lab measurement) y(1:death 0:alive) I perform logistic regression to find out the factors that influence y. result<-glm(y ~ factor(age) + factor(gender) + CD4,family = binomial) >From the result,I can get OR(Odds Ratio) of gender via exp(Estimate of Female,

Hi All, I am using "glm" function to build logistic regression. I noticed that glm function glm function is computing many other statistics which are not required for our analysis. As our dataset is very big and we have to run logistic regression on several samples the run time drastically increases if all those statistics are computed. Is these any way to skip computation in glm

I would like to know how R computes the probability of an event in a logistic model (P(y=1)) from the score s, linear combination of x and beta. I noticed that there are differences (small, less than e-16) between the fitting values automatically computed in the glm procedure by R, and the values "manually" computed by me applying the reverse formula p=e^s/(1+e^s); moreover I noticed

Dear List, With the R 2.12.0 addition of table methods for points(), dotchart() struggles with tables. I found several possible solutions, but it is beyond my skill to decide what is "best". Here is a small example: ############################################# x <- table(infert$education) y <- 1:3L dotchart(x) # error about incorrect plot type ## moving closer to the cause, dotchart() essentially does plot(x = c(0, 150), y = c(0, 5), type = "n") points(x, y) # error, because arg #2 is the type ## One option would be to explicitly name the arguments in dotc...

I have a case-control study that I'm analysing using the conditional logistic regression function clogit from the survival package. I would like to calculate the AICc of the models I fit using clogit. I have a variety of scripts that can calculate AICc for models with a logLik method, but clogit does not appear to use this method. Is there a way I can calculate AICc from clogit in R? Many

Dear R-help, I'm doing conditional logit regression for a discrete choice model. I want to know whether there's a way to get the fitted probabilities. In Stata, "predict" works for clogit, but it seems that in R "predict" does not. Thank you very much! Best wishes. Sincerely, Min -- Min Chen Graduate Student Department of Agricultural,

...an Cities faithful Old Faithful Geyser Data fdeaths (UKLungDeaths) Monthly Deaths from Lung Diseases in the UK freeny Freeny's Revenue Data freeny.x (freeny) Freeny's Revenue Data freeny.y (freeny) Freeny's Revenue Data infert Infertility after Spontaneous and Induced Abortion iris Edgar Anderson's Iris Data iris3 Edgar Anderson's Iris Data islands Areas of the World's Major Landmasses ldeaths (UKLungDeaths)...

I have a factor vector of subject races (Asian, Black, Hispanic, White; n=30) that I want to plot with a Cleveland dotplot or dotchart. I tried the following in R2.12.1 : > dotchart(table(school$Race)) Error in plot.xy(xy.coords(x, y), type = type, ...) : invalid plot type Using the same data set in R2.11.1 the operation succeeded (I tried several variations to be sure): >

i have a data frame(dat) which has many variables.and i use the following script to get the crosstable. >danx2<-c("x1.1","x1.2","x1.3","x1.4","x1.5","x2","x4","x5","x6","x7","x8.1","x8.2","x8.3","x8.4","x11",

...x. Dimnames are propagated in a sensible fashion. o New chisq.test(.) and prop.test() from Kurt Hornik. o New read.fwf(.) for reading fixed width format (KH). o New str(.) [alternative to summary(.) for programmers] (MM). o New example data sets "esoph", "infert" and "anscombe" (TL), "iris3" (KH) and "stackloss" (MM). o source(.) has several new arguments, notably ``echo = FALSE''. This is applied in the new function demo(.) which runs all the code in demos/ (but dynload). o strheight(.) is new,...

...x. Dimnames are propagated in a sensible fashion. o New chisq.test(.) and prop.test() from Kurt Hornik. o New read.fwf(.) for reading fixed width format (KH). o New str(.) [alternative to summary(.) for programmers] (MM). o New example data sets "esoph", "infert" and "anscombe" (TL), "iris3" (KH) and "stackloss" (MM). o source(.) has several new arguments, notably ``echo = FALSE''. This is applied in the new function demo(.) which runs all the code in demos/ (but dynload). o strheight(.) is new,...

...b/R/library/base/R-ex/hist.R action: create D: file: /usr/local/lib/R/library/base/R-ex/hsv.R action: create D: file: /usr/local/lib/R/library/base/R-ex/ifelse.R action: create D: file: /usr/local/lib/R/library/base/R-ex/image.R action: create D: file: /usr/local/lib/R/library/base/R-ex/infert.R action: create D: file: /usr/local/lib/R/library/base/R-ex/influence.measures.R action: create D: file: /usr/local/lib/R/library/base/R-ex/interaction.R action: create D: file: /usr/local/lib/R/library/base/R-ex/interactive.R action: create D: file: /usr/local/lib/R/library/base/R-ex/...

...text html latex example identify text html latex ifelse text html latex example image text html latex example index.search text html latex infert text html latex example influence.measures text html latex example integer text html latex integrate text html latex example interaction text ht...