similar to: R-alpha: fitted = 0 of 1 in logistic regression

Dear All I am attempting some power analyses, based on simulated data. My experimental set up is thus: Bleach: main effect, three levels (control, med, high), Fixed. Temp: main effect, two levels (cold, hot), Fixed. Main effect interactions, six levels (fixed) For each main-effect combination I have three replicates. Within each replicate I can take varying numbers of measurements (response

1) hist() does not take NA's. Incompatible with Splus, probably just a bug? 2) I do wish we could somehow get rid of the misfeatures of indexing with logical NA's: > table(juul$menarche,juul$tanner) I II III IV V No 221 43 32 14 2 Yes 1 1 5 26 202 > juul$menarche=="Yes"&juul$tanner=="I",] ...and you find yourself with a listing of 477

Dear R-users, After running a logistic regression, I need to calculate OR by exponentiating the coefficient, and then I need the 95% CL for the OR as well. For the following example (taken from P. Dalaagard's book), what would be the most straightforward method of getting what I need? Could anyone enlight me please? Thank you! Lucho > summary(glm(menarche~age,binomial)) Call:

Sorry if this is a duplicate: This is a re-post because the pdf's mentioned below did not go through. Hello, I'm new'ish to R, and very new to glm. I've read a lot about my issue: Warning message: glm.fit: fitted probabilities numerically 0 or 1 occurred ...including: http://tolstoy.newcastle.edu.au/R/help/05/07/7759.html

Hello all, A really noddy question for you all: I''m trying without success to do some subhypothesis testing. Using simple anova model, with a toy dataset from a book. I have four factors A,B,C,D, and wish to test mu_C = mu_D. This is what I have tried: > contrasts(infants$group,how.many=1) <- c(0,0,1,-1) > contrasts(infants$group) [,1] A 0 B 0 C 1

Hi Flac developers! I'm using the seekable_stream_decoder API of libFLAC 1.03 and I'm stuck with a mysterious problem: my FLAC plugin works nicely as long as the seek_absolute() function is never used, but as soon as I seek around in a sample (ie. call seek_absolute() with an offset != 0), I cannot detect the end of a file anymore. decoder_get_state() never returns

I simply put, > NEVER=subset(infants$bwt,ISNO1) > UNTILPREGNANT=subset(infants$bwt, ISNO2) > ONCENOTNOW=subset(infants$bwt, ISNO3) and I wanna combine those three. I do it like ISNO=NEVER&UNTILPREGNANT&ONCENOTNOW and R tells me 1: In NEVER & UNTILPREGNANT : longer object length is not a multiple of shorter object length 2: In NEVER & UNTILPREGNANT & ONCENOTNOW

Hi, r-gurus: I happened to have a question in my work: I have a dataset, which has only one dimention, like 0.99037297527605 0.991179836732708 0.995635340631367 0.997186769599305 0.991632565640424 0.984047197106486 0.99225943762649 1.00555642128421 0.993725402926564 .... the data is saved in a file called f392.txt. I used the following codes to play around :)

The x11() window can be a nuisance to have popping up at startup (esp. on small screens) when you're not working with graphics. However, currently you can't get rid of it without modifying the systemwide Rprofile. Current logic is: Run $RHOME/library/Rprofile if ./.Rprofile exists run it else if $HOME/.Rprofile exists run that endif I think it should be Run $RHOME/library/Rsetup if

Hello, In a study exploring transgenerational transmission of anxiety disorder we investigate whether infants react to experimentally induced mood changes of their mothers. We measured the time that an infant needed to cross a cliff (=crossing time) depending on whether his mother had previously undergone a mood induction (treatment) or not (control). The treatment is thus a

Ben Bolstad experienced something odd while running R CMD check over a package of ours using r-devel (yesterday's fresh). He tracked down the problem to a vignette calling 'library(MASS)'. Please kindly ignore if this is a known (and temporary) issue. Laurent ----- Forwarded message from Ben Bolstad <bolstad@stat.berkeley.edu> ----- Return-Path:

http://bugzilla.mindrot.org/show_bug.cgi?id=946 dtucker at zip.com.au changed: What |Removed |Added ---------------------------------------------------------------------------- Status|NEW |RESOLVED Resolution| |WORKSFORME ------- Additional Comments From dtucker at zip.com.au

Hi, I noticed that mount_ext2 isn't "connected" to the normal "buildworld" but sys/*/ext2fs seems to. Building and using it from recent RELENG_7 sources works for me (having ext2fs loaded as module in GENERIC). I mounted an ext3 rw, edited grub's menu.lst, wrote and umounted without any problems. It seems recent bugs got pathched in January

Hi. I used Samba 2.2.x on my server with LDAP backend. I also had linux users in LDAP (pam_ldap). I had a separate tree for Samba and for linux users. Because I had problems browsing the user lists when trying to set permisions on a share on Win2k and WinXP members of domain and also for other reasons (printer drivers,...) I started an upgrade process to Samba 3.0 . After playing a lot with

Johnny, Thought I'd try these out for you. Looks like a prerequisite is missing? yum --enablerepo=c5-testing install java-1.6.0-openjdk\* <snip the usual gobbledygook> 330 packages excluded due to repository priority protections Parsing package install arguments Resolving Dependencies --> Populating transaction set with selected packages. Please wait. ---> Package

Hello! I'm using the seekable stream decoder API of libFLAC 1.02 and I think that I found a possible bug in process_metadata(). The problem is as follows: I have a file which isn't a FLAC sample (it's actually an ACE archive) where process_metadata() returns TRUE. And even worse, the metadataCallback() is never called (which means process_metadata() succeeds, although no metadata is

https://bugzilla.mindrot.org/show_bug.cgi?id=2483 Bug ID: 2483 Summary: use AI_ADDRCONFIG flag for getaddrinfo() hints on Solaris Product: Portable OpenSSH Version: 6.8p1 Hardware: All OS: Solaris Status: NEW Severity: enhancement Priority: P5 Component:

Hello! I would like to send some bugs and bugfix(es): 1. There is an encoding bug in the beta 3 (downloaded exe from the vorbis.com) and in the latest CVS (jan. 10.) too (compiled with MSVC 6.0): If the file contains some non-audible noise at the begin or at the end, the encoder/decoder makes some large clips in the output file, if the bitrate is 350kbit/s . If I use 256kbit/s or

I have a logistic model fitted with the following R function: glmfit<-glm(formula, data, family=binomial) A reasonable cutoff value in order to get a good data classification (or confusion matrix) with the fitted model is 0.2 instead of the mostly used 0.5. And I want to use the `cv.glm` function with the fitted model: cv.glm(data, glmfit, cost, K) Since the response in the fitted

I am fitting a logistic binomial model of the form glm(y ~ a*x,family=binomial) where a is a factor (with 5 levels) and x is a continuous predictor. To assess how much ``impact'' x has, I want to compare the fitted success probability when x = its maximum value with the fitted probability when x = its mean value. (The mean and the max are to be taken by level of the factor