# Analys av fältförsök

Agronomer använder försöksanalys för att bedöma hur olika sorter, odlingsmetoder och insatser presterar, inklusive resultaten av variabel giva i precisionsjordbruk. Genom att samla in, analysera och tolka data från fältförsök får forskare insikter i samspelet mellan genetik, miljö och odlingsåtgärder. Denna kunskap ligger till grund för utvecklingen av växtodlingsstrategier som optimerar skördepotentialen samtidigt som insatsmedlen minimeras. Dessutom gör försöksanalys inte bara att man kan utvärdera hur effektiva precisionsodlingsmetoderna är, utan hjälper också till att identifiera tåliga sorter som kan trivas under olika och utmanande förhållanden, vilket i sin tur bidrar till livsmedelsförsörjningen.

{% hint style="info" %}
GeoPard stödjer också split-plot-försök för tvåfaktordesign, till exempel samma givor över olika hybrider. Utvärdera både huvudverkan och interaktionen giva × hybrid på samma georefererade försökslager.
{% endhint %}

## Databeredning

För effektiv försöksanalys krävs några viktiga dataset:

1. **Skördedataset**:\
   Detta dataset innehåller skördedatan.\
   Vi kan importera detta från [JohnDeere Operation Center](https://docs.geopard.tech/geopard-tutorials/swe/produktgenomgang-webbapp/import/importera-fran-myjohndeere) eller ladda upp det manuellt som [shapefile](https://docs.geopard.tech/geopard-tutorials/swe/produktgenomgang-webbapp/import/import-av-skordedata) eller som [maskinens proprietära format](https://docs.geopard.tech/geopard-tutorials/swe/produktgenomgang-webbapp/import/maskintillverkares-egna-format).
2. **Applikationsdataset**:\
   Detta är avgörande för att förstå den faktiskt utförda körningen på fältet. Minst innehåller det attribut som TargetRate, AppliedRate och vissa maskinrelaterade mått.\
   Precis som med skördedatasetet har vi alternativ att importera det från [JohnDeere Operation Center](https://docs.geopard.tech/geopard-tutorials/swe/produktgenomgang-webbapp/import/importera-fran-myjohndeere) eller ladda upp det manuellt som [shapefile](https://docs.geopard.tech/geopard-tutorials/swe/produktgenomgang-webbapp/import/import-av-data-for-kord-giva-sadd) eller som [maskinens proprietära format](https://docs.geopard.tech/geopard-tutorials/swe/produktgenomgang-webbapp/import/maskintillverkares-egna-format).
3. **Zoner/Fältparceller med försök/experiment**:\
   Dessa visar de planerade givorna för våra försök och ger insikt i försöksupplägget.\
   Om ett sådant datalager finns laddar vi upp det som [shapefile](https://docs.geopard.tech/geopard-tutorials/swe/produktgenomgang-webbapp/import/import-av-data-for-kord-giva-sadd) i AsApplied/AsPlanted eller Yield-kontrollen. Detta säkerställer kompatibilitet när du bygger EquationMaps och effektiviserar din upplevelse av försöksanalys.\
   Detta kan vara en enkelfaktorslayout eller en split-plot-layout med en andra behandlingsdimension, till exempel hybrid eller sort.\
   Om ett sådant datalager inte finns kan attributet TargetRate från applikationsdatasetet användas som ersättning för försöksevalueringar.
4. **Historiska fältpotentialzoner:**\
   Dessa zoner genereras av GeoPard (detaljer finns [HÄR](https://docs.geopard.tech/geopard-tutorials/swe/produktgenomgang-webbapp/zonkartor-och-analyser/flerariga-zoner)). De är användbara för att analysera försök med jämn historisk produktivitet. Detta är särskilt värdefullt när försöken är fördelade över regioner med varierande historisk produktivitet.

När vi har samlat in dessa dataset är nästa steg att starta utvärderingen av försöket.

## Dataöversikt

Följande data finns för jordbrukssäsongen 2023 för höstvete:

* Skördedataset som visar fördelningen av våtvikt *(Fig.1)*

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2Fg9Iv1AMpACkykujxvJ9s%2Fimage.png?alt=media&#x26;token=c2a61d7b-12e9-40b5-ad88-e6f820372638" alt=""><figcaption><p>Fig.1 Originalt skördedataset</p></figcaption></figure>

* Kväve-(N34)-VRA-plan (150 kg/ha) med 2 försöksrutor (120 kg/ha och 180 kg/ha)*(Fig.2)*

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FwkGO1Tp0Ji2KzUPhfQAE%2Fimage.png?alt=media&#x26;token=28a94d78-7360-43e2-865f-4a330028ab20" alt=""><figcaption><p>Fig.2 Kväve-(N34)-VRA-plan med försöksrutor</p></figcaption></figure>

* Applikationsdataset som visar tillämpad statistik *(Fig.3)*

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FF1YKbPqDw2HBDdnTqazX%2Fimage.png?alt=media&#x26;token=0f83b5b0-9220-4cea-984c-735338716048" alt=""><figcaption><p>Fig.3 Applikationsdataset</p></figcaption></figure>

* Historisk fältproduktivitet (*Fig.4*)

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FCrFaBaZxtjd0xEXZ9TUd%2Fimage.png?alt=media&#x26;token=9cf68b49-a17c-4f2f-be15-1b0bfee4993e" alt=""><figcaption><p>Fig.4 Historisk fältproduktivitet</p></figcaption></figure>

{% hint style="warning" %}
Skördedatasetet har inte kalibrerats: flera skördetröskor arbetade där, det finns vändtegs- och saknade dataspår, och störningar är tydliga. Det rekommenderas att köra operationerna Kalibrera skörd och Rensa ovanpå för bästa resultat.\
En steg-för-steg-guide finns på [LÄNK](https://docs.geopard.tech/geopard-tutorials/swe/agronomi/skordkalibrering-och-rensning).
{% endhint %}

Skördedatasetet, efter kalibrering och rensning, visas i *Fig.5*, tillsammans med uppdaterad statistik. Detta dataset kommer att användas i de följande stegen.

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2F62J6z5wt7ndZ4lWxZq58%2Fimage.png?alt=media&#x26;token=9e6ea416-81a0-4dfc-a94a-aa19f69c0ff3" alt=""><figcaption><p>Fig.5 Kalibrerat och rensat skördedataset</p></figcaption></figure>

## Koncept

Här är målet med försöksanalysen att fastställa den mest effektiva kvävegivan (N34) för fältet. Det finns utpekade områden med kvävegivor på 120 kg/ha, 150 kg/ha och 180 kg/ha. Denna data kommer dels från applikationsdatasetet och dels från det kalibrerade skördedatasetet.

Vi fokuserar vår analys på tre tydliga zoner:

* 120 kg/ha (utpekad som försökszon)
* 150 kg/ha (betraktad som huvudzon)
* 180 kg/ha (en annan försökszon)

Vårt arbetssätt kommer att omfatta följande utvärderingar:

1. **Planbaserat:** med den planerade variabla givan (VRA) kopplad till den kalibrerade skörden.
2. **Baserat på faktisk körning:** jämföra de faktiska tillämpade datamängderna mot den kalibrerade skörden.
3. **Baserat på faktisk körning och historisk produktivitet:** jämföra de faktiska tillämpade datamängderna mot den kalibrerade skörden överlagrad med historiska fältpotentialzoner.

Detta metodiska arbetssätt möjliggör en heltäckande bedömning av kvävets påverkan på skörden, baserat på både planerad och faktiskt utförd tillförsel.

## Planbaserat

Påverkan av ~~tillförd~~ planerad kväve-(N34) på skördefördelningen syns tydligt i de följande skärmbilderna *(Fig.6, Fig.7, Fig.8)*. Här är en kort sammanfattning av resultaten:

* <mark style="color:blå;">Huvudzonen, med en kvävegiva på 150 kg/ha, omfattar 45,8 ha och har en genomsnittlig skörd på 4,99 t/ha</mark> (*Fig.6*).
* <mark style="color:blå;">Den första försökszonen, med en kvävegiva på 180 kg/ha, omfattar 1,76 ha och ger en genomsnittlig skörd på 6,5 t/ha</mark> (*Fig.7*).
* <mark style="color:blå;">Den andra försökszonen, med 120 kg/ha kväve, omfattar 1,86 ha och ger en genomsnittlig skörd på 6,39 t/ha</mark> (*Fig.8*).

Resultaten väcker en viktig fråga: <mark style="color:orange;">Varför verkar den lägre givan vara mer effektiv än den högre?</mark> För att få djupare insikter innebär nästa steg att[ utvärdera försöken med de faktiska tillämpade datan](#applied-based-evaluation).

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FraVRCgxuEVTGfTKK4c0s%2Fimage.png?alt=media&#x26;token=57602ad8-bd9f-46b4-adcf-5de53f99fe7b" alt=""><figcaption><p>Fig.6 Huvudzon med N34 150 kg/ha</p></figcaption></figure>

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FNM80FFM2TQULwV54LuqW%2Fimage.png?alt=media&#x26;token=1542336e-20a2-4a77-a343-9a92f67c90eb" alt=""><figcaption><p>Fig.7 Försökszon med N34 180 kg/ha</p></figcaption></figure>

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FITg4CWAJb27X28jquqyO%2Fimage.png?alt=media&#x26;token=cff8d0b3-5896-4f34-bbdc-6a34a7352982" alt=""><figcaption><p>Fig.8 Försökszon med N34 120 kg/ha</p></figcaption></figure>

Längre ned hittar du en fördjupad genomgång av formlerna och konfigurationerna som användes under utvärderingen.

{% hint style="info" %}
För att gå djupare in i ekvationsmetoden och hur den körs, se våra guider för både [användargränssnittet](https://docs.geopard.tech/geopard-tutorials/swe/produktgenomgang-webbapp/ekvationsbaserad-analys) och [API:et](https://docs.geopard.tech/geopard-tutorials/swe/api-dokumentation/diagram-med-grundlaggande-floden/5.-kor-ekvationer).
{% endhint %}

Här är ekvationerna som ska köras för att återskapa beräkningarna.

1. Huvudzon med 150 kg/ha:\
   `Yield_Main = np.where(Zone==1, Yield_WetMass, np.nan)`
2. Försök med 120 kg/ha:\
   `Yield_Zone = np.where(Zone==3, Yield_WetMass, np.nan)`
3. Försök med 180 kg/ha:\
   `Yield_Zone = np.where(Zone==2, Yield_WetMass, np.nan)`

Det är viktigt att aktivera *Numpy* *(Fig.9)* och stänga av *Interpolering* *(Fig.10)*.

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FixfesbJQxzVC93JW3H7b%2Fimage.png?alt=media&#x26;token=dcab1611-8a33-4450-9b06-5ba24760e8c3" alt=""><figcaption><p>Fig.9 Aktivera "Numpy"</p></figcaption></figure>

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2F0V4AxWhE0zraNfH1axHQ%2Fimage.png?alt=media&#x26;token=d86b3921-c313-4946-9455-35fbef29cc13" alt=""><figcaption><p>Fig.10 Stäng av användning av "interpolerad" data</p></figcaption></figure>

## Baserat på faktisk körning

En tydlig iakttagelse är att den faktiska tillämpade givan under försöket inte konsekvent överensstämmer med den planerade (mål-)givan. Fördelningen sträcker sig specifikt från 120 kg/ha upp till 189 kg/ha *(Fig.11)*. Med tanke på denna variation blev det viktigt att sätta en toleransnivå för fel. Därför bedömdes en noggrannhet på ±5 % vara en acceptabel gräns för att anse att försöket lämpar sig för utvärdering.

I de följande skärmbilderna (*Fig.12, Fig.13, Fig.14)* visas skördens statistiska fördelning med fokus på de faktiskt tillämpade kväve-(N34)-mängderna. Här är sammanfattad statistik, med hänsyn till acceptansgränsen på ±5 % noggrannhet:

* <mark style="color:blå;">Huvudzonen vid 150 kg/ha hade en tillämpad yta på 43,5 ha och gav en genomsnittlig skörd på 4,9 t/ha</mark> (*Fig.12*).
* <mark style="color:blå;">Den första försökszonen vid 180 kg/ha omfattade en yta på 1,47 ha och gav en genomsnittlig skörd på 6,5 t/ha</mark> (*Fig.13*).
* <mark style="color:blå;">Den andra försökszonen vid 120 kg/ha omfattade en yta på 1,44 ha och gav en genomsnittlig skörd på 6,3 t/ha</mark> (*Fig.14*).

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FQ5Aawg78wAVNsybfq2RB%2Fimage.png?alt=media&#x26;token=293f63d1-414b-4e23-be09-b54844f21988" alt=""><figcaption><p>Fig.11 Faktiskt tillämpade givor i försök</p></figcaption></figure>

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FyBHR1SLHA9ujXr5wk8xA%2Fimage.png?alt=media&#x26;token=78f67150-560b-4495-974e-30407e25aca9" alt=""><figcaption><p>Fig.12 Huvudzon med N34 150 kg/ha ±5%</p></figcaption></figure>

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FeOMxE76sp3ugO1NffWPN%2Fimage.png?alt=media&#x26;token=9e3f47ae-250a-41dd-b9ce-f3d4bcaab85e" alt=""><figcaption><p>Fig.13 Försökszon med N34 180 kg/ha ±5%</p></figcaption></figure>

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FJUBlkogPAF2fTSYa0EOW%2Fimage.png?alt=media&#x26;token=48dd3f24-df15-40da-905c-9012b70971d2" alt=""><figcaption><p>Fig.14 Försökszon med N34 120 kg/ha ±5%</p></figcaption></figure>

För att bättre förstå metodiken och detaljerna i dessa resultat finns de använda ekvationerna nedan:

1. Faktiskt tillförd kväve i försök:\
   `Applied_Trial = np.where((Zone == 3) | (Zone == 2), Applied_Value, np.nan)`
2. Huvudzon med 150 kg/ha med 5 % acceptans:\
   `Yield_Main = np.where((Zone == 1) & (Applied_Value >= 142.5) & (Applied_Value <= 157.5), Yield_WetMass, np.nan)`
3. Försök med 120 kg/ha med 5 % acceptans:\
   `Yield_Trial = np.where((Zone == 3) & (Applied_Value >= 114.0) & (Applied_Value <= 126.0), Yield_WetMass, np.nan)`
4. Försök med 180 kg/ha med 5 % acceptans:\
   `Yield_Trial = np.where((Zone == 2) & (Applied_Value >= 171.0) & (Applied_Value <= 189.0), Yield_WetMass, np.nan)`

## **Baserat på faktisk körning och historisk produktivitet**

Skördesiffrorna från försöken ligger konsekvent över den genomsnittliga skörden över hela fältet. En viktig faktor bakom denna skillnad verkar vara den historiskt högproduktiva zon där försöken genomfördes, som visas i *Fig.15* och *Fig.16*. För en mer nyanserad bedömning av försöken är det viktigt att ta hänsyn till produktivitetszonerna när resultaten analyseras.

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FH0aA9ZAe2IkSFtqcJe5B%2Fimage.png?alt=media&#x26;token=c58e81d6-b2df-4d6d-8682-3c750f38f2a1" alt=""><figcaption><p>Fig.15 Historiska fältpotentialzoner</p></figcaption></figure>

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2Fi2agD6onOT4aVaW9LvHh%2Fimage.png?alt=media&#x26;token=feecbfd5-2433-4061-85b9-8233446572d8" alt=""><figcaption><p>Fig.16 Historiska fältpotentialzoner som skördedataset</p></figcaption></figure>

I de följande skärmbilderna (*Fig.17, Fig.18, Fig.19)* visas skördens statistiska fördelning med fokus på de faktiskt tillämpade kväve-(N34)-mängderna överlagrade med historiska produktivitetszoner (skapade i GeoPard). Här är sammanfattad statistik, med hänsyn till ±5 % acceptans för de tillämpade mängderna:

* <mark style="color:blå;">Huvudzonen vid 150 kg/ha hade en tillämpad yta på 2,65 ha och gav en genomsnittlig skörd på 6,34 t/ha</mark> (*Fig.17*).
* <mark style="color:blå;">Den första försökszonen vid 180 kg/ha omfattade en yta på 1,08 ha och gav en genomsnittlig skörd på 6,41 t/ha</mark> (*Fig.18*).
* <mark style="color:blå;">Den andra försökszonen vid 120 kg/ha omfattade en yta på 1,78 ha och gav en genomsnittlig skörd på 6,33 t/ha</mark> (*Fig.19*).

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FPqw7u62bJAezDn9UW5qs%2Fimage.png?alt=media&#x26;token=ff101b42-7882-4dc3-a2ad-b57053900768" alt=""><figcaption><p>Fig.17 Huvudzon med N34 150 kg/ha överlagrad med historisk produktivitet</p></figcaption></figure>

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FmzkQRoU3HMMOm1fLLjNY%2Fimage.png?alt=media&#x26;token=ba8105da-877b-4234-9e1d-f59826949bf0" alt=""><figcaption><p>Fig.18 Försökszon med N34 180 kg/ha ±5% överlagrad med historisk produktivitet</p></figcaption></figure>

<figure><img src="https://1016759462-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2FYICBELdyAXXebKAzfLOR%2Fuploads%2FzTOI3v14HdZLANjM3o36%2Fimage.png?alt=media&#x26;token=f83a50dd-585d-4b66-937d-eb0d2f1c78f0" alt=""><figcaption><p>Fig.19 Försökszon med N34 120 kg/ha ±5% överlagrad med historisk produktivitet</p></figcaption></figure>

För att bättre förstå metodiken och detaljerna i dessa resultat finns de använda ekvationerna nedan:

1. Huvudzon med 150 kg/ha med 5 % acceptans överlagrad med historisk produktivitet:\
   `Yield_Main = np.where((Application_Zone == 1) & (Productivity_SubZone == 51) & (Applied_Value >= 142.5) & (Applied_Value <= 157.5), Yield_WetMass, np.nan)`
2. Försök med 120 kg/ha med 5 % acceptans överlagrad med historisk produktivitet:\
   `Yield_Trial = np.where((Application_Zone == 3) & (Productivity_SubZone == 51) & (Applied_Value >= 114.0) & (Applied_Value <= 126.0), Yield_WetMass, np.nan)`
3. Försök med 180 kg/ha med 5 % acceptans överlagrad med historisk produktivitet:\
   `Yield_Trial = np.where((Application_Zone == 2) & (Productivity_SubZone == 51) & (Applied_Value >= 171.0) & (Applied_Value <= 189.0), Yield_WetMass, np.nan)`

där

* delen `Productivity_SubZone == 51` återspeglar de högproduktiva zonerna med de tillämpade försöken,
* delarna `(Applied_Value >= 142.5) & (Applied_Value <= 157.5)` , `(Applied_Value >= 114.0) & (Applied_Value <= 126.0)`, `(Applied_Value >= 171.0) & (Applied_Value <= 189.0)` innefattar ±5 % noggrannhet från givorna `150`, `120`, `180` kg/ha.

## Sammanfattning

Skörderesultaten från försöken ligger nära den genomsnittliga skörd som observerats i fältets högproduktiva historiska zon. Med andra ord gav den experimentella tillförseln av N34-produkt i givorna <mark style="color:blå;">120 kg/ha - 150 kg/ha - 180 kg/ha</mark>, genomsnittliga skördar på <mark style="color:blå;">6,33 t/ha - 6,34 t/ha - 6,41 t/ha</mark> respektive, vilket inte har någon betydande påverkan på den skördade avkastningen inom den högproduktiva zonen.
