tsur<-function(fokszam){
  
  library(stats)
  library(graphics)
  library(utils)
  
  xtengely<-seq(-3,3,0.01)
  
  for(i in 1:fokszam){
    ytengely<-xtengely
    ytengelynorm<-xtengely
    
    for(j in 1:length(xtengely)){
      ytengely[j]<-dt(xtengely[j],i)
      ytengelynorm[j]<-dnorm(xtengely[j],0,1)
    }
    
    seged<-paste("tsuruseg",toString(i),".jpg",sep="")
    
    jpeg(seged)
    
    plot(xtengely,ytengely,ylim=c(0,0.4))
    lines(xtengely,ytengelynorm,col="GREEN")
    dev.off()
  }  
  
}

tsur(8)


koleszterinregi<-read.table("cholestr.txt",header=TRUE)

View(koleszterinregi)

colnames(koleszterinregi)

str(koleszterinregi)

t.test(koleszterinregi$Day2,mu=280)

t.test(koleszterinregi$Day2,mu=280, alternative="less")

?t.test

koleszterin<-read.table("cholestg.txt",header=TRUE)

head(koleszterin)


koleszterin2<-subset(koleszterin,day==2)
koleszterin2

koleszterin2<-koleszterin2[c("patient","group","cholest")]
koleszterin2

colnames(koleszterin2)<-c("patient","group","cholest2")



koleszterin14<-subset(koleszterin,day==14)

koleszterin14<-koleszterin14[c("patient","group","cholest")]


colnames(koleszterin14)<-c("patient","group","cholest14")

koleszterin14

?merge

kol<-merge(koleszterin2,koleszterin14,by=c("patient","group"))

View(kol)



koleszterinNA<-subset(koleszterin,is.na(day))

koleszterinNA<-koleszterinNA[c(1,2,4)]

koleszterinNA

colnames(koleszterinNA)<-c("patient","group","cholest2")

koleszterinNA$cholest14<-koleszterinNA$cholest2

koleszterinNA

kolegyes<-rbind(kol,koleszterinNA)

View(kolegyes)

save(kolegyes,file="kolegyes.RData")

load("kolegyes.RData")

View(kolegyes)