similar to: Fwd: Re: Finney's fiducial confidence intervals of LD50

I am working with Probit regression (I cannot switch to logit) can anybody help me in finding out how to obtain with R Finney's fiducial confidence intervals for the levels of the predictor (Dose) needed to produce a proportion of 50% of responses(LD50, ED50 etc.)? If the Pearson chi-square goodness-of-fit test is significant (by default), a heterogeneity factor should be used to calculate

It appears that the chromatic adaptation feature of `grDevices::convertColor`is broken, and likely has been for many years.? While a little surprising, it is an obscure enough feature that there is some possibility this is actually broken, as opposed to user error on my part.? If it turns out to the latter, I apologize in advance for spamming this list. Consider: ??? rgb.in <-

Classification: UNCLASSIFIED Caveats: NONE A similar question has been posted in the past but never answered. My question is this: for probit analysis, how do you program a 95% confidence interval for the LD50 (or LC50, ec50, etc.), including a heterogeneity factor as written about in "Probit Analysis" by Finney(1971)? The heterogeneity factor comes into play through the chi-squared

Would it be possible to set the ringtone based on the number that was dialed? Example of what the goal is: Dial Denver number Incoming calls ring with ringtone 1 Dial main number Incoming calls ring with ringtone 2 We are currently using Digium D40, D50, D70 phones. -------------- next part -------------- An HTML attachment was scrubbed... URL:

Dear R users, I want to creat a group of data frames, such as: d1<-data.frame(x=factor(),y=numeric(),...) d2<-data.frame(x=factor(),y=numeric(),...) d3<-data.frame(x=factor(),y=numeric(),...) : : d50<-data.frame(x=factor(),y=numeric()) How can I use the loop to simply my procedure? Thanks advance Ted -- View this message in context:

R users, I have df like this a <- data.frame( indx = 1:20, var1 = rep(c("I20", "I40", "A50", "B60"), each=5), var1_lab= rep(c("cat", "dog", "mouse", "horse"), each=5), var2 = rep(c("B20", "X40", "D50",

Hi, I just got a T1 interface for a AT&T (became Lucent) System75 (uses same cards as Definity). I would take a crack at making a cable but can't determine the pinout for a cable and it is not apparent from the board. Asking you guys makes sense as one of you may have one of these systems. The cable has a amphenol male D50 connector on one end and probably a rj45 on the other. I also

I'm trying to solve a mystery for the last couple of days. I have a mix of D70, D50 and D40 behind NAT. Server is in a colocation, not a VPS. For several years, everything was working fine, no issues. A few days ago I started having problems at one particular site. NO CHANGES have been made to this office network - same router, switch and internet provider. No new equipment added or

Hi, Can anybody advise me how to debug this crash? it happens when I try to start a 32 bit kernel on my 64 bit Xen system, the kernel is my own using opensuse 11.1 rc1 2.6.27.7 xen patches, the kernel works as 32 and 64 bit dom0, and as 64 bit domU PV, is there a technique to find out where in the kernel source the crash happens? i.e what line in the source code? (XEN) mm.c:707:d49 Error getting

Hi, I'm not a programmer, so I have a question about R functions, I run the Random Forest regression models, but I would like to run the random forest model 1000 times with different random sample set. to check the uncertainty of the regression model estimates. exemple of data: ################################# table= all Y: all$AGB X variables: Variables=as.matrix(all[, c( "min",

Dear R-list, I have a data set from 20 pigs, each of which is tested at crossed 9 doses (logdose -4:4) and 3 skin treatment substances when exposed to a standard polluted environment. So there are 27 patches on each pig. The response is irritation=yes/no. I want to determine "equally effective 50% doses" (similar to old LD50), and to test the treatments against each other. I am looking

Hi, I'm not yet familar with GLM and still learning. How can I perform a BNL (to estimate LDp values) with matrixes like this (N indicates the observed objects): data from: Kerr D, Meador J. 1996. "Modeling Dose Response Using Generalized Linear Models". Environ Toxicol Chem 15(3) conc respond n 0 0 20 0 1 19 0 1 20 0 0 20 13.5 5 23 13.5 2 20 29 9 20 29 6 20 53 12 20 53 15 20 85

I'm writing a function and keep getting the following error message. myfunc <- function(lst) { lst <- list(roots = c("car insurance", "auto insurance"), roots2 = c("insurance"), prefix = c("cheap", "budget"), prefix2 = c("low cost"), suffix = c("quote", "quotes"), suffix2 = c("rate",

Hi All! I am desperately needing some help figuring out how to calculate LD50 with a GLMM (probit link) or, more importantly, the standard error of the LD50. I conducted a cold temperature experiment and am trying to assess after how long 50% of the insects had died (I had 3 different instars (non significant fixed effect) and several different blocks (I did 4 replicates at a time)=

Hi, I have recently been attempting to find the LD50 from two predicted fits (For male and females) in a Generalised linear model which models the effect of both sex + logdose (and sex*logdose interaction) on proportion survival (formula = y ~ ldose * sex, family = "binomial", data = dat (y is the survival data)). I can obtain the LD50 for females using the dose.p() command in the MASS

Quiero comparar varias dosis letales 50% (LD50) usando análisis probit. He seguido un ejemplo que viene en paquete DRC, pero no obtengo el resultado esperado. Lo que quiero es saber si las LD50s, son diferentes y si la diferencias son estadísticamente significativas. Gracias de antemano. José Arturo e-mail. jafarfan@uady.mx <grejon@uady.mx> e-mail alterno. jafarfan@gmail.com

Hello, I am attempting to write a script that adds error bars to a barchart. I basing my attempt heavily on the following thread: http://tolstoy.newcastle.edu.au/R/e2/help/06/10/2791.html I can't seem to get around the problem that was discussed in the thread. The following example should illustrate my problem. Sorry about the messy example but I am 1) trying to make it as close as possible

R users, My intention is to take factors out of DF, create list of tables and export these tables separately using write.table or sink function. write.table writes tables out as DF:s, should I use sink instead? Here is my example: a <- data.frame( indx = 1:20, var1 = rep(c("I20", "I40", "A50", "B60"), each=5),

Today I attempted to upgrade to S10_U2 and migrate some mirrored UFS SVM partitions to ZFS. I used Live Upgrade to migrate from U1 to U2 and that went without a hitch on my SunBlade 2000. And the initial conversion of one side of the UFS mirrors to a ZFS pool and subsequent data migration went fine. However, when I attempted to attach the second side mirrors as a mirror of the ZFS pool, all

Dear All, I have a binomial response with one continuous predictor (d) and one factor (g) (8 levels dummy-coded). glm(resp~d*g, data, family=binomial) Y=b0+b1*X1+b2*X2 ... b7*X7 how can I get the inflection point per group, e.g., P(d)=.5 I would be grateful for any help. Thanks in advance, Ren?