Experimental Design and Process Optimization with R

Gerhard Krennrich

1 Introduction

The present document is a short and elementary course on the Design of Experiments (DoE) and empirical process optimization with the open-source Software R . The course is self-contained and does not assume any preknowledge in statistics or mathematics beyond high school level. Statistical concepts will be introduced on an elementary level and made tangible with R-code and R-graphics based on simulated and real world data. So, then, what is DoE and why should the reader become familiar with the concepts of DoE? Very briefly, DoE is the science of varying many experimental parameters in a systematic way to gain insight on how to further improve and optimize these parameters. Chapter 2 will show how and why multidimensional DoE techniques are superiour to the classical “one-dimensional” optimization approach. Chapter 6 will demonstrate why and how DoE can be combined with optimization. Finally, the use of DoE and optimization will be practically demonstrated in chapter 7 for improving the performance of a catalytic system. Historically, Experimental Design started as a branch of statistics in the early years of the 20 th century and has meanwhile grown into a mature method with a plethora of applications in the experimental sciences. Consequently, there are many good and comprehensive books available about DoE, some of which we will make frequent reference to in the present text, namely (George E.P. Box, Norman R. Draper 1987 ) , (D.C. Montgomery 2013 ) and (G.E.P. Box, W.G. Hunter, J.S. Hunter 2005 ) . A more recent text with emphasis on the use of R in conjuction with DoE is (John Lawson 2015 ) . Linear models are comprehensively covered, e.g., by the text book (A. Sen, M. Srivastava 1990 ) . A general, however fairly technical text on linear and nonlinear statistical model building is the excellent book (T. Hastie, R. Tibshirani, J. Friedman 2009 ) . (J.G. Kalbfleisch 1985 ) is a smooth introduction into statistics, probability and statistical inference. The present text draws on these books and on many years of experience as a statistical consultant in the chemical industry. Most examples in this course are therefore taken from applications and optimization projects in the chemical sciences. The primarily intended readers of this document are chemists and engineers entrusted with empirical optimization in research and development. However, the presented methods and concepts are fairly generic and scientist working in other areas such as biology or the medical sciences might benefit from the text. As to software, R, probably together with Phyton, is the only open-source software which combines the whole spectrum of DoE and optimization with the flexibility of a powerful script language that allows any kind of data pre- and postprocessing within one software environment. That makes, in my opinion, R superior to many commercial GUI based tools which often buy userfriendlyness at the expense of flexibility.

1.1 How to install R

The R-software can be downloaded free of charge from the R repository CRAN

An IDE ( I ntegrated D evelopment E nvironment) is reqired for smoothly working with R. An IDE allows editing, running and debugging of R code and managing programm in- and output. In principle any IDE can be used but we recommend R-Studio as the de-facto standard.

Get R-Studio IDE

The R-introduction at CRAN is a concise introduction into the R-language. A short R-introduction

1.2 Some remarks on how to read the present text

This document is not an introduction into the R language, rather the document follows the philosophy of “learning by doing”. In this spirit the above mentioned text R-introduction is recommended as a first reference together with the present R examples on DoE and optimization. As it is usually easier to modify existing code than writing code from scratch, it is hoped that the R-examples in this course will help learning both R and DoE more rapidly. The course is divided into seven chapters. There is, however, one stand-alone chapter, chapter 5, which can be skipped by those readers not explicitly dealing with mixture problems. The final chapter 7 is a published, (Siebert M., Krennrich G., Seibicke M., Siegle A.F., Trapp O. 2019 ) , real-world example combining many elements of DoE and optimization for improving the performance of a catalytic system. This application should encourage readers to use these powerful methods for the sake of their own projects.

A. Sen, M. Srivastava. 1990. Regression Analysis, Theory, Methods and Applications . 1st ed. Springer-Verlag, New York.

D.C. Montgomery. 2013. Design and Analysis of Experiments . 8th ed. John Wiley & Sons Inc.

G.E.P. Box, W.G. Hunter, J.S. Hunter. 2005. Statistics for Experimenters: Design, Innovation, and Discovery . 2nd ed. John Wiley & Sons, Hoboken.

George E.P. Box, Norman R. Draper. 1987. Empirical Model-Building and Response Surfaces . 1st ed. John Wiley & Sons.

J.G. Kalbfleisch. 1985. Probability and Statistical Inference, Vol 1&2 . 2nd ed. Springer.

John Lawson. 2015. Design and Analysis of Experiments with R . 1st ed. Chapman & Hall.

Siebert M., Krennrich G., Seibicke M., Siegle A.F., Trapp O. 2019. “Identifying High-Performance Catalytic Conditions for Carbon Dioxide Reduction to Dimethoxymethane by Multivariate Modelling.” Chemical Science 10:45. https://pubs.rsc.org/en/content/articlelanding/2019/sc/c9sc04591k#!divAbstract .

T. Hastie, R. Tibshirani, J. Friedman. 2009. The Elements of Statistical Learning . 2nd ed. Springer-Verlag.

Breadcrumbs Section. Click here to navigate to respective pages.

Design and Analysis of Experiments with R

Design and Analysis of Experiments with R

DOI link for Design and Analysis of Experiments with R

Get Citation

Design and Analysis of Experiments with R presents a unified treatment of experimental designs and design concepts commonly used in practice. It connects the objectives of research to the type of experimental design required, describes the process of creating the design and collecting the data, shows how to perform the proper analysis of the data,

TABLE OF CONTENTS

Chapter chapter 1 | 15  pages, introduction, chapter chapter 2 | 37  pages, completely randomized designs with one factor, chapter chapter 3 | 55  pages, factorial designs, chapter chapter 4 | 28  pages, randomized block designs, chapter chapter 5 | 52  pages, designs to study variances, chapter chapter 6 | 67  pages, fractional factorial designs, chapter 7 | 45  pages, incomplete and confounded block designs, chapter 8 | 44  pages, split-plot designs, chapter chapter 9 | 32  pages, crossover and repeated measures designs, chapter chapter 10 | 63  pages, response surface designs, chapter chapter 11 | 56  pages, mixture experiments, chapter chapter 12 | 53  pages, robust parameter design experiments, chapter chapter 13 | 14  pages, experimental strategies for increasing knowledge.

  • Privacy Policy
  • Terms & Conditions
  • Cookie Policy
  • Taylor & Francis Online
  • Taylor & Francis Group
  • Students/Researchers
  • Librarians/Institutions

Connect with us

Registered in England & Wales No. 3099067 5 Howick Place | London | SW1P 1WG © 2024 Informa UK Limited

R-bloggers

R news and tutorials contributed by hundreds of R bloggers

Some basic concepts about design of experiments and how to perform their analysis in r.

Posted on January 15, 2022 by R in the Lab in R bloggers | 0 Comments

Basic design of experiments in R for one factor and two factors designs.

You can find all the code, data and results in the GitHub repository for this post: Basic design of experiments .

There is no signal without noise

It never hurts to go back to basics before tackling more complex things. The purpose of this post is to give a brief overview of the basics of design of experiments, their analysis and how to present results using R and packages like ggplot2 and agricolae . Included are one- and two-factor experiments.

What is the design of experiments?

The design of experiments (DOE) deals with the planning and performance of tests with the objective of generating data. Statistical analysis of these data will provide objective evidence that will allow the researcher to resolve questions about a given situation, process or phenomenon.

DOE can be applied to scientific research and in industry. Its goal is to generate knowledge and learning in an efficient manner. Ideally, this process can be considered as a cycle where initial hypotheses are refined as more information is obtained. Thus, we start with an initial hypothesis that we test through experimentation. If the data (our evidence) does not agree with the hypothesis, a new hypothesis is sought to explain the observed discrepancy.

design and analysis of experiments in r

Inductive-deductive loop – Pamela Toman

Hypothesis?

Hypotheses are tentative explanations of the research phenomenon that are stated as propositions or assertions. There are different types of hypotheses, such as:

  • Research hypotheses. Tentative propositions about possible relationships between two or more variables.
  • Null hypothesis. Propositions that deny or refute the relationship between variables.
  • Alternative hypotheses. These are different or “alternate” possibilities to the research and null hypotheses.

design and analysis of experiments in r

Null hypothesis

Data and measurements

A datum is a number that is the product of a measurement. A measurement is a process that links abstract concepts with empirical indicators. Of course, measurements could not be made without a measurement instrument, which can be defined as the resource used by the researcher to record information or data.

design and analysis of experiments in r

Measuring absorbance

Other basic definitions in the design of experiments

It is a change in the operating conditions of a system or process, which is made with the objective of measuring the effect of the change on one or more properties of the product or result.

design and analysis of experiments in r

Schematic view of experiment design

design and analysis of experiments in r

Controlled experiment

Experimental unit

Piece(s) or sample(s) used to generate a value that is representative of the test result.

design and analysis of experiments in r

Response variable (Depedent variable)

Through these variable(s) we know the effect or results of each experimental test. They are the product of our measurements once we make changes in the system.

Controllable factors (Independent variable)

They are process variables that can be modified and fixed at a given level.

design and analysis of experiments in r

Response and independent variables

Non-controllable factors

These are variables that cannot be controlled during the experiment or the normal operation of the process.

Factors studied

These are the variables that are investigated in the experiment to observe how they affect or influence the response variables.

Levels and treatments

The different values assigned to each factor studied in an experimental design are called levels. A combination of levels of all the factors studied is called a treatment or design point. In the case of experimenting with a single factor, each level is a treatment.

Design matrix

It is the arrangement formed by the treatments that will be carried out, including the repetitions. In this same array are usually included the results of each response variable for each treatment and repetition. For me this is a concept of great importance, since from this matrix it is possible to perform directly the different statistical tests in R.

design and analysis of experiments in r

Random error

It is the observed variability that cannot be explained by the factors studied. It is the result of the effect of factors not studied and experimental error.

Experimental error

A component of random error that reflects the experimenter’s errors in the planning and execution of the experiment.

Basic principles of experimental design

  • Randomization. Performing the experimental runs (treatments) in a randomized manner is essential to ensure, as far as possible, the assumption of independence of errors.
  • Repetition. It is running a treatment or combination of factors more than once. This makes it possible to estimate random error and in turn calculate realistic statistics in data analysis.
  • Blocking. This consists of nullifying the effect of factors in which we are not really interested and which may affect the response variables. Not taking into account their effect can lead us to erroneous conclusions about the factors in which we are interested.

Block example: Male and females

Some concepts of statistical inference

  • Population. In a very general way, it can be defined as the totality of possible individuals, specimens, objects or measurements of interest on which a study is made. Populations can be finite or infinite.
  • Parameters. Characteristics that, by means of their numerical value, describe the entire population.
  • Representative sample. An appropriately selected part of a population that retains key aspects of a population.
  • Statistic. A quantity obtained from the data of a sample that summarizes the characteristics of the sample.
  • Statistical inference. These are statistical statements about the population or process based on the information contained in the sample.
  • Probability distribution. It relates the set of values of a characteristic in a given population to the probability associated with each of these values.

design and analysis of experiments in r

Statistical inference

design and analysis of experiments in r

Probability distribution

Some hypothesis testing concepts

  • Statistical hypothesis. An assertion about the values of the parameters of a population or process, which can be tested from the information in a sample.
  • Test statistic. Formula with which a number is calculated from the data, the magnitude of which makes it possible to determine whether or not the null hypothesis is rejected.
  • Acceptance region. These are the possible values of the test statistic where the null hypothesis is not rejected.
  • Rejection region. It is the set of possible values of the test statistic that lead to rejecting the null hypothesis.
  • Type I error. It is when a null hypothesis that is true is rejected. It can also be called false positive.
  • Type II error. This is when a null hypothesis that is false is accepted. It is basically a false negative.
  • Power of the test. It is the probability of rejecting the null hypothesis when it is false.
  • Predefined significance. It is the maximum risk that the experimenter is willing to take with respect to the type I error. It is denoted by α.
  • Observed significance. It is the area under the reference distribution beyond the value of the test statistic. It is called p-value.

design and analysis of experiments in r

Hypothesis testing

Experiments with a single factor

When the objective is to study the effect of a single factor on the mean of a response variable taking into account more than two values of the factor under consideration, it is advisable to perform a completely randomized experiment and analyze the results by analysis of variance (ANOVA). The reason for comparing means by ANOVA and not by Student’s t-tests is that the latter test increases false positives as the number of means to be compared increases, i.e., we will find differences between means where there are none.

The analysis of variance consists of separating the total variation observed in each of the sources that contribute to it.

Design and results

Obtaining a single factor design is very simple using R commands. Let’s say we want to compare the effect of three baking times (35, 45 and 60 min) in the average diameter of cookie batches. For each time we baked ten cookie batches so we have ten replicates. Note that each time and replicate have to be run in aleatory order:

Subsequently, the design can be exported in some other format such as CSV:

The results can be integrated into the design matrix directly, let’s simulate the average diameter with a small function:

Or, once the results have been recorded externally, they can be imported from a CSV file:

The aov function is used for the analysis of factorial designs. Previously, it is recommended to convert the factor values into the factor class and then perform the analysis. Here I going to analyze our first design ( one_fct_dgn ):

Note the aov requires a formula that explicitly specifies the relationship between the response and the factor(s) in the design. Once this has been done, it is possible to obtain an ANOVA table using the “anova” function:

And export it in the format of your choice:

Multiple range comparisons or tests

Once significant differences are established between at least two of the means, the next step is to perform a multiple comparison test, which will allow us to define which means are significantly different. There are several methods that can be used as the LSD method, Tukey’s method, or Duncan’s method. Nice functions for these methods can be found in agricolae package.

The results are consistent among the three methods. The section groups indicates which means are different, different letters mean significant differences. You can export the results of this kind of methods, as a TXT file, with the function capture.output . For example, to export the results of the Tukey method:

Verification of model assumptions

Something that is not very widespread in the scientific literature is the verification of the normality of the residuals, the equality between variances for the means of each treatment and the independence between each data recorded. Verification of these assumptions helps to validate the differences established as significant.

Shapiro-Wilks test to verify normality

This test is easy to perform:

Here, if p-value > 0.05, we can say that data come, approximately, from a normal distribution.

The results of this analysis can be exported with the capture.output function:

Bartlett’s test for equality of variances

To perform this test it is possible to use the bartlett.test function:

A p-value > 0.05 mean that treatment means have equal variances. We can export the results in the same way:

Independence

To verify independence between observations, residuals are plotted as a function of run order. If the residuals follow a defined pattern it will be a clear indication of lack of independence. For this plot I used the ggplot2 package:

design and analysis of experiments in r

Tests such as the Durbin-Watson test are also available to verify this assumption. Use the function durbinWatsonTest in the car package:

A p-value > 0.05 means that data observations are not auto correlated. You can also save the results in the usual way:

Display of results

To visualize the results of this type of experiments, bar graphs are usually used, whose height represents the magnitude of each mean, together with error bars representing the standard deviation and letters showing the significant differences established by the multiple comparisons test:

design and analysis of experiments in r

Alternative

Although bar graphs are widely used in the scientific literature of various types, their use hides the true dispersion of the data for each treatment. One thing I have also noticed is that people tend to judge differences between means by the overlap or separation of standard deviations, which seems to me to be more of a bias than an objective criterion.

As an alternative to the bar chart, it is possible to make a graph showing each observation, the mean, and the letters obtained in the multiple comparisons test:

design and analysis of experiments in r

Two-factor designs

In this type of design the objective is to verify the individual and interaction effect of two factors. Before showing how to obtain and analyze them, it is necessary to recall some concepts:

  • Qualitative factor. Its levels take discrete or nominal values.
  • Quantitative factor. Its levels can take any value within a certain interval. The scale is continuous.
  • Effect of a factor. It is the change observed in the response variable due to a change in the level of the factor.
  • Interaction effect. Two factors interact significantly on the response variable when the effect of one depends on the level of the other.

To obtain this kind of designs we can use the function fac.design from DoE.base package. For this we need specify the number of factors and the number of levels of each factor:

For more details about this function, just type ?fac.design in your console. Also note that there is an extra column with the name “Blocks”, this is not big deal since there is just one and we won’t taking it into account further in the analysis. In future posts I will address block designs.

This design can also be exported in the desired format:

The results can also be integrated directly or imported once recorded in the design matrix. Here, I’m going to import and work with data previously recorded:

To analyze this type of design, the “aov” function is used, specifying each main effect and the interaction effect in the formula:

You can also use the short form Y ~ A*B , which will be displayed in full once you get the ANOVA table:

Note that previously I also converted the values of each factor in the factor class.

Multiple range comparisons tests

Once a significant interaction effect is established, the next step is to establish which means differ. For this, a Tukey test can be used, for which the interaction effect must first be explicitly specified in a separated aov model:

The results of Tukey’s test can be exported directly:

For the verification of the assumptions it is possible to proceed in a manner similar to that of the single-factor designs.

For Bartlett test I previously made an extra column specifying the combination between factors (treatments). This will indicate the groups to bartlett.test function:

design and analysis of experiments in r

The realization of this type of graphs requires keeping in mind the individual means of each treatment and the interaction effects, if there is some. To calculate them together with the standard deviations we used the “group_by” function and then “summarise” from the “dplyr” package:

design and analysis of experiments in r

To make the graph showing each observation, the mean, a confidence interval and the letters obtained in the multiple comparisons test, it is necessary to use the data.frame with the observations next to the one with the means:

design and analysis of experiments in r

A good way to visualize the interaction effect is through interaction.plot function:

design and analysis of experiments in r

That’s it! Thank you very much for visiting this site, I hope you find the content of this post useful. See you soon!

Juan Pablo Carreón Hidalgo 🤓

[email protected] https://twitter.com/JuanPa2601

The text and code on this tutorial is under Creative Commons Attribution 4.0 International License .

CC BY 4.0

Copyright © 2022 | MH Corporate basic by MH Themes

Navigation Menu

Search code, repositories, users, issues, pull requests..., provide feedback.

We read every piece of feedback, and take your input very seriously.

Saved searches

Use saved searches to filter your results more quickly.

To see all available qualifiers, see our documentation .

  • Notifications You must be signed in to change notification settings

An R package for Design and Analysis of Experiment

ehassler/MontgomeryDAE

Folders and files.

NameName
6 Commits

Repository files navigation

R package for montgomery's design and analysis of experiments 10th ed..

This repository is for Wiley's companion R package to Douglas C. Montgomery's "Design and Analysis of Experiments" 10th edition. The most current version of the PDF Guide will walk you through how to use R to work examples and exercises in the text.

The PDF guide also shows how to install the most current package . The package contains the content of the PDF in vignettes and also all of the data sets used in the book. It makes doing the exercises much easier.

One of the great features of github is that you can report bugs through the "Issues" tab. Please look there to see if any error you find has been reported, and if it has not then post a report of the error there.

Design and Analysis of Experiments and Observational Studies using R

A volume in the chapman & hall/crc texts in statistical science series.

design and analysis of experiments in r

Sample Course Documents

Sample course documents (i.e., slides, assignments, etc.) are available here .

scidesignR is an R package that contains the data used in the book.

Creative Commons License

design and analysis of experiments in r

1st Edition

Design and Analysis of Experiments with R

VitalSource Logo

  • Taylor & Francis eBooks (Institutional Purchase) Opens in new tab or window

Description

Design and Analysis of Experiments with R presents a unified treatment of experimental designs and design concepts commonly used in practice. It connects the objectives of research to the type of experimental design required, describes the process of creating the design and collecting the data, shows how to perform the proper analysis of the data, and illustrates the interpretation of results. Drawing on his many years of working in the pharmaceutical, agricultural, industrial chemicals, and machinery industries, the author teaches students how to: Make an appropriate design choice based on the objectives of a research project Create a design and perform an experiment Interpret the results of computer data analysis The book emphasizes the connection among the experimental units, the way treatments are randomized to experimental units, and the proper error term for data analysis. R code is used to create and analyze all the example experiments. The code examples from the text are available for download on the author’s website, enabling students to duplicate all the designs and data analysis. Intended for a one-semester or two-quarter course on experimental design, this text covers classical ideas in experimental design as well as the latest research topics. It gives students practical guidance on using R to analyze experimental data.

Table of Contents

John Lawson is a professor in the Department of Statistics at Brigham Young University.

Critics' Reviews

"This is an excellent but demanding text. … This book should be mandatory reading for anyone teaching a course in the statistical design of experiments. … reading this text is likely to influence their course for the better." — MAA Reviews , March 2015 "In my opinion, this is a very valuable book. It covers the topics that I judge should be in such a book including what might be called the standard designs and more … it has become my go to text on experimental design." David E. Booth, Technometrics

About VitalSource eBooks

VitalSource is a leading provider of eBooks.

  • Access your materials anywhere, at anytime.
  • Customer preferences like text size, font type, page color and more.
  • Take annotations in line as you read.

Multiple eBook Copies

This eBook is already in your shopping cart. If you would like to replace it with a different purchasing option please remove the current eBook option from your cart.

Book Preview

design and analysis of experiments in r

The country you have selected will result in the following:

  • Product pricing will be adjusted to match the corresponding currency.
  • The title Perception will be removed from your cart because it is not available in this region.

Ask the publishers to restore access to 500,000+ books.

Internet Archive Audio

