[R] Probit Analysis

Wed Jan 1 13:09:03 CET 2003

A better option for a confidence interval for LD50
might be Fiellers method. 

To see how to get the numbers from the object returned from dose.p, use for
example
str(object) or names(object).

Kjetil Halvorsen

-------Original Message-------

From: Ngayee J Law
Date: Tuesday, December 31, 2002 7:17:39 PM
To: apjaworski at mmm.com
Cc: r-help at stat.math.ethz.ch
Subject: Re: [R] Probit Analysis

Hi Andy,

Thanks! glm and dose.p can do exaclty what I want to do. Just a quick
question,
how can I extract the numbers in the glm.dose object created by the dose.p?
My
goal is to calculate a 95% CI for the 'LD50'.

- Jacqueline

apjaworski at mm 
m.com To: "Ngayee J Law" <Ngayee.Law at celeradiagnostics.com> 
cc: r-help at stat.math.ethz.ch 
12/31/2002 Subject: Re: [R] Probit Analysis 
02:47 PM 

Jacqueline,

This should be a simple application of GLM. Check the help pages for glm,
family and dose.p in library MASS. Look at the example at the bottom of
the dose.p help page - it does almost exactly what you need.

Hope this helps,

Andy

__________________________________
Andy Jaworski
Engineering Systems Technology Center
3M Center, 518-1-01
St. Paul, MN 55144-1000
-----
E-mail: apjaworski at mmm.com
Tel: (651) 733-6092
Fax: (651) 736-3122

|---------+---------------------------------->
| | "Ngayee J Law" |
| | <Ngayee.Law at celeradiagn|
| | ostics.com> |
| | Sent by: |
| | r-help-admin at stat.math.|
| | ethz.ch |
| | |
| | |
| | 12/31/2002 16:17 |
| | |
|---------+---------------------------------->
>
-----------------------------------------------------------------------------
-----------------------------------------------|

|
|
| To: r-help at stat.math.ethz.ch
|
| cc:
|
| Subject: [R] Probit Analysis
|
>
-----------------------------------------------------------------------------
-----------------------------------------------|

Hello all,

I have a very simple set of data and I would like to analyze them with
probit analysis.
The data are:

X Event Trial
100 8 8
75 8 8
50 6 8
25 4 8
10 2 8
0 0 8

I want to estimate the value of X that will give a 95% hit rate
(Event/Trial) and the corresponding
95% CI. Anyone can offer some help? Thanks!!

- Jacqueline

______________________________________________
R-help at stat.math.ethz.ch mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help

______________________________________________
R-help at stat.math.ethz.ch mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help
.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: https://stat.ethz.ch/pipermail/r-help/attachments/20030101/a3a7fab6/attachment.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: image/gif
Size: 494 bytes
Desc: not available
Url : https://stat.ethz.ch/pipermail/r-help/attachments/20030101/a3a7fab6/attachment.gif
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: image/jpeg
Size: 1431 bytes
Desc: not available
Url : https://stat.ethz.ch/pipermail/r-help/attachments/20030101/a3a7fab6/attachment.jpe