[캐나다 농업 일기] 2022년 궬프 지역 옥수수 수확
2022년 옥수수 필드 연구를 마무리 짓고 수확을 마무리 하였습니다. 파종일: 2022년 5월 9일개화일: 2022년 7월 27일등숙일: 2022년 10월 28일수확일: 2022년 11월 14일
Stories about crops and statistics (plus coding)
2022년 옥수수 필드 연구를 마무리 짓고 수확을 마무리 하였습니다. 파종일: 2022년 5월 9일개화일: 2022년 7월 27일등숙일: 2022년 10월 28일수확일: 2022년 11월 14일
I attended the 2022 ASA, CSSA and SSSA International Annual Meeting at Baltimore, MD, and presented my Ph.D research topic.
geom_text(aes(family=”serif”,fontface=6), x=65, y=0.02, label=paste(“Pistolo\n”,pistolo_CT_MEAN, “±”, round(pistolo_CT_SE,digits=2)),size=7, col=”Black”)
Let’s upload one data to R. Now, I’d like to change the name of column as field → locationgenotype → varietyblock → repstreatment → experimentshoot → branchgrain_number → GNgrain_weight → GW I introduce two ways to change column names. 1) using colnames() 2) using rename() in dplyr package In this time, I’ll use dplyr package.
I uploaded one .csv file to GitHub. Now I want to analyze this data in R. Simply I can download this file and upload to R. But let’s directly upload this file from GitHub to R. First, we need to know the URL address of this file. If you click your file name in GitHub, … Read more
I’ll generate one data. Then, I’ll make a bar graph about this data. To make a bar graph, data should be summarized. Now, I want to change legend name from N0 to 0kg N/ha, and N1 to 200kg N/ha. Simply we can add more code like this; scale_fill_manual(label=c(“0kg N/ha”,”200kg N/ha”), values=c(“grey75″,”grey25”)) What if we want … Read more
label=paste(round(mean(na.rm=T,long_low$KW),digit=1),”mg”)
The reason why we use t-test, not z-test is because we don’t know the variance (σ2) of the population. The sample mean (x̅) is an unbiased estimator for the population mean (μ), and therefore we can estimate μ from x̅ (E(x̅) = μ). How about variance? If we know σ2, sample variance could be estimated … Read more
The genotype is dependent on environmental changes. One genotype may strongly respond to certain environmental conditions, while another genotype may weakly respond to the same conditions. If some genotypes strongly respond under better conditions, they would be adaptable to the environment. Adaptability refers to the flexibility of a genotype in its response to improved environments. … Read more
One tomato farmer is growing tomato seedlings, and all of sudden he wants to investigate the amount of calcium in leaves. So, he selected four tomato seedlings, and he randomly chose three leaves in each seedling and investigated the amount of calcium. He measured twice in each leaf. This experimental design would be explained by … Read more
How much fertilizer do farmers lose when they apply nitrogen in season? Is volatilization a significant problem? Could 4R nutrient management and urease inhibitors help? To answer those questions, our research team had an interview with RealAgriculture in Canada. ArticleCorn School: How much nitrogen is lost when applying in season?
After 18 – 20 days from R1 stage, the kernels become mostly yellow and contain “milky” white fluid (i.e. starch accumulates). However, moisture content of kernel at R3 is around 80%. Starch accumulation continuously occurs in the endosperm, and cell division of endosperm is almost complete, indicating that growth of kernels since R3 would be … Read more
Here is an example of missing data. There are four different cultivars and I would like to determine if there is a difference in yield among them. I have five replicates as blocks, so the experimental design is a Randomized Complete Block Design (RCBD), and we can analyze the data using one-way ANOVA with blocks. … Read more
This is our corn field for split N experiment in University of Guelph. Video recording by drone @ Guelph, Ontario in Canada.
New equation I suggest!! You might be already familiar with how to calculate pooled variance. This story is about pooled variance when blocks exist. If you run statistics programs, you’ll simply obtain pooled variance (also known as MSE), but you’ll never understand the concept of pooled variance if you just run software programs. Here is … Read more
I harvested some corn plants at R1 stage 12 days ago. Around R1 stage, final kernel number would be determined. In corn, kernel row number (girth) and number per row (length) in an ear are main yield components. Girth determination of the uppermost ear (the ear harvested) begins shortly after the ear shoot is initiated … Read more
When we set up our experimental design, it is not easy to decide the sample size because we don’t know exactly how many samples are required for our experiments. Of course the more, the better. However, eventually, we need to decide appropriate sample size according to our time and resources. For example, if we want … Read more
In corn, R1 stage (silking) is one of the most critical period to determine yield potential. The emergence of silk is regarded as the initiation of reproductive period. Each ovule (potential kernel) on the ear develops its own silk (as the stigma of the female flower). Silks initiate to elongate around the V12. It stars … Read more
Sowing date: 12/05/2022 Today: 28/07/2022 (78 days from sowing) Current stage: R1 Site: Ontario, Canada Finally corn reached to R1. Now we can see corn silk will appear. Also, we regard corn developmental stage as R1 (silking) when tassels appears. It’s very important silking and tasseling are synchronized timely for successful pollination and fertilization (It … Read more
When we use WordPress, sometimes we want to increase width of page. Of course, we can do that using several plugin, but through simple CSS codes, we can adjust page width and font size. Let’s go to Dashboard > Appearance > Customize Then, choose Additional CSS Then write the below code in the window. After … Read more
This post is the summary of the article, ‘Ear Initiation & Size Determination in Corn’ by prof. Bob Nielsen in Purdue University, and therefore all storytelling in this post is based on the article. In a corn plant, kernel number of per ear (KN) is an important contributor to the grain yield potential. Therefore, developmental … Read more
Why plants at the front are always smaller than others inside plot? If there is a border effect, it’s likely to be taller. This mainly results from more soil compaction and also might be due to competition against grasses. Furthermore, N fertilizers would be less supplied at the front where the applicator just starts to … Read more
Sowing date: 12/05/2022 Today: 19/07/2022 (69 days from sowing) Current stage: V13-V14 Site: Guelph, Ontario, Canada Now corn is reaching to VT. Sooner or later, tassels will come out and also corn silk will appear. Now it’s so drought. When corn goes through drought stress, they roll up leaves to protect water loss. So, when … Read more
Sowing date: 30/04/2022 Today: 19/07/2022 (81 days from sowing) Current stage: R1 Site: Guelph, Ontario, Canada Corn variety which was planted at the end of April reached R1 stage.
Previous post□ What is ANCOVA (1/3)? The basic concept In previous post, I explained how to interpret ANCOVA table (red box in below tables). In this post, I’ll explain how to interpret Parameter Estimates (blue box in below table) in ANCOVA analysis. Let’s check the ‘Parameter Estimates’ table. Most statistical programs set up one level of … Read more
Today, I will explain Analysis of Covariance (ANCOVA). ANCOVA is a statistical technique that involves including covariates, which are additional variables that may impact the dependent variable (y) in addition to the independent variable (x). I have a dataset as shown below, and I would like to analyze crop yield based on different fertilizer types … Read more
Previous post□ What is ANCOVA (1/3)? The basic concept□ What is ANCOVA (2/3)? How to interpret Parameter Estimates In the previous posts, I explained the basic concept of ANCOVA and how to interpret statistical results. Now, I will discuss the most important concept that is not commonly mentioned. The statistical program provided the following model in … Read more
Here is one data Then I make a regression graph. Now, I’d like to reverse the order of x-axis by descending (60 to 0). So I delete the code, scale_x_continuous(breaks = seq(-0, 60, 10), limits = c(-0, 60)) and add a new code, scale_x_reverse(limits = c(60,0)) So the whole code is below. Now, you can … Read more
Sowing date: 12/05/2022 Today: 28/06/2022 (61 days from sowing) Current stage: V7-V8 Site: Guelph, Ontario, Canada Today is 61 days from sowing corn in Guelph. The current stage of corn is between V7 and V8. From now on, the demand of nitrogen in corn will be increased. Therefore, it would be important to track down … Read more
The current growth stage of corn in Guelph is between V8 and V9.