################################################### ### 块1号:setup1 ################################################### 图书馆(“RNAither ") ################################################### ### 块2号:################################################### # 基因名字plateLayout1 < - c(“test1”、“空”、“test3”,“test4”、“test5”,“test6”、“test7”、“空”、“test9”,“test10”、“test11”、“test12”)plateLayout2 < - c(“test1”、“test2”,“test3”、“test4”,“test5”,“test6”、“test7”、“test8”、“test9”、“test10”、“test11”、“test12”)plateLayout < - cbind (plateLayout1 plateLayout2) emptyWells < -列表(c (2, 8), NA_integer_) #第一板有两个空井在位置2和8,#第二板没有任何空井poorWells <-NA_integer_ #没有质量差的井controlCoordsOutput <- list(list(NA_integer_, NA_integer_), list(NA_integer_, c(9,10))) #第一个板没有任何控制siRNAs, #第二板在位置9和10有两个负控制backgroundValOutput<-NA_integer_ #没有背景信号强度可用sigPlate1<-c(2578, NA_integer_, 3784, 3784, 2578, 5555, 5555, NA_integer_, 8154, 2578, 3784, 2578) sigPlate2<-c(8154, NA_integer_)3784, 5555, 3784,11969,2578,55555,2578) #板上的信号强度meanSignalOutput<-list(sigPlate1, sigPlate2) sdmesignal <- na_integer# no standard deviation available objnumOutput<- na_integer# no cell count available cellnumOutput<-NA_integer_ generateDatasetFile("第一个测试屏幕","RNAi in virus-infected cells", NA_character_, "testscreen_output.txt", plateLayout, plateLayout, 3,4,1, emptyWells, poorWells, controlCoordsOutput,backgroundValOutput, meanSignalOutput, sdmeanssignal, objnumOutput, cellnumOutput) #将数据集加载到R: header<-readLines("testscreen_output.txt",3) dataset<-read.table("testscreen_output.txt", skip=3, colClasses=c(NA, NA, NA, NA, "factor", NA, NA, NA, NA, NA, NA, NA, NA), stringsAsFactors=FALSE) ################################################### ### chunk number 3:################################################### 数据(exampleDataset包=“RNAither”)doubledataset < - joinDatasets(列表(数据集,数据集 )) ################################################### ### 块数量4:################################################### data(exampleHeader, package="RNAither") data(exampleDataset, package="RNAither") saveDataset(header, dataset, "save_testfile1.txt") header[[1]] <- "external_experiment_name,Test screen" header[[2]] <- "comments,包含两次screen Nb. 1" joinDatasetFiles(list("save_testfile1.txt", "save_testfile1.txt"), 3, header, "concatenated_testfile.txt")