It is possible to create a vector by typing data directly into R using the combine function ‘c’

x

same as

x

creates the vector x with the numbers between 1 and 5.

You can see what is in an object at any time by typing its name;

x

will produce the output ‘[1] 1 2 3 4 5′

Note that names need to be quoted

daysofweek ← c(‘Monday’, ‘Tuesday’, ‘Wednesday’, ‘Thursday’, ‘Friday’);

Usually however you want to input from a file. We have touched on the ‘read.table’ function already.

mydata

Now mydata is a data frame with multiple vectors

each vector can be identified by the default syntax

#if any of these are typed it will print to screen

mydata$V1 mydata$V2 mydata$V3

By default the function assumes certain things from the file

• The file is a plain text file (there are function to read excel files: not covered here)
• columns are separated by any number of tabs or spaces
• there is the same number of data points in each column
• there is no header row (labels for the columns)
• there is no column with names for the rows** [I’ll explain].

If any of these are false, we need to tell that to the function

If it has a header column

mydata header=T also works

Note that there is a comma between different parts of the functions arguments

If there is one less column in the header row, then R assumes that the 1^st column of data after the header are the row names

Now the vectors (columns) are identified by their name

#if any of these are typed it will print to screen

mydata$A mydata$B mydata$C

# Summary about the whole data frame

summary(mydata)

# Summary information of column A

summary(mydata$A)

We can shortcut having to type the data frame each time by attaching it

attach(mydata)

# summary of column B as ‘mydata’ is attached

summary(B)

Two other important options for read.table

If is is separated only by tabs and has a header

mydata

Really useful if you have spaces in the contents of some columns, so R does not mess up reading the columns . However if the columns or of an uneven length it will tell you.

If you know that the file has uneven columns

mydata

This causes R to fill empty spaces in a columns with ‘NA’ .

The last two examples will still work with our file and give the same result as with only headers=T

Graphs

to get an idea of what R is capable of type

demo(graphics)

steps through the examples, and the code is printed to the screen

We will work with simpler examples that have immediate use to biologists.

Remember to get more information about the options to a function type ‘?function’

Histogram of A

hist(mydata$A)

If there was more data we could increase the number of vertical columns with the option, breaks=50 (or another relevant number).

boxplot(mydata)

We can get rid of the need to type the data frame each time by using the attach function

# if not already done so

attach(mydata)

boxplot(mydata$A, mydata$B, name=c(“Value A”, “Value B”) , ylab=“Count of Something”)

same as

boxplot(A, B, name=c(“Value A”, “Value B”) , ylab=“Count of Something”)

Scatter plot

# if not already done so

attach(mydata)

plot(A,B) # or plot(mydata$A, mydata$B)

SAVING an image

Windows users (Rgui) RIGHT click on image and select which you want.

These instructions work for everyone.

You need to create a new device of the type of file you need, then send the data to that device

