Lab 16 my_result.csv

Blinda Sui (PID: A17117043)

b <- read.table("my_results.tsv")


colnames(b) <- c("qseqid", "sseqid", "pident", "length", "mismatch", "gapopen", "qstart", "qend", "sstart", "send", "evalue", "bitscore")

library(ggplot2)
ggplot(b, aes(pident, bitscore)) + geom_point(alpha=0.1) 

ggplot(b, aes((b$pident * (b$qend - b$qstart)), bitscore)) + geom_point(alpha=0.1) + geom_smooth()

Warning: Use of `b$pident` is discouraged.
ℹ Use `pident` instead.

Warning: Use of `b$qend` is discouraged.
ℹ Use `qend` instead.

Warning: Use of `b$qstart` is discouraged.
ℹ Use `qstart` instead.

Warning: Use of `b$pident` is discouraged.
ℹ Use `pident` instead.

Warning: Use of `b$qend` is discouraged.
ℹ Use `qend` instead.

Warning: Use of `b$qstart` is discouraged.
ℹ Use `qstart` instead.

`geom_smooth()` using method = 'gam' and formula = 'y ~ s(x, bs = "cs")'