Skip to main content

Java Template code

package com.cbrc.pipeline2;

import java.io.BufferedReader;
import java.io.BufferedWriter;
import java.io.DataInputStream;
import java.io.FileInputStream;
import java.io.FileWriter;
import java.io.InputStreamReader;

public class Test {


    String foldIn;
    String fnmFasta;

    String foldOut;
    String fnmOut;
   
    void init(String rootIn, String rootOut, String fnmInFasta, String fOut) {
       
        this.foldIn = rootIn;
        this.fnmFasta = this.foldIn + fnmInFasta;
   
        this.foldOut = rootOut;
        this.foldOut = this.foldOut + fOut;

    }


    void loadFasta(){

        try {
            FileInputStream fstream = new FileInputStream(fnmFasta);
            DataInputStream in = new DataInputStream(fstream);
            BufferedReader br = new BufferedReader(new InputStreamReader(in));

            String befTrim , strLine;
            while ((befTrim = br.readLine()) != null) {
                strLine = befTrim.trim();


            }

            br.close();
            in.close();
            fstream.close();


        } catch (Exception e) {

            e.printStackTrace();
        }


    }



    void writeSelectedRegion(){


        try {

            BufferedWriter bwr = new BufferedWriter(new FileWriter(this.fnmOut));
            StringBuffer buf = new StringBuffer();




            bwr.write(buf+"");

            bwr.close();
        } catch (Exception e) {
            e.printStackTrace();
        }


    }

    void doProcessing() {
       
        loadFasta();
        writeSelectedRegion();
    }

    public static void main(String[] args) {

        Test obj = new Test();

        obj.init(args[0], args[1] , args[2] , args[3]);
        //        obj.init("./", "test.fa", "20", "72" , "seq.prf");

        obj.doProcessing();


    }

}

Comments

Popular posts from this blog

MATLAB cross validation

// use built-in function samplesize = size( matrix , 1); c = cvpartition(samplesize,  'kfold' , k); % return the indexes on each fold ///// output in matlab console K-fold cross validation partition              N: 10    NumTestSets: 4      TrainSize: 8  7  7  8       TestSize: 2  3  3  2 ////////////////////// for i=1 : k    trainIdxs = find(training(c,i) ); %training(c,i);  // 1 means in train , 0 means in test    testInxs  = find(test(c,i)       ); % test(c,i);       // 1 means in test , 0 means in train    trainMatrix = matrix (  matrix(trainIdxs ), : );    testMatrix  = matrix (  matrix(testIdxs  ), : ); end //// now calculate performance %%  calculate performance of a partition     selectedKfoldSen=[];selectedKfoldSpe=[];selectedKfoldAcc=[];     indexSen=1;indexSpe=1;indexAcc=1;     if ( kfold == (P+N) )% leave one out         sensitivity = sum(cvtp) /( sum(cvtp) + sum(cvfn) )         specificity = sum(cvtn) /( sum(cvfp) + sum(cvtn) )         acc

R tutorial

Install R in linux ============ In CRAN home page, the latest version is not available. So, in fedora, Open the terminal yum list R  --> To check the latest available version of r yum install R --> install R version yum update R --> update current version to latest one 0 find help ============ ?exact topic name (  i.e.   ?mean ) 0.0 INSTALL 3rd party package  ==================== install.packages('mvtnorm' , dependencies = TRUE , lib='/home/alamt/myRlibrary/')   #  install new package BED file parsing (Always use read.delim it is the best) library(MASS) #library(ggplot2) dirRoot="D:/research/F5shortRNA/TestRIKEN/Rscripts/" dirData="D:/research/F5shortRNA/TestRIKEN/" setwd(dirRoot) getwd() myBed="test.bed" fnmBed=paste(dirData, myBed, sep="") # ccdsHh19.bed   tmp.bed ## Read bed use read.delim - it is the  best mybed=read.delim(fnmBed, header = FALSE, sep = "\t", quote = &q