Ich habe eine data.frame
df
, die aus 4 Standorten besteht. Jeder Standort hat vier Parameter (A, B, C und D), und ich habe den Wert für jeden Parameter 2011-2014Trend-Analyse (steigende, abnehmende und keine Tendenz) in einem Tabellenformat
siteID <- c(rep("site1", 16), rep("site2", 16), rep("site3", 16),rep("site4", 16))
YEAR <- rep(c("2011", "2012", "2013", "2014"), 16)
parameter <- c(rep("A", 4), rep("B", 4), rep("C", 4), rep("D", 4),
rep("A", 4), rep("B", 4), rep("C", 4), rep("D", 4),
rep("A", 4), rep("B", 4), rep("C", 4), rep("D", 4),
rep("A", 4), rep("B", 4), rep("C", 4), rep("D", 4))
value <- c(seq(1, 4, by=1), seq(10, 40, by=10), seq(12, 18, by=2), seq(5, 20, by=5),
seq(3, 12, by=3), sample(13:18, 4), sample(15:22, 4), sample(10:18, 4),
seq(7, 1, by=-2), sample(15:22, 4), sample(15:19, 4), sample(10:20, 4),
seq(8, 5, by=-1), seq(50, 20, by=-10), seq(16, 10, by=-2), seq(20, 5, by=-5))
df <- data.frame(siteID, YEAR, parameter, value)
head(df, 20)
# siteID YEAR parameter value
#1 site1 2011 A 1
#2 site1 2012 A 2
#3 site1 2013 A 3
#4 site1 2014 A 4
#5 site1 2011 B 10
#6 site1 2012 B 20
#7 site1 2013 B 30
#8 site1 2014 B 40
#9 site1 2011 C 12
#10 site1 2012 C 14
#11 site1 2013 C 16
#12 site1 2014 C 18
#13 site1 2011 D 5
#14 site1 2012 D 10
#15 site1 2013 D 15
#16 site1 2014 D 20
#17 site2 2011 A 3
#18 site2 2012 A 6
#19 site2 2013 A 9
#20 site2 2014 A 12
ich eine Idee über die Entwicklung der Daten über das Grundstück bekommen kann unter
library(ggplot2)
ggplot(df, aes(x = YEAR, y= value, col= parameter))+
geom_line(aes(group=parameter))+
facet_wrap(~siteID)
wenn ich jedoch viele Websites (zB 70 Seiten), wird es mehr Zeit nehmen, um die Entwicklung visuell für alle Parameter zu analysieren. Gibt es trotzdem ich, wie sie der Trend in Daten in einer Tabelle erhalten könnte unter
siteID parameter trend
site1 A increasing
site1 B decreasing
site1 C no trend
site1 D increasing
.
site4 A decreasing
site4 B decreasing
site4 C increasing
site4 D decreasing
Irgendwelche Vorschläge wäre willkommen?
Lesen Sie mehr über 'lm()'. – zx8754
@ zx8754 Vielen Dank für Ihren Vorschlag – aelwan