design and analysis of experiments in r

  • This Just In
  • Grateful Dead
  • Old Time Radio
  • 78 RPMs and Cylinder Recordings
  • Audio Books & Poetry
  • Computers, Technology and Science
  • Music, Arts & Culture
  • News & Public Affairs
  • Spirituality & Religion
  • Radio News Archive

design and analysis of experiments in r

  • Flickr Commons
  • Occupy Wall Street Flickr
  • NASA Images
  • Solar System Collection
  • Ames Research Center

design and analysis of experiments in r

  • All Software
  • Old School Emulation
  • MS-DOS Games
  • Historical Software
  • Classic PC Games
  • Software Library
  • Kodi Archive and Support File
  • Vintage Software
  • CD-ROM Software
  • CD-ROM Software Library
  • Software Sites
  • Tucows Software Library
  • Shareware CD-ROMs
  • Software Capsules Compilation
  • CD-ROM Images
  • ZX Spectrum
  • DOOM Level CD

design and analysis of experiments in r

  • Smithsonian Libraries
  • FEDLINK (US)
  • Lincoln Collection
  • American Libraries
  • Canadian Libraries
  • Universal Library
  • Project Gutenberg
  • Children's Library
  • Biodiversity Heritage Library
  • Books by Language
  • Additional Collections

design and analysis of experiments in r

  • Prelinger Archives
  • Democracy Now!
  • Occupy Wall Street
  • TV NSA Clip Library
  • Animation & Cartoons
  • Arts & Music
  • Computers & Technology
  • Cultural & Academic Films
  • Ephemeral Films
  • Sports Videos
  • Videogame Videos
  • Youth Media

Search the history of over 866 billion web pages on the Internet.

Mobile Apps

  • Wayback Machine (iOS)
  • Wayback Machine (Android)

Browser Extensions

Archive-it subscription.

  • Explore the Collections
  • Build Collections

Save Page Now

Capture a web page as it appears now for use as a trusted citation in the future.

Please enter a valid web address

  • Donate Donate icon An illustration of a heart shape

Design and analysis of experiments with R

Bookreader item preview, share or embed this item, flag this item for.

  • Graphic Violence
  • Explicit Sexual Content
  • Hate Speech
  • Misinformation/Disinformation
  • Marketing/Phishing/Advertising
  • Misleading/Inaccurate/Missing Metadata

[WorldCat (this item)]

plus-circle Add Review comment Reviews

2 Favorites

Better World Books

DOWNLOAD OPTIONS

No suitable files to display here.

IN COLLECTIONS

Uploaded by station42.cebu on June 27, 2023

SIMILAR ITEMS (based on metadata)

Design and Analysis of Experiments with R

[Up]

Documentation for package ‘daewr’ version 1.2-11

  • DESCRIPTION file .

A B C D E F G H I L M N O P Q R S T V W

Data frames and functions for Design and Analysis of Experiments with R
Alternate 16 run screening designs
Two-period crossover study of antifungal agent
apolipoprotein survey varaince component study
Confounded apple slice browning experiment
2^{(7-3)} arsenic removal experiment
2^{(7-3)} arsenic removal experiment augmented with mirror image
Confounded Block Dishwashing Experiment
Confounded block fractional mouse growth experiment
mouse liver enzyme experiment
Balanced incomplete blocksize
Extra-period crossover bioequivalence study
Latin Square bioequivalence experiment
Variance component study of calcium in blood serum
Box and Meyer's unreplicated 2^4 from Chapter 3
blood pressure monitor experiment
Bread rise experiment data from Chapter 2
Split-Plot response surface for cake baking experiment
CCD design for cement workability experiment
Chemical process experiment data from Chapter 3
Williams' crossover design for sprinting experiment
CO emmisions experiment data from Chapter 3
This function makes a colormap of correlations in a design matrix
Table 12.21 Experiment with Elastometric Connector
Control factor array and summary statistics for controller circuit design experiment
Split-plot response surface for ceramic pipe experiment
paecilomyces variotii culture experiment
Data frames and functions for Design and Analysis of Experiments with R
Repeated measures study with dairy cow diets
Definitive Screening Designs
Data from rat behavior experiment in Chapter 4
D-efficient Estimation Equivalent Response Surface Designs
D-efficient Estimation Equivalent Response Surface Designs
D-efficient Estimation Equivalent Response Surface Designs
D-efficient Estimation Equivalent Response Surface Designs
D-efficient Estimation Equivalent Response Surface Designs
D-efficient Estimation Equivalent Response Surface Designs
D-efficient Estimation Equivalent Response Surface Designs
Single array and raw response for silicon layer growth experiment
Control array and variance of response for silicon layer growth experiment
Control array and mean response for silicon layer growth experiment
F-Distribution critical values
Subsequent steps in a forward stepwise regression that preserves model hierarchy
Forward Stepwise modeling taking into account special structure of Definitive Screening Design
Find first term to enter forward stepwise regression that preserves model hierarchy
An Effective Design Based Model Fitting Method for Definitive Screening Designs
Find first term to enter forward stepwise regression that preserves model hierarchy
Find first term to enter forward stepwise regression that preserves model hierarchy
F-Distribution Power Calculation
F-Distribution Power Calculation
F-Distribution Power Calculation
This function makes a full normal plot of the elements of the vector called effects
Gauge R&R Study
This function computes the gap statistic which is used to test for an outlier using Daniels method
This function uses Daniel's Method to find an outlier in an unreplicated 2^{(k-p)} design.
Unreplicated split-plot fractional-factorial experiment on geometric distortion of drive gears
This function makes a half normal plot of the elements of the vector called effects
low grade hardwood conjoint study
RSM forward regression keeping model hierarchy
First step in a forward stepwise regression that preserves model hierarchy
Single array for injection molding experiment
Interleave vectors
Lenth's Plot of Effects
This function uses the LGB Method to detect significant effects in unreplicated fractional factorials.
This function does the calculations for the LGB Method to detect significant effects in unreplicated fractional factorials.
Mod function
Model Robust Factorial Designs
mixture process variable experiment with mayonnaise
Yields of naphthalene black
Optimum Plackett-Burman Designs
Blocked response surface design for pastry dough experiment
Plackett-Burman Designs
Pesticide formulation experiment
pesticide application experiment
Unreplicated split-plot 2^5 experiment on plasma treatment of paper
Polvoron mixture experiment
polymerization strength variability study
Complete control factor array and noise factor array for connector experiment
Library of substituted hydroxyphenylurea compounds
Cattle rations design experiment data from Table 10.16
generalized RCB golf driving experiment
Herbicide degradation experiment
Rubber Elasticity data
Split-plot experiment on sausage casing with RCB in whole plot
Single array for starting motor experiment
dry mix soup experiment
dry soup mix variance component study
Split-plot cookie baking experiment
Split-plot mixture process variable experiment with vinyl
Standard Order
Repeated measures study with dairy cow diets
Strung out control factor array and raw response data for Ina tile experiment
Sugarbeet data from Chapter 2
taste test panel experiment
Teaching experiment data from Chapter 2
Tetracycline concentration in plasma
Control factor array and summary statistics for Ina tile experiment
Box-Behnken design for trebuchet experiment
This function performs Tukey's single degree of freedom test for interaction in an unreplicated two-factor design
confidence limits for method of moments estimators of variance components
Vinysl plasticizer formulations experiment data
Assay of Viral Contamination experiment data from Chapter 3
Volt meter experiment data from Chapter 3
Web page design experiment data from Chapter 3
Table 12.24 Experiment with Weld Tensile Strength

Note that this "residual" for the within plot \(subplot\) part of the analysis is actually the sum of squares for the interaction of rows \(w\ hole plots\) with varieties \(subplot treatments\)---as in an RCBD.

- r_k\(i\) ~ N\(0, sigma^2_r\)

- e_ijk ~ N\(0, sigma^2_e\)

design and analysis of experiments in r

  • Science & Math
  • Mathematics

Sorry, there was a problem.

Kindle app logo image

Download the free Kindle app and start reading Kindle books instantly on your smartphone, tablet, or computer - no Kindle device required .

Read instantly on your browser with Kindle for Web.

Using your mobile phone camera - scan the code below and download the Kindle app.

QR code to download the Kindle App

Image Unavailable

Design and Analysis of Experiments with R (Chapman & Hall/CRC Texts in Statistical Science)

  • To view this video download Flash Player

design and analysis of experiments in r

Follow the author

John Lawson

Design and Analysis of Experiments with R (Chapman & Hall/CRC Texts in Statistical Science) 1st Edition

Design and Analysis of Experiments with R presents a unified treatment of experimental designs and design concepts commonly used in practice. It connects the objectives of research to the type of experimental design required, describes the process of creating the design and collecting the data, shows how to perform the proper analysis of the data, and illustrates the interpretation of results.

Drawing on his many years of working in the pharmaceutical, agricultural, industrial chemicals, and machinery industries, the author teaches students how to:

  • Make an appropriate design choice based on the objectives of a research project
  • Create a design and perform an experiment
  • Interpret the results of computer data analysis

The book emphasizes the connection among the experimental units, the way treatments are randomized to experimental units, and the proper error term for data analysis. R code is used to create and analyze all the example experiments. The code examples from the text are available for download on the author’s website, enabling students to duplicate all the designs and data analysis.

Intended for a one-semester or two-quarter course on experimental design, this text covers classical ideas in experimental design as well as the latest research topics. It gives students practical guidance on using R to analyze experimental data.

  • ISBN-10 9781439868133
  • ISBN-13 978-1439868133
  • Edition 1st
  • Publisher Chapman and Hall/CRC
  • Publication date December 17, 2014
  • Part of series Chapman & Hall/CRC Texts in Statistical Science
  • Language English
  • Dimensions 6.75 x 1.25 x 9.75 inches
  • Print length 628 pages
  • See all details

Editorial Reviews

"This is an excellent but demanding text. … This book should be mandatory reading for anyone teaching a course in the statistical design of experiments. … reading this text is likely to influence their course for the better." ― MAA Reviews , March 2015

"Thank you for writing your phenomenal book "Design and Analysis of Experiments with R". I'm teaching a new course this spring on experimental design and reinforcement learning. The students are graduate bioengineers, so I was having difficulty finding a text that blends theory, practice, and computation. Your book excels at all three. The first chapter I read clarified several topics and improved both my teaching and research. After testing a dozen DOE and RSM books, yours is the clear winner. I understand the enormous time that goes into a well-constructed textbook. I hope this message conveys my deep appreciation for your effort." ― Paul Jensen , Ph.D., Assistant Professor , Department of Bioengineering and Carl R. Woese Institute for Genomic Biology, University of Illinois at Urbana-Champaign

"In my opinion, this is a very valuable book. It covers the topics that I judge should be in such a book including what might be called the standard designs and more … it has become my go to text on experimental design."

David E. Booth, Technometrics

About the Author

John Lawson is a professor in the Department of Statistics at Brigham Young University.

Product details

  • ASIN ‏ : ‎ 1439868131
  • Publisher ‏ : ‎ Chapman and Hall/CRC; 1st edition (December 17, 2014)
  • Language ‏ : ‎ English
  • Hardcover ‏ : ‎ 628 pages
  • ISBN-10 ‏ : ‎ 9781439868133
  • ISBN-13 ‏ : ‎ 978-1439868133
  • Item Weight ‏ : ‎ 2.25 pounds
  • Dimensions ‏ : ‎ 6.75 x 1.25 x 9.75 inches
  • #748 in Statistics (Books)
  • #2,159 in Probability & Statistics (Books)

About the author

John lawson.

John Lawson is a Professor Emeritus from the Statistics Department at Brigham Young University where he taught from 1986 to 2019. He is an ASQ-CQE and he has a Masters Degree in Statistics from Rutgers University and a PhD in Applied Statistics from the Polytechnic Institute of N.Y. He worked as a statistician for Johnson & Johnson Corporation from 1971 to 1976, and he worked at FMC Corporation Chemical Division from 1976 to 1986 where he was the Manager of Statistical Services. In industry he used designed experiments and statistical analysis to help engineers and chemists on product development and manufacturing process improvements. At BYU he taught courses on experimental design and quality control and consults with faculty and graduate students involved in research projects through the BYU Center for Statistical Consultation and Collaborative Research. He is the the co-author (with John Erjavec) of Basic Experimental Strategies and Data Analysis for Science and Engineering, CRC Press, the author of Design and Analysis of Experiments with R, CRC Press, and the author of An Introduction to Acceptance Sampling and SPC with R, CRC Press. Additional resources for these books, such as electronic versions of computer code in the books, lecture slides, etc. can be downloaded from: https://lawsonjsl7.netlify.app/webbook/. Additional resources for an earlier book Design and Analysis of Experiments with SAS, CRC Press can be downloaded from: https://sasbook.netlify.com

Customer reviews

  • 5 star 4 star 3 star 2 star 1 star 5 star 77% 10% 4% 4% 6% 77%
  • 5 star 4 star 3 star 2 star 1 star 4 star 77% 10% 4% 4% 6% 10%
  • 5 star 4 star 3 star 2 star 1 star 3 star 77% 10% 4% 4% 6% 4%
  • 5 star 4 star 3 star 2 star 1 star 2 star 77% 10% 4% 4% 6% 4%
  • 5 star 4 star 3 star 2 star 1 star 1 star 77% 10% 4% 4% 6% 6%

Customer Reviews, including Product Star Ratings help customers to learn more about the product and decide whether it is the right product for them.

To calculate the overall star rating and percentage breakdown by star, we don’t use a simple average. Instead, our system considers things like how recent a review is and if the reviewer bought the item on Amazon. It also analyzed reviews to verify trustworthiness.

  • Sort reviews by Top reviews Most recent Top reviews

Top reviews from the United States

There was a problem filtering reviews right now. please try again later..

design and analysis of experiments in r

Top reviews from other countries

design and analysis of experiments in r

  • About Amazon
  • Investor Relations
  • Amazon Devices
  • Amazon Science
  • Sell products on Amazon
  • Sell on Amazon Business
  • Sell apps on Amazon
  • Become an Affiliate
  • Advertise Your Products
  • Self-Publish with Us
  • Host an Amazon Hub
  • › See More Make Money with Us
  • Amazon Business Card
  • Shop with Points
  • Reload Your Balance
  • Amazon Currency Converter
  • Amazon and COVID-19
  • Your Account
  • Your Orders
  • Shipping Rates & Policies
  • Returns & Replacements
  • Manage Your Content and Devices
 
 
 
 
  • Conditions of Use
  • Privacy Notice
  • Consumer Health Data Privacy Disclosure
  • Your Ads Privacy Choices

design and analysis of experiments in r

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • View all journals
  • Explore content
  • About the journal
  • Publish with us
  • Sign up for alerts
  • Open access
  • Published: 17 September 2024

Improving rigor and reproducibility in western blot experiments with the blotRig analysis

  • Cleopa Omondi 1 ,
  • Austin Chou 1 ,
  • Kenneth A. Fond 1 ,
  • Kazuhito Morioka 1 ,
  • Nadine R. Joseph 1 ,
  • Jeffrey A. Sacramento 1 ,
  • Emma Iorio 1 ,
  • Abel Torres-Espin 1 , 3 , 4 ,
  • Hannah L. Radabaugh 1 ,
  • Jacob A. Davis 1 ,
  • Jason H. Gumbel 1 ,
  • J. Russell Huie 1 , 2 &
  • Adam R. Ferguson 1 , 2  

Scientific Reports volume  14 , Article number:  21644 ( 2024 ) Cite this article

Metrics details

  • Biochemistry
  • Biological techniques
  • Computational biology and bioinformatics
  • Neuroscience

Western blot is a popular biomolecular analysis method for measuring the relative quantities of independent proteins in complex biological samples. However, variability in quantitative western blot data analysis poses a challenge in designing reproducible experiments. The lack of rigorous quantitative approaches in current western blot statistical methodology may result in irreproducible inferences. Here we describe best practices for the design and analysis of western blot experiments, with examples and demonstrations of how different analytical approaches can lead to widely varying outcomes. To facilitate best practices, we have developed the blotRig tool for designing and analyzing western blot experiments to improve their rigor and reproducibility. The blotRig application includes functions for counterbalancing experimental design by lane position, batch management across gels, and analytics with covariates and random effects.

Similar content being viewed by others

design and analysis of experiments in r

A critical path to producing high quality, reproducible data from quantitative western blot experiments

design and analysis of experiments in r

A data analysis framework for combining multiple batches increases the power of isobaric proteomics experiments

design and analysis of experiments in r

OmicScope unravels systems-level insights from quantitative proteomics data

Introduction.

Proteomic technologies such as protein measurement with folin phenol reagent were first introduced by Lowry et al. in 1951 1 . The resulting qualitative data are typically confirmed by a second, independent method such as western blot (WB) 2 , 3 . The WB method, first described by Towbin et al. 4 and Burnette 5 in 1979 and 1981, respectively, uses specific antibody-antigen interactions to confirm the protein present in the sample mixture. Quantitative WB (qWB assay) is a technique to measure protein concentrations in biological samples with four main steps: (1) protein separation by size, (2) protein transfer to a solid support, (3) marking a target protein using proper primary and secondary antibodies for visualization, and (4) semi -quantitative analysis 6 . Importantly, qWB data is considered semi -quantitative because methods to control for experimental variability ultimately yield relative comparisons of protein levels rather than absolute protein concentrations 2 , 3 , 7 , 8 . Similarly, western blotting applying ECL (enhanced chemiluminescence) is considered a semi-quantitative method because it lacks cumulative luminescence linearity and offers limited quantitative reproducibility 9 . However, the emergence of highly sensitive fluorescent labeling techniques, which exhibit a wider quantifiable linear range, greater sensitivity, and improved stability when compared to the conventional ECL detection method, now permits the legitimate characterization of protein expression as linearly quantitative 10 . Current methodologies do not sufficiently account for diverse sources of variability, producing highly variable results between different laboratories and even within the same lab 11 , 12 , 13 . Indeed, qWB data exhibits more variability compared to other experimental techniques such as enzyme linked immunosorbent assay (ELISA) 14 . For example, results have shown that qWB can produce significant variability in detecting host cell proteins and lead to researchers missing or overestimating true biological effects 15 . This in turn results in publication of irreproducible qWB interpretations, which leads to loss of its credibility 13 . In the serious cases, qWB results may even provide clinical misdiagnosis 16 that could impact on a larger public health concern due to the prevalence of WB in biomedical research, such as diagnosis of SARS-CoV2 infection 17 .

The process of recognizing and accounting for variability in WB analyses will ultimately improve reproducibility between experiments. A growing body of studies has shown that this requires a fundamental shift in the experimental methodology across data acquisition, analysis, and interpretation to achieve precise and accurate results 2 , 3 , 11 , 12 , 13 .

Here we highlight experimental design practices that enable a statistics-driven approach to improve the reproducibility of qWBs. Specifically, we discuss major sources of variability in qWB including the non-linearity in antibody signal 2 , 3 ; imbalanced experimental design 13 ; lack of standardization in the treatment of technical replicates 3 , 18 ; and variability between protein loading, lanes, and blots 2 , 7 , 19 . To address these issues, we provide new comprehensive suggestions for quantitative evaluation of protein expression by combining linear range characterization for antibodies, appropriate counterbalancing during gel loading, running technical replicates across multiple gels, and by taking careful consideration of the analysis method. By applying these experimental practices, we can then account for more sources of variability by running analysis of covariance (ANCOVA) or generalized linear mixed models (LMM). Such approaches have been shown to successfully improve reproducibility compared to other methods 13 .

Good options for qWB protein bands analysis using free, downloadable tools are available for researchers. Amongst others, LI-COR Image Studio Lite can be used to measure the intensity of protein bands in western blots and calculate their relative abundance . Likewise, ThermoFisher ImageQuant Lite offers features such as the ability to perform background subtraction and normalization. However, to date, no specific tools are freely available to provide a map to counterbalance samples, which overcome imperfect uniform protein electrophoresis/transfer and perform statistical analysis. Here, we present blotRig, a tool for researchers with functionalities to counterbalance samples and perform statistical analysis.

To help improve WB rigor we developed the blotRig protocol and application harnessing a database of 6000 + western blots from N = 281 subjects (rats and mice) collected by multiple UCSF labs on core equipment. To demonstrate blotRig best practices in a real-world experiment, we carried out prospective multiplexed WB analysis of protein lysate from lumbar cord in rodent models of spinal cord injury (SCI) (N = 29 rats) in 2 groups (experimental group & control group). In order to show that these experimental suggestions could improve qWB reproducibility, we compared different statistical approaches to handling loading controls and technical replicates. Specifically, we applied two strategies to integrate loading controls: (i) normalizing the target protein levels by dividing by the loading control or (ii) treating the loading control as a covariate in a LMM. Additionally, we analyzed technical replicates in four ways: (1) assume each sample was only run once without replication, (2) treat each technical replicate as an independent sample, (3) use the mean of the three technical replicate values, and 4) treat the replicate as a random effect in a LMM. Altogether, we found that the statistical power of the experiment was significantly increased when we used loading control as a covariate with technical replicates as a random effect during analysis. In addition, the effect size was increased, and the p-value of our analysis decreased when using this LMM, suggesting the potential for greater sensitivity in our WB experiment when using this approach 20 . Through rigorous experimental design and statistical analysis we show that we can account for greater variability in the data and more clearly identify underlying biological effects.

