我自己的電腦是win10系統1803版本，如果是其他版本系統執行的話有可能會出現系統庫缺失等問題，所以最好自己執行一遍。

step2:開啟終端，
win+R鍵，輸入cmd，確定，進入終端

step3：開啟R程式（下面是預設安裝位置的進入方法）
在終端輸入
cd C:\Program Files\R\R-3.5.1
cd bin
R.exe

step4:執行演算法指令碼
把NIPT演算法壓縮包解壓到c盤根目錄下
在終端輸入
setwd(“c:/”)
source(“test.R”)
等一會，程式執行完，就會在c盤根目錄生成三個檔案，就是所求的結果
example.bf.txt
example.zscore.txt
example.hh.txt

題目2：
step1: 使用管理員模式開啟R程式（預設位置是C:\Program Files\R\R-3.5.1\bin\R.exe 右鍵點選R.exe，以管理員模式執行）

step2:在C盤根目錄新建一個txt檔案,改名為test2.R,用記事本開啟，寫入以下內容並儲存。

library(reshape2);
library(ggplot2);
#資料讀入
data <- read.table("example.bf.txt", sep=";", header=T, row.names=1, quote="");
data<-as.matrix(data);
data<-as.data
 
.frame(data);
#is.matrix(data);

# 增加行的名字
rownames(data)<- c("chr1","chr2","chr3","chr4","chr5","chr6","chr7","chr8","chr9","chr10","chr11","chr12","chr13","chr14","chr15","chr16","chr17","chr18","chr19","chr20","chr21","chr22")

# 增加列的名字
colnames(data) <- paste("ID1","ID2","ID3","ID4","ID5","ID6"
 
,"ID7","ID8","ID9","ID10","ID11","ID12","ID13","ID14","ID15","ID16","ID17","ID18","ID19","ID20")

#增加一列ID列，儲存行名字
data$ID <- rownames(data)

# melt：把正常矩陣轉換為長表格模式的函式。工作原理是把全部的非id列的數值列轉為1列，命名為value；所有字元列轉為variable列。
# id.vars 列用於指定哪些列為id列；這些列不會被merge，會保留為完整一列。
data_m <- melt(data, id.vars=c("ID"))
head(data_m)

p <- ggplot(data_m, aes(x=variable,y=ID)) + xlab("samples") + theme_bw() + theme(panel.grid.major = element_blank()) + theme(legend.key=element_blank())  + theme(axis.text.x=element_text(angle=45,hjust=1, vjust=1)) + theme(legend.position="top") +  geom_tile(aes(fill=value)) + scale_fill_gradient(low = "white", high = "red")

# 可以跟輸出檔案不同的字尾，以獲得不同的輸出格式
# colormode支援srgb （螢幕）和cmyk （列印，部分雜誌需要，看上去有點褪色的感覺）格式
ggsave(p, filename="heatmap.pdf", width=10, height=15, units=c("cm"),colormodel="srgb")

step3： 在開啟R.exe後的終端上輸入以下內容：

setwd("c:/")
source("test2.R")

等程式執行完成後就可以在c盤根目錄下看到輸出的圖片，就是所求bf的結果

step4：

將第二步驟中程式碼段第四行的

data <- read.table("example.bf.txt", sep=";", header=T, row.names=1, quote="");

中的example.bf.txt分別改成
example.zscore.txt
example.hh.txt
，重複step2，step3就可以得到zscore,hh的結果

題目3：

#include<stdio.h>
#include<math.h>

float calculate(float bf,float pror)
{
    float poor=bf*pror,ppv=poor/(poor+1);

    return ppv;
}
int main()
{
    float bf,pror,ppv;

    //從終端輸入bf和pror值
    printf("input the value of BF:\n");
    scanf("%f",&bf);
    printf("input the value of pror:\n");
    scanf("%f",&pror);


    ppv = calculate(bf,pror);
    printf("The value of ppv is %f\n",ppv);
    return 0;

}