View tutorialsXLSTAT-Dose is an Excel add-in which has been developed for dose analysis in the chemical and pharmaceutical industries. XLSTAT-Dose is an extension of the XLSTAT-Pro logistic regression for these industries.
- Dose effects:
- Models: Logit, Probit, Gompertz, complementary Log-log
- Computes effect doses from ED01 to ED99 (ED50 et ED90 included)
- Option to take into account the natural mortality
- Option to automatically take the Log of the doses
- Parallel 4 parameter fit:
- Fitting of model a-(d-a)/(1+(x/c)^b)
- Comparison of the model for a standard sample and a sample of interest
- Fitting of the model on the sample of interest, taking account constraints from the standard sample
- Optional removal of outliers (Dixon's test)
- Parallelism test (Fisher)
- Plot of parallel curves
You can find tutorials that explains how XLSTAT-Dose works here.
A trial version of XLSTAT-Dose is included in the main XLSTAT-Pro download.
Prices and ordering
For prices, on-line ordering and other purchasing information please go to our ordering page.
View a tutorialWhat is dose effect analysis
Logistic regression (Logit, Probit, complementary Log-log, Gompertz models) is used to model the impact of doses of chemical components (for example a medicine or phytosanitary product) on a binary phenomenon (healing, death).
Natural mortality in dose effect analysis
Natural mortality should be taken into account in order to model the phenomenon studied more accurately. Indeed, if we consider an experiment carried out on insects, certain will die because of the dose injected, and others from other phenomenon. None of these associated phenomena are relevant to the experiment concerning the effects of the dose but may be taken into account. If p is the probability from a logistic regression model corresponding only to the effect of the dose and if m is natural mortality, then the observed probability that the insect will succumb is:
P(obs) = m + (1- m) * p
Abbot's formula (Finney, 1971) is written as:
p = (P(obs) – m) / (1 – m)
The natural mortality m may be entered by the user as it is known from previous experiments, or is determined by XLSTAT.
Compute effect dose from from ED01 to ED99 including ED 50, ED 90
XLSTAT calculates ED 50 (or median dose), ED 90 and ED 99 doses which correspond to doses leading to an effect respectively on 50%, 90% and 99% of the population.
Results for dose effect analysis in XLSTAT
- Goodness of fit coefficients:
This table displays a series of statistics for the independent model (corresponding to the case where the linear combination of explanatory variables reduces to a constant) and for the adjusted model.
- Observations: The total number of observations taken into account (sum of the weights of the observations);
- Sum of weights: The total number of observations taken into account (sum of the weights of the observations multiplied by the weights in the regression);
- DF: Degrees of freedom;
- -2 Log(Like.): The logarithm of the likelihood function associated with the model;
- R² (McFadden): Coefficient, like the R2, between 0 and 1 which measures how well the model is adjusted. This coefficient is equal to 1 minus the ratio of the likelihood of the adjusted model to the likelihood of the independent model;
- R²(Cox and Snell): Coefficient, like the R2, between 0 and 1 which measures how well the model is adjusted. This coefficient is equal to 1 minus the ratio of the likelihood of the adjusted model to the likelihood of the independent model raised to the power 2/Sw, where Sw is the sum of weights.
- R²(Nagelkerke): Coefficient, like the R2, between 0 and 1 which measures how well the model is adjusted. This coefficient is equal to ratio of the R² of Cox and Snell, divided by 1 minus the likelihood of the independent model raised to the power 2/Sw;
- AIC: Akaike’s Information Criterion;
- SBC: Schwarz’s Bayesian Criterion.
- Test of the null hypothesis H0: Y=p0:
The H0 hypothesis corresponds to the independent model which gives probability p0 whatever the values of the explanatory variables. We seek to check if the adjusted model is significantly more powerful than this model. Three tests are available: the likelihood ratio test (-2 Log(Like.)), the Score test and the Wald test. The three statistics follow a Chi2 distribution whose degrees of freedom are shown.
- Type III analysis:
This table is only useful if there is more than one explanatory variable. Here, the adjusted model is tested against a test model where the variable in the row of the table in question has been removed. If the probability Pr > LR is less than a significance threshold which has been set (typically 0.05), then the contribution of the variable to the adjustment of the model is significant. Otherwise, it can be removed from the model.
- Model parameters:
The parameter estimate, corresponding standard deviation, Wald's Chi2, the corresponding p-value and the confidence interval are displayed for the constant and each variable of the model. If the corresponding option has been activated, the "profile likelihood" intervals are also displayed.
- Model equation:
The equation of the model is then displayed to make it easier to read or re-use the model.
- Standardized coefficients table:
The table of standardized coefficients (also called beta coefficients) is used to compare the relative weights of the variables. The higher the absolute value of a coefficient, the more important the weight of the corresponding variable. When the confidence interval around standardized coefficients has value 0 (this can be easily seen on the chart of normalized coefficients), the weight of a variable in the model is not significant.
- Predictions and residuals table:
The predictions and residuals table shows, for each observation, its weight, the value of the qualitative explanatory variable, if there is only one, the observed value of the dependent variable, the model's prediction, the same values divided by the weights, the standardized residuals and a confidence interval.
- Probability analysis table:
If only one quantitative variable has been selected, the probability analysis table allows to see to which value of the explanatory variable corresponds a given probability of success.
View a tutorialWhat is four-parameter parallel lines logistic regression
The four parameter logistic model writes:
y = a + (d -a) / [1 + (x / c)b] model (1)
where a, b, c, d are the parameters of the model, and where x corresponds to the explanatory variable and y to the response variable. a and d are parameters that respectively represent the lower and upper asymptotes, and b is the slope parameter. c is the abscissa of the mid-height point which ordinate is (a+b)/2. When a is lower than d, the curve decreases from d to a, and when a is greater than d, the curve increases from a to d.
The parallel lines four parameter logistic model writes:
y = a + (d -a) / [1 + (st * x / c1 + sp * x / c2)b] model (2)
where st is 1 if the observation comes from the standard sample (STD), and 0 if not, and where sp is 1 if the observation is from the sample of interest (SOI), and 0 if not. This is a constrained model because the observations corresponding to the standard sample influence the optimization of the values of a, b, and d. From the above writing of the model, one can understand that this model generates two parallel curves, which only difference is the positioning of the curve, the shift being given by (c2-c1). If c2 is greater than c1, the curve corresponding to the sample of interest is shifted to the right of the curve corresponding to the standard sample, and vice-versa.
XLSTAT allows you to fit:
- model (1) to a standard sample or to the sample of interest,
- model (2) to both STD and SOI at the same time.
If the Dixon’s test option is activated, XLSTAT can test for each sample (STD and SOI) if some outliers influence too much the fit of the model. In the (a) case, a Dixon’s test is performed once the model (1) is fitted. If an outlier is detected, it is removed, and the model is fitted again, and so on, until no outlier is detected. In the (b) case, we first perform a Dixon’s test on the STD, than on the SOI, and then, the model (2) is fitted on the merge of the two samples, without the outliers.
In the (b) case, and if the sum of the sample sizes is greater than 9, a Fisher’s F test is performed to detect if the a, b and d parameters obtained for both models with the model (1) are not significantly different from those obtained with model (2).
Results for four-parameter parallel lines logistic regressionin XLSTAT
- Fisher's test assessing parallelism between curves:
The Fisher’s F test is used to determine if one can consider that the models corresponding to the standard sample and the sample of interest are significantly different or not. If the probability corresponding to the F value is lower than the significance level, then one can consider that the difference is significant.
- Goodness of fit coefficients:
This table shows the following statistics:
- The number of observations;
- The number of degrees of freedom (DF);
- The determination coefficient R2;
- The sum of squares of the errors (or residuals) of the model (SSE or SSR respectively);
- The means of the squares of the errors (or residuals) of the model (MSE or MSR);
- The root mean squares of the errors (or residuals) of the model (RMSE or RMSR);
- Model parameters:
This table displays the estimator and the standard error of the estimator for each parameter of the model. It is followed by the equation of the model.
- Predictions and residuals:
This table displays giving for each observation the input data and corresponding prediction and residual. The outliers detected by the Dixon’s test, if any, are displayed at the bottom of the table.
Charts for four-parameter parallel lines logistic regression is four-parameter parallel lines logistic regression in XLSTAT
- Parallel logistics cruves
The Parallel logistics cruves is plotted displaying in blue color, the data and the curve corresponding to the standard sample, and in red color, the data and the curve corresponding to the sample of interest. This makes it easy for you to compare the distribution.
- Plots to assess the model quality:
- Predictions and observed values
- Bar chart of the residuals
Copyright © 2013 Kovach Computing Services, Anglesey, Wales. All Rights Reserved. Portions copyright Addinsoft and Provalis Research.
Last modified 5 April, 2013