Recently, I set up my database workflow using a Virtual Private Server (VPS) and developed R code to easily access and analyze this database. I’ll introduce them one by one. 1) Data in local PC to MySQL Server I want to upload my data in my PC to a MySQL server hosted on a Virtual Private Server (VPS) managed by Digital Ocean. To facilitate this, I have developed a Python script, which you can find in my Github; https://github.com/agronomy4future/python_sql_autoload. To…
In my previous post, I explained What logistic regression is and also explained odds, odds ratio and model equation. Today, I’ll introduce Firth’s Logistic Regression. Common logistic regression is mostly used when the sample size is sufficiently large, and the outcome variable (0 and 1) is well-balanced. Also, it is used when there is no perfect separation between predictors and outcomes. However, common logistic regression fails or becomes unstable in certain situations. Firth’s Logistic Regression is usually used when there…
In a modern workflow, we constantly bounce between various web-based note-taking tools. For a developer or researcher, the act of opening a browser, logging in, and navigating to a page is a significant context switch. When you are deep in a coding session or analyzing data in the terminal, grabbing the mouse to record a quick thought or a draft email breaks your ‘Flow State.’ Every second spent away from the terminal is a second of lost concentration. To solve…
On Linux, OneDrive cannot be downloaded and synchronized as seamlessly as it is in a Windows environment. Today, I will provide a step-by-step guide on how to install and sync OneDrive on your Linux PC. 1) Set up a new remote connection. First, let’s open the terminal using the shortcut Ctrl + Alt + T and type rclone config Once the output appears, type n for a ‘New remote’ connection. Next, you’ll need to name your new remote. For this…
When analyzing grain size, we’ve used high-throughput image scanning machines. However, if we can detect grains using R code, estimating grain size becomes possible. In my previous R function, colorcapture(), we were able to detect fruit size and estimate its surface area in 2D. Recently, I developed a new R function called phenokio() specifically designed to estimate the grain area of cereals. 1. Install the function Before installing, please download Rtools (https://cran.r-project.org/bin/windows/Rtools) If the function is not able to be…
1) Nitrogen Immobilization Nitrogen immobilization is the process by which soil microorganisms take up inorganic nitrogen (mainly NH₄⁺ and NO₃⁻) from the soil solution and convert it into organic forms inside their own cells. In other words, the N is temporarily “locked up” in microbial cells and is not available for plants. Microbes need nitrogen to build proteins, enzymes, and DNA. When residues with high carbon-to-nitrogen (C:N) ratios (e.g., straw) are added to soil, microbes require extra N to decompose…
Meta-analysis is a quantitative method that synthesizes results from multiple independent studies to identify overall patterns, effect sizes, and sources of variability for specific treatments or research questions. It is particularly powerful for summarizing existing evidence and placing new findings within the context of current scientific trends through transparent and reproducible statistical integration of prior studies. A rigorous meta-analysis begins with a comprehensive and systematic literature search using predefined search terms and Boolean combinations across multiple databases. Study selection and…
When data is arranged, it can be structured either vertically (row-based) or horizontally (column-based). The choice depends on your preference for organizing data. However, when running statistics, data should be arranged row-based, as variables need to be in the same column. On the other hand, when calculating per variable, it is much easier to organize data column-based, allowing for simpler calculations. Regardless of the approach, well-organized data is essential, and the ability to restructure data is a valuable skill. Today,…
When conducting a meta-analysis, results from multiple studies must be synthesized despite differences in measurement scales, experimental designs, and reported units. To enable meaningful comparisons across studies, effect sizes are commonly used. One of the most widely applied effect size metrics in ecological and agricultural research is the Log Response Ratio (LRR). The Log Response Ratio is calculated as: LRR = ln(Xt / Xc)where Xt represents the mean response of the treatment group and Xc represents the mean response of…
When performing a t-test, a confidence interval can be obtained. Below, I describe how to calculate the confidence interval manually, step by step. 1) Difference in means 2) Pooled variance 3) Standard error 4) Degrees of freedom 5) Confidence interval and let’s calculate the confidence interval We aim to develop open-source code for agronomy ([email protected]) © 2022 – 2025 https://agronomy4future.com – All Rights Reserved. Last Updated: 12/12/2025
Recently, I’ve been focusing on Bayesian statistics. To organize the concepts for myself as well, I’m going to explain Bayes’ theorem as simply as possible. Let’s import a dataset from Kaggle. https://www.kaggle.com/datasets/cameronseamons/electronic-sales-sep2023-sep2024 This dataset contains information used to analyze customer purchasing behavior at an electronics store. You can download it from Kaggle after creating an account. I’ll load the data directly using Python code. I’m using Google Colab. This is how I imported the dataset directly from Kaggle. For more…
I don’t manage my data with Excel, but primarily with SQL. I also don’t store any data on my local PC—not even SQL files. Instead, I keep everything on the server and access it whenever I need it. Previously, I shared how to import data to a cloud MySQL server and how to import data from the MySQL server into R. □ How to Import Data to a Cloud MySQL Server (DigitalOcean) Using Python from Command Prompt□ How to Import…
In my previous post, I explained how to calculate Cook’s Distance step by step, and noted that in R you can simply use the function cooks.distance(). However, this simple function only provides the Cook’s Distance values. In my previous post, I explained the formula for Cook’s Distance step by step, including how to compute residuals, leverage, and internal studentized residuals. To make it easier to calculate these values, I recently developed an R package called datacooks(). □ datacooks() Before installing…
When analyzing yield, we usually use bar or box plots based on the mean. However, sometimes we need to visualize the spatial variation of yield. In those cases, GIS maps are commonly used, often created with powerful software such as ArcGIS, which I have used and found to be excellent. The drawback is that ArcGIS is not free. If you want to create a simple yield map using R, it is possible, although the map quality may be slightly lower….
1) to replace text with an empty string 2) to replace text with other texts
Previous post□ How to analyze linear plateau model in R Studio? In my previous post, I explained how to analyze the linear plateau model. I simulated yield data for five different crop varieties with varying sulphur applications and suggested that the optimum sulphur application would be 23.3 kg/ha based on the linear plateau model. In this post, I’ll explain how to analyze the quadratic plateau model using the same data in R Studio. 1) Data upload If you run the…
The Gamma distribution is a flexible family of continuous probability distributions defined only for non-negative values (x>0). It’s commonly used to model quantities that represent time, size, or waiting periods—anything that can’t go below zero and often shows right-skewed behavior (a long tail toward larger values). In essence, the Gamma distribution describes how likely different positive values are to occur, determined by two key parameters: the shape (α) and the scale (θ). Together, these parameters control the curve’s form and…
#guides(fill=”none”) +
A Virtual Private Server (VPS) is a virtualized computer within a larger physical server. It acts like an independent server, offering dedicated resources and control at a lower cost than a full physical machine. VPS hosting uses this setup to give users private, customizable environments for web hosting or applications. Some reputable VPS providers include DigitalOcean, Linode, and Vultr — I currently use DigitalOcean. A Droplet in DigitalOcean is basically a virtual private server (VPS) — a self-contained Linux machine…
In my previous post, I introduced how to import data into a Cloud MySQL database using Python from the Command Prompt. By typing the code below in your Command Prompt, you can automatically import data into your MySQL server. This is the next step. After importing data into the MySQL server, what if I want to bring that data into R? Is that possible? Yes—it is. 1) install and load the package and library. 2) Create a connection to the…
This colorcapture() R function provides easy image analysis to estimate fruit surface area. 1. Install the pacakge Before installing colorcapture(), please download Rtools (https://cran.r-project.org/bin/windows/Rtools), and install the following package. □ colorcapture() 2. Basic code If you want to change other criteria, type ?colorcapture to see detailed information about the colorcapture() function. The default code is set to capture the yellow color. If you run the code, it will detect the yellow surface area. The frame size was 24 cm × 24 cm. and the…
In my previous post, I introduced how to import data into SQL. That process was done manually to import data into a local SQL server. I use MySQL. A local SQL server means you can only access the data from your own computer. If you use another computer elsewhere, you won’t be able to access that SQL data. However, I want to access my data universally—from any device, at any location—to manage it more efficiently. To achieve this, we can…
When measuring leaf area or fruit surface area, these have usually been measured manually, which is time-consuming and often inaccurate, especially when the shape is not perpendicular. For more reliable measurement, image analysis provides a good alternative. To facilitate this process, I developed an R function, greencapture() that captures the green area of leaves or fruits in the lab using a fixed frame. 1. Install the greencapture() package 2. Set up folders for input and output images I will save…
When analyzing data, we sometimes need to analyze cumulative data. When calculating cumulative data, grouping is important, and it takes time to perform the grouping and calculations. To simplify this process, I developed an R package called datacume(). Let’s upload the dataset. This dataset contains biomass measurements across different treatments over time, recorded for various plants and their branches. I want to calculate the cumulative biomass over time. To do this accurately, I’ll first calculate the average biomass of the…
In my previous post, I suggested different statistical models for agrivoltaics studies to explain why we should consider Linear Mixed Models in this field. In many agrivoltaics studies, researchers overlook the actual experimental layout and analyze data using split-plot or RCBD models, focusing only on treatment variables (e.g., inside vs. outside the solar panel array). However, this approach does not accurately reflect real field conditions. ■ [STAT Article] Statistical Models in Agrivoltaics: Linear Mixed Models Across Different Field Layouts Randomization…
https://github.com/agronomy4future/r_code/blob/main/geom_mark_ellipse.ipynb We aim to develop open-source code for agronomy ([email protected]) © 2022 – 2025 https://agronomy4future.com – All Rights Reserved. Last Updated: 19/07/2025
Agrivoltaics is the study and practice of combining agriculture and solar energy production on the same land. The core idea is to install solar panels above or among crops, allowing for simultaneous food and energy production. Generally, an agrivoltaics study investigates how this dual-use approach affects crop growth and yield (due to changes in light, temperature, and moisture), microclimate conditions under the panels, solar panel efficiency influenced by vegetation, as well as land-use efficiency, economic outcomes, and sustainability metrics. Today,…
This post is part of a series introducing R packages that present methods for calculating phenotypic plasticity. □ [R package] Finlay-Wilkinson Regression model (feat. fwrmodel)□ [R package] Calculate the responsiveness of each treatment relative to a control (Feat. deltactrl)□ [R package] Quantifying Reaction Norm Plasticity from Slopes to Individual Responses (Feat. nrmodel) In my previous post, I explained how to quantify phenotypic plasticity in crops and described three different approaches: Responsiveness, Reaction Norm, and the Finlay-Wilkinson Regression Model. □ Quantifying…
□ Quantifying Phenotypic Plasticity of Crops In my previous post, I explained how to quantify phenotypic plasticity in crops and described three different approaches: 1) Responsiveness, 2) Reaction Norm, and the 3) Finlay-Wilkinson Regression Model. Responsiveness is calculated as (Treatment − Control) / Control. It indicates how the dependent variable (e.g., yield) responds to a given treatment relative to the control. The responsiveness value can range from −1 (complete reduction) to values greater than 0, depending on the magnitude of…