Materials and methods

All experiments protocol were approved by the University Laboratory Animal Care Committee at University of California, San Francisco (UCSF, CA, USA) and followed the animal guidelines of the National Institutes of Health Guide for the Care and Use of Laboratory animals (National Research Council (US) Committee for the Update of the Guide for the Care and Use of Laboratory Animals, 2011). We followed The ARRIVE guidelines (Animal Research: Reporting In Vivo Experiments) to describe our in vivo experiments.

Male Simonsen Long Evans rats (188–385 g; Gilroy (Santa Clara, CA, USA), (N = 29) aged 3 weeks were housed under standard conditions with a 12-h light–dark cycle (6:30 am to 6:30 pm) and were given food and water ad libitum. The animals were housed mostly in pairs in 30 × 30 × 19-cm isolator cages with solid floors covered with a 3 cm layer of wood chip bedding. The experimenters were blind to the identity of treatments and experimental conditions, and all experiments were designed to minimize suffering and limit the number of animals required.

Anesthesia and surgery

We performed non-survival spinal cord injury and spared nerve injury surgeries on animals. Specifically, 3 week old female rats were anesthetized with continuous inhalation of isoflurane (1–5% mg/kg) while on oxygen (0.6–1 mg/kg) in accordance with the IACUC surgical and anesthesia guidelines. Preoperative 0.5% lidocaine local infiltration was applied once at surgical site, avoiding injection into muscle. Fur over the T7–T9 thoracic level was shaved. The dorsal skin was aseptically prepared with surgical iodine or chlorhexidine and 70% ethanol. A small longitudinal incision was made along the spine through the skin, fascia, and muscle to expose the T7-T9 vertebrae. Animals undergoing sham procedure did not undergo laminectomy and immediately proceeded to wound closure. Overlying muscle and subcutaneous tissue was sutured closed using an absorbable suture in a layered fashion. External skin was reinforced using monofilament suture or tissue glue as needed. Animals were euthanized after 30 min to extract spinal cord tissue through fluid expulsion.

Experimental methodology

In accordance with established quality standards for preclinical neurological research 21 , experimenters were kept blind to experimental group conditions throughout the entire study. Western blot loading order was determined a priori by a third-party coder, who ensured that a representative sample from each condition was included on each gel in a randomized block design. The number of subjects per condition was kept consistent across groups for each experiment to ensure that proper counterbalancing could be achieved across independent western runs. All representative western images presented in the figures represent lanes from the same gel. Sometimes, the analytical comparisons of interest were not available on adjacent lanes even though they come from the same gel because of our randomized counterbalancing procedure.

  • Western blot

The example western blot data used in this paper are taken from a model of spared nerve injury in animals with spinal cord injury. The nerve injury model used is based on models from pain literature 22 , where two of the three branches of the sciatic nerve are transected, sparing the sural nerve (SNI) 23 . Two surgeons perform the procedure simultaneously, with injuries occurring 5 min apart. The spinal cord of animals was obtained based on fluid expulsion model 24 and a 1 cm section of the lumbar region was excised at the lumbar enlargement section. The tissue was then preserved in a -80 degree freezer until it was needed for an experiment, at which point it was thawed and used to run a Western blot. We conducted a Western blot analysis on 29 samples from animals using standard biochemical methods. We measured the protein levels of the AMPA receptor subunit GluA2 and used beta-actin as a loading control. The data from these experiments was then aggregated and used for statistical analysis.

Protein assay

We assayed sample protein concentration using a bicinchoninic acid (BCA assay (Pierce) for reliable quantification of total protein using a plate reader (Tecan; GeNios) with triplicate samples (technical replicates) detected against a Bradford Assay (BSA) standard curve. Technical replicates are multiple measurements that are performed under the same conditions in order to quantify and correct for technical variability and improve the accuracy and precision of the results (48). We ran the same WB loading scheme three times (technical replicates of the entire gel) and measured the protein levels of AMPA receptors.

Polyacrylamide gel electrophoresis and multiplexed near-infrared immunoblotting

The approach involved performing serial 1:2 dilutions with cold Laemmli sample buffer in room temperature; 15 μg of total protein per sample was loaded into separate lanes on a precast 10–20% electrophoresis gel (Tris–HCl polyacrylamide, BioRad) to establish linear range (Fig.  1 ). The blotRig software helps counterbalance sample positions across the gel by treatment condition. (Fig.  2 ). A kaleidoscope ladder was loaded on the first lane of each gel to confirm molecular weight (Fig.  2 ). The gel was electrophoresed for 30 min at 200 V in SDS buffer (25 mm Tris, 192 mm glycine, 0.1% SDS, pH 8.3; BioRad). Protein was transferred to a nitrocellulose membrane in cold transfer buffer (25 mm Tris, 192 mm glycine, 20% ethanol, pH 8.3). Membrane transfer was confirmed using Ponceau S stain (67) followed by a quick rinse and blocking in Odyssey blocking buffer (Li-Cor) containing Tween-20.

figure 1

Determining linear range of antibodies to optimize parametric analysis of Western blot data. When small or large protein concentrations are loaded, there is often a possibility that their representation on western blot band density may become non-linear. If there is a disconnect between the observed and expected protein concentrations, results may be inaccurate. Thus determining the linear range wherein, a one-unit increase in protein is reflected in a linear increase in band density for each western blot antibody is a crucial initial step to ensure confidence in reproducibility of the linear models commonly applied to western blot data analysis.

figure 2

Counterbalancing to reduce bias. ( A) Experimental design. A simple hypothetical experimental design for illustrating counterbalancing. Two experimental groups (Wild Type vs Transgenic), with two treatments (Drug vs Vehicle) analyzed within each individual. This 2 (Experimental Condition) by 2 (Tissue Area) design yields four groups. ( B) Counter-balanced Gel Loading. The goal of appropriate counterbalancing is to optimize the sequence in which samples are loaded such that groups are represented equally across the gel. Those with red X have with the experimental groups and treatment condition grouped in the same area of the gel, and thus variability across the gel may be conflated with group differences. In contrast, those with the green check are organized so that experimental condition and treatment condition are better placed to reduce the possibility of any single group being over-represented in a particular area of the gel.

The membrane was blocked for 1 h in Odyssey Blocking Buffer (Li-Cor) containing 0.1% Tween-20, followed by an overnight incubation in primary antibody solution at 4 °C. Membrane incubation was done in a primary antibody solution containing Odyssey blocking buffer, Tween-20, appropriate primary antibody receptor targeting1:2000 mouse PSD-95 (cat # MA1-046,Thermofisher), 1:200 rabbit GluA1 (cat # AB1504, Millipore), 1:200 rabbit GluA2 (cat # AB1766, Millipore), 1:200 rabbit pS831(cat # 04–823, Millipore), 1:200 p880 (cat#07–294, Millipore) or 1:1,500 mouse actin loading control (cat # 612,857, BD Transduction)]. Following incubation, the membrane was washed 4 × 5 min with Tris-buffered saline containing 0.1% Tween 20 (TTBS) and incubated in fluorescent-labeled secondary antibody (1:30 K LiCor IRdye appropriate goat anti-rabbit in Odyssey blocking buffer plus 0.2% Tween 20) for 1 h in the dark. Subsequent to 4 × 5 min washes in TTBS, followed by a 5 min wash in TBS.

Membrane incubation was used to detect the presence of a specific protein or antigen on a membrane. In this case, the membrane was incubated with a fluorescently labeled secondary antibody solution that was specifically tuned to the emission spectra of the laser lines used by the Li-Cor Odyssey quantitative near-infrared molecular imaging system instrument. This allows for specific detection of the protein of interest on the membrane. The sample is then imaged using an infrared imaging system that is optimized for detecting the specific wavelengths of light emitted by the fluorescent label. Additional rounds of incubation and imaging are performed to detect additional proteins using the multiplexing functionality of the Li-Cor instrument, with each round adding new bands at different molecular weight ranges. This allows for the detection of multiple proteins in the same sample, maximizing the proteomic detection.

Quantitative near-IR densitometric analysis

Using techniques optimized in the our lab 25 , 26 , we established near-infrared labeling and detection techniques (Odyssey Infrared Imaging System, Li-Cor) to quantify linear intensity detection of fluorescently labeled protein bands. The biochemistry is performed in a blinded, counterbalanced fashion, and three independent replications of the assay are run on different days 27 . Fluorescent Western blotting utilizes fluorescent-labeled secondary antibodies to detect the target protein, which allows for more sensitive and specific detection compared to chemiluminescence 11 , 28 , 29 . Additionally, fluorescence imaging allows multiple detection of a target protein and internal loading control in the same blot, which enables more accurate correction of sample-to-sample and lane-to-lane variation 11 , 30 , 31 . This provides a more accurate and reliable quantification of the target protein, making it a popular choice for quantitative analysis of WB data.

It is good practice for the pipetting experimenter to remain blind to experimental conditions during gel loading, transfer, and densitometric quantification. We achieved this using de-identified tube codes and a priori gel loading sequences that were developed by an outside experimenter using the method implemented in the blotRig software.

Statistical analyses

Statistical analyses were performed using the R statistical software. Our WB data was analyzed using parametric statistics. The WB was run using three independent replications and covariance corrected by beta-actin loading control, with replication statistically controlled as a random factor. Significance was assessed at p < 0.05 25 , 26 , 32 , 33 , 34 . We report estimated statistical power and standardized regression coefficient effect sizes in the results section.

All ANOVAs were run using the stats R package; standardized effect size was calculated using the parameters R package 35 . Linear mixed models were run using the lme4 R package. Observed power was calculated by Monte Carlo simulation (1000x) run on the fitted model (either ANOVA or LMM) using the simR package 36 . For the development of the blotRig interface, the R packages used included: shiny, tidyverse, DT, shinythemes, shinyjs, and sortable ) 37 , 38 , 39 , 40 , 41 , 42 . You can access the blotRig analysis software, which includes code for inputting experimental parameters for all Western blot analysis, through the following link: https://atpspin.shinyapps.io/BlotRig/.

Designing reproducible western blot experiments

Determining linear range for each primary antibody.

Most WB analyses assume semi -quantitatively that the relationship between qWB assay optical density data (i.e. western band signal) and protein abundance is linear 2 , 3 , 11 , 18 . Accordingly, most qWB analyses use statistical tests (t-test; ANOVA) that assume a linear effect. However, recent studies have shown that the relationship can potentially be highly non-linear 19 As Fig.  1 illustrates, the WB band signal can become non-linearly correlated with protein concentrations at low and high values. This may result in inaccurate quantification of relative target protein amount in the experiment and violates the assumptions for linear model which can lead to false inferences. To address the assumption of linearity, it is important to first determine the optimal linear range for each protein of interest so that one can be confident that a unit change in band density reflects a linear change in protein concentration. This enables an experimenter to accurately quantify the protein of interest and apply linear statistical methods appropriately for hypothesis testing.

Counterbalancing during experimental design

Counterbalancing is the practice of having each experimental condition represented on each gel and evenly distributing them to prevent overrepresentation of the same experimental groups in consecutive lanes. For example, imagine an experimental design in which we are studying two experimental groups (wild type and transgenic animals) and are also looking at two treatment conditions (Drug and Vehicle). The best way to determine the effects and interactions between our experimental and treatment groups would be to create a balanced factorial design. A factorial design is one in which all combinations of levels across factors are represented. For the current example, a balanced factorial design would produce four groups, covering each possible combination (Drug-treated Wild Type, Vehicle-treated Wild Type, Drug-Treated Transgenic and Vehicle-treated Transgenic) (Fig.  2 A). During WB gel loading, experimenters often distribute their samples unevenly such that certain experimental conditions may be missing on some gels or samples from the same experimental condition are loaded adjacently on a gel. This is problematic because we know that polyacrylamide gel electrophoresis (PAGE) gels are not perfectly uniform, reflecting a source of technical variability 43 ; in the worst case, if we have only loaded a single experimental group on a gel and found a significant effect of the group, we cannot conclude if the effect is due to the experimental condition or a technical problem of the gel. At minimum, experimenters should ensure that every group in a factorial design is represented on each gel to avoid confounding technical gel effects with experimental differences. If the number of combinations is too large to represent on a single gel because of the number of factors or the number of levels of the factors, then a smaller "fractional factorial" design will provide maximal counterbalancing to ensure unbiased estimates of all factor effects and the most important interactions.

In addition, experimenters can further counter technical variability by arranging experimental groups on each gel to ensure adequately counterbalanced design assuming the uniformed protein concentration and fluid volume of all samples. This importantly addresses the variability due to physical effects within an individual gel. In our example, this means alternating the tissue areas and experimental conditions as much as possible to minimize similar samples from being loaded next to one another (Fig.  2 B). By spreading the possibility of technical variability across all samples by counterbalancing across and within gels, we can mitigate potential technical effects that can bias our results. Proper counterbalancing also enables us to implement more rigorous statistical analysis to account for and remove more technical variability 25 , 26 , 32 , 33 . Overall, this will help to ensure that experimenters can find the same result in the future and improve reproducibility.

Technical replication

Technical replicates are used to measure the precision of an assay or method by repeating the measurement of the same sample multiple times. The results of these replicates can then be used to calculate the variability and error of the assay or method 13 . This is important to establish the reliability and accuracy of the results. Most experimenters acknowledge the importance of running technical replicates to avoid false positives and negatives due to technical error 13 . Even beyond extreme results, technical replicates can account for the differences in gel makeup, human variability in gel loading, and potential procedural discrepancies. In fact, most studies run at least duplicates; however, the experimental implementation of replicates (e.g., running replicates on the same gel or separate gels) as well as the statistical analysis of replicates (e.g., dropping “odd-man-out” or taking the mean or standard deviation) can differ greatly 44 , 45 . This experimental variability ultimately impedes our ability to meaningfully compare results. For experimenters to establish accuracy and advance reproducibility in WB experiments, it is important to implement standardized and rigorous protocols to handle technical replicates 11 , 13 . In doing so, we can further reduce the technical variability with statistical methods during analysis.

As underscored previously, we recommend that technical replicates are counterbalanced on separate gels to mitigate any possible gel effect. Additionally, by running triplicates, we can treat replicates as a random effect in a LMM during statistical analysis. Importantly, triplicates provide more values to measure the distribution of technical variance to ensure the robustness of the LMM than only running duplicates. This approach isolates and removes technical variance from biological variation which ultimately improves our sensitivity for true experimental effects 46 .

In the following demonstration of statistical methods, we replicated all WB analyses in triplicate with a randomized counterbalanced design. We then explore how the way in which technical replicates and loading controls are incorporated into analysis can have a significant impact on both the sensitivity of our results and the interpretation of the findings. An example mockup of a dataset illustrating the various ways in which western blot data are typically prepared for analysis can be found in Fig.  3 .

figure 3

Western Blot Gel and Replication Strategies. ( A) Illustration of Western Blot Gel. This depiction of a typical multiplexed western blot gel highlights the antibody-labeled target protein bands of interest (green/yellow) and housekeeping protein loading control that is always run and quantified in the same sample and lane as the target of interest. Total protein stain (fluorescent ponceau stain) is shown in red can can be used as an alternative loading control. Specific, quantification is typically executed on a single antibody-labeled channel for the target protein and housekeeping protein loading control (gray scale image). ( B) Balanced Factorial Technical Replicate Strategy. Here we show the western blot data for the first 3 subjects from an example dataset. In a balanced factorial design, an equal number of samples from all possible experimental groups are represented on each gel. This table shows the subject number, the technical replicate, experimental group, and the band quantifications for both the target protein and the loading control. A ratio of target protein and loading control is also calculated. ( C) Other Common Technical Replicate Strategies. In this example table are two of the other ways western blot data are typically formatted. Some experimenters choose to not include technical replicates, with only one sample from each subject quantified. In another replication strategy, technical replicates are averaged. Averaging may bias or skew the data. We recommend running technical replicates on separate gels or batches, and using gel/batch as a random factor when analyzing western blot data.

Statistical methodology to improve western blot analysis

Loading control as a covariate.

Most qWB assay studies use loading controls (either a housekeeping protein or total protein within lane) to ensure that there are no biases in total protein loaded in a particular lane 2 , 11 , 27 . The most common way that loading controls are used to account for variability between lanes is by normalizing the target protein expression values by dividing it by the loading control values (Fig.  3 ) , resulting in a ratio between target protein to loading control 2 , 47 , 48 . However, ratios may violate assumptions of common statistical test used to analyze qWB (e.g., t-test, ANOVA, etc.) 49 This ultimately hinders the ability to statistically account for the variance in qWB outcomes and have a reliable estimate of the statistics. An alternative approach to improve the parametric properties would be to include loading control values as a covariate—a variable that is not our experimental factors but that may affect the outcome of interest and presents a source of variance that we may account for 50 . For instance, we know the amount of protein loaded is a source of variability in WB quantification, so we can use the loading control as a covariate to adjust for that variance. In doing so, we extend the method of ANOVA into that of ANCOVA 51 . This approach accounts for the technical variability present between lanes while meeting the necessary assumptions for parametric statistics which helps curb bias and averts false discoveries.

Replication and subject as a random effect

Most WB studies use ANOVA, a test that allows comparison of the means of three or more independent samples, for quantitative analysis of WB data 49 . One of the assumptions in ANOVA is the independence of observations 49 . This is problematic because we often collect multiple observations from the same analytical unit, for example different tissue samples from a single subject, or technical replicates. As a result, those observations don’t qualify as independent and should be analyzed using models controlling for variability within units of observations (e.g., the animal) to mitigate inferential errors (false positives and negatives) 52 caused by what is known as pseudoreplication. This arises when the quantity of measured values or data points surpasses the number of actual replicates, and the statistical analysis treats all data points as independent, resulting in their full contribution to the final result 53 .

In addition, when conducting experiments, it is important to consider the randomness of the conditions being observed. Treating both subjects and conditions as fixed effects can lead to inaccurate p-values. Instead, subjects/ animals should be treated as random effects and the conditions should be considered as a sample from a larger population 54 . This is especially important when collecting data from different replicates or gels, as the separate technical replicate runs should be considered as random.

In Fig.  4 we use a simple experimental design comparing the difference in a target protein between two experimental groups to demonstrate four of the most common ways researchers tend to analyze western blot data: (1) running each sample once without replication, (2) treating each technical replicate as an independent sample, (3) taking the mean of technical replicate values, and (4) treating subject and replication as a random effect (Fig.  4 ). We then tested how effect size, power, and p value are affected by each of these strategies to get a sense of how much these estimates vary between analyses. For each of these strategies, we also tested the difference between using the ratio of target protein to loading controls versus using loading control as a statistical covariate. For further exploration of the way these data are prepared and analyzed, see the data workup in Supplementary Figs.  1 and 2 .

figure 4

Effect of different replication and loading control strategies on statistical outcomes. Eight possible strategies are shown, representing the most common ways in which replication and loading controls are treated in a typical Western blot analysis. Four replication strategies: either no replication at all, 3 technical replicate gels treated as independent, mean of three replicates, or replicate treated as a random effect in a linear mixed model. These are crossed with two loading control strategies: either target protein is divided by loading control, or loading control is treated as a covariate in a linear mixed model. ( A) Effect Size: Standardized effect size coefficient is generally improved when loading control is treated as a covariate, compared to using a ratio of the target protein and loading control values. ( B) Power: By treating each replication as independent the statistical power is increased (due to the inaccurate assumption that technical replicates are not related, thus artificially tripling the n). Conversely, including the variability inherent in technical replicates as a part of the statistical model, we work to identify and account for a major source of variability, thus improving power in a more appropriate way. ( C) P value: As expected, when each replication is inaccurately treated as independent the p value is low (due to artificially inflated n). We found that using the mean of replications and loading controls as covariates also resulted in a p value below 0.05. The smallest p value was found when including replication as a random factor. Across each of these statistical measures, only when replication is included as a random factor and loading control as a covariate do we see a strong effect size, high power, and low p value.

In the first scenario, we imagined that no technical replication was run at all (by using only the first replication). With this strategy, we found that standardized effect size is weak, power is low, and the p value was high (Fig.  4 ). Second, we demonstrate how analytical output would be different if we did run three technical replicates, but treated each as independent. As discussed above, this strategy does not take into account the fact that each sample is being run three times, and consequently the overall n of your experiment is artificially tripled! As one might expect, observed power is quite high, and our p value is low (< 0.05). Power is increased by an increase in sample size, so it is not surprising that the power is much higher if we erroneously report that we have a 3X larger sample size (i.e., pseudoreplication) 53 . In this case, the observed power is inflated and an artifact of inappropriate statistics, and the probability of a false positive is considerably increased with respect to the expected 5%.

So, what would be a more appropriate way to handle technical replicates? One method that researchers often use is to take the mean of their technical replicates. This does ensure that we are not artificially inflating our sample size, which is certainly an improvement over the previous strategy. With this strategy, we do find that our p value is less than 0.05 (when loading control is treated as a covariate). But we also see that our power is still low. We have effectively taken our replicates into account by collapsing across them within each sample, but this can be dangerous. If there is wide variation across replicates of a particular sample, then taking the mean of three replicates could produce an inaccurate estimate of the ‘true’ sample value. Ideally, we want to find a solution where instead of collapsing this variation, we add it to our statistical model so that we can better understand what amount of variation is randomly coming from within technical replicates, and in turn what amount of variation is actually due to potential differences in our experimental groups.

To achieve this, we need to model both the fixed effect of all groups in a full factorial design, and the random effect of replication across western blot gels. When we use both fixed and random effects, this is referred to as a linear mixed model (LMM). When using this strategy, we find that our effect size remains strong, and our p value is low. But importantly, we now have strong observed power (Fig.  4 ). This suggests that we can achieve greater sensitivity in our WB experiment when using this approach . Specifically, if we implement careful counterbalancing while designing our experiments, then we can use the variability between gels to our advantage during analysis using linear mixed effects model 55 .

LMM is recommended because it takes into account both the multiple observations within a single subject/animal in a given condition and differences across subjects observed in multiple conditions. This reduces chances of inaccurate p-values and improves reliability 56 . Further, treating both subjects and replication as random effects generalizes the results to the population of subjects and also to the population of conditions 57 .

Real world application of blotRig software for western blot experimental design, technical replication, and statistical analysis

We have designed a user interface that is designed to facilitate appropriate counterbalancing and technical replication for western blot experimental design. The ‘blotRig’ application is run through RStudio, and can be found here: https://atpspin.shinyapps.io/BlotRig/ Upon starting the blotRig application, the user is prompted to upload a comma separated values (CSV) spreadsheet. This spreadsheet should include separate columns for subject ID and experimental group. The user is then prompted to enter the total number of lanes that are available on their particular western blot gel apparatus. The blotRig software will first run a quality check to confirm that each subject ID (unique sample or subject) is only found in one experimental group. If duplicates are found, a warning will be shown that specifies which subjects are repeated across groups. If no errors are found, a centered gel map will be generated that illustrates the western blot gel lanes into which each subject should be loaded (Fig.  5 A). The decision for each lane loading is based on two main principles outlined above: (1) each western blot gel should hold a representative sample of each experimental group (2) samples from the same experimental group are not loaded in adjacent lanes whenever possible. This ensures that proper counterbalancing is achieved so that we can limit the chances that the inherent variability within and across western blot gels is confounded with the experimental groups that we are interested in experimentally testing.

figure 5

Example of the blotRig Gel Creator interface. ( A ) Illustration of the blotRig interface. User has entered their sample IDs, experimental groups, and the number of lanes per western blot gel. ( B) The blotRig system then creates a counterbalanced gel map that ensures each gel contains a representative from each experimental group. This illustration shows the exact lane for each gel in which each sample should be run.

Once the gel map has been generated, the user can then select to export this gel map to a CSV spreadsheet. This sheet is designed to clearly show which gel each sample is on, which lane on each gel a sample is found, what experimental group each sample belongs to, and importantly, a repetition of each of these values for three technical replicates (Fig.  5 B). User will also see columns for Target Protein and Loading Control. These are the cells where the user can then input their densitometry values upon completing their western blot runs. Once this spreadsheet is filled out, it is then ready to go for blotRig analysis.

To analyze western blot data, users can upload the completed template that was exported in the blotRig experimental design phase or their own CSV file under the ‘Analysis’ tab (Fig.  6 ). The blotRig software will first ask the user to identify which columns from the spreadsheet represent Subject/SampleID, Experimental Group, Protein Target, Loading Control, and Replication. The blotRig software will again run a quality check to confirm that there are no subject/sample IDs that are duplicated across experimental groups. If no errors are found, the data will then be ready to analyze. The blotRig analysis will then be run, using the principles discussed above. Specifically, a linear mixed-model runs using the lmer R package, with Experimental Group as a fixed effect, Loading Control as a covariate, and Replication (nested within Subject/Sample ID) as a random factor. Analytical output is then displayed, giving a variety of statistical results from the linear mixed model output table, including fixed and random effects and associated p values (Fig.  6 ). A bar graph of group means and 95% confidence interval error bars will also be generated, along with a summary of the group means, standard error of the mean, and upper/lower 95% confidence intervals. These outputs can be directly reported in the results sections of papers, improve the statistical rigor of published WB reports. In addition, since the entire pipeline is opensource, the blotRig code itself can be reported to support transparency and reproducibility.

figure 6

Workflow for running statistical analysis of replicate western blot data using blotRig. First, fill out spreadsheet with subject ID, experimental group assignment, number of technical replication, the densitometry values for your target proteins and loading controls. After saving this spreadsheet as a.csv file, the file can be uploaded to blotRig. Tell blotRig the exact names of each of your variables, then click ‘Run Analysis’. This will produce a statistical output using linear mixed model testing for group differences using loading control as a covariate and replication as a random effect. Bar graph with error bars and summary statistics can then be exported.

Although the western blot technique has proven to be a workhorse for biological research, the need to enhance its reproducibility is critical 13 , 19 , 27 . Current qWB assay methods are still lacking for reproducibly identifying true biological effects 13 . We provide a systematic approach to generate quantitative data from western blot experiments that incorporates key technical and statistical recommendations which minimize sources of error and variability throughout the western blot process. First, our study shows that experimenters can improve the reproducibility of western blots by applying the experimental recommendations of determining the linear range for each primary antibody, counterbalancing during experimental design, and running technical triplicates 13 , 27 . Furthermore, these experimental implementations allow for application of the statistical recommendations of incorporating loading controls as covariates and analyzing gel and subject as random effects 58 , 59 . Altogether, these enable more rigorous statistical analysis that accounts for more technical variability which can improve the effect size, observed power, and p-value of our experiments and ultimately better identify true biological effects.

Biomedical research has continued to rely on p-values for determining and reporting differences between experimental groups, despite calls to retire the p-value 60 . Power (sensitivity) calculations have also become increasingly common. In brief, p-values and the related alpha value are associated with Type I error rate—the probability of rejecting the null hypothesis (i.e., claiming there is an effect) when there is no true effect 61 . On the other hand, power effectively measures the probability of rejecting the null hypothesis (i.e. stating there is not effect) when there is indeed a true underlying effect—a concept that is closely related to reducing the Type II error rate 59 , 62 . Critically, empirical evidence estimates that the median statistical power of studies in neuroscience is between ∼ 8% and ∼ 31%, yet best practices suggest that an experimenter should aim to achieve a power of 80% with an alpha of 0.05 20 . By being underpowered, experiments are at higher likelihood of producing a false inference. If an underpowered experiment is seeking to reproduce a previous observation, the resulting false negative may throw into question the original findings and directly exacerbate the reproducibility crisis 59 . Even more alarmingly, a low power also increases the likelihood that a statistically significant result is actually a false positive due to small sample size problems 61 . In our analyses, we show that our technical and statistical recommendations lower the p-value (indicating that the observed relationship between variables is less likely to be due to chance) as well as observed power of our experiments. This translates into the ability to better avoid false negatives when there is a true effect as well as reduce the likelihood of false positives when there is not a true experimental effect, both of which will ultimately improve the reproducibility of qWB assay experiments.

Another useful component of statistical analyses that is not as commonly reported but is critically related to p-value and power is effect size. Effect size is a statistical measure that describes the magnitude of the difference between two groups in an experiment 63 . It is used to quantify the strength of the relationship between the variables being studied 63 . The estimated effect size is important because it answers the most frequent question that researchers ask: how big is the difference between experimental groups, or how strong is the relationship or association? 63 . The combination of standardized effect size, p-value and power reflect crucial experimental results that can be broadly understood and compared with findings from other studies 62 , thus improving comparability of qWB experiments 49 , 63 . In particular, studies with large effect sizes have more power: we are more likely to detect a true positive experimental effect and avoid the false negative if the underlying difference between experimental groups is large 46 . In some cases, the calculated effect size is greatly influenced by how sources of variance are handled during analysis 13 . Our results demonstrate that by reducing the residual variance (by modeling the random effect of replication) the estimated effect size of our experiment increases. This could mean that the magnitude of the difference between the groups in our experiment is larger than it was originally thought to be. This could be due to a variety of factors such as improving the experimental design, sample size, or the measurement of the variables 13 . Likewise, conducting a power analysis is an essential step in experimental design that should be done before collecting data to ensure that the study is adequately powered to detect an effect of a certain size 64 .

Increasingly, power analysis is becoming a requirement for publications and grant proposals 65 . This is because a study with low statistical power is more likely to produce false negative results, which means that the study may fail to detect a real effect that actually exists. This can lead to the rejection of true hypotheses, wasted resources, and potentially harmful conclusions. In brief, given an experimental effect size and variance, we can calculate the sample size needed to achieve an alpha of 0.05 and power of 0.8; an increased sample size reduces the standard error of mean (SEM), which is the measured spread of sample means and consequently increases the power of the experiment 66 . We have demonstrated that our experimental and statistical recommendations lead to a lower p value (Fig.  3 C) and effect size (Fig.  3 B) without changing the sample size. This may be of greatest interest to researchers: more rigorous analytics ultimately improves experimental sensitivity without relying solely on increasing the sample size.

Reducing the sample size of an experiment can be beneficial for several reasons, one of which is cost-effectiveness. A smaller sample size can lead to a reduction in the number of animals or other resources that are needed for the study, which can result in lower costs. Additionally, it can also save time and reduce the duration of the experiment, as fewer subjects need to be recruited, and the data collection process can be completed more quickly. However, it is important to note that reducing the sample size can also lead to decreased statistical power. As a result, reducing sample size too much can increase the risk of a type II error, failing to detect significance when there is a true effect 62 .Therefore, it is important to consider the trade-off between sample size and power when designing an experiment, and to use statistical techniques like power analysis to ensure that the sample size is sufficient to detect an effect of a certain size. Moreover, when using animals in research, it's always important to consider the ethical aspect and the 3Rs principles of reduction, refinement, and replacement 55 .

Despite our best efforts in creating a balanced, full factorial experimental design, there will always be random variation in biological experiments. Fixed effects such as experimental group differences are expected to be generalizable if the experiment is replicated. Random effects (such as gel variation) on the other hand are unpredictable across experiments. Western blot analyses are particularly susceptible to this random gel variation, as different values may be observed for technical replicates run on different gels. By using a linear mixed model paired with rigorous full factorial design, we can ensure that we account for as much of that random variation as possible. When we acknowledge, identify, and model random effects we enhance the possibility of discovering our fixed effect of experimental treatment, if one exists.

The linear mixed model framework discussed above assumes that our western blot outcome measures are on a linear scale. As described above, parametric work to identify the linear range of a protein of interest is critical for ensuring that the results of a LMM (or ANOVA and t-test) are accurate and interpretable. While we recommend using loading control (or total protein control) as a covariate in a linear mixed model, many bench researchers may prefer to use the within-lane loading control (or total protein) to normalize target protein values. It is important to consider that in doing so, one creates a ratio value that is multiplicative instead of linear. This property has the side effect of artificially distorting the variance. To account for this non-linearity, we recommend that one uses semi-parametric mixed models such as generalized estimating equations with a gamma distribution link function that appropriately represents ratio data.

There has been recent recognition that an appropriate study design can be achieved by balancing sample size (n), effect size, and power 31 . The experimental and statistical approach presented in this study provide insight into how more rigorous planning for western blot experimental design and corresponding statistical analysis without depending on p-values only can acquire precise data resulting in true biological effects. Using blotRig as a standardized, integrated western blot methodology, quantitative western blot may become highly reproducible, reliable, and a less controversial protein measurement technique 18 , 28 , 67 .

Study reporting

This study is reported in accordance with ARRIVE guidelines.

Supporting information

This article contains supporting information. You can access the blotRig analysis software, which includes code for inputting experimental parameters for all Western blot analysis, through the following link: https://atpspin.shinyapps.io/BlotRig/

Data availability

The datasets and computer code generated or used in this study are accessible in a public, open-access repository at https://doi.org/10.34945/F51C7B and https://github.com/ucsf-ferguson-lab/blotRig/ respectively.

Abbreviations

American association for accreditation of laboratory animal care

Animal research reporting of in vivo experiments

American veterinary medical association

Institutional animal care and use committee

Quantitative western blot

Enzyme linked immunosorbent assay

Severe acute respiratory syndrome coronavirus 2

Analysis of covariance

Analysis of variance

Spinal cord injury

Spared nerve injury

α-Amino-3-hydroxy-5-methyl-4-isoxazoleproprionic acid

Glutamate receptor 1

Glutamate receptor 2

Linear mixed models

Tris-buffered saline containing 0.1% Tween 20

Polyacrylamide gel electrophoresis

Standard error of mean

Lowry, O., Rosebrough, N., Farr, A. L. & Randall, R. Protein measurement with the Folin phenol reagent. J. Biol. Chem. 193 , 265–275. https://doi.org/10.1016/S0021-9258(19)52451-6 (1951).

PubMed   Google Scholar  

Aldridge, G. M., Podrebarac, D. M., Greenough, W. T. & Weiler, I. J. The use of total protein stains as loading controls: An alternative to high-abundance single protein controls in semi-quantitative immunoblotting. J. Neurosci. Methods 172 , 250–254. https://doi.org/10.1016/j.jneumeth.2008.05.00 (2008).

PubMed   PubMed Central   Google Scholar  

McDonough, A. A., Veiras, L. C., Minas, J. N. & Ralph, D. L. Considerations when quantitating protein abundance by immunoblot. Am. J. Physiol. Cell Physiol. 308 , C426-433. https://doi.org/10.1152/ajpcell.00400.2014 (2015).

Towbin, H., Staehelin, T. & Gordon, J. Electrophoretic transfer of proteins from polyacrylamide gels to nitrocellulose sheets: Procedure and some applications. PNAS 76 , 4350–4354. https://doi.org/10.1073/pnas.76.9.4350 (1979).

ADS   PubMed   PubMed Central   Google Scholar  

Burnette, W. N. “Western blotting”: Electrophoretic transfer of proteins from sodium dodecyl sulfate-polyacrylamide gels to unmodified nitrocellulose and radiographic detection with antibody and radioiodinated protein A. Anal. Biochem. 112 , 195–203. https://doi.org/10.1016/0003-2697(81)90281-5 (1981).

Mahmood, T. & Yang, P.-C. Western blot: Technique, theory, and trouble shooting. N. Am. J. Med. Sci. 4 , 429–434. https://doi.org/10.4103/1947-2714.100998 (2012).

Alegria-Schaffer, A., Lodge, A. & Vattem, K. Performing and optimizing Western blots with an emphasis on chemiluminescent detection. Methods Enzymol. 463 , 573–599. https://doi.org/10.1016/S0076-6879(09)63033-0 (2009).

Khoury, M. K., Parker, I. & Aswad, D. W. Acquisition of chemiluminescent signals from immunoblots with a digital SLR camera. Anal. Biochem. 397 , 129–131. https://doi.org/10.1016/j.ab.2009.09.041 (2010).

Zellner, M. et al. Fluorescence-based western blotting for quantitation of protein biomarkers in clinical samples. Electrophoresis 29 , 3621–3627. https://doi.org/10.1002/elps.200700935 (2008).

Gingrich, J. C., Davis, D. R. & Nguyen, Q. Multiplex detection and quantitation of proteins on western blots using fluorescent probes. Biotechniques 29 , 636–642. https://doi.org/10.2144/00293pf02 (2000).

Janes, K. A. An analysis of critical factors for quantitative immunoblotting. Sci. Signal 8 , rs2. https://doi.org/10.1126/scisignal.2005966 (2015).

Mollica, J. P., Oakhill, J. S., Lamb, G. D. & Murphy, R. M. Are genuine changes in protein expression being overlooked? Reassessing western blotting. Anal. Biochem. 386 , 270–275. https://doi.org/10.1016/j.ab.2008.12.029 (2009).

Pillai-Kastoori, L., Schutz-Geschwender, A. R. & Harford, J. A. A systematic approach to quantitative western blot analysis. Anal. Biochem. 593 , 113608. https://doi.org/10.1016/j.ab.2020.113608 (2020).

Aydin, S. A short history, principles, and types of ELISA, and our laboratory experience with peptide/protein analyses using ELISA. Peptides 72 , 4–15. https://doi.org/10.1016/j.peptides.2015.04.012 (2015).

Seisenberger, C. et al. Questioning coverage values determined by 2D western blots: A critical study on the characterization of anti-HCP ELISA reagents. Biotechnol. Bioeng. 118 , 1116–1126. https://doi.org/10.1002/bit.27635 (2021).

Edwards, V. M. & Mosley, J. W. Reproducibility in quality control of protein (western) immunoblot assay for antibodies to human immunodeficiency virus. Am. J. Clin. Pathol. 91 , 75–78. https://doi.org/10.1093/ajcp/91.1.75 (1989).

Matschke, J. et al. Neuropathology of patients with COVID-19 in Germany: A post-mortem case series. Lancet Neurol. 19 , 919–929. https://doi.org/10.1016/S1474-4422(20)30308-2 (2020).

Murphy, R. M. & Lamb, G. D. Important considerations for protein analyses using antibody based techniques: Down-sizing western blotting up-sizes outcomes. J. Physiol. 591 , 5823–5831. https://doi.org/10.1113/jphysiol.2013.263251 (2013).

Butler, T. A. J., Paul, J. W., Chan, E.-C., Smith, R. & Tolosa, J. M. Misleading westerns: Common quantification mistakes in western blot densitometry and proposed corrective measures. Biomed. Res. Int. 2019 , 5214821. https://doi.org/10.1155/2019/5214821 (2019).

Button, K. S. et al. Power failure: Why small sample size undermines the reliability of neuroscience. Nat. Rev. Neurosci. 14 , 365–376. https://doi.org/10.1038/nrn3475 (2013).

Landis, S. C. et al. A call for transparent reporting to optimize the predictive value of preclinical research. Nature 490 , 187–191. https://doi.org/10.1038/nature11556 (2012).

Shields, S. D., Eckert, W. A. & Basbaum, A. I. Spared nerve injury model of neuropathic pain in the mouse: A behavioral and anatomic analysis. J. Pain 4 , 465–470. https://doi.org/10.1067/s1526-5900(03)00781-8 (2003).

Decosterd, I. & Woolf, C. Spared nerve injury: An animal model of persistent peripheral neuropathic pain. Pain 87 , 149–158. https://doi.org/10.1016/S0304-3959(00)00276-1 (2000).

Richner, M., Jager, S. B., Siupka, P. & Vaegter, C. B. Hydraulic extrusion of the spinal cord and isolation of dorsal root ganglia in rodents. J. Vis. Exp. https://doi.org/10.3791/55226 (2017).

Ferguson, A. R. et al. Cell death after spinal cord injury is exacerbated by rapid TNFα-induced trafficking of GluR2-lacking AMPARS to the plasma membrane. J Neurosci 28 , 11391–11400. https://doi.org/10.1523/JNEUROSCI.3708-08.2008 (2008).

Ferguson, A. R., Huie, J. R., Crown, E. D. & Grau, J. W. Central nociceptive sensitization vs. spinal cord training: Opposing forms of plasticity that dictate function after complete spinal cord injury. Front. Physiol. 3 , 1. https://doi.org/10.3389/fphys.2012.00396 (2012).

Google Scholar  

Taylor, S. C., Berkelman, T., Yadav, G. & Hammond, M. A defined methodology for reliable quantification of western blot data. Mol. Biotechnol. 55 , 217–226. https://doi.org/10.1007/s12033-013-9672-6 (2013).

Bakkenist, C. J. et al. A quasi-quantitative dual multiplexed immunoblot method to simultaneously analyze ATM and H2AX phosphorylation in human peripheral blood mononuclear cells. Oncoscience 2 , 542–554. https://doi.org/10.18632/oncoscience.162 (2015).

Wang, Y. V. et al. Quantitative analyses reveal the importance of regulated Hdmx degradation for p53 activation. Proc. Natl. Acad. Sci. USA 104 , 12365–12370. https://doi.org/10.1073/pnas.0701497104 (2007).

Bass, J. et al. An overview of technical considerations for western blotting applications to physiological research. Scand. J. Med. Sci. Sports 27 , 4–25. https://doi.org/10.1111/sms.12702 (2017).

Lazzeroni, L. C. & Ray, A. The cost of large numbers of hypothesis tests on power, effect size and sample size. Mol. Psychiatry 17 , 108–114. https://doi.org/10.1038/mp.2010.117 (2012).

Huie, J. R. et al. AMPA receptor phosphorylation and synaptic colocalization on motor neurons drive maladaptive plasticity below complete spinal cord injury. eNeuro https://doi.org/10.1523/ENEURO.0091-15.2015 (2015).

Stück, E. D. et al. Tumor necrosis factor alpha mediates GABAA receptor trafficking to the plasma membrane of spinal cord neurons in vivo. Neural Plast https://doi.org/10.1155/2012/261345 (2012).

Krzywinski, M. & Altman, N. Points of significance: Power and sample size. Nat. Method. 10 , 1139–1140. https://doi.org/10.1038/nmeth.2738 (2013).

R Core Team R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/ (2021).

Green, P. & MacLeod C. J. “simr: An R package for power analysis of generalised linear mixed models by simulation.” Meth. Ecol. Evolut. 7 (4), 493–498. https://doi.org/10.1111/2041-210X.12504 , https://CRAN.R-project.org/package=simr (2016).

Attali, D. shinyjs: Easily Improve the User Experience of Your Shiny Apps in Seconds. R package version 2.1.0, https://deanattali.com/shinyjs/ (2022).

Chang, W. et al. shiny: Web Application Framework for R. R package version 1.9.1.9000, https://github.com/rstudio/shiny , https://shiny.posit.co/ (2024).

Chang, W. shinythemes: Themes for Shiny. R package version 1.2.0, https://github.com/rstudio/shinythemes (2024).

de Vries, A., Schloerke, B., Russell, K. sortable: Drag-and-Drop in ‘shiny’ Apps with ‘SortableJS’. R package version 0.5.0, https://github.com/rstudio/sortable (2024).

Wickham, H. et al. Welcome to the tidyverse. JOSS 4 (43), 1686. https://doi.org/10.21105/joss.01686 (2019).

ADS   Google Scholar  

Xie, Y., Cheng, J., Tan, X. DT: A Wrapper of the JavaScript Library ‘DataTables’. R package version 0.33.1, dt. https://github.com/rstudio/ (2024).

Krzywinski, M. & Altman, N. Points of significance: Analysis of variance and blocking. Nat Methods 11 , 699–700. https://doi.org/10.1038/nmeth.3005 (2014).

Heidebrecht, F., Heidebrecht, A., Schulz, I., Behrens, S.-E. & Bader, A. Improved semiquantitative western blot technique with increased quantification range. J. Immunol. Methods 345 , 40–48. https://doi.org/10.1016/j.jim.2009.03.018 (2009).

Huang, Y.-T. et al. Robust comparison of protein levels across tissues and throughout development using standardized quantitative western blotting. J. Vis. Exp. https://doi.org/10.3791/59438 (2019).

Krzywinski, M. & Altman, N. Points of view: Designing comparative experiments. Nat. Methods 11 , 597–598. https://doi.org/10.1038/nmeth.2974 (2014).

Thacker, J. S., Yeung, D. H., Staines, W. R. & Mielke, J. G. Total protein or high-abundance protein: Which offers the best loading control for western blotting?. Anal. Biochem. 496 , 76–78. https://doi.org/10.1016/j.ab.2015.11.022 (2016).

Zeng, L. et al. Direct blue 71 staining as a destaining-free alternative loading control method for western blotting. Electrophoresis 34 , 2234–2239. https://doi.org/10.1002/elps.201300140 (2013).

Jaeger, T. F. Categorical data analysis: Away from ANOVAs (transformation or not) and towards logit mixed models. J. Mem. Lang. 59 , 434–446. https://doi.org/10.1016/j.jml.2007.11.007 (2008).

Mefford, J. & Witte, J. S. The covariate’s dilemma. PLoS Genet. 8 , e1003096. https://doi.org/10.1371/journal.pgen.1003096 (2012).

Schneider, B. A., Avivi-Reich, M. & Mozuraitis, M. A cautionary note on the use of the analysis of covariance (ANCOVA) in classification designs with and without within-subject factors. Front. Psychol. 6 , 474. https://doi.org/10.3389/fpsyg.2015.00474 (2015).

Nieuwenhuis, S., Forstmann, B. U. & Wagenmakers, E.-J. Erroneous analyses of interactions in neuroscience: A problem of significance. Nat. Neurosci. 14 , 1105–1107. https://doi.org/10.1038/nn.2886 (2011).

Freeberg, T. M. & Lucas, J. R. Pseudoreplication is (still) a problem. J. Com. Psychol. 123 , 450–451. https://doi.org/10.1037/a0017031 (2009).

Judd, C. M., Westfall, J. & Kenny, D. A. Treating stimuli as a random factor in social psychology: A new and comprehensive solution to a pervasive but largely ignored problem. J. Pers. Soc. Psychol. 103 , 54–69. https://doi.org/10.1037/a0028347 (2012).

Lee, O. E. & Braun, T. M. Permutation tests for random effects in linear mixed models. Biometrics 68 , 486–493. https://doi.org/10.1111/j.1541-0420.2011.01675.x (2012).

MathSciNet   PubMed   Google Scholar  

Baayen, R. H., Davidson, D. J. & Bates, D. M. Mixed-effects modeling with crossed random effects for subjects and items. J. Mem. Lang 59 , 390–412. https://doi.org/10.1016/j.jml.2007.12.005 (2008).

Barr, D. J., Levy, R., Scheepers, C. & Tily, H. J. Random effects structure for confirmatory hypothesis testing: Keep it maximal. J. Mem. Lang. https://doi.org/10.1016/j.jml.2012.11.001 (2013).

Blainey, P., Krzywinski, M. & Altman, N. Points of significance: Replication. Nat. Methods 11 , 879–880. https://doi.org/10.1038/nmeth.3091 (2014).

Drubin, D. G. Great science inspires us to tackle the issue of data reproducibility. Mol. Biol. Cell 26 , 3679–3680. https://doi.org/10.1091/mbc.E15-09-0643 (2015).

Amrhein, V., Greenland, S. & McShane, B. Scientists rise up against statistical significance. Nature 567 , 305–307. https://doi.org/10.1038/d41586-019-00857-9 (2019).

ADS   PubMed   Google Scholar  

Cohen, J. The earth is round (p <.05). Am. Psychol. 49 , 997–1003. https://doi.org/10.1037/0003-066X.49.12.997 (1994).

Ioannidis, J. P. A., Tarone, R. & McLaughlin, J. K. The false-positive to false-negative ratio in epidemiologic studies. Epidemiology 22 , 450–456. https://doi.org/10.1097/EDE.0b013e31821b506e (2011).

Sullivan, G. M. & Feinn, R. Using effect size-or why the P value Is not enough. J. Grad. Med. Educ. 4 , 279–282. https://doi.org/10.4300/JGME-D-12-00156.1 (2012).

Brysbaert, M. & Stevens, M. Power analysis and effect size in mixed effects models: A tutorial. J. Cogn. 1 , 9. https://doi.org/10.5334/joc.10 (2018).

Kline, R. B. Beyond significance testing: Reforming data analysis methods in behavioral research. Am. Psychol. Associat . https://doi.org/10.1037/10693-000 (2024).

Rosner, Bernard (Bernard A.). Fundamentals of biostatistics. (Boston, Brooks/Cole, Cengage Learning, 2011).

Bromage, E., Carpenter, L., Kaattari, S. & Patterson, M. Quantification of coral heat shock proteins from individual coral polyps. Mar. Ecol. Progress Ser. 376 , 123–132 (2009).

Download references

Acknowledgements

The authors would like to thank Alexys Maliga Davis for data librarian services.

This work was supported by a National Institutes of Health/National Institute of Neurological Disorders and Stroke grant (R01NS088475) to A. R. F. NIH NINDS: R01NS122888, UH3NS106899, U24NS122732, US Veterans Affairs (VA): I01RX002245, I01RX002787, I50BX005878, Wings for Life Foundation, Craig H. Neilsen Foundation. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. Correspondence and requests for materials should be addressed to A.R.F.

Author information

Authors and affiliations.

Weill Institute for Neurosciences, University of California, San Francisco, CA, USA

Cleopa Omondi, Austin Chou, Kenneth A. Fond, Kazuhito Morioka, Nadine R. Joseph, Jeffrey A. Sacramento, Emma Iorio, Abel Torres-Espin, Hannah L. Radabaugh, Jacob A. Davis, Jason H. Gumbel, J. Russell Huie & Adam R. Ferguson

San Francisco Veterans Affairs Medical Center, San Francisco, CA, USA

J. Russell Huie & Adam R. Ferguson

School of Public Health Sciences, Faculty of Health Sciences, University of Waterloo, Waterloo, ON, Canada

Abel Torres-Espin

Department of Physical Therapy, Faculty of Rehabilitation Medicine, University of Alberta, Edmonton, AB, Canada

You can also search for this author in PubMed   Google Scholar

Contributions

C.O: Writing-original draft preparation, Investigation, Validation, Data Curation, Visualization, Formal Analysis, Writing-Review & Editing; A. C: Formal analysis, Writing-Review & Editing; K. A. F: Software, Writing-Review & Editing; K. M: Methodology, Writing-Review & Editing; N. R. J: Writing—Review & Editing; J. A. S: Investigation, Project Administration, Writing—Review & Editing; E. I: Resources, Writing-Review & Editing; A.T.E: Software, Writing—Review & Editing; H. L. R: Software, Writing—Review & Editing; J. A. D: Investigation, Writing—Review & Editing; J. H. G: Investigation, Writing – Review & Editing; J. R. H: Conceptualization, Methodology, Validation, Formal Analysis, Investigation, Data Curation, Writing- Review & Editing, Visualization, Supervision; A. R. F:Conceptualization, Methodology, Validation, Formal Analysis, Resources, Investigation, Data Curation, Writing-Review & Editing, Visualization, Supervision, Project Administration, Funding Acquisition.

Corresponding authors

Correspondence to J. Russell Huie or Adam R. Ferguson .

Ethics declarations

Competing interests.

The authors declare no competing interests.

Additional information

Publisher's note.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Supplementary figure 1., supplementary figure 2., rights and permissions.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/ .

Reprints and permissions

About this article

Cite this article.

Omondi, C., Chou, A., Fond, K.A. et al. Improving rigor and reproducibility in western blot experiments with the blotRig analysis. Sci Rep 14 , 21644 (2024). https://doi.org/10.1038/s41598-024-70096-0

Download citation

Received : 19 December 2023

Accepted : 13 August 2024

Published : 17 September 2024

DOI : https://doi.org/10.1038/s41598-024-70096-0

Share this article

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

  • Analytical chemistry
  • Biostatistics
  • Computational biology
  • Computational chemistry

By submitting a comment you agree to abide by our Terms and Community Guidelines . If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Quick links

  • Explore articles by subject
  • Guide to authors
  • Editorial policies

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

design and analysis of experiments in r

{{errorText}}

Visit www.aminer.cn directly for more AI services.

Go to AMiner

error code: {{errorCode}}

design and analysis of experiments in r

Logging off

Search for peer-reviewed journals with full access.

Discover the SciOpen Platform and Achieve Your Research Goals with Ease.

Experimental and computational analysis of the structure-activity relationship of ionic gel electrolytes based on bistrifluoromethanesulfonimide salts for supercapacitors

Peer review under responsibility of The Chinese Ceramic Society.

Graphical Abstract

design and analysis of experiments in r

Ionic gel (IG) electrolytes are emerging as promising components for the development of next-generation supercapacitors (SCs), offering benefits in terms of safety, cost-effectiveness, and flexibility. The ionic conductivity, stability, and mechanical properties of the gel electrolyte are relevant factors to be considered and the key to improving the performance of the SC. However, the structure–activity relationship between the internal structure of IGs and their SC properties is not fully understood. In the current study, the intuitive and regular structure–activity relationship between the structure and properties of IGs was revealed via combining computational simulation and experiment. In terms of conductivity, the ionic liquid (IL) ([EMIM][TFSI]) in the IG has a high self-diffusion coefficient calculated by molecular dynamics simulation (MDS), which is conductive to transfer and then improves the conductivity. The radial distribution function of the MDS shows that the larger the g ( r ) between the particles in the polymer network, the stronger the interaction. For stability, IGs based on [EMIM][TFSI] and [EOMIM][TFSI] ILs have higher density functional theory calculated binding energy, which is reflected in the excellent thermal stability and excellent capacitor cycle stability. Based on the internal pore size distribution and stress-strain characterization of the gel network ([ME3MePy][TFSI] and [BMIM][TFSI] as additives), the highly crosslinked aggregate network significantly reduces the internal mesoporous distribution and plays a leading role in improving the mechanical properties of the network. By using this strategy, it will be possible to design the ideal structure of the IG and achieve excellent performance.

No abstract is available for this article. Click the button above to view the PDF directly.

Hu O, Lu J, Chen G, Chen K, Gu J, Weng S, et al. An antifreezing, tough, rehydratable, and thermoplastic poly(vinyl alcohol)/sodium alginate/poly(ethylene glycol) organohydrogel electrolyte for flexible supercapacitors. ACS Sustainable Chem Eng 2021;9:9833–45.

Chen Z, Shen H, Zhu Y, Hua M, Pan X, Liu Y, et al. Advanced low-flammable pyrrole ionic liquid electrolytes for high safety lithium-ion batteries. J Energy Storage 2023;72:108289.

Li J, Qiao J, Lian K. Hydroxide ion conducting polymer electrolytes and their applications in solid supercapacitors: a review. Energy Storage Mater 2020;24:6–21.

Bhat TS, Patil PS, Rakhi RB. Recent trends in electrolytes for supercapacitors. J Energy Storage 2022;50:104222.

Qian J, Jin B, Li Y, Zhan X, Hou Y, Zhang Q. Research progress on gel polymer electrolytes for lithium-sulfur batteries. J Energy Chem 2021;56:420–37.

Cheng Y, Zhu H, Li S, Xu M, Li T, Yang X, et al. Stretchable, low-hysteresis, and recyclable ionogel by ionic liquid catalyst and mixed ionic liquid-induced phase separation. ACS Sustainable Chem Eng 2023;11:15031–42.

Barbosa JC, Pinto RS, Correia DM, Tubio CR, Gonçalves R, Costa CM, et al. Solid polymer electrolytes based on a high dielectric polymer and ionic liquids for lithium batteries. J Power Sources 2023;585:233630.

Liu J, Mei X, Peng F. Lignin derived porous carbon with favorable mesoporous contributions for highly efficient ionic liquid-based supercapacitors. Chin Chem Lett 2023;34:108187.

Zaccagnini P, Serrapede M, Armandi M, Bianco S, Carminati S, Zampato M, et al. A high-temperature high-pressure supercapacitor based on ionic liquids for harsh environment applications. Electrochim Acta 2023;447:142124.

Zhao R, Xu X, Zhou Y, Wang Z, Zhou Y. Study on the structural characteristics and interaction mechanisms of ionic liquid mixtures with a common imidazolium cation. J Mol Liq 2023;380:121743.

Zhang Q, Liu D, Li Q, Zhang X, Wei Y, Lang X. Thermodynamic properties, excess properties, and molecular interactions of ionic liquids 1-cyanopropyl-3-methyl-imidazolium bis(fluorosulfonyl)imide/trifluoromethanesulfonate and binary systems containing acetonitrile. J Mol Liq 2018;268:770–80.

Marzouki M, Zarrougui R, Ghodbane O. Application of aprotic ionic liquids based on bis(trifluoromethylsulfonyl)imide anion as polymer gel electrolytes for cobalt oxide symmetric supercapacitors. J Energy Storage 2021;40:102761.

Wei Y, Chen W, Ge X, Liang J, Xing Z, Zhang Q, et al. A flexible, highly conductive, tough ionogel electrolyte containing LiTFSI salt and ionic liquid [EMIM][TFSI] based on PVDF-HFP for high-performance supercapacitors. Polymer 2023;289:126501.

Sang J, Tang B, Qiu Y, Fang Y, Pan K, Zhou Z. How does stacking pressure affect the performance of solid electrolytes and all-solid-state lithium metal batteries? Energy Environ Mater 2023:e12670.

Rao A, Bhat S, De S, Cyriac V, AR S. Study of [EMIM] [EtSO 4 ] ionic liquid-based gel polymer electrolyte mediated with hydroquinone redox additive for flexible supercapacitors. J Energy Storage 2023;68:107716.

Tiwari RK, Mishra R, Patel A, Tiwari A, Meghnani D, Singh RK. Polysulfide rejection strategy in lithium−sulfur batteries using an ion-conducting gel-polymer interlayer membrane. ACS Appl Mater Interfaces 2023;15:33957–71.

Qian Y, Yu Y, Wu W, Fan Q, Chai C, Hao J. Wide-temperature flexible supercapacitor from an organohydrogel electrolyte and its combined electrode. Chem Eur J 2023;29:e202300123.

Xu Y, Gao L, Wu X, Zhang S, Wang X, Gu C, et al. Porous composite gel polymer electrolyte with interfacial transport pathways for flexible quasi solid lithium-ion batteries. ACS Appl Mater Interfaces 2021;13:23743–50.

Castillo J, Santiago A, Judez X, Garbayo I, Clemente JAC, Morant-Miňana MC, et al. Safe, flexible, and high-performing gel-polymer electrolyte for rechargeable lithium metal batteries. Chem Mater 2021;33:8812–21.

Obeidat AM, Rastogi AC. Graphene and poly (3,4-ethylenedioxy-thiophene) (PEDOT) based hybrid supercapacitors with ionic liquid gel electrolyte in solid state design and their electrochemical performance in storage of solar photovoltaic generated electricity. MRS Adv 2016;1:3553–64.

Pandey GP, Liu T, Hancock C, Li Y, Sun XS, Li J. Thermostable gel polymer electrolyte based on succinonitrile and ionic liquid for high-performance solid-state supercapacitors. J Power Sources 2016;328:510–9.

Bhat Y, Yadav N, Hashmi SA. Pinecone-derived porous activated carbon for high performance all-solid-state electrical double layer capacitors fabricated with flexible gel polymer electrolytes. Electrochim Acta 2019;304:94–108.

Gu J, Wang H, Li S, Muhammad SR, Ning J, Pu X, et al. Tuning pyridinic-N and graphitic-N doping with 4,40-bipyridine in honeycomb-like porous carbon and distinct electrochemical roles in aqueous and ionic liquid gel electrolytes for symmetric supercapacitors. J Colloid Interface Sci 2023;635:254–64.

Ebadi M, Eriksson T, Mandal P, Costa LT, Araujo CM, Mindemark J, et al. Restricted ion transport by plasticizing side chains in polycarbonatebased solid electrolytes. Macromolecules 2020;53:764–74.

Pal P, Ghosh A. Solid-state gel polymer electrolytes based on ionic liquids containing imidazolium cations and tetrafluoroborate anions for electrochemical double layer capacitors: influence of cations size and viscosity of ionic liquids. J Power Sources 2018;406:128–40.

Shalu Chaurasia SK, Singh RK, Thermal Stability Chandra S, Behavior Complexing, Transport Ionic. Of polymeric gel membranes based on polymer PVdF-HFP and ionic liquid, [BMIM][BF4]. J Phys Chem B 2013;117:897–906.

Jeschke S, Mutke M, Jiang Z, Alt B, Wiemhöfer HD. Study of carbamate-modified disiloxane in porous PVDF-HFP membranes: new electrolytes/separators for lithiumIon batteries. ChemPhysChem 2014;15:1761–71.

Gan H, Li S, Zhang Y, Yu L, Wang J, Xue Z. Mechanically strong and electrochemically stable single-ion conducting polymer electrolytes constructed from hydrogen bonding. Langmuir 2021;37:8270–80.

Sha Y, Yu T, Dong T, Wu X, Ha Tao, Zhang H. In situ network electrolyte based on a functional polymerized ionic liquid with high conductivity toward lithium metal batteries. ACS Appl Energy Mater 2021;4:14755–65.

Yan T, Zou Y, Zhang X, Li D, Guo X, Yang D. Hydrogen bond interpenetrated Agarose/PVA network: a highly ionic conductive and flame-retardant gel polymer electrolyte. ACS Appl Mater Interfaces 2021;13:9856–64.

Guo P, Su A, Wei Y, Liu X, Li Y, Guo F, et al. Healable, highly conductive, flexible and nonflammable supramolecular ionogel electrolytes for lithium ion batteries. ACS Appl Mater Interfaces 2019;11:19413–20.

Bhat MY, Hashmi SA. Mixture of non-ionic and organic ionic plastic crystals immobilized in poly (vinylidene fluoride-co-hexafluoropropylene): a flexible gel polymer electrolyte composition for high performance carbon supercapacitors. J Energy Storage 2022;51:104514.

Li W, Li L, Liu Z, Zheng S, Li Q, Yan F. Supramolecular ionogels tougher than metals. Adv Mater 2023;35:2301383.

Qin Y, Song Z, Miao L, Hu C, Chen Y, Liu P, et al. Hydrogen-bond-mediated micelle aggregating self-assembly towards carbon nanofiber networks for high-energy and long-life zinc ion capacitors. Chem Eng J 2023;470:144256.

Zhang Y, Song Z, Miao L, Lv Y, Gan L, Liu M. All-round enhancement in zn-ion storage enabled by solvent-guided lewis acid–base self-assembly of heterodiatomic carbon nanotubes. ACS Appl Mater Interfaces 2023;15:35380–90.

Song Z, Miao L, Lv Y, Gan L, Liu M. NH4+ charge carrier coordinated h-bonded organic small molecule for fast and superstable rechargeable zinc batteries. Angew Chem, Int Ed 2023;62:e202309446.

France-Lanord A, Grossman JC. Correlations from ion pairing and the Nernst-Einstein equation. Phys Rev Lett 2019;122:136001.

Peng Y, Lai L, Tai Y, Zhang K, Xu X, Cheng X. Diffusion of ellipsoids in bacterial suspensions. Phys Rev Lett 2016;116:068303.

Van Der Spoel D, Lindahl E, Hess B, Groenhof G, Mark AE, Berendsen HJC. GROMACS: fast, flexible, and free. J Comput Chem 2005;26:1701–18.

Martínez L, Andrade R, Birgin EG, Martínez JM. PACKMOL: a package for building initial configurations for molecular dynamics simulations. J Comput Chem 2009;30:2157–64.

Rabet S, Raabe G. Comparison of the GAFF, OPLSAA and CHARMM27 force field for the reproduction of the thermodynamics properties of furfural, 2-methylfuran, 2,5-dimethylfuran and 5-hydroxymethylfurfural. Fluid Phase Equil 2022;554:113331.

Guo G, Ji C, Mi H, Yang C, Li M, Sun C, et al. Zincophilic anionic hydrogel electrolyte with interfacial specific adsorption of solvation structures for durable zinc ion hybrid supercapacitors. Adv Funct Mater 2023:2308405.

Lourenco TC, Ebadi M, Brandell D, Da Silva JLF, Costa LT. Interfacial structures in ionic liquid-based ternary electrolytes for lithium-metal batteries: a molecular dynamics study. J Phys Chem B 2020;124:9648–57.

Wang Y, Xie T, France-Lanord A, Berkley A, Johnson JA, Shao-Horn Y, et al. Toward designing highly conductive polymer electrolytes by machine learning assisted coarse-grained molecular dynamics. Chem Mater 2020;32:4144–51.

Chen W, Xing Z, Wei Y, Zhang X, Zhang Q. High thermal safety and conductivity gel polymer electrolyte composed of ionic liquid [EMIM][BF 4 ] and PVDF-HFP for EDLCs. Polymer 2023;268:125727.

Alexandre SA, Silva GG, Santamaría R, Trigueiro JPC, Lavall RL. A highly adhesive PIL/IL gel polymer electrolyte for use in flexible solid state supercapacitors. Electrochim Acta 2019;299:789–99.

Silva FCA, Ortega PFRdos Reis RA, Lavall RL, Costa LT. Polymer-ion interactions in PVDF@ionic liquid polymer electrolytes: a combined experimental and computational study. Electrochim Acta 2022;427:140831.

Hu P, Duan Y, Hu D, Qin B, Zhang J, Wang Q, et al. A rigid-flexible coupling high ionic conductivity polymer electrolyte for an enhanced performance of LiMn 2 O 4 /graphite battery at elevated temperature. ACS Appl Mater Interfaces 2015;7:4720–7.

Akbar ZA, Malik YT, Kim DH, Cho S, Jang SY, Jeon JW. Self-healable and stretchable ionic-liquid-based thermoelectric composites with high ionic seebeck coefficient. Small 2021;18:2106937.

Zhang Q, Feng S, Zhang X, Wei Y. Thermodynamic properties and intermolecular interactions of ionic liquids [DEME][BF 4 ] or [DEME][TFSI] and their binary mixture systems with GBL. J Mol Liq 2021;328:115373.

Hu Z, Chen J, Guo Y, Zhu J, Qu X, Niu W, et al. Fire-resistant, high-performance gel polymer electrolytes derived from poly(ionic liquid)/P(VDF-HFP) composite membranes for lithium ion batteries. J Membr Sci 2020;599:117827.

Lukatskaya MR, Dunn B, Gogotsi Y. Multidimensional materials and device architectures for future hybrid energy storage. Nat Commun 2016;7:12647–60.

Rana HH, Park JH, Gund GS, Park HS. Highly conducting, extremely durable, phosphorylated cellulose-based ionogels for renewable flexible supercapacitors. Energy Storage Mater 2020;25:70–5.

Yang K, Luo M, Zhang D, Liu C, Li Z, Wang L, et al. Ti 3 C 2 T x /carbon nanotube/porous carbon film for flexible supercapacitor. Chem Eng J 2022;427:132002.

Park JH, Rana HH, Kim JS, Hong JW, Lee SJ, Park HS. Inorganic−organic double network ionogels based on silica nanoparticles for high-temperature flexible supercapacitors. ACS Appl Mater Interfaces 2023;15:37344–53.

Feng M, Zhang Y, Zhu X, Chen W, Lu W, Wu G. Interface-Anchored covalent organic frameworks@amino modified Ti 3 C 2 T x MXene on Nylon 6 film for high-performance deformable supercapacitors. Angew Chem Int Ed 2023;62:e202307195.

Feng L, Wang K, Zhang X, Sun X, Li C, Ge X, et al. Flexible solid-state supercapacitors with enhanced performance from hierarchically graphene nanocomposite electrodes and ionic liquid incorporated gel polymer electrolyte. Adv Funct Mater 2017;28:1704463.

Bai Y, Yang C, Yuan B, Li H, Chen W, Yin H, et al. A UV cross-linked gel polymer electrolyte enabling high-rate and high voltage window for quasi-solid-state supercapacitors. J Energy Chem 2023;6:41–50.

Pandey GP, Hashmi SA. Ionic liquid 1-ethyl-3-methylimidazolium tetracyanoborate-based gel polymer electrolyte for electrochemical capacitors. J Mater Chem A 2013;1:3372–8.

Na R, Huo G, Zhang S, Huo P, Du Y, Luan J, et al. A novel poly(ethylene glycol)-grafted poly(arylene ether ketone) blend micro-porous polymer electrolyte for solid-state electric double layer capacitors formed by incorporating a chitosanm-based LiClO 4 gel electrolyte. J Mater Chem A 2016;4:18116–27.

Zhang Y, Li M, Qin B, Chen L, Liu Y, Zhang X, et al. Highly transparent, underwater self-healing, and ionic conductive elastomer based on multivalent ion-dipole interactions. Chem Mater 2020;32:6310–7.

Yadav N, Singh MK, Yadav N, Hashmi SA. High performance quasi-solid-state supercapacitors with peanut-shell-derived porous carbon. J Power Sources 2018;402:133–46.

Suleman M, Kumar Y, Hashmi SA. Solid-state electric double layer capacitors fabricated with plastic crystal based flexible gel polymer electrolytes: effective role of electrolyte anions. Mater Chem Phys 2015;163:161–71.

Mohit Hashmi SA. Biodegradable poly- ε -caprolactone based porous polymer electrolytes for high performance supercapacitors with carbon electrodes. J Power Sources 2023;557:232548.

Jia Z, Wang Y, Chen J, Cao Z, Pan S, Zhou Y, et al. Metal-organic frameworks derived low-crystalline NiCo 2 S 4 /Co 3 S 4 nanocages with dual heterogeneous interfaces for high-performance supercapacitors. Chin Chem Lett 2023;34:107137.

Wang X, Zhang Q, Zhao L, Hadi MK, Sambasivam S, Zhou Q, et al. A renewable hydrogel electrolyte membrane prepared by carboxylated chitosan and polyacrylamide for solid-state supercapacitors with wide working temperature range. J Power Sources 2023;560:232704.

Liu C, Wu H, Wang X, Fan J, Su H, Yang D, et al. Flexible solid-state supercapacitor integrated by methanesulfonic acid/polyvinyl acetate hydrogel and Ti 3 C 2 T x . Energy Storage Mater 2023;54:164–71.

Lashkari S, Valappil MO, Pala R, Pope MA. Electrolyte-mediated assembly of graphene-based supercapacitors using adsorbed ionic liquid/nonionic surfactant complexes. J Mater Chem A 2023;11:11222.

design and analysis of experiments in r

Web of Science

This is an open access article under the CC BY-NC-ND license ( http://creativecommons.org/licenses/by-nc-nd/4.0/ ).

design and analysis of experiments in r

京ICP备 10035462号-42

design and analysis of experiments in r

Development and analysis of Hastelloy-X alloy butt joint made by laser beam welding

  • Published: 18 September 2024
  • Volume 49 , article number  262 , ( 2024 )

Cite this article

design and analysis of experiments in r

  • G Sathishkumar 1 ,
  • S Senthil Murugan 2 &
  • P Sathiya 1  

The investigation describes the processing and analysis of laser welded Hastelloy-X (HX) alloy joints through this paper. HX alloys are employed in aerospace and nuclear industries, especially for high-temperature applications. In this research, a series of laser beam welding (LBW) experiments denoted as E1 to E9 were conducted on the HX base metal (BM), adhering to the L9-orthogonal array (L9-OA) design matrix. CO 2 laser technology was employed to fabricate HX butt joints. The laser power, focal length and welding speed were the variables. Then, the character of each joint was analyzed by sophisticated testing methods as per ASTM standards. The results showed that the character of each sample was varied depending on the selection of parameters. The E5 sample had a maximum tensile strength (TS) and ductility with 93 % joint efficiency. The grain elongation and refinement in the weld zone (WZ) were confirmed through microstructures and electron back scattered diffraction (EBSD) studies. The corrosion character of each joint (E1 to E9) was analyzed using the potentiostatic polarization method. The E1 sample had the highest corrosion resistance. The corrosion rate was in the range of 7.4E−03 to 8.6E−05 mm/yr. The dry sliding wear test (A1 to A9) was carried out as per L9-OA on the E1 sample, since this weld parameter had good corrosion resistance, by varying applied load, sliding distances and sliding velocities. Wear test results showed an increase in wear rate with the increase in load and sliding distance. A wear map was also drawn using the results to find out the association between wear rates and wear parameters. Weld speed had influenced the strength of the joint and laser power had shown an impact on the corrosion and wear of the HX alloy joints. The weld zone, corroded sample and the worn-out surfaces of weld joints were further analyzed using an optical microscope (OM) and Field Emission Scanning Electron Microscopy (FESEM).

This is a preview of subscription content, log in via an institution to check access.

Access this article

Subscribe and save.

  • Get 10 units per month
  • Download Article/Chapter or eBook
  • 1 Unit = 1 Article or 1 Chapter
  • Cancel anytime

Price includes VAT (Russian Federation)

Instant access to the full article PDF.

Rent this article via DeepDyve

Institutional subscriptions

design and analysis of experiments in r

Technical Data Sheet, 2005, Inconel Alloy HX. Special Metals. Special Metal Corporation, New Hartford, NY. No: SMC-007. pp.1–8

Sindo K 2003 Welding Metallurgy . 3rd edn. Wiley, Hoboken, pp. 1–688

Google Scholar  

DuPont J N, Lippold J C and Kiser D S 2009 Welding Metallurgy and Weldability of Nickel-Base Alloys . 1st edn. Wiley, Hoboken, pp. 1–440. ISBN: 9780470087145. 1-440. https://doi.org/10.1002/9780470500262

Book   Google Scholar  

Radhakrishnan B and Thompson R G 1993 The effect of weld heat-affected zone (HAZ) liquation kinetics on the hot cracking susceptibility of alloy 718. Metall. Mater. Trans. A. 24(6): 1409–1422. https://doi.org/10.1007/BF02668209

Article   Google Scholar  

Ren W, Lu F, Yang R, Liu X, Li Z and Hosseini S R E 2015 A comparative study on fiber laser and CO 2 laser welding of inconel 617. Mater. Design. 76: 207–214. https://doi.org/10.1016/j.matdes.2015.03.033

Manikandan M, Arivazhagan N, Nageswara Rao M and Reddy G M 2014 Microstructure and mechanical properties of alloy C-276 weldments fabricated by continuous and pulsed current gas tungsten arc welding techniques. J. Manuf. Process. 16(4): 563–572. https://doi.org/10.1016/j.jmapro.2014.08.002

Zhang X, Xu H, Li Z, Dong A, Du D and Lei L et al . 2021 Effect of the scanning strategy on microstructure and mechanical anisotropy of Hastelloy X superalloy produced by laser powder bed fusion. Mater. Charact. 173: 110951. https://doi.org/10.1016/j.matchar.2021.110951

Pakniat M, Ghaini F M and Torkamany M J 2016 Hot Cracking in laser welding of Hastelloy X with pulsed Nd: YAG and continuous fiber lasers. Mater. Design. 106: 177–183. https://doi.org/10.1016/j.matdes.2016.05.124

Manikandan S G K, Sivakumar D, Rao K P and Kamaraj M 2014 Effect of weld cooling rate on laves phase formation in inconel 718 fusion zone. J. Mater. Process. Technol. 214(2): 358–364. https://doi.org/10.1016/j.jmatprotec.2013.09.006

Reza Abedi M, Sabet H and Razavi H 2016 The effect of repair welding number on microstructure of Hastelloy X fabricated via TIG process. Int. J. Mater. Sci. Appl. 5(2): 43–48. https://doi.org/10.11648/j.ijmsa.20160502.12

Shimizu S and Mutoh Y 1984 Weldability and weld performance of a special grade Hastelloy X modified for high temperature gas-cooled reactors. Nucl. Technol. 66(1): 44–53. https://doi.org/10.13182/NT84-A33453

Lertora E, Mandolfino C and Gambaro C 2014 Mechanical behaviour of inconel 718 thin walled laser welded components for aircraft engines. Int. J. Aerosp. Eng. 2014: 721680. https://doi.org/10.1155/2014/721680

Thirugnanasambantham K G and Natarajan S 2016 Mechanistic studies on degradation in sliding wear behavior of IN718 and Hastelloy X superalloys at 500°C. Tribol. Int. 101: 324–330. https://doi.org/10.1016/j.triboint.2016.04.016

Khan Adam M, Sundarrajan S, Duraiselvam M, Natarajan S and Kumar Senthil A 2017 Sliding wear behaviour of plasma sprayed coatings on nickel based superalloy. Surf. Eng. 33(1): 35–41. https://doi.org/10.1179/1743294415Y.0000000087

Stott F H 1998 The role of oxidation in the wear of alloys. Tribol. Int. 31(1–3): 61–71. https://doi.org/10.1016/S0301-679X(98)00008-5

Stott F H and Jordan M P 2001 The effects of load and substrate hardness on the development and maintenance of wear-protective layers during sliding at elevated temperatures. Wear. 250(1–12): 391–400. https://doi.org/10.1016/S0043-1648(01)00601-9

Zhen J, Cheng J, Zhu S, Hao J, Qiao Z and Yang J et al . 2017 High-temperature tribological behavior of a nickel alloy matrix solid-lubricating composite under vacuum. Tribol. Int. 110: 52–56. https://doi.org/10.1016/j.triboint.2017.02.011

Zimmerman C 2013, Boriding (Boronizing) of metals, steel heat treating fundamental and processes. In: Dossett J L and George E (eds.) ASM Handbook, Totten, ASM International. ISBN: 978-1-62708-165-8. pp. 709–724. https://doi.org/10.31399/asm.hb.v04a.a0005772

Ekambaram P 2019 Study of mechanical and metallurgical properties of Hastelloy X at cryogenic condition. J. Mater. Res. Technol. 8(6): 6413–6419. https://doi.org/10.1016/j.jmrt.2019.10.048

Beamer C, Denlinger D, Rao S and Christina D 2021, High Pressure heat treatment for L-PBF Hastelloy X. In: Proceedings of the HT2021. Heat Treat 2021: Proceedings from the 31st Heat Treating Society Conference and Exposition . St. Louis, Missouri, pp. 44–50. https://doi.org/10.31399/asm.cp.ht2021p0044

Lee Y S and Sung J H 2023 Microstructure and mechanical properties of Hastelloy X fabricated using directed energy deposition. Metals. 13(5): 885. https://doi.org/10.3390/met13050885

Datta S, Raza Mohammad S, Saha P and Pratihar D K 2019 Effects of process parameters on the quality aspects of weld-bead in laser welding of NiTinol sheets. Mater. Manuf. Process. 34(6): 648–659. https://doi.org/10.1080/10426914.2019.1566608

Download references

No funding was obtained for this study.

Author information

Authors and affiliations.

Department of Production Engineering, National Institute of Technology, Trichy, Tamil Nadu, 620 015, India

G Sathishkumar & P Sathiya

Department of Mechanical Engineering, Rajalakshmi Engineering College, Chennai, Tamil Nadu, 602 105, India

S Senthil Murugan

You can also search for this author in PubMed   Google Scholar

Corresponding author

Correspondence to G Sathishkumar .

Ethics declarations

Conflict of interest.

The authors declare that there they no conflict of interest in publishing this article.

Rights and permissions

Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.

Reprints and permissions

About this article

Sathishkumar, G., Senthil Murugan, S. & Sathiya, P. Development and analysis of Hastelloy-X alloy butt joint made by laser beam welding. Sādhanā 49 , 262 (2024). https://doi.org/10.1007/s12046-024-02603-y

Download citation

Received : 22 February 2024

Revised : 13 May 2024

Accepted : 27 June 2024

Published : 18 September 2024

DOI : https://doi.org/10.1007/s12046-024-02603-y

Share this article

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

  • Hastelloy-X
  • CO 2 laser welding
  • potentiostatic polarization
  • Find a journal
  • Publish with us
  • Track your research

Information

  • Author Services

Initiatives

You are accessing a machine-readable page. In order to be human-readable, please install an RSS reader.

All articles published by MDPI are made immediately available worldwide under an open access license. No special permission is required to reuse all or part of the article published by MDPI, including figures and tables. For articles published under an open access Creative Common CC BY license, any part of the article may be reused without permission provided that the original article is clearly cited. For more information, please refer to https://www.mdpi.com/openaccess .

Feature papers represent the most advanced research with significant potential for high impact in the field. A Feature Paper should be a substantial original Article that involves several techniques or approaches, provides an outlook for future research directions and describes possible research applications.

Feature papers are submitted upon individual invitation or recommendation by the scientific editors and must receive positive feedback from the reviewers.

Editor’s Choice articles are based on recommendations by the scientific editors of MDPI journals from around the world. Editors select a small number of articles recently published in the journal that they believe will be particularly interesting to readers, or important in the respective research area. The aim is to provide a snapshot of some of the most exciting work published in the various research areas of the journal.

Original Submission Date Received: .

  • Active Journals
  • Find a Journal
  • Journal Proposal
  • Proceedings Series
  • For Authors
  • For Reviewers
  • For Editors
  • For Librarians
  • For Publishers
  • For Societies
  • For Conference Organizers
  • Open Access Policy
  • Institutional Open Access Program
  • Special Issues Guidelines
  • Editorial Process
  • Research and Publication Ethics
  • Article Processing Charges
  • Testimonials
  • Preprints.org
  • SciProfiles
  • Encyclopedia

sensors-logo

Article Menu

design and analysis of experiments in r

  • Subscribe SciFeed
  • Recommended Articles
  • Google Scholar
  • on Google Scholar
  • Table of Contents

Find support for a specific problem in the support section of our website.

Please let us know what you think of our products and services.

Visit our dedicated information section to learn more about MDPI.

JSmol Viewer

Bioelectrical impedance vector analysis (biva) for assessment of hydration status: a comparison between endurance and strength university athletes.

design and analysis of experiments in r

1. Introduction

2. materials and methods, 2.1. study design and sampling, 2.1.1. athletic population, 2.1.2. reference population, 2.1.3. selection criteria, 2.2. recruitment, 2.3. data collection, 2.3.1. pre-hydration assessment questionnaire, 2.3.2. anthropometric measurements and body composition, 2.3.3. urine collection, 2.3.4. sweat rate, 2.4. ethical considerations, 2.5. statistical analysis, 3.1. urine color and biva parameters, 3.2. agreement between biva and (usg and sr), 3.2.1. bland–altman plots.

Click here to enlarge figure

3.2.2. Canonical Correlation Analysis

3.2.3. examination of data characteristics and comparative analysis, 3.3. pre–post measurements, 4. discussion, 4.1. strengths, 4.2. limitations, 4.3. future studies, 5. conclusions, supplementary materials, author contributions, institutional review board statement, informed consent statement, data availability statement, acknowledgments, conflicts of interest, appendix a. urine color chart, appendix b. urine specific gravity equivalent hydration status, appendix c. eligibility questionnaire.

DISEASEPlease tick (✓) as appropriateDuration of disease (years)Do you still suffer from any of these conditions? (YES/NO)Are you following any treatment? (YES/NO)
Diabetes
Heart Disease
Hypertension
High Blood Cholesterol
Thyroid Disease
Sleep Apnea
Other (Asthma, Claustrophobia)
Cancer
Other Health Problems

Appendix D. Pre-Hydration Assessment Questionnaire

  • Other, please specify: _____________
  • How much time do you usually spend doing vigorous physical activities? _____ hours per day       _____ minutes per day _____Don’t know/Not sure
  • Medical/Biochemical/Clinical Assessment:
  • How often do you have a bowel movement? _____________________________
  • If you take laxatives, what type/brand and how often? _____________________
  • Heartburn, Bloating, Gas, Constipation, Diarrhea, Nausea, Vomiting:
  • When was the last time you had Blood Test? ____________
  • Recent Blood test results:
  • Water Retention: Yes No
  • Dietary Assessment:
  • If yes, what is the reason?
  • Are you taking any drugs or medications? If yes, specify name, quantity and duration:______________________________________________________ ▢ Yes, please specify: _______________ ▢ No ▢ Don’t know
  • Do you consume any of the following foods? ▢ Beets ▢ Blackberries ▢ Carrots ▢ Fava beans ▢ Rhubarb
  • Which meals do you eat regularly, check all that apply: ▢ Breakfast ▢ AM Snack ▢ Lunch ▢ Afternoon Snack ▢ Dinner ▢ Late Night Snack
  • Are you vegetarian? ▢ Yes ▢ No
  • Have you ever followed a dietary regime? ▢ Yes, now ▢ Yes, previously ▢ No ▢ If yes, please give details: ____________________________________________________________________________________________________________________
Supplement Type and BrandHow much? When? (Before or after Training)Reason for Consumption
  • How many fruits and vegetables in all do you consume daily on average? ____________ Fruits _________Vegetables
  • How often do you consume fast foods? ▢ Never ▢ Two to three times monthly ▢ One or two times weekly ▢ Daily
  • How many liters of water do you consume daily? ▢ Less than one liter ▢ One liter to two liters ▢ Two to three liters ▢ More than three liters
Tea  None to occasionally
Coffee  One to four
Cola/Pepsi  Four to ten
Redbull  More than ten
  • Do you consume alcohol? If yes, specify ________________ ▢ Never or rarely ▢ Only in weekends ▢ Once or twice weekly ▢ Once or twice daily ▢ More than twice daily
  • Do you smoke? ▢ Yes ▢ No If yes, how much on average? __________________________
  • On average how many hours do you sleep per night? ▢ More than 10 h ▢ Between 8–10 h ▢ Between 6–8 h ▢ Between 4–6 h ▢ Less than 4 h
  • How important do you think good nutrition is to sports performance? ▢ Very important ▢ Important ▢ Moderately important ▢ Of little importance ▢ Unimportant
  • How important do you think hydration status is to sports performance? ▢ Very important ▢ Important ▢ Moderately important ▢ Of little importance ▢ Unimportant
  • “The more supplements I take, the better I will perform”: ▢ Agree ▢ Disagree ▢ Strongly disagree ▢ Neither agree nor disagree
  • Have you received any previous nutritional advice? ▢ Yes ▢ No
  • Do you have access to a sports nutritionist/dietitian? ▢ Yes, through CHDC ▢ Yes, outside AUB ▢ No Anthropometric Assessment: Usual body weight: _________
  • During these past 6 months, did your weight change? ▢ Yes, I gained weight ▢ Yes, I lost weight ▢ No ▢ Not sure How many kg(s)_____________________?
  • Would you like to receive the hydration tests and body composition results by phone/mail? ▢ Yes, by phone ▢ Yes, by mail ▢ No, I don’t want to receive the results Instructions to prepare you for the hydration tests: - Avoid caffeine and alcohol consumption prior to the experiment. - Do not exercise intensely at least 12 h prior to the experiment. - Avoid food in the 4 h prior to the experiment.

Appendix E. Data Collection Sheet

Gender   
Date of birth/Age   
Sport type   
Eligibility questionnaire   
Pre-hydration assessment questionnaire   
Empty urine cup 1 weight   
Empty urine cup 2 weight   
Empty urine cup 3 weight   
Urine collection 1   
Urine collection 2   
Urine collection 3   
Height   
Body weight 1   
Body weight 2   
Water bottle weight A   
Water bottle weight B   
BIA test 1   
BIA test 2   
BIA test 3   
BIA test 4   
  • Fink, H.H.; Mikesky, A.E. Practical Applications in Sports Nutrition ; Jones & Bartlett Learning: Burlington, MA, USA, 2015. [ Google Scholar ]
  • Meyer, F.; Volterman, K.A.; Timmons, B.W.; Wilk, B. Fluid Balance and Dehydration in the Young Athlete: Assessment Considerations and Effects on Health and Performance. Am. J. Lifestyle Med. 2012 , 6 , 489–501. [ Google Scholar ] [ CrossRef ]
  • Chumlea, W.C.; Guo, S.S.; Zeller, C.M.; Reo, N.V.; Baumgartner, R.N.; Garry, P.J.; Wang, J.; Pierson, R.N., Jr.; Heymsfield, S.B.; Siervogel, R.M. Total body water reference values and prediction equations for adults. Kidney Int. 2001 , 59 , 2250–2258. [ Google Scholar ] [ CrossRef ] [ PubMed ]
  • Barley, O.R.; Chapman, D.W.; Abbiss, C.R. Reviewing the current methods of assessing hydration in athletes. J. Int. Soc. Sports Nutr. 2020 , 17 , 52. [ Google Scholar ] [ CrossRef ]
  • Campa, F.; Toselli, S.; Mazzilli, M.; Gobbo, L.A.; Coratella, G. Assessment of body composition in athletes: A narrative review of available methods with special reference to quantitative and qualitative bioimpedance analysis. Nutrients 2021 , 13 , 1620. [ Google Scholar ] [ CrossRef ]
  • Casa, D.J. Performing in extreme environments. J. Athl. Train. 2000 , 35 , 111. [ Google Scholar ]
  • Zubac, D.; Marusic, U.; Karninčič, H. Hydration Status Assessment Techniques and Their Applicability among Olympic Combat Sports Athletes: Literature Review. Strength Cond. J. 2016 , 38 , 80–89. [ Google Scholar ] [ CrossRef ]
  • Achamrah, N.; Colange, G.; Delay, J.; Rimbert, A.; Folope, V.; Petit, A.; Grigioni, S.; Déchelotte, P.; Coëffier, M. Comparison of body composition assessment by DXA and BIA according to the body mass index: A retrospective study on 3655 measures. PLoS ONE 2018 , 13 , e0200465. [ Google Scholar ] [ CrossRef ] [ PubMed ]
  • Kyle, U.G.; Bosaeus, I.; De Lorenzo, A.D.; Deurenberg, P.; Elia, M.; Gómez, J.M.; Heitmann, B.L.; Kent-Smith, L.; Melchior, J.-C.; Pirlich, M. Bioelectrical impedance analysis—Part I: Review of principles and methods. Clin. Nutr. 2004 , 23 , 1226–1243. [ Google Scholar ] [ CrossRef ]
  • Nwosu, A.C.; Mayland, C.R.; Mason, S.; Cox, T.F.; Varro, A.; Stanley, S.; Ellershaw, J. Bioelectrical impedance vector analysis (BIVA) as a method to compare body composition differences according to cancer stage and type. Clin. Nutr. ESPEN 2019 , 30 , 59–66. [ Google Scholar ] [ CrossRef ]
  • Tinsley, G.M.; Moore, M.L.; Silva, A.M.; Sardinha, L.B. Cross-sectional and longitudinal agreement between two multifrequency bioimpedance devices for resistance, reactance, and phase angle values. Eur. J. Clin. Nutr. 2020 , 74 , 900–911. [ Google Scholar ] [ CrossRef ]
  • Stagi, S.; Silva, A.M.; Jesus, F.; Campa, F.; Cabras, S.; Earthman, C.P.; Marini, E. Usability of classic and specific bioelectrical impedance vector analysis in measuring body composition of children. Clin. Nutr. (Edinb. Scotl.) 2022 , 41 , 673–679. [ Google Scholar ] [ CrossRef ] [ PubMed ]
  • De la Cruz Marcos, S.; Redondo Del Río, M.P.; de Mateo Silleras, B. Applications of bioelectrical impedance vector analysis (Biva) in the study of body composition in athletes. Appl. Sci. 2021 , 11 , 9781. [ Google Scholar ] [ CrossRef ]
  • Castizo-Olier, J.; Irurtia, A.; Jemni, M.; Carrasco-Marginet, M.; Fernández-García, R.; Rodríguez, F.A. Bioelectrical impedance vector analysis (BIVA) in sport and exercise: Systematic review and future perspectives. PLoS ONE 2018 , 13 , e0197957. [ Google Scholar ] [ CrossRef ]
  • Di Vincenzo, O.; Marra, M.; Scalfi, L. Bioelectrical impedance phase angle in sport: A systematic review. J. Int. Soc. Sports Nutr. 2019 , 16 , 1–49. [ Google Scholar ] [ CrossRef ]
  • Maughan, R.J.; Shirreffs, S.M. Development of hydration strategies to optimize performance for athletes in high-intensity sports and in sports with repeated intense efforts. Scand. J. Med. Sci. Sports 2010 , 20 , 59–69. [ Google Scholar ] [ CrossRef ]
  • Suppiah, H.T.; Ng, E.L.; Wee, J.; Taim, B.C.; Huynh, M.; Gastin, P.B.; Chia, M.; Low, C.Y.; Lee, J.K.W. Hydration status and fluid replacement strategies of high-performance adolescent athletes: An application of machine learning to distinguish hydration characteristics. Nutrients 2021 , 13 , 4073. [ Google Scholar ] [ CrossRef ] [ PubMed ]
  • Martins, P.C.; Junior, C.A.S.A.; Silva, A.M.; Silva, D.A.S. Phase angle and body composition: A scoping review. Clin. Nutr. ESPEN 2023 , 56 , 237–250. [ Google Scholar ] [ CrossRef ]
  • Martins, P.C.; Gobbo, L.A.; Silva, D.A.S. Bioelectrical impedance vector analysis (BIVA) in university athletes. J. Int. Soc. Sports Nutr. 2021 , 18 , 7. [ Google Scholar ] [ CrossRef ]
  • Carrasco-Marginet, M.; Castizo-Olier, J.; Rodríguez-Zamora, L.; Iglesias, X.; Rodríguez, F.A.; Chaverri, D.; Brotons, D.; Irurtia, A. Bioelectrical impedance vector analysis (BIVA) for measuring the hydration status in young elite synchronized swimmers. PLoS ONE 2017 , 12 , e0178819. [ Google Scholar ] [ CrossRef ]
  • Thorpe, B.R. Validation of Urinary Biomarkers of Hydration Status in College Athletes. Master’s Thesis, Virginia Tech, Blacksburg, VA, USA, 2018. [ Google Scholar ]
  • Richardson, A.J. The Physiological and Renal Responses to Hydration Status in Hypoxia. Ph.D. Thesis, Brighton University, Brighton and HoveBrighton, UK, 2010. [ Google Scholar ]
  • Baker, L.B. Sweating rate and sweat sodium concentration in athletes: A review of methodology and intra/interindividual variability. Sports Med. 2017 , 47 , 111–128. [ Google Scholar ] [ CrossRef ]
  • Piccoli, A.; Pastori, G. BIVA Software ; Department of Medical and Surgical Sciences, University of Padova: Padova, Italy, 2002. [ Google Scholar ]
  • Evans, W.D.; McClagish, H.; Trudgett, C. Factors affecting the in vivo precision of bioelectrical impedance analysis. Appl. Radiat. Isot. 1998 , 49 , 485–487. [ Google Scholar ] [ CrossRef ] [ PubMed ]
  • Micheli, M.L.; Pagani, L.; Marella, M.; Gulisano, M.; Piccoli, A.; Angelini, F.; Burtscher, M.; Gatterer, H. Bioimpedance and impedance vector patterns as predictors of league level in male soccer players. Int. J. Sports Physiol. Perform. 2014 , 9 , 532–539. [ Google Scholar ] [ CrossRef ] [ PubMed ]
  • Caton, J.R.; Molé, P.A.; Adams, W.C.; Heustis, D.S. Body composition analysis by bioelectrical impedance: Effect of skin temperature. Med. Sci. Sports Exerc. 1988 , 20 , 489–491. [ Google Scholar ] [ CrossRef ] [ PubMed ]
  • Silleras, B.d.M.; Ares, G.C.; Marcos, S.d.l.C.; Enciso, L.C.; Fernández, E.Q.; Río, P.R.d. Bioelectrical Impedance Vector Analysis (BIVA) and Somatotype in Female Rugby Players. Appl. Sci. 2023 , 13 , 5242. [ Google Scholar ] [ CrossRef ]
  • Koury, J.C.; Trugo, N.M.F.; Torres, A.G. Phase angle and bioelectrical impedance vectors in adolescent and adult male athletes. Int. J. Sports Physiol. Perform. 2014 , 9 , 798–804. [ Google Scholar ] [ CrossRef ]
  • Gatterer, H.; Schenk, K.; Laninschegg, L.; Schlemmer, P.; Lukaski, H.; Burtscher, M. Bioimpedance identifies body fluid loss after exercise in the heat: A pilot study with body cooling. PLoS ONE 2014 , 9 , e109729. [ Google Scholar ] [ CrossRef ]
  • Andreoli, A.; Monteleone, M.; Van Loan, M.; Promenzio, L.; Tarantino, U.; De Lorenzo, A. Effects of different sports on bone density and muscle mass in highly trained athletes. Med. Sci. Sports Exerc. 2001 , 33 , 507–511. [ Google Scholar ] [ CrossRef ]
  • Tinsley, G.M.; Stratton, M.T.; Harty, P.S.; Williams, A.D.; White, S.J.; Rodriguez, C.; Dellinger, J.R.; Johnson, B.A.; Smith, R.W.; Trexler, E.T. Influence of acute water ingestion and prolonged standing on raw bioimpedance and subsequent body fluid and composition estimates. J. Electr. Bioimpedance 2022 , 13 , 10–20. [ Google Scholar ] [ CrossRef ]
  • Cheuvront, S.N.; Kenefick, R.W.; Montain, S.J.; Sawka, M.N. Mechanisms of aerobic performance impairment with heat stress and dehydration. J. Appl. Physiol. 2010 , 109 , 1989–1995. [ Google Scholar ] [ CrossRef ]
  • Sagayama, H.; Yamada, Y.; Ichikawa, M.; Kondo, E.; Yasukata, J.; Tanabe, Y.; Higaki, Y.; Takahashi, H. Evaluation of fat-free mass hydration in athletes and non-athletes. Eur. J. Appl. Physiol. 2020 , 120 , 1179–1188. [ Google Scholar ] [ CrossRef ]
  • Azmy, U.; Rahmaniah, N.; Renzytha, A.R.; Atmaka, D.R.; Pratiwi, R.; Rizal, M.; Adiningsih, S.; Herawati, L. Comparison of Body Compositions among Endurance, Strength, and Team Sports Athletes. Sport Mont. 2023 , 21 , 45–50. [ Google Scholar ] [ CrossRef ]
  • Heavens, K.R.; Charkoudian, N.; O’Brien, C.; Kenefick, R.W.; Cheuvront, S.N. Noninvasive assessment of extracellular and intracellular dehydration in healthy humans using the resistance-reactance-score graph method. Am. J. Clin. Nutr. 2016 , 103 , 724–729. [ Google Scholar ] [ CrossRef ]
  • Sawka, M.N. Physiological consequences of hypohydration: Exercise performance and thermoregulation. Med. Sci. Sports Exerc. 1992 , 24 , 657–670. [ Google Scholar ] [ CrossRef ] [ PubMed ]
  • Brener, A.; Waksman, Y.; Rosenfeld, T.; Levy, S.; Peleg, I.; Raviv, A.; Interator, H.; Lebenthal, Y. The heritability of body composition. BMC Pediatr. 2021 , 21 , 225. [ Google Scholar ] [ CrossRef ]
  • Segal, K.R.; Burastero, S.; Chun, A.; Coronel, P.; Pierson, R.N.; Wang, J. Estimation of extracellular and total body water by multiple-frequency bioelectrical-impedance measurement. Am. J. Clin. Nutr. 1991 , 54 , 26–29. [ Google Scholar ] [ CrossRef ] [ PubMed ]
  • De Lorenzo, A.; Andreoli, A.; Matthie, J.; Withers, P. Predicting body cell mass with bioimpedance by using theoretical methods: A technological review. J. Appl. Physiol. 1997 , 82 , 1542–1558. [ Google Scholar ] [ CrossRef ]
  • Castizo-Olier, J.; Carrasco-Marginet, M.; Roy, A.; Chaverri, D.; Iglesias, X.; Pérez-Chirinos, C.; Rodríguez, F.; Irurtia, A. Bioelectrical impedance vector analysis (BIVA) and body mass changes in an ultra-endurance triathlon event. J. Sports Sci. Med. 2018 , 17 , 571–579. [ Google Scholar ]
  • Shirreffs, S.M.; Aragon-Vargas, L.F.; Chamorro, M.; Maughan, R.J.; Serratosa, L.; Zachwieja, J.J. The Sweating Response of Elite Professional Soccer Players to Training in the Heat. Int. J. Sports Med. 2005 , 26 , 90–95. [ Google Scholar ] [ CrossRef ]
  • Campa, F.; Toselli, S. Bioimpedance vector analysis of elite, subelite, and low-level male volleyball players. Int. J. Sports Physiol. Perform. 2018 , 13 , 1250–1253. [ Google Scholar ] [ CrossRef ]
  • Giorgi, A.; Vicini, M.; Pollastri, L.; Lombardi, E.; Magni, E.; Andreazzoli, A.; Orsini, M.; Bonifazi, M.; Lukaski, H.; Gatterer, H. Bioimpedance patterns and bioelectrical impedance vector analysis (BIVA) of road cyclists. J. Sports Sci. 2018 , 36 , 2608–2613. [ Google Scholar ] [ CrossRef ]
  • Cavedon, V.; Milanese, C.; Marchi, A.; Zancanaro, C. Different amount of training affects body composition and performance in high-intensity functional training participants. PLoS ONE 2020 , 15 , e0237887. [ Google Scholar ] [ CrossRef ] [ PubMed ]
  • Alvero-Cruz, J.R.; Carnero, E.A.; García, M.A.G.; Cárceles, F.A.; Correas-Gómez, L.; Rosemann, T.; Nikolaidis, P.T.; Knechtle, B. Predictive performance models in long-distance runners: A narrative review. Int. J. Environ. Res. Public Health 2020 , 17 , 8289. [ Google Scholar ] [ CrossRef ] [ PubMed ]
Correlations
Urine Color PREAverage Xc PREAverage R PREAverage Z PRE
Spearman’s rhoUrine color PRECorrelation Coefficient1.0000.0120.0580.057
Sig. (2-tailed).0.8880.4940.498
N142142142142
Average Xc PRECorrelation Coefficient0.0121.0000.746 **0.749 **
Sig. (2-tailed)0.888.0.0000.000
N142148148148
Average R PRECorrelation Coefficient0.0580.746 **1.0001.000 **
Sig. (2-tailed)0.4940.000.0.000
N142148148148
Average Z PRECorrelation Coefficient0.0570.749 **1.000 **1.000
Sig. (2-tailed)0.4980.0000.000.
N142148148148
( ) Canonical correlation analysis coefficients
Canonical FunctionBIVA VariablesCoefficientOther methodsCoefficient
Function 1PA0.7Urine Color0.4
Xc0.4Sweat Rate−2.3
R5.8USG−36.3
Z−5.8Body Mass Loss0.4
Function 2PA7.1Urine Color0.3
Xc−0.5Sweat Rate1.4
R6.7USG62.3
Z−6.6Body Mass Loss−1.8
Function 3PA−6.6Urine Color0.4
Xc−2.0Sweat Rate0.3
R−27.1USG21.0
Z27.2Body Mass Loss1.1
Function 4PA1.1Urine Color−0.2
Xc−0.9Sweat Rate−1.9
R−8.7USG83.0
Z8.7Body Mass Loss0.8
( ) Correlation table for BIVA observed and BIVA dimensions
BIVA VariablesDimension 1Dimension 2Dimension 3Dimension 4
PA−0.9−0.1−0.30.3
XC0.5−0.4−0.20.7
R0.9−0.10.10.3
Z0.9−0.10.10.3
( ) Correlation table for other methods observed and BIVA dimensions (cross-loadings).
Other MethodsDimension 1Dimension 2Dimension 3Dimension 4
Urine Color0.20.10.1−0.0
Sweat Rate−0.3−0.00.1−0.0
USG−0.10.10.00.0
Body Mass Loss−0.1−0.10.10.0
( ) Correlation table for BIVA observed and other methods dimensions (cross-loadings).
BIVA VariablesDimension 1Dimension 2Dimension 3Dimension 4
PA−0.4−0.0−0.00.0
XC0.2−0.1−0.00.0
R0.4−0.00.00.0
Z0.4−0.00.00.0
( ) Correlation table for other methods and other methods dimensions.
Other MethodsDimension 1Dimension 2Dimension 3Dimension 4
Urine Color0.60.40.7−0.3
Sweat Rate−0.8−0.10.5−0.4
USG−0.30.50.30.8
Body Mass Loss−0.4−0.60.70.1
( ) Statistical test results for canonical dimensions.
Canonical DimensionsHotelling–Lawley TraceF-Approximation StatisticApproximate df1Approximate df2p-Value
1 to 40.24702.0766165380.0082
2 to 40.04210.638195460.7647
3 to 40.01660.574545540.6812
4 to 40.00050.065515620.7981
Statistics
GenderSportUSG_MUSG_PREUSG_POST
Male1NValid545353
Missing011
Mean1.0208651.0171591.016575
Median1.0211001.0168001.016200
Std. Deviation0.00788260.00828620.0091917
2NValid363436
Missing020
Mean1.0194691.0165961.017628
Median1.0199001.0155001.016300
Std. Deviation0.00769170.00927460.0081198
Female1NValid292729
Missing131
Mean1.0174661.0148871.014226
Median1.0174001.0150001.011600
Std. Deviation0.00889950.00856370.0097921
2NValid282828
Missing000
Mean1.0193001.0169481.016348
Median1.0183501.0188001.018350
Std. Deviation0.00684510.00894060.0091557
Multivariate Tests
EffectValueFHypothesis dfError dfSig.Partial Eta Squared
USGPillai’s Trace0.1292.371 2.00032.0000.1100.129
Wilks’ Lambda0.8712.371 2.00032.0000.1100.129
Hotelling’s Trace0.1482.371 2.00032.0000.1100.129
Roy’s Largest Root0.1482.371 2.00032.0000.1100.129
Delta ValueDelta Value/hPaired t-Testp-ValueCohen’s d
Mean ± Error Term
Gender = 1 (Male)
R/h−4.50 ± 1.20 −1.32 ± 0.347.381.1117 × 10 ***0.799
Xc/h−0.37 ± 0.20−1.00 ± 0.553.530.000549 ***0.46
Z−4.51 ± 1.21−1.32 ± 0.347.331.4247 × 10 ***0.7957
PA0.02 ± 0.02 0.32 ± 0.35−1.820.07 *0.116
Gender = 2 (Female)
R/h−4.84 ± 1.27−1.57 ± 0.427.543.995 × 10 ***0.437
Xc/h−0.49 ± 0.22−1.36 ± 0.624.343.794 × 10 ***0.118
Z−4.85 ± 1.28−1.57 ± 0.427.514.665 × 10 ***0.433
PA0.01 ± 0.030.23 ± 0.42−1.040.270.202
Sport = 1 (Weightlifting)
R/h−3.99 ± 2.40−0.94 ± 0.573.330.001527 ***0.791
Xc/h−0.17 ± 0.39−0.45 ± 1.010.90.370.616
Z−3.98 ± 2.42−0.94 ± 0.573.30.001677 ***0.792
PA0.03 ± 0.030.46 ± 0.63−1.540.130.179
Sport = 2 (Endurance)
R/h−4.64 ± 1.27−1.40 ± 0.387.251.956 × 10 ***0.475
Xc/h−0.62 ± 0.22−1.70 ± 0.605.652.202 × 10 ***0.0212
Z−4.67 ± 1.28−1.39 ± 0.397.261.896 × 10 ***0.468
PA−0.02 ± 0.02−0.32 ± 0.401.640.110.541
Xc/H + R/H+ Z/H + PA (PRE)Two-Sample Hotelling’s T TestPaired One-Sample Hotelling’s T-Test
All athletes and reference populationT = 18.2p-value = 8.79 × 10 ***
Male athletes and male reference populationT = 18.2p-value = 0.001775 ***
Female athletes and female reference populationT = 4.8p-value = 0.3208
Strength athletes and general athletic populationT = 24.1p-value = 0.000179 ***
Endurance athletes and general athletic populationT = 11.5p-value = 0.02721 ***
Female endurance athletes and female strength athletes pre–post training Endurance females: T = 17.3; p = 0.77
Strength females: T = 11.8; p = 0.64
Male endurance athletes and male strength athletes pre–post training Endurance males: T = 22.7; p = 0.79
Strength males: T = 43.7; p = 0.9
Welch’s Two-Sample T-Test
GroupsZ/HR/HXc/HPA
Endurance females vs. reference femalesT = 94.707
df = 100.05
p < 2.2 × 10 ***
T = 1.1654
df = 72.47
p = 0.2477
T = −0.6267
df = 65.93
p = 0.533
T = −1.8205
df = 42.21
p = 0.0758
95% CI:
[417.46, 435.32]
95% CI:
[−7.38, 28.18]
95% CI:
[−2.40, 1.25]
95% CI:
[−0.45, 0.02]
: 430.55
: 4.16
: 424.95
: 414.55
: 39.34
: 39.91
: 5.31
: 5.52
Strength females vs. reference females T = 94.642
df = 100.11
p < 2.2 × 10 ***
T = −1.0636
df = 46.66
p = 0.293
T = −1.9459
df = 55.08
p = 0.0568
T = −1.0687
df = 43.61
p = 0.2911
95% CI:
[417.22, 435.09]
95% CI:
[−36.58, 11.28]
95% CI:
[−4.02, 0.06]
95% CI:
[−0.35, 0.11]
: 430.55
: 4.39
: 424.94
: 437.60
: 39.34
: 41.32
: 5.31
: 5.43
Endurance males vs. reference males T = −14.444
df = 76.35
p < 2.2 × 10 ***
T = 7.3219
df = 131.25
p = 2.198 × 10 ***
T = 1.4205
df = 73.37
p = 0.1597
T = −8.3279
df = 108.39
p = 2.757 × 10 ***
95% CI:
[−194.14, −147.10]
95% CI:
[49.41, 85.98]
95% CI:
[−0.57, 3.40]
95% CI:
[−1.28, −0.79]
: 376.25
: 546.87
: 371.39
: 303.69
: 36.43
: 35.01
: 5.65
: 6.69
Strength males vs. reference males T = 52.851
df = 98.008
p < 2.2 × 10 ***
T = 8.9766
df = 148.03
p = 1.168 × 10 ***
T = 2.4131
df = 149.38
p = 0.01703***
T = −9.3513
df = 148.25
p < 2.2 × 10 ***
95% CI:
[359.29, 387.32]
95% CI:
[61.40, 96.07]
95% CI:
[0.33, 3.35]
95%CI:
[−1.35, −0.88]
: 376.25
: 2.95
: 371.39
: 292.65
: 36.43
: 34.58
: 5.65
: 6.77
11747.957.1750.14.4
12742.856.47454.3
13743.155.9745.24.3
SD 2.90.62.90.1
2.01797.267.1800.14.8
2.02790.165.4792.84.7
2.03794.866.5797.64.8
SD 3.60.93.70.1
31788.075.4791.75.5
32784.574.0788.05.4
33786.774.1790.25.4
SD 1.80.81.90.1
41471.359.9475.27.2
42466.858.6470.57.2
43469.058.8472.77.1
SD 2.30.72.40.1
51503.360.2506.96.8
52499.760.1503.46.9
53499.560.8503.26.9
SD 2.10.42.10.1
61580.562.9583.96.2
62576.762.9580.26.2
63574.162.7577.66.2
SD 3.20.13.20.0
SD 2.60.62.70.0
Reference Males (n = 100)Endurance Males (n = 36)
Parameter PREPOSTΔ-Pre–Post (%)Δ-Value
Ref-Group
BM (kg)76.7 ± 12.976.1 ± 13.0−0.9 ± 0.8_
Bioelectrical
R (Ω)631.6 ± 91.1543.1 ± 56.7535.5 ± 55.9−1.4 ± 2.388.6
Xc (Ω)62.2 ± 6.562.6 ± 8.363.4 ± 7.02.8 ± 18.7−0.4
R/h (Ω/m)304.5 ± 37.2303.7 ± 32.0299.4 ± 31.9−1.4 ± 2.30.8
Xc/h (Ω/m)45.0 ± 4.135.0 ± 4.935.5 ± 4.22.8 ± 18.710
PA (°)5.7 ± 0.86.7 ± 0.56.8 ± 0.51.3 ± 1.9−1
Z (Ω/m)634.8 ± 91.0546.9 ± 56.8539.3 ± 56.1−1.3 ± 2.387.9
r (R/h, Xc/h)0.670.470.75_—_—
= 100) = 54)

BM (kg)77.8 ± 11.377.86 ± 11.30.1 ± 0.7
Bioelectrical
R (Ω)631.6 ± 91.1515.8 ± 53.9507 ± 50.0−1.6 ± 1.9124.6
Xc (Ω)62.2 ± 6.560.9 ± 5.759.6 ± 5.42.18±5.51.3
R/h (Ω/m)304.5 ± 37.2292.6 ± 32.8287.7 ± 30.6−1.6 ± 1.816.8
Xc/h (Ω/m)45.0 ± 4.134.6 ± 3.633.8 ± 3.4−2.0 ± 2.811.2
PA (°)5.7 ± 0.86.6 ± 0.46.7 ± 0.6−0.3 ± 1.8−1
Z (Ω/m)634.8 ± 91.0519.5 ± 53.9510.6 ± 50.1−1.6 ± 1.8124.2
r (R/h, Xc/h)0.670.670.62_—_—
Reference Females (n = 100) Endurance Females (n = 29)
Parameter PREPOSTΔ-Pre–Post (%)Δ-Value
Ref-Group
BM (kg)62.2 ± 7.262.0 ± 7.30.3 ± 1.37_
Bioelectrical
R (Ω)697.8 ± 69.0680.4 ± 60.1674.1 ± 59.30.93 ± 1.3517.4
Xc (Ω)64.6 ± 8.765.4 ± 5.665.5 ± 6.2−0.15 ± 9.7−0.8
R/h (Ω/m)428.7 ± 45.2414.6 ± 36.6410.69 ± 35.70.95 ± 0.0218
Xc/h (Ω/m)39.7 ± 4.639.9 ± 3.939.96 ± 4.2−0.15 ± 7.14−0.2
PA (Ω)5.3 ± 0.56.7 ± 0.56.8 ± 0.5−1.47 ± 0−1.4
Z (Ω/m)700.9 ± 69.1683.6 ± 59.9677.4 ± 59.20.91 ± 0.0217.3
r (R/h, Xc/h)0.600.390.27_—_—
= 100) = 29)

