Ask History Expert

Problem 1: From ourMoodel site, Run the lines provided at the very bottom of the code working with the skull data and compare it to the output on pages 65 and 66. Run the command two more times using Euclidean and Mahalanobis distances. Put all three distance matrices into your lab report where you rounding to three significant figures. (1a.) Provide a paragraph of commentary on the distances with respect to the data. (2b.)Carefully look at the provided R code. Tell me what two commands are new to you. Research these two new commands and in your own words explain to me what they do.

Distances based on proportions

Ex 1. When there is NO overlap of classes:

 

Type 1

Type 2

Type 3

Colony 1 proportions

0.7

0.3

0.0

Colony 2 proportions

0

0

1.0

1115_dissimilarity index.pngThis is called the dissimilarity index.

833_similarity index.png, and this is called the similarity index.

Another index:

2042_similarity_index_1.png

 

Ex 2. When there is a complete overlap of classes:

 

Type 1

Type 2

Type 3

Colony 1 proportions

0.7

0.3

0.0

Colony 2 proportions

0.7

0.3

0.0

2337_distance measure.png

* In other cases these indices have values between 0 and 1.

* Similarity measures are often constructed as 1/D or 1/(1+D), where D = distance measure.

Problem 2:(2a.) Create a function that returns both the results of equation 5.5 and 5.6. (Use the list( ) command as done in the distmatrix( ) function.) Input should be two vectors p1 and p2 which represent the proportions for each species and each sum to 1. Comment your code.(2b.) Show your data for three different scenarios: No overlap, complete overlap, and partial overlap.

For problem 2a, use the following code to help you.
dissimilarity<- function(p1,p2)

{

Your code here

return( list(d1=d1, d2=d2) ) # returning list of both indices

}

Presence-absence data

(Table 5.5 on text p. 68)

Site

1

2

3

4

5

6

7

8

9

10

Species 1

0

0

1

1

1

0

1

1

1

0

Species 2

1

1

1

1

0

0

0

0

1

1

Summary:

 

Species 2

 

 

Species 1

present

absent

TOTAL

present

a (=3)

b (=3)

a + b

absent

c (=3)

d (=1)

c + d

TOTAL

a + c

b + d

n

Similarity measures(all vary between 0 (=no similarity) and 1 (=complete similarity)):

• Simple matching index: (a+d)/n

• Ochiai index:

936_Ochiai index.png

• Dice-Sorensen index: 2a/(2a+b+c)

• Jaccard index: a/(a+b+c)

Problem 3:(3a.) Create a function that returns a list containing all four indices shown above (and on page 68). The input should be four values: PP,AP,PA,AA which match with a, b, c, and d. Comment your code.(3b.) Run your function for the three scenarios: (1) 0's for AP and PA, but non-zero values for PP and AA; (2) 0's for PP and AA, but non-zero values for PA and AP; (3) some blend between scenarios (1) and (2).

The following commands may help you:

sp1 <- c(0,1,0,1,1,0,1,1,0,0,0,1)

sp2 <- c(1,1,1,1,0,0,0,0,1,1,1,1)

table(sp1,sp2)

PresAbsIndex<- function(PP,PA,AP,AA)

{

[ Your code goes here ]

}