to save as a png file (easy to load into the likes of powerpoint, also great for web applications.

png(‘filename’)

boxplot(A, B, name=c(“Value A”, “Value B”) , ylab=“Count of Something”)

or to save as a pdf

pdf(‘filename’)

boxplot(A, B, name=c(“Value A”, “Value B”) , ylab=“Count of Something”)

Note

• Nothing will appear on screen, the output is going to the file
• Also it may not be saved immediately but will once the device (or R) is turned quit.

To quit R type

q() # If you save your session, next time you start R, you will have your data preloaded.

Or if you want to remain in R

dev.off() #turns of the png (or pdf etc) device, thus forces the data to save

PRINCIPAL SCIENTIST Pay Band: Minimum pay of `43,000 in the PB-4 of `37400-67000/- + RGP of `10,000/-.

Age: The candidates must not have attained the age of 52 years as on 24.03.2015. There shall be no age limit for the Council’s employees.

ICAR-NATIONAL INSTITUTE OF BIOTIC STRESS MANAGEMENT, RAIPUR (CHHATTISGARH)

57. Principal Scientist (Agricultural Entomology) (Two post)

Qualifications Essential:

(i) Doctoral degree in Agricultural Entomology including relevant basic sciences.

(ii) 10 years experience in the relevant subject out of which at least 8 years should be as Scientist/ Lecturer/Extension Specialist or in an equivalent position in the Pay Band- 3 of `15600-39100 with Grade Pay of `5400/`6000/`7000/`8000 and 2 years as a Senior Scientist or in an equivalent position in the Pay Band- 4 of ` 37400-67000 with Grade Pay of ` 8700/ ` 9000.

(iii) The candidate should have made contribution to research/teaching/extension education as evidenced by published work/innovations and impact.

Desirable:

(i) Experience of using frontiers research tools in management of insect pests of crop plants.

(ii) Evidence of contributions to relevant field through publications/ patents/citation index to suggest a vision/perspective in biotic stress research.

61. Principal Scientist (Bioinformatics) (One post)

Qualifications Essential:

(i) Doctoral degree in Bioinformatics including relevant basic sciences. (ii) & (iii) As in item no. 57 above.

Desirable:

(i) Experience of using bioinformatics for advancement of knowledge and for research on biotic stress management.

(ii) Evidence of contributions to relevant field through publications/patents/citation index to suggest a vision/perspective in biotic stress research.

http://asrb.org.in/administrator/uploads_dir/1424859407english.pdf

ICAR-NPTC: Fibre development in flax/linseed.

(Job # 1) Research Associate (one) (Bioinformatics)

Rs.24000+ 30% HRA) for Ph.D. and for M. Sc Rs.23000/‐ (+ 30% HRA)

Ph.D. Degree in Bioinformatics/Molecular Biology/Biotechnology/ Genetics/allied sciences; or M. Sc in Bioinformatics/ Biotechnology/Life Sciences/ allied sciences with 1st division or 60% marks or equivalent overall grade point average with at least two years of research experience as evidenced from Fellowship/ Associate ship. 2 years research experience in bioinformatic data analysis/molecular biology techniques, and high throughput DNA/RNA sequencing, and transcriptome data analysis. Research paper with IF>1 will be desirable

ICAR-NPTC: Shade avoidance/low-light tolerance in rice.

General Terms & Conditions applicable to all the positions:
Age Limit: 35 years max. (5 years relaxation for SC/ST/OBC and woman candidates as per ICAR rules).
The positions are purely temporary, on a contractual basis and are initially offered for one year.
Originals must be shown for verification. 7. Research experience (Experience certificate from previous employer to be attached): I hereby declare that the information provided above is true to the best of my knowledge. Date: Signature

Advertisement:

www.nrcpb.org/sites/default/files/ICAR-NPTC%20DBT%20RA%20SRF%20interview%2024th%20March.pdf

A postdoctoral training position is currently available in the Computational and Statistical Genomics Branch (CSGB) of the National Human Genome Research Institute (NHGRI). The position is located in the laboratory of Andy Baxevanis, Ph.D., whose research group uses comparative genomics approaches to better-understand the molecular innovations that drove the surge of diversity in early animal evolution. The overarching theme of Dr. Baxevanis’ research program is focused on how non-traditional animal models convey critical insights into human disease research.

Candidates should have or be close to obtaining a Ph.D. or equivalent degree in bioinformatics, computational biology, computer science, molecular biology, or a closely related field. Candidates with a background in evolutionary biology are particularly encouraged to apply. Programming skills and experience in the application of computational methods to genomic data are highly desirable. Applicants must possess good communication skills and be fluent in both spoken and written English. The ability to learn how to use new software and quickly become expert in its use, critical thinking, problem-solving abilities, and the ability to work semi-independently are required.
How to Apply

Interested applicants should submit a curriculum vitae, a detailed letter of interest, and the names of three potential referees to Dr. Baxevanis at andy@mail.nih.gov.
About Our Organization

The NIH Intramural Research Program is on the Bethesda, Maryland campus and offers a wide array of training opportunities for scientists early in their careers. The funding for this position is stable and offers the trainee wide latitude in the design and pursuit of their research project. The successful candidate will have access to NHGRI’s established and robust bioinformatics infrastructure, as well as resources made available through NIH’s Center for Information Technology (CIT) and the National Center for Biotechnology Information (NCBI).

For more information on CSGB and NHGRI’s Intramural Research Program, please see http://genome.gov/DIR/.

http://www.thestreet.com/story/13151062/1/affymetrix-acquires-assets-of-eureka-genomics-corporation-to-provide-high-throughput-and-economical-crop-and-animal-genotyping.html

Our main research themes are:

Domain annotation and metagenomics
Transcriptomics and sequence analysis
Protein evolution and interactions
Protein conformational dynamics

More at http://www.lcqb.upmc.fr/AnalGenom/home.html

Bioinformatics
Cancer
Epigenetics

More at http://liulab.dfci.harvard.edu/

#Find all occurrences of a pattern in a string.
#Given: Strings Pattern and Genome.
#Return: All starting positions in Genome where Pattern appears as a substring. Use 0-based indexing.

use strict;
use warnings;

my $string="GATATATGCATATACTT";
my $subStr="ATAT";
my $kmer=length($subStr);

kmerMatch ($string, $subStr, $kmer);

sub kmerMatch { #Check the exact matching kmers with sliding window
my ($string, $myStr, $kmer)=@_;
for (my $aa=0; $aa<=(length($string)-$kmer); $aa++) {
    my $myWin=substr  $string, $aa,$kmer;
    if ($myWin eq $myStr) {
        #print "$myWin eq $myStr\n";
        print $aa;
    }
}
}

• Expansion of Protists and Fungi with hundreds of annotated genomes
• Variation data for bread wheat, rice, Aedes aegypti, and Ixodes scapularis
• Whole genome alignments for O. longistaminata and T. cacao
• Non-coding RNA gene models in Bacteria
• New assembly of tomato (version 2.50)
• Full support for UCSC Track Hub format for hosting your own data in Ensembl

More at http://www.ensembl.info/blog/2015/06/16/ensembl-genomes-release-27-is-out/