BM (kg)57.9 ± 8.157.9 ± 8.10
Bioelectrical
R (Ω)697.8 ± 69.0705.0 ± 85.5698.3 ± 87.90.95 ± 2.73−7.2
Xc (Ω)64.6 ± 8.766.5 ± 6.165.8 ± 6.91.06 ± 11.6−1.9
R/h (Ω/m)428.7 ± 45.2437.6 ± 56.0433.5 ± 57.70.95 ± 2.95−8.9
Xc/h (Ω/m)39.7 ± 4.641.3 ± 4.640.9 ± 5.10.98 ± 0.09−1.6
PA (Ω)5.3 ± 0.55.4 ± 0.55.4 ± 0.60−0.1
Z (Ω/m)700.9 ± 69.1708.2 ± 85.4701.5 ± 87.90.95 ± 2.91−7.3
r (R/h, Xc/h)0.600.630.69_—_—
The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

Abdelnour, M.; Berkachy, R.; Nasreddine, L.; Fares, E.-J. Bioelectrical Impedance Vector Analysis (BIVA) for Assessment of Hydration Status: A Comparison between Endurance and Strength University Athletes. Sensors 2024 , 24 , 6024. https://doi.org/10.3390/s24186024

Abdelnour M, Berkachy R, Nasreddine L, Fares E-J. Bioelectrical Impedance Vector Analysis (BIVA) for Assessment of Hydration Status: A Comparison between Endurance and Strength University Athletes. Sensors . 2024; 24(18):6024. https://doi.org/10.3390/s24186024

