[#6498] Model matrix constructed different with link = identity and link = logit

Date:
2017-03-23 08:54 Priority:
3

State:
Closed

Submitted by:
Joris Meys (jorismeys) Assigned to:
Nobody (None)

Hardware:
None Product:
None

Operating System:
None Component:
None

Version:
None Severity:
None

Resolution:
Fixed

URL:

Summary:
Model matrix constructed different with link = identity and link = logit

Detailed description
Ik moest een PIM fitten met het PIM package en denk een inconsistentie gevonden te hebben. Het was een model met een 0/1 dummy. Als ik een logit link functie gebruik is de regressor Z=X^-X (notatie uit PIM paper), maar met een identity link is het Z=X-X^.

Comments:

Message
Date: 2017-03-23 09:42 Sender: Joris Meys Problem is the fact that using a simple formula with the identity link doesn't add an intercept to the model matrix. This needs to be solved.
Date: 2017-03-23 09:24 Sender: Joris Meys Code : > m<-pim(FEV~Smoke,link="identity",data=FEVSubset) > summary(m) pim.summary of following model : FEV ~ Smoke Type: difference Link: identity Estimate Std. Error z value Pr(>\|z\|) Smoke -0.24064 0.09116 -2.64 0.0083 --- Signif. codes: 0 '' 0.001 '' 0.01 '' 0.05 '.' 0.1 ' ' 1 Null hypothesis: b = 0 > m<-pim(FEV~Smoke,link="logit",data=FEVSubset) > summary(m) pim.summary of following model : FEV ~ Smoke Type: difference Link: logit Estimate Std. Error z value Pr(>\|z\|) Smoke 1.1492 0.4989 2.303 0.0213 *

Field	Old Value	Date	By
Resolution	None	2018-02-08 13:05	jorismeys
status_id	Open	2018-02-08 13:05	jorismeys
close_date	None	2018-02-08 13:05	jorismeys