This is a follow-up to Automate Query.
Dr. Undómiel is really testing your kindness now. She’s seen that all of the average hind foot length and average weights are pretty similar within species, but she wants to support the idea with some statistics.
Write a function that selects the sex
, hindfoot_length
and weight
data
for a given species.
Write a function that uses the data from the first function and returns the results of a
t.test()
of each measurement factored by sex. (HINT: A t.test()
takes two primary
arguments: a vector or data frame of values for group a and the same class of data for group
b. The arguments can also be expressed in formula form similar to aov()
. Be sure to
assign the results to an object (test <- t.test()
), so you can extract the *p-value from the
object (test$p.value
) to include in your results.*)
Loop through the rodent species at Portal and generate an output data frame
with the species_id
, p_hindfoot_length
, p_weight
. (HINT: You are
likely to run into some issues of insufficient data that you should skip in your loop.)