Abdelnour, Maria, Rédina Berkachy, Lara Nasreddine, and Elie-Jacques Fares. 2024. "Bioelectrical Impedance Vector Analysis (BIVA) for Assessment of Hydration Status: A Comparison between Endurance and Strength University Athletes" Sensors 24, no. 18: 6024. https://doi.org/10.3390/s24186024

Article Metrics

Article access statistics, supplementary material.

ZIP-Document (ZIP, 135 KiB)

Further Information

Mdpi initiatives, follow mdpi.

MDPI

Subscribe to receive issue release notifications and newsletters from MDPI journals

IMAGES

  1. (PDF) Design and Analysis of Experiments with Examples in R

    design and analysis of experiments in r

  2. Design and Analysis of Experiments, 8th Edition, Douglas C. Montgomery

    design and analysis of experiments in r

  3. 4 Design of Experiments (DoE)

    design and analysis of experiments in r

  4. Design and Analysis of Experiments and Observational Studies Using R

    design and analysis of experiments in r

  5. (PDF) Design and Analysis of Experiments

    design and analysis of experiments in r

  6. Design and Analysis of Experiments with R

    design and analysis of experiments in r

VIDEO

  1. salt analysis experiments

  2. 典型金属材料金相组织的观察与分析实验Observation and analysis experiments on the metallographic structure of typical

  3. Крутые эксперименты!

  4. Design of Experiments, Online Statistical software's

  5. To Enhance Experiments, R&D, and Optimize FaaS Services

  6. Completely Randomized Design (CRD)

