Ik ben nieuw voor r en heb geen programmering gedaan voor …
Wanneer ik probeer een vakje met standaardfoutstaven te maken, krijg ik het foutbericht dat in de titel wordt vermeld.
Ik heb een script gebruikt dat ik heb gevonden op R-kookboek die ik een beetje heb aangepast:
ggplot(GVW, aes(x="variable",y="value",fill="Genotype")) +
geom_bar(position=position_dodge(),stat="identity",colour="black", size=.3)+
geom_errorbar(data=GVW[1:64,3],aes(ymin=value-seSKO, ymax=value+seSKO), size=.3, width=.2, position=position_dodge(.9))+
geom_errorbar(data=GVW[65:131,3],aes(ymin=value-seSWT, ymax=value+seSWT), size=.3, width=.2, position=position_dodge(.9))+
geom_errorbar(data=GVW[132:195,3],aes(ymin=value-seEKO, ymax=value+seEKO), size=.3, width=.2, position=position_dodge(.9))+
geom_errorbar(data=GVW[196:262,3],aes(ymin=value-seEWT, ymax=value+seEWT), size=.3, width=.2, position=position_dodge(.9))+
xlab("Time")+
ylab("Weight [g]")+
scale_fill_hue(name="Genotype", breaks=c("KO", "WT"), labels=c("Knock-out", "Wild type"))+
ggtitle("Effect of genotype on weight-gain")+
scale_y_continuous(breaks=0:20*4) +
theme_bw()
Data<- data.frame(
Genotype<- sample(c("KO","WT"), 262, replace=T),
variable<- sample(c("Start","End"), 262, replace=T),
value<- runif(262,20,40)
)
names(Data)[1] <- "Genotype"
names(Data)[2] <- "variable"
names(Data)[3] <- "value"
Antwoord 1, Autoriteit 100%
De fout gebeurt omdat u probeert een numerieke vector in kaart te brengen voor data
in geom_errorbar
: GVW[1:64,3]
. ggplot
Werkt alleen met data.frame
.
Over het algemeen moet u geen subsets maken binnen ggplot
-aanroepen. U doet dit omdat uw standaardfouten in vier afzonderlijke objecten worden opgeslagen. Voeg ze toe aan je originele data.frame
en je kunt alles in één gesprek plotten.
Hier met een dplyr
-oplossing om de gegevens samen te vatten en vooraf de standaardfout te berekenen.
library(dplyr)
d <- GVW %>% group_by(Genotype,variable) %>%
summarise(mean = mean(value),se = sd(value) / sqrt(n()))
ggplot(d, aes(x = variable, y = mean, fill = Genotype)) +
geom_bar(position = position_dodge(), stat = "identity",
colour="black", size=.3) +
geom_errorbar(aes(ymin = mean - se, ymax = mean + se),
size=.3, width=.2, position=position_dodge(.9)) +
xlab("Time") +
ylab("Weight [g]") +
scale_fill_hue(name = "Genotype", breaks = c("KO", "WT"),
labels = c("Knock-out", "Wild type")) +
ggtitle("Effect of genotype on weight-gain") +
scale_y_continuous(breaks = 0:20*4) +
theme_bw()