ggplot2の物理実験からデータを整理する

Question

現在、簡単な電流 - 電圧実験の結果を視覚化するためにggplot2を使用しようとしています。私はもちろん、1セットのデータに対して良い結果を達成することができました。

しかし、私は（最小限のコードを参照してください）再帰的に以下の組織を得るためにRでI入力電流 - 電圧データセットの数を、持っている：

data.frame(cbind(batch(string list), sample(string list), dataset(data.frame list))) 

編集：私のデータはテキストファイル名で保存されていますbatchname_samplenumber.txtには電圧と電流の列があります。

require(plyr) 
require(ggplot2) 

#VARIABLES 
regex <- "([[:alnum:]_]+).([[:alpha:]]+)" 
regex2 <- "G5_([[:alnum:]]+)_([[:alnum:]]+).([[:alpha:]]+)" 

#FUNCTIONS 
getJ <- function(list, k) llply(list, function(i) llply(i, function(i, indix) getElement(i,indix), indix = k)) 

#FILES 
files <- list.files("Data/",full.names= T) 

#NAMES FOR FILES 
paths <- llply(llply(files, basename),function(i) regmatches(i,regexec(regex,i))) 
paths2 <- llply(llply(files, basename),function(i) regmatches(i,regexec(regex2,i))) 
names <- llply(llply(getJ(paths, 2)),unlist) 
batches <- llply(llply(getJ(paths2, 2)),unlist) 
samples <- llply(llply(getJ(paths2, 3)),unlist) 

#SETS OF DATA, NAMED 
sets <- llply(files,function(i) read.table(i,skip = 0, header = F)) 
names(sets) <- names 
for (i in as.list(names)) names(sets[[i]]) <- c("voltage","current") 

df<-data.frame(cbind(batches,samples,sets))  

、最小限のデータを介して生成することができます：私はそれらをインポートするために使用するコードです

require(plyr) 

batch <- list("A","A","B","B") 
sample <- list(1,2,1,2) 
set <- list(data.frame(voltage = runif(10), current = runif(10)),data.frame(voltage = runif(10), current = runif(10)),data.frame(voltage = runif(10), current = runif(10)),data.frame(voltage = runif(10), current = runif(10))) 

df<-data.frame(cbind(batch,sample,set)) 

私の質問です：コードを使用してプロットすることであるとして、それの可能なデータを使用することです次のような（動作しません）？

ggplot(data, aes(x = dataset$current, y = dataset$voltage, colour = sample)) + facet_wrap(~batch) 

より一般的なバージョンであろうが：（ダイヤモンド、車など）の統計データを離散とは対照的に、生の物理的データをhandelingのggplot2ことができますか？

Answer 1

、私は次のような戦略を示唆している：

read_custom <- function(f, ...) { 
d <- read.table(f, ...) 
names(d) <- c("V", "I") 
## extract sample and batch from the base filename 
ids <- strsplit(gsub(".txt", "", f), "_") 
d$batch <- ids[[1]][1] 
d$sample <- ids[[1]][2] 
d 
} 

## list files to read 
files <- list.files(pattern=".txt") 
## read them all in a single data.frame 
m <- ldply(files, read_custom) 

Answer 2

サンプル名がデータセットに対してどのように定義されているかは不明です。 ggplot2の一般的な考え方は、すべてのデータを溶けた（長い形式の）data.frame形式でグループ化することです。

新しく定義された問題（「batchname_samplenumber.txt」という名前の2つの列のファイル）で

library(ggplot2) 
library(plyr) 
library(reshape2) 

l1 <- list(batch="b1", sample=paste("s", 1:4, sep=""), 
      dataset=data.frame(current=rnorm(10*4), voltage=rnorm(10*4))) 
l2 <- list(batch="b2", sample=paste("s", 1:4, sep=""), 
      dataset=data.frame(current=rnorm(10*4), voltage=rnorm(10*4))) 
l3 <- list(batch="b3", sample=paste("s", 1:4, sep=""), 
      dataset=data.frame(current=rnorm(10*4), voltage=rnorm(10*4))) 

list_to_df <- function(l, n=10){ 

    m <- l[["dataset"]] 
    m$batch <- l[["batch"]] 
    m$sample <- rep(l[["sample"]], each=n) 
    m 
} 

## list_to_df(l1) 

m <- ldply(list(l1, l2, l3), list_to_df) 

ggplot(m) + facet_wrap(~batch)+ 
    geom_path(aes(current, voltage, colour=sample))