Problem 4:Install the package "ade4" followed by data(butterfly) to get the butterfly dataset described on pages 7 and 8 of our text. Our question is: "Is genetic similarity correlated with geographic distances between butterfly colonies?" Repeat the below code and fill in the genetic distance matrix using the distance measure we called d1 in problem 2. (Don't worry about using percentages rather than proportions, the relative distances end up just being 0-100 instead of 0-1.) Perform a Mantel Randomization test and interpret the results. Repeat the test using the measure of distance d2. Did it make a difference? Start off with the below R code to help you. (Show your for-loop code as part of your answer as well as the scatter plots and randomization plots from the Mantel Randomization tests.)

library(ade4)

data(butterfly)

names(butterfly)

help(butterfly)

help(pch)

help(polygon)

help(segments)

plot(butterfly$contour[,1:2], pch=16, cex=.4) # set up x and y limits for graph

#segments(x0=butterfly$contour[,1],y0=butterfly$contour[,2],

# x1=butterfly$contour[,3],y1=butterfly$contour[,4])

polygon(butterfly$contour[,1:2], lty=2) # does same as segments()

points(butterfly$xy, pch=7)

nrow(butterfly$xy)

text(butterfly$xy, labels=1:16, pos=2, cex=.8)

apply(butterfly$genet,1,sum) # confirming rows add to 100

(Ddist<- dist(butterfly$xy))

dissimilarity(butterfly$genet[1,],butterfly$genet[2,])$d2 # testing it out

Dgenet<- matrix(NA,nrow=16,ncol=16)

[Your code, consisting of two nested for-loops, to fill in Dgenetmatrix . Use the function you created in problem 2.]

Dgenet # to see what matrix now looks like

Dgenet<- as.dist(Dgenet)

plot(Ddist,Dgenet)

cor(Ddist, Dgenet)

(mantel.out<- mantel.rtest(Dgenet,Ddist, nrepet=10000))

plot(mantel.out).

History, Academics

  • Category:- History
  • Reference No.:- M91295558
  • Price:- $20

Guranteed 24 Hours Delivery, In Price:- $20

Have any Question?


Related Questions in History

Complete an 8-page report with a strong thesis statement on

Complete an 8-page report with a Strong Thesis statement on the history of St. Augustine. The oldest European settlement in the U.S., St. Augustine was discovered in 1565 - 42 years before the Jamestown colony in Virgini ...

Question what is history why should we study historywithin

Question : What is history? Why should we study history? Within the context of our stories concerning Christopher Columbus, Native Americans, the Pilgrims, or slavery, discuss how historians (Zinn and Loewen) have dealt ...

Clarifying some things on the revolutioni am going to say

Clarifying some things on the Revolution I am going to say something, and I want you to hear me. I am a scholar of the Revolution. That's the topic of my dissertation. Please believe me when I say that I know a lot about ...

Discussion for assignment the protestant

Discussion for Assignment : The Protestant Reformation Sections A: What if anything did you know about the Protestant Reformation prior to this class? Were you aware of how much control the Roman Catholic Church had over ...

Paper assignmentthen each faction member should choose one

PAPER ASSIGNMENT Then each faction member should choose one of the characters listed below (it's OK if some characters are shared by more than one person) and drawing on the textbook readings and the links provided on yo ...

1 discussion forumcompare the various motivations for

1. Discussion Forum Compare the various motivations for establishing the colonies. How well did the colonies achieve their original purpose or have to adapt to unexpected situations in the New World? Give specific exampl ...

The purpose of this discussion assignment is to reflect

The purpose of this discussion assignment is to reflect upon your understanding of different historical and philosophical approaches to the self. After completing your readings, write a 2-3-page paper explaining the avoc ...

Answer the following question what is a human beings true

Answer the following Question : What is a human being's true identity? To what degree does society form our identity? Is this our "true" identity, or is it something else? Must be 500 words, APA format , must incorporate ...

Consider the contemporary expressions of these faiths and

Consider the contemporary expressions of these faiths and to reflect upon their relevance to your own and our collective societal experience--To consider how central teachings, principals and practices of these faith tra ...

This exam is open-book open note you may use any sources

This exam is open-book, open note. You may use any sources you wish. However, you must document (with end notes or in-text references) all sources that you use. Write in your own words, quote (from primary sources) spari ...

  • 4,153,160 Questions Asked
  • 13,132 Experts
  • 2,558,936 Questions Answered

Ask Experts for help!!

Looking for Assignment Help?

Start excelling in your Courses, Get help with Assignment

Write us your full requirement for evaluation and you will receive response within 20 minutes turnaround time.

Ask Now Help with Problems, Get a Best Answer

Why might a bank avoid the use of interest rate swaps even

Why might a bank avoid the use of interest rate swaps, even when the institution is exposed to significant interest rate

Describe the difference between zero coupon bonds and

Describe the difference between zero coupon bonds and coupon bonds. Under what conditions will a coupon bond sell at a p

Compute the present value of an annuity of 880 per year

Compute the present value of an annuity of $ 880 per year for 16 years, given a discount rate of 6 percent per annum. As

Compute the present value of an 1150 payment made in ten

Compute the present value of an $1,150 payment made in ten years when the discount rate is 12 percent. (Do not round int

Compute the present value of an annuity of 699 per year

Compute the present value of an annuity of $ 699 per year for 19 years, given a discount rate of 6 percent per annum. As