similar to: Off Topic: Re: Calculating NOEL using R and logistic regression - Toxicology

Hello, I used the glm function in R to fit a dose-response relationship and then have been using dose.p to calculate the LC50, however I would like to calculate the NOEL (no observed effect level), ie the lowest dose above which responses start occurring. Does anyone know how to do this? [[alternative HTML version deleted]]

Hi all, I'm having a strange error that prevents me from installing new packages, or updating packages after reinstalling. The error message is Warning: unable to move temporary installation ‘C:\Program Files\R\R-2.14.2\library\file15045004ac2\sandwich’ to ‘C:\Program Files\R\R-2.14.2\library\sandwich’ for one of the packages that is failing to install/update. This error started happening

I have a comma delimited text file in which many columns of numbers are also quoted and have commas as well as decimals. I was surprised to find read.csv() didn't import this seamlessly, even after messing around with the colClasses argument. I did find a solution to convert the character strings after reading them in, but wonder if there isn't a better one I overlooked. test =

Greetings, I have a question that I'd like to get input on. I have a classic toxicology study where I artificially fertilized and exposed embryos to a chemical and counted defects. In addition, I kept track of male-female pairs that I used to artificially fertilize and generate embryos with. I need to use logistic regression to model the response, but also check that the genetics of the

Dear all, I've been having problems running update() to re-fit a cpglm model inside a function (as in the code below). The solution is probably simple, but I'm stuck. If anyone could help, I'd greatly appreciate it. Regards, Rubem ## R code library(cplm) ## Data simulation period<-factor(1:4)                        herd<-factor(1:50)  

Hello, I would like to use a permanova to analyze my repeated measures, multivariate data set using the package vegan (adonis function). I have several explanatory variables and many clinical biochemistry/hematological response variables (all continuous, non-normally distributed). My explanatory variables are: ID - 29 levels - individual subjects Treatment - 2 levels Sex - 2 levels Time -

Hi, I appreciate you're busy so I just wanted to follow up on the email I sent you the other day, copy included below for reference. On Sat, Nov 12, 2016 at 12:40 PM, Danielle James <danielle.james at whoishostingthismail.com> wrote: Hi, I just wanted to let you know about a link that seems to be broken on this page http://www.hdt-project.org/. Here's the link

I am trying to write a wrapper for the last example in help(cv.glm) that deals with leave-one-out-cross-validation (LOOCV) for a logistic model. This wrapper will be used as part of a bigger program. Here is my wrapper funtion : logistic.LOOCV.err <- function( formu=NULL, data=NULL ){ cost.fn <- function(cl, pred) mean( abs(cl-pred) > 0.5 ) glmfit <- glm(

I have been using R the past couple of years to run models on data we are collecting. I recently got a new computer and updated to a new version of R (2.60 -> 2.90). Since the update, I cannot get my syntax to run. I have tried copying the file it is looking for into many different directories to try and run it. In the last version, I found that it was easiest if the file was copied into the R

Hi all, I've now got a problem with some modified lmer code (function lmer1 pasted at end) - I've made only three changes to the lmer code (marked), and I'm not really looking for comments on this function, but would like to know why execution of the following commands that use it almost invariably (but not quite predictably) leads to the R session terminating. Here's the command

Hi, New to Wine but have been using Ubuntu (Lucid) for a year now with Wine 1.3.0. My problem is getting a program called URL TYRE (http://www.janboersma.nl/gett/index.php) to run. This is a nice program that allows you to create GPS road routes using Google Maps. The first time I try to start it TYRE there is some HDD activity but it doesn't run. If I try a second time to launch it there

Dear all, I would like to Boostrap the stepAIC() procedure from package MASS for variety of model objects, i.e., fn <- function(object, data, B = 2){ n <- nrow(data) res <- vector(mode = "list", length = B) index <- sample(n, n * B, replace = TRUE) dim(index) <- c(n, B) for (i in 1:B) { up.obj <- update(object, data = data[index[, i], ])

Sorry, spoke too quickly! The problem was caused by my editing the MkRules file with an editor that replaces tabs with spaces. Make evidently requires that all lines start with a tab (see http://www.delorie.com/djgpp/v2faq/faq22_17.html). Perhaps worth warning windoze users like myself of such perils in readme.packages or somewhere similar! Now I have a new problem, but it looks fixable.

Hello all- I am doing a survival analysis for two species of invasive plants I outplanted to edges and interiors of island and mainland sites in a local reservoir. I am using the KM estimate and had no problem doing survdiff for my data using the following code: S4<-Surv(outplant$SurvTime, outplant$StatusD6) diff4=survdiff(S4 ~ outplant$Species+outplant$SiteType+outplant$EdgInt) diff4

Running dovecot 1.2.12 on ubuntu 10.10. Here is the problem: I have mysql backend for expire. I am using the dovecot-postfix package which creates /etc/dovecot/conf.d/01-mail-stack-delivery.conf which looks like this: # Some general options protocols = imap imaps managesieve disable_plaintext_auth = no ssl = yes ssl_cert_file = /etc/ssl/certs/ssl-mail.pem ssl_key_file =

On Fri, 3 May 2002 drew.tyre@csiro.au wrote: > Full_Name: Drew Tyre > Version: 1.4.0 > OS: Win 2000 > Submission from: (NULL) (140.253.74.225) > > > When creating a rather complex plot, the resulting figures 'shrink' when either > copied to the clipboard or saved as a metafile using the GUI menu. Also saving > as postscript or pdf causes the same problem, as

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 model to binary data. The response variable is a factor (0 or 1) and all predictors are continuous variables. The main predictor is LT (I expect a logistic relation between LT and the probability of being mature) and the other are variables I expect to modify this relation. I want to test if all predictors contribute significantly for the fit or not I fit the full

Hello, I am new to R. I often collect data at multiple sites and need to create separate graphs (such as scatterplots or histograms) of specific variables for each site. I have tried to do this by splitting the data frame and then using lapply, but it seems that the graphing commands cannot be called as functions. Here is a sample of my data, called "seeddist2": site

Hi All, I am trying to build a couple of packages from source, one of my own and others. I have had this working under R1.5.0, now trying to set up on a new machine under R1.6.0. Starting from various places, I am getting an error in MkRules that I do not understand. C:\PROGRA~1\R\rw1060\src\gnuwin32>make libR.a libRblas.a MkRules:91: *** missing separator. Stop. The example is from the