COMMENTS

  1. CRAN Task View: Design of Experiments (DoE) & Analysis of Experimental Data

    This task view collects information on R packages for experimental design and analysis of data from experiments. Packages that focus on analysis only and do not make relevant contributions for design creation are not considered in the scope of this task view. Please feel free to suggest enhancements, and please send information on new packages or major package updates if you think they belong ...

  2. PDF Design of Experiments in R

    „Design of Experiments and Analysis of Experimental Data" (or brief: Experimental Design) started February 2008 currently contains 37 . R. packages related to Design of Experiments Main purposes Pointer to existing functionality support synergies, avoid double work

  3. Experimental Design and Process Optimization with R

    The present document is a short and elementary course on the Design of Experiments (DoE) and empirical process optimization with the open-source Software R. The course is self-contained and does not assume any preknowledge in statistics or mathematics beyond high school level. Statistical concepts will be introduced on an elementary level and ...

  4. Design and Analysis of Experiments with R

    ABSTRACT. Design and Analysis of Experiments with R presents a unified treatment of experimental designs and design concepts commonly used in practice. It connects the objectives of research to the type of experimental design required, describes the process of creating the design and collecting the data, shows how to perform the proper analysis ...

  5. PDF daewr: Design and Analysis of Experiments with R

    Title Design and Analysis of Experiments with R Version 1.2-11 Date 2023-09-04 Maintainer John Lawson <[email protected]> Description Contains Data frames and functions used in the book ``Design and Analysis of Experi-ments with R'', Lawson(2015) ISBN-13:978-1-4398-6813-3. License GPL-2 Depends R (>= 3.5.0)

  6. Some Basic Concepts about Design of Experiments and How to ...

    The design of experiments (DOE) deals with the planning and performance of tests with the objective of generating data. Statistical analysis of these data will provide objective evidence that will allow the researcher to resolve questions about a given situation, process or phenomenon.

  7. PDF Design and Analysis of Experiments with R

    1.2 Beginnings of Statistically Planned Experiments 2 1.3 De nitions and Preliminaries 2 1.4 Purposes of Experimental Design 5 1.5 Types of Experimental Designs 6 1.6 Planning Experiments 7 1.7 Performing the Experiments 9 1.8 Use of R Software 12 1.9 Review of Important Concepts 12 1.10 Exercises 15 2 Completely Randomized Designs with One ...

  8. Experimental Design in R Course

    Experimental design is a crucial part of data analysis in any field, whether you work in business, health or tech. If you want to use data to answer a question, you need to design an experiment! In this course you will learn about basic experimental design, including block and factorial designs, and commonly used statistical tests, such as the ...

  9. Design and Analysis of Experiments with R

    John Lawson. CRC Press, Dec 17, 2014 - Mathematics - 620 pages. Design and Analysis of Experiments with R presents a unified treatment of experimental designs and design concepts commonly used in practice. It connects the objectives of research to the type of experimental design required, describes the process of creating the design and ...

  10. Design and Analysis of Experiments

    Daniel Voss, Danel Draguljić. Second edition includes new material on screening experiments and analysis of mixed models, a new chapter on computer experiments, added "Using R" sections, updated SAS output, and use of SAS Proc Mixed. Presents a step-by-step guide to design, including a planning checklist that emphasizes practical ...

  11. PDF Design of Experiments in R

    Design of Experiments is also used in many other fields of application, e.g. medicine, psychology, sensometrics or quality improvement. The talk gives an overview of what functionality is available in R for experimental design and analysis of experimental data in general (cf. also the Task View on Experimental Design and Analysis of ...

  12. PDF R Companion to Montgomery's Design and Analysis of Experiments (2005)

    The 6th edition of Montgomery's book, Design and Analysis of Experiments, has many more to do with the various kind of experimental setups commonly used in biomedical research or industrial engineering, and how to reach significant conclusions from the observed results. This is an art and it is called the Design of Experiment (doe). The

  13. PDF Design of Experiments in R

    Make existing analysis capabilities accessible through RcmdrPlugin.DoE package: linear model functions in general (are in R-commander already) simple plotting facilities for orthogonal 2-level experiments from package FrF2 analysis facilities for response surface designs from package rsm Later

  14. An R package for Design and Analysis of Experiment

    This repository is for Wiley's companion R package to Douglas C. Montgomery's "Design and Analysis of Experiments" 10th edition. The most current version of the PDF Guide will walk you through how to use R to work examples and exercises in the text. The PDF guide also shows how to install the most current package.

  15. Extended R Examples for A First Course in Design and Analysis of

    1.3 Packages. Base R can be extended with packages (libraries) that contain additional functions and data. There are literally tens of thousands of R packages available; fortunately, we won't need all of them. Start up RStudio. I believe the default will bring up a window with four panes.

  16. Design and Analysis of Experiments and Observational Studies using R

    Design and Analysis of Experiments and Observational Studies using R A Volume in the Chapman & Hall/CRC Texts in Statistical Science Series. Click on the cover to read. Sample Course Documents. Sample course documents (i.e., slides, assignments, etc.) are available here. R package. scidesignR is an R package that contains the data used in the book.

  17. PDF A First Course in Design and Analysis of Experiments

    Contents Preface xvii 1 Introduction 1 1.1 Why Experiment? . . . . . . . . . . . . . . . . . . . . . . 1 1.2 Components of an Experiment . . . . . . . . . . . . . . . 4

  18. Design and Analysis of Experiments with R

    Description. Design and Analysis of Experiments with R presents a unified treatment of experimental designs and design concepts commonly used in practice. It connects the objectives of research to the type of experimental design required, describes the process of creating the design and collecting the data, shows how to perform the proper ...

  19. Design and Analysis of Experiments with Examples in R

    This introductory text on Design and analysis of experiments is developed to satisfy the curriculum of the course at the polytechnics and undergraduate level in the universities in the physical ...

  20. Design and analysis of experiments with R

    Design and analysis of experiments with R by Lawson, John, 1947-Publication date 2015 Topics Sampling (Statistics), Statistics -- Data processing, R (Computer program language) Publisher Boca Raton : CRC Press, Taylor & Francis Group ... Robust parameter design experiments -- Experimental strategies for increasing knowledge ...

  21. Design and Analysis of Experiments with R

    This function uses Daniel's Method to find an outlier in an unreplicated 2^ { (k-p)} design. gear. Unreplicated split-plot fractional-factorial experiment on geometric distortion of drive gears. -- H --. halfnorm. This function makes a half normal plot of the elements of the vector called effects. hardwood.

  22. PDF Design and Analysis of Experiments

    The Three R's of Experimental Design. Randomization. Replication. Stratify (block) Randomization - It is important to randomize because it averages out the effect of all other lurking variables - it doesn't remove their effects, but makes, on average, their effects equal in all groups.

  23. Design and Analysis of Experiments with R (Chapman & Hall/CRC Texts in

    Design and Analysis of Experiments with R presents a unified treatment of experimental designs and design concepts commonly used in practice. It connects the objectives of research to the type of experimental design required, describes the process of creating the design and collecting the data, shows how to perform the proper analysis of the ...

  24. Improving rigor and reproducibility in western blot experiments with

    The experimental and statistical approach presented in this study provide insight into how more rigorous planning for western blot experimental design and corresponding statistical analysis ...

  25. Structure design and vibration control of vibration reduction boring

    2.2 Dynamic model analysis of vibration reduction boring bar with variable damping. As shown in Fig. 2, the dynamic model of the boring bar can be regarded as a two degree of freedom mass spring system.The boring bar body can be equivalent to a mass block M 1, the boring bar stiffness can be equivalent to a spring, and the stiffness coefficient is represented by K 1.

  26. Multi-experimental seismic analysis of low-rise thin reinforced

    Likewise, the experimental response of the TRCW-M3R10 prototype (F m T R C W) was acquired by a horizontal cell load and processed by the number of structural specimens located at the first level of the structure, which comprised the feedback force F f e e d T R C W to feedback the NS, completing the second experimental dynamic loop, as shown ...

  27. Experimental and computational analysis of the structure-activity

    By using this strategy, it will be possible to design the ideal structure of the IG and achieve excellent performance. ... et al. Experimental and computational analysis of the structure-activity relationship of ionic gel electrolytes based on bistrifluoromethanesulfonimide salts for supercapacitors. Journal of Materiomics, 2024, 10(6): 1222 ...

  28. Development and analysis of Hastelloy-X alloy butt joint ...

    The investigation describes the processing and analysis of laser welded Hastelloy-X (HX) alloy joints through this paper. HX alloys are employed in aerospace and nuclear industries, especially for high-temperature applications. In this research, a series of laser beam welding (LBW) experiments denoted as E1 to E9 were conducted on the HX base metal (BM), adhering to the L9-orthogonal array (L9 ...

  29. The impact of food packaging on measured food intake: A systematic

    Chronic, excess energy intake contributes to overweight and obesity. Solutions are needed to guide consumers towards portion control, especially for high energy density (HED), palatable foods. Food packaging, a key element of the eating environment, offers a potential solution. To investigate whether packaging design influences measured food intake, a systematic search was undertaken in four ...

  30. Bioelectrical Impedance Vector Analysis (BIVA) for Assessment of ...

    The study has some limitations including the quasi-experimental design with no blinding of the procedure. This could have biased athletes to drink more water before and during the experiment. Moreover, some participants fasted for 4 h whereas others fasted for 8 h and this could have affected the accuracy of BIA measurements and group analysis.