library(samr)
    library(RColorBrewer)
    library(gplots) #���������heatmap.2�����������ʹ�ý���ɫ��
    library(survival)
    library(randomForest)
    library(e1071)
    library(pROC)

    setwd(“D:\����\GSE35713”)
    linchuang <- read.csv(“GSE35713.csv”,stringsAsFactors = F)
    data <- read.table(“GSE35713_lnc.txt”)
    data <- data[,linchuang[,1]]
    w<-rep(2,20) #184��type 1 diabetes�ĸ�ֵ1
    e<-rep(1,42) #14��healthy�ĸ�ֵ0
    data <- data[,linchuang[,1]]
    o<-rownames(data)
    data <- as.matrix(data)
    rownames(data) <- o
    set.seed(1) #�趨�������
    q<-SAM(data[,c(1:62)],c(w,e),resp.type=”Two class unpaired”, genenames=o) #x�Ǿ���ͨ��cbind�ϲ��˸����Ͳ����������������԰���������˳�������ľ���
    results.sam <- rbind(q差异表达 - 图1genes.up,q差异表达 - 图2genes.lo)
    genename <- as.character(results.sam[,1])
    results.sam <- data.frame(results.sam)
    rownames(results.sam) <- genename
    for(i in 1:ncol(results.sam)) {
    results.sam[,i] <- as.numeric(as.character(results.sam[,i]))
    } #��results.sam��ÿһ�н���numeric��

    results.sam <- results.sam[(results.sam$q.value…)<1,-1] #qֵȡ0.01�����DZ����е��ǰٷ���
    results.sam