Subset Datenrahmen und zeichnen alle Untermengen mit einer Schleife [R]

Question

Im arbeiten mit einem Datenrahmen mit 8 nützliche Variablen, die Idee des Codes ist es, 4 Variablen (3 auf y-Achse und eine gemeinsame X-Achse) zu plotten. Der Datenrahmen sieht wie folgt aus:

It has like 6500 rows

ich die data.frame aus der Datei Spalte der Teilmenge werden soll, und dann plotten LogP als x-Achse und Temperatur, relativer Luftfeuchtigkeit und Ozon als y-Achse. Ich habe versucht, Teilmenge innerhalb der Plot-Funktion, aber ging nicht gut. Ich habe diesen Code für das Grundstück mit einem der ursprünglichen Dateien, aber keine Ahnung, wie die Teilmenge

> plot(DataOzono$LogP, DataOzono$Temperature, axes= F,type="l",col="red", ylab = NULL, xlab = 'LogP',xaxt="n",yaxt="n") 
axis(2,ylim(c(min(DataOzono$Temperature),max(DataOzono$Temperature)), layout.widths(2))) 
mtext(text = 'T',line = 2,side = 2) 
par(new=TRUE) 
plot(DataOzono$LogP, DataOzono$RH,type="l",col="blue",xaxt="n",yaxt="n",xlab="",ylab="") 
axis(4) 
mtext("RH",side=4,line=2) 
par(new=TRUE) 
plot(DataOzono$LogP, DataOzono$Ozone,type="l",col="green",xaxt="n",yaxt="n",xlab="",ylab="") 
mtext("O3",side=5,line=3) 
axis(2, line = 4) 

jede Beratung wird es sein, umfasst sehr hilfreich.

Answer 1

So können Sie die Diagramme in einem loop plotten. In dem Beispiel, das Sie angegeben haben, haben wir nur eine Dateinummer. Es sollte jedoch für jede Zahl in der Dateispalte ein Diagramm erstellt werden. Unter Windows können Sie savePlot auf Ihrem Laufwerk speichern. Ich habe Ihr Beispiel vereinfacht, weil ich Fehler bekommen habe.

DataOzono <- read.table(text="pressure height Temperature RH Ozone file LogP 
753.6 2541 16.8 76 0 80131 0.3475673 
748.0 2604 17.7 32 0 80131 0.347959 
743.5 2656 15.9 38 0 80131 0.3482766 
739.8 2697 15.4 39 0 80131 0.3485396 
736.6 2734 15.0 41 0 80131 0.3487685 
731.8 2790 14.5 42 0 80131 0.3491142", header=TRUE, stringsAsFactors=FALSE) 

original_par <- par() 
par(mar=c(5.1, 8.1, 4.1, 3.1)) 

for (i in unique(DataOzono$file)){ 
DataOzono_subset <- DataOzono[DataOzono$file==i,] #keep only rows for that file number 

plot(DataOzono_subset$LogP, DataOzono_subset$Temperature, axes= F,type="l",col="red", ylab = "", xlab = 'LogP',xaxt="n",yaxt="n") 
axis(2,col="red",col.axis="red") 
mtext(text = 'T',line = 2,side = 2,col="red",col.lab="red") 
par(new=TRUE) 
plot(DataOzono_subset$LogP, DataOzono_subset$RH,type="l",col="blue",xaxt="n",yaxt="n",xlab="",ylab="") 
axis(4,col="blue",col.axis="blue") 
mtext("RH",side=4,line=2,col="blue",col.lab="blue") 
par(new=TRUE) 
plot(DataOzono_subset$LogP, DataOzono_subset$Ozone,type="l",col="darkgreen",xaxt="n",yaxt="n",xlab="",ylab="") 
mtext("O3",side=2,line=6,,col="darkgreen",col.lab="darkgreen") 
axis(2, line = 4,col="darkgreen",col.axis="darkgreen") 

savePlot(filename=paste0("c:/temp/",i,".png"),type="png") 
} 

par() <- original_par #restore par to initial value.