This workshop is aimed at bringing together ~40 diverse speciation researchers:

- to collaborate on populating a database of published reproductive barriers based on a standardized RIO framework (https://ecoevorxiv.org/repository/view/10083/). This will involve working through papers during the workshop to extract RI measures and other metadata and entering them into a draft database (some preparatory work before the workshop may be requested to facilitate these steps during the workshop)

- to start working towards a manuscript using this database to answer an outstanding question in speciation

- to network, learn about reproductive isolation, and have fun!

If you are interested in applying to participate in the workshop, please fill out the form in the link below by ** May 20th **. Room & board will be covered by the organizers; all other travel costs are the responsibility of the attendee.

Application link:

https://docs.google.com/forms/d/e/1FAIpQLSenAMqSdZjeRmKEDtBNa5tpjPn7IukPyT5BzfZ4JMpIk2YOEw/viewform

The previous IOS workshop was held in Finland in 2023 (see more details at https://speciation-network.pages.ist.ac.at/workshops) and resulted in new interactions between speciation researchers and the successful publication of an Integration of Speciation research manuscript (https://academic.oup.com/evolinnean/article/3/1/kzae001/7609448).

We look forward to seeing you in Glasgow!

BioCoder features:

• Novel therapeutic discovery strategies
• Hardware such as low-cost lab equipment or diagnostics
• Open or low­-cost bioinformatics tools
• Engineered organisms for the production of small molecules, biologics, or other products
• Research projects at a community labspace or projects for science education or public engagement
• Hardware or software for lab automation
• Citizen science or DIY research projects
• Science policy
• Tools to increase reproducibility in research, or anything related

Address of the bookmark: https://www.oreilly.com/biocoder/

During the first phase of ENCODE, published in 2007, microarray-based technologies were used to detect regions associated with transcription factors, certain histone modifications and open chromatin within a pre-specified 1% of the human genome.

ENCODE’s second phase saw a switch to sequencing-based technologies, the addition of new assay types and the analysis of functional elements genome-wide, described in a collection of research articles in 2012.

The Encyclopedia paper of ENCODE 3, published in Nature, gives an overview of the various assays that were performed in human and mouse cell lines and tissues and describes a Registry of human and mouse candidate cis-regulatory elements (cCREs).

More at https://www.nature.com/immersive/d42859-020-00027-2/index.html

https://bioalgorithms.ucsd.edu/research4.html

For example ... contig-stats.pl is a Perl script that will automatically describe features of a sequence assembly.

http://milkweedgenome.org/?q=scripts

Address of the bookmark: http://milkweedgenome.org/?q=scripts

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

More at http://www.binc.co.in/College/Index_New.aspx

BINC notification http://www.binc.co.in/PdfDocuments/Notification.pdf

Few dates to remember:

Starting of online submission of application: March 9, 2015
Last date for submission of application: April 30,2015
Examination consists of two parts:
Part I (Paper I) : June 7, 2015 (10 AM-12 PM)
Part II ( Paper II & III) :June 28, 2015 (9 AM-12 PM & 2 PM-4 PM)