cancer<-read.table("group2.txt", header=TRUE)
noncancer<-read.table("group1.txt", header=TRUE)

#meanc<-mean(cancer)
#meannc<-mean(noncancer)
#varc<-var(cancer)
#varnc<-var(noncancer)
#Sx1x2<-sqrt((varc+varnc)/2)
#n<-4
#t<-(noncancer-cancer)/(Sx1x2*(sqrt(2/n)))
#df<-(2*n)-2
#pt(t,df)*2


answer = matrix(,1733,1);
alpha = 0.05;


#for(i=1:1733)
#{
#  meanc = mean(cancer[i,]);
#  meannc = mean(noncancer[i,]);
#  
#  varc = var(cancer[i,]);
#  varnc = var(noncancer[i,]); 
#
#  Sx1x2 = sqrt((varc+varnc)/2);
#  n = 4;
#  t = (noncancer-cancer)/(Sx1x2*(sqrt(2/n)));
#  df = 2*n - 2;
#
#
#  answer[i,1] = pt(t,df)*2;
#
#  if (answer[i,1] = alpha)
#  {
#    sing.
#  }
#  else
#  {
#    nonsing.
#  }
#}

for(i in 1:1733)
{
  meanc = mean(cancer[i,]);
  meannc = mean(noncancer[i,]);
  
  varc = var(cancer[i,]);
  varnc = var(noncancer[i,]); 

  Sx1x2 = sqrt((varc+varnc)/2);
  n = 4;
  t = (noncancer-cancer)/(Sx1x2*(sqrt(2/n)));
  df = 2*n - 2;

  answer[i,1] = pt(t,df)*2 ;

  if (answer[i,1] < alpha)
  {
    banswer[i,1] = 1;
  }
  else
  {
    banswer[i,1] = 0;
  }

  #ifelse(answer[i,1] < alpha, banswer[i,1] = 1, banswer[i,1] = 0);



}