turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- JMP User Community
- :
- Discussions
- :
- Dose Response Probit analysis

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Oct 28, 2014 6:20 AM
(2988 views)

Hi group,

I am trying to reproduce a simple dose response analysis that I had performed in SPSS. So in JMP i am using the "Fit model" module, selecting for Generalized linear models, binomial distribution and probit link.

I have three variables, number of assays performed (10 assays), number of assays with positive result and concentration (from 0 to 36)

0 10 2

1 10 4

2 10 9

6 10 16

8 10 25

9 10 36

0 10 0

I select number of assays performed as "Y" variable, number of assays with positive result as "Freq" and concentration in logs as variable to model effects on. But when I hit run I get this message in the results:

"Warning: Hessian not positive definite. Singularity is likely due to a poor fit of the model, or linear dependencies among the model covariates. Fit and result are of questionable value"

And the analysis is not performed.

Can anyone share some comments?

Thanks,

Dave

1 ACCEPTED SOLUTION

Accepted Solutions

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Oct 28, 2014 9:47 AM
(3604 views)

Solution

Dave,

Your "number of assays with positive result" needs to go in as a "Y" variable also. The Frequency role is for indicating how many times a row in the data table actually occurred. In this case, the GLM model expects the number of assays performed and the number of positive results as "separate" Y values. The analysis will link them together.

One thing you did not mention that is often done for this type of work is to take the log of the concentration. I did this in the model dialog shown above. Highlight the concentration model effect, click the transform button and choose log. Results are posted here.

Dan

Dan Obermiller

6 REPLIES

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Oct 28, 2014 9:47 AM
(3605 views)

Dave,

Your "number of assays with positive result" needs to go in as a "Y" variable also. The Frequency role is for indicating how many times a row in the data table actually occurred. In this case, the GLM model expects the number of assays performed and the number of positive results as "separate" Y values. The analysis will link them together.

One thing you did not mention that is often done for this type of work is to take the log of the concentration. I did this in the model dialog shown above. Highlight the concentration model effect, click the transform button and choose log. Results are posted here.

Dan

Dan Obermiller

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Oct 29, 2014 12:57 AM
(1802 views)

Thanks a lot DanO, that is a very clear answer. It is interesting that JMP knows they are linked.

Regarding the log transformation, yes I will take that into account (it was indeed mentioned in my original post, although pretty vague I must admit!)

Best,

Dave

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Oct 29, 2014 9:35 AM
(1802 views)

Dear DanO,

sorry for bothering you again, but changing the numbers to the demo, I seem to be getting wrong results. Please have a look at this dataset, with same variables but different numbers. For the three highest concentrations I get 100% positives

Number of positive results Number of assays performed Concentration

0 1140 0.01

241 1140 0.1

590 1140 0.5

963 1140 1

1134 1140 5

1140 1140 10

1140 1140 50

1140 1140 100

I run the Probit as you suggested

An after running I get a strange result which doesn't make sense to me. Instead of doing percentages of corrects on the Y, it seems to be doing something different. What am I doing wrong?

Thanks ever so much for your help

Dave

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Oct 29, 2014 1:11 PM
(1802 views)

Been following this thread with interest.

From the JMP docs:

If your data is summarized as frequencies of events and trials, specify two continuous

columns** in this order**: a count of the number of successes, and a count of the number of

trials. Alternatively, you can specify the number of failures instead of successes.

(my emphasis).

So if you switch the order the 2 y variables are entered I think you will see more of what you expected especially if you change the x axis to be logarithmic.

That said, I can't explain why the order doesn't seem to matter with your original data.

Michael

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Oct 29, 2014 1:31 PM
(1802 views)

The 0 10 0 last row from the original data for some reason insulates the analysis from needing the correct order. If you leave it out then order does matter.

Michael

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Nov 2, 2014 11:34 PM
(1802 views)

That is indeed the right order to feed in the variables. Thanks for pointing that out, Michael.

Dave