У меня есть полурасплавленный фрейм данных, который выглядит так:

head(final_melt)

   Group       Source variable   value
 Control Whole Kidney     MZF1 0.23879
 Control Whole Kidney     MZF1 0.49381
 Control Whole Kidney     MZF1 0.40827
 Control Whole Kidney     MZF1 0.55548
 Control Whole Kidney     MZF1 0.34664
 Control Whole Kidney     MZF1 0.68102

Группа имеет два уровня (контроль и заболевание), источник имеет 4 уровня (цельная почка, клубочек, тубулоинтерстиций и HK-2 + TGF-B). Переменная также имеет четыре уровня (TFAP2A, MZF1, YY1, SP1). Я хотел бы сделать что-то вроде следующего в цикле

d = subset(final_melt, final_melt$Source=="Whole Kidney")
qplot(data=d, Group, value, facets=.~variable, geom="boxplot")
pdf("Whole Kideny.pdf")
dev.off()

Хотя я знаю, что могу сказатьfacets=Source~variableотдельные участки оказываются слишком маленькими, чтобы быть информативными. Поэтому мне нужно наносить один уровень фактора источника за раз.
Проблема в том, что я не могу даже получитьpdf() функция для работы. Он создает файл с правильным именем, но когда я пытаюсь открыть его, Adobe говорит, что при открытии файла произошла ошибка, и он уже открыт в другом приложении (почему я добавилdev.off(), но это, похоже, ничего не делает).

Любая помощь приветствуется.
Ура, Дэви.