Skip to content

summary(m) does not work with error: Error in solve.default(object$fit$hessian) : Lapack routine dgesv: system is exactly singular: U[1,1] = 0 #36

@lijuanyanbasf

Description

@lijuanyanbasf

dataframe <- data.frame(conc=c(0, 0.944, 2.18, 4.14, 8.37, 16.1), total=c(160, 80, 80, 80, 80, 80), response=c(3, 1, 1, 5, 8, 80))

1) summary(m) does not work as below (increasing trend)

m <- drm(response/total ~ conc,
weights=total,
data=dataframe,
fct = W1.4(),
type = "binomial")
plot(m)
summary(m)

**# Error in solve.default(object$fit$hessian) :

Lapack routine dgesv: system is exactly singular: U[1,1] = 0**

2) However, summary(m) works if I model the non-responding fractions (decreasing trend)

m <- drm((total-response)/total ~ conc,
weights=total,
data=dataframe,
fct = W1.4(),
type = "binomial")
plot(m)
summary(m)

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions