2/1/2016. Species Abundance Curves Plot of rank abundance (x-axis) vs abundance or P i (yaxis).

Size: px
Start display at page:

Download "2/1/2016. Species Abundance Curves Plot of rank abundance (x-axis) vs abundance or P i (yaxis)."

Transcription

1 Specie Abundance Curve Plot of rank abundance (x-axi) v abundance or P i (yaxi). More divere communitie lack numerically dominant pecie, flatter line. Proportion abundance A C DB F E A B C D E F abundance pecie rank pecie rank Rank # Compare rank abundance among creek or year # THIS FUNCTION WILL WAIT FOR YOU TO CLICK ON THE GRAPH TO PLACE THE LEGEND rank_data<-rankabuncomp(community,y=ample_factor,factor="creek") # Compare rank abundance among creek or year # THIS FUNCTION WILL WAIT FOR YOU TO CLICK ON THE GRAPH TO PLACE THE LEGEND rank_data<-rankabuncomp(community,y=ample_factor,factor="year") Specie Abundance Curve Function: rankabundance, rankabuncomp in the BiodiverityR package rankabundance(community,y=environmental,factor) Can ue either raw abundance or proportion data rankabuncomp allow for comparion among factor rankabunplot will plot the reult Baic Plotting Function plot doe a baic catter plot. plot(loblolly$age, Loblolly$height) abundance Specie2 Specie10 Specie4 Specie5 Specie7 rank_data<-rankabundance(community) rankabunplot(rank_data) pecie rank 1

2 Baic Plotting You can cutomize: Axi label with xlab or ylab variable Symbol with pch variable Symbol color with col variable Symbol ize with cex variable Heatmap Heatmap provide a quick way to viualize data. Function heatmap will plot raw value in a matrix: heatmap(a.matrix(community),rowv=na,colv=na) Can only plot a matrix, not a data frame. Ue a.matrix to convert, auming everything in your data frame i numeric. plot(loblolly$age,loblolly$height,xlab="age",ylab="height",pch=15,col="blue",cex=1.5) Baic Plotting Plot function will behave differently if one of the variable in a factor. Note pacing on a-axi. plot(a.factor(loblolly$age),loblolly$height,xlab="age",ylab="height",col="blue",cex=1. 5, cex.lab=1.5, cex.axi=1.5) Community Similarity or Diimilarity Community imilarity indice quantify imilarity among two ample. For a full community matrix do all poible pairwie comparion among communitie. Specie Sample Sample Sample 2

3 Quantitative Indice of Similarity (0-1.0) Symmetrical v. aymmetrical metric Are hared zero indicator of actual imilarity? Community data typically ue aymmetric metric Qualitative v. quantitative Ordinal v categorical, bionmical etc. Q mode v. R mode Q: Are actual object being compared (quetion i how imilar are A and B) E.g.: how imilar i community A and B R: Are relationhip or dependence of meaure of interet (quetion i if A i correlated with B) E.g. I pecie X abundance correlated with temperature Ruzicka( PSI) min Bray Curti Both are bound Typically log tranform data P i Where Pi = the proportion of the community compoed of pecie i. ( xij x ) ik x x ij Where x ij = i the abundance of pecie i in community j ik Community Similarity or Diimilarity Function: vegdit (vegan package), dvdi (labdv package), daiy (cluter package), deigndit(vegan) Very often raw abundance data can be ued Variable depending on propertie of metric, look at how each i calculated Mot are bound (0-1 range) Converion to imilarity or diimilarity Sim_matrix<-vegdit(community) Dim_matrix<-1-Sim_matrix ruzicka Bray-Curti v PSI ditance for local fih community dataet bray_curt 3

4 Quantitative Indice of Similarity (unbound) Qualitative Indice of Similarity (0-1.0) Euclidian ( xij xik ) 2 Manhattan= x Typically done without tranformation. Some ue preence/abence matrix with thee metric. No upper limit. ij x ik Where a = number of pecie in both communitie, b= number of pecie unique to community 1, c = number of pecie unique to community 2 Both convert data to preence/abence. Can be converted to diimilarity by 1-Steinhau or Sorenen. Steinhau a a b c 2a Sorenen 1 2( a b c) Euclidian v Manhattan ditance for local fih community dataet Steinhau v Sorenen ditance for local fih community dataet manhattan teinhau euclidian orenen 4

5 Bray-Curti v Sorenen ditance for local fih community dataet orenen Genetic Ditance A variety of package will work with genetic data and calculate ditance matrice. Package adegenet ha a function dit.genpop to calculate five common meaure (?dit.genpop for detail) # genetic ditance library(adegenet) #load a microat dataet that come with the package data(microatt) # table of allele frequencie (18 population x 112 allel) microatt$tab dim(microatt$tab) # euclidean ditance mat<-euclidean<-dit(microatt$tab) # Nei' ditance nei<-dit.genpop(a.genpop(microatt$tab), method=1) plot(mat,nei) INRA5.133 INRA5.137 INRA5.139 INRA5.141 INRA5.143 INRA Baoule Borgou BPN Charolai Holtein Jerey Lagunaire Limouin MaineAnjou Mtbeliard NDama Normand Parthenai Somba Vogienne ZChoa ZMbororo Zpeul bray_curt Defining your own function Function deigndit in the vegan package allow you to define any imilarity index. mat<-euclidean<-dit(microatt$tab) # Nei' ditance nei<-dit.genpop(a.genpop(microatt$tab), method=1) plot(mat,nei) heatmap(a.matrix(nei),rowv=t,colv=t) deigndit(community,method = (B+C)/2, abcd=true) review of 24 meaure of beta diverity 5

6 Heat Map of imilarity indice Unordered diimilarity (braycurti) heatmap(x=a.matrix(braycurt),rowv=na,colv=na) Ordered diimilarity (braycurti) with dendrogram in the margin heatmap(x=a.matrix(braycurt), Rowv=T,Colv=T) New function thi week pecaccum diverityreult vegdi, dvdi, deigndi rankabundance, rankabuncomp, rankabunplot plot heatmap dit.genepop RColorBrewer RColorBreer i a package that ha nicer palate of color for thing like heat map. library(rcolorbrewer) diplay.brewer.all(type="eq") heatmap(x=a.matrix(braycurt),col = brewer.pal(9,"purd"), Rowv=T,Colv=T) Aignment Ue the ame dataet a lat week (Fun_08_09_community.cv) Load data Eliminate rare pecie (any pecie with le than three occurrence) Eliminate ample with le than 50 individual. Calculate the etimated rarefied richne (and tandard deviation of the etimate) at 20 ample. Produce a rank abundance plot for the community (all pecie and all ample) With rare pecie and low abundance ample eliminated, calculate Bray Curti and Steinhau imilarity indice. Preent the Bray Curti imilarity a a heatmap uing any color cheme except the default. 6

Rarefaction Example. Consider this dataset: Original matrix:

Rarefaction Example. Consider this dataset: Original matrix: Rarefaction Example Conider thi dataet: Where i diverity highet? S 4 6 6 6 Shannon 1.0375911 0.9176461 0.9908044 1.0397044 What about rarefied diverity? rarefy(community,ample=10) 3.175905 2.576947 2.889674

More information

Week 3 Statistics for bioinformatics and escience

Week 3 Statistics for bioinformatics and escience Week 3 Statitic for bioinformatic and escience Line Skotte 28. november 2008 2.9.3-4) In thi eercie we conider microrna data from Human and Moue. The data et repreent 685 independent realiation of the

More information

Clustering Methods without Given Number of Clusters

Clustering Methods without Given Number of Clusters Clutering Method without Given Number of Cluter Peng Xu, Fei Liu Introduction A we now, mean method i a very effective algorithm of clutering. It mot powerful feature i the calability and implicity. However,

More information

Confusion matrices. True / False positives / negatives. INF 4300 Classification III Anne Solberg The agenda today: E.g., testing for cancer

Confusion matrices. True / False positives / negatives. INF 4300 Classification III Anne Solberg The agenda today: E.g., testing for cancer INF 4300 Claification III Anne Solberg 29.10.14 The agenda today: More on etimating claifier accuracy Cure of dimenionality knn-claification K-mean clutering x i feature vector for pixel i i- The cla label

More information

Multivariate analysis of genetic data: exploring group diversity

Multivariate analysis of genetic data: exploring group diversity Practical course using the software Multivariate analysis of genetic data: exploring group diversity Thibaut Jombart Abstract This practical course tackles the question of group diversity in genetic data

More information

Comparing Means: t-tests for Two Independent Samples

Comparing Means: t-tests for Two Independent Samples Comparing ean: t-tet for Two Independent Sample Independent-eaure Deign t-tet for Two Independent Sample Allow reearcher to evaluate the mean difference between two population uing data from two eparate

More information

Social Studies 201 Notes for March 18, 2005

Social Studies 201 Notes for March 18, 2005 1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the

More information

Source slideplayer.com/fundamentals of Analytical Chemistry, F.J. Holler, S.R.Crouch. Chapter 6: Random Errors in Chemical Analysis

Source slideplayer.com/fundamentals of Analytical Chemistry, F.J. Holler, S.R.Crouch. Chapter 6: Random Errors in Chemical Analysis Source lideplayer.com/fundamental of Analytical Chemitry, F.J. Holler, S.R.Crouch Chapter 6: Random Error in Chemical Analyi Random error are preent in every meaurement no matter how careful the experimenter.

More information

Suggested Answers To Exercises. estimates variability in a sampling distribution of random means. About 68% of means fall

Suggested Answers To Exercises. estimates variability in a sampling distribution of random means. About 68% of means fall Beyond Significance Teting ( nd Edition), Rex B. Kline Suggeted Anwer To Exercie Chapter. The tatitic meaure variability among core at the cae level. In a normal ditribution, about 68% of the core fall

More information

Social Studies 201 Notes for November 14, 2003

Social Studies 201 Notes for November 14, 2003 1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the

More information

Avoiding Forbidden Submatrices by Row Deletions

Avoiding Forbidden Submatrices by Row Deletions Avoiding Forbidden Submatrice by Row Deletion Sebatian Wernicke, Jochen Alber, Jen Gramm, Jiong Guo, and Rolf Niedermeier Wilhelm-Schickard-Intitut für Informatik, niverität Tübingen, Sand 13, D-72076

More information

MINITAB Stat Lab 3

MINITAB Stat Lab 3 MINITAB Stat 20080 Lab 3. Statitical Inference In the previou lab we explained how to make prediction from a imple linear regreion model and alo examined the relationhip between the repone and predictor

More information

Multiple Sequence Alignment. Progressive Alignment Iterative Pairwise Guide Tree ClustalW Co-linearity Multiple Sequence Alignment Editors

Multiple Sequence Alignment. Progressive Alignment Iterative Pairwise Guide Tree ClustalW Co-linearity Multiple Sequence Alignment Editors Yverdon Le Bain Introduction to Bioinformatic Introduction to Bioinformatic ami Khuri Department of Computer cience an Joé tate Univerity an Joé, California, UA khuri@c.ju.edu www.c.ju.edu/faculty/khuri

More information

Advanced Digital Signal Processing. Stationary/nonstationary signals. Time-Frequency Analysis... Some nonstationary signals. Time-Frequency Analysis

Advanced Digital Signal Processing. Stationary/nonstationary signals. Time-Frequency Analysis... Some nonstationary signals. Time-Frequency Analysis Advanced Digital ignal Proceing Prof. Nizamettin AYDIN naydin@yildiz.edu.tr Time-Frequency Analyi http://www.yildiz.edu.tr/~naydin 2 tationary/nontationary ignal Time-Frequency Analyi Fourier Tranform

More information

SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD

SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD S.P. Teeuwen, I. Erlich U. Bachmann Univerity of Duiburg, Germany Department of Electrical Power Sytem

More information

Multivariate analysis of genetic data exploring group diversity

Multivariate analysis of genetic data exploring group diversity Multivariate analysis of genetic data exploring group diversity Thibaut Jombart, Marie-Pauline Beugin MRC Centre for Outbreak Analysis and Modelling Imperial College London Genetic data analysis with PR

More information

Z a>2 s 1n = X L - m. X L = m + Z a>2 s 1n X L = The decision rule for this one-tail test is

Z a>2 s 1n = X L - m. X L = m + Z a>2 s 1n X L = The decision rule for this one-tail test is M09_BERE8380_12_OM_C09.QD 2/21/11 3:44 PM Page 1 9.6 The Power of a Tet 9.6 The Power of a Tet 1 Section 9.1 defined Type I and Type II error and their aociated rik. Recall that a repreent the probability

More information

Suggestions - Problem Set (a) Show the discriminant condition (1) takes the form. ln ln, # # R R

Suggestions - Problem Set (a) Show the discriminant condition (1) takes the form. ln ln, # # R R Suggetion - Problem Set 3 4.2 (a) Show the dicriminant condition (1) take the form x D Ð.. Ñ. D.. D. ln ln, a deired. We then replace the quantitie. 3ß D3 by their etimate to get the proper form for thi

More information

Lecture 7: Testing Distributions

Lecture 7: Testing Distributions CSE 5: Sublinear (and Streaming) Algorithm Spring 014 Lecture 7: Teting Ditribution April 1, 014 Lecturer: Paul Beame Scribe: Paul Beame 1 Teting Uniformity of Ditribution We return today to property teting

More information

Exercise set 1 Solution

Exercise set 1 Solution Exercie et 1 Solution 1. What are the main poible advantage of DNA-baed computing? How about diadvantage? Advantage: - More parallel computation - Able to tore tremendou amount of information in very mall

More information

LINEAR ALGEBRA METHOD IN COMBINATORICS. Theorem 1.1 (Oddtown theorem). In a town of n citizens, no more than n clubs can be formed under the rules

LINEAR ALGEBRA METHOD IN COMBINATORICS. Theorem 1.1 (Oddtown theorem). In a town of n citizens, no more than n clubs can be formed under the rules LINEAR ALGEBRA METHOD IN COMBINATORICS 1 Warming-up example Theorem 11 (Oddtown theorem) In a town of n citizen, no more tha club can be formed under the rule each club have an odd number of member each

More information

Publication V by authors

Publication V by authors Publication Kontantin S. Kotov and Jorma J. Kyyrä. 008. nertion lo and network parameter in the analyi of power filter. n: Proceeding of the 008 Nordic Workhop on Power and ndutrial Electronic (NORPE 008).

More information

1 Routh Array: 15 points

1 Routh Array: 15 points EE C28 / ME34 Problem Set 3 Solution Fall 2 Routh Array: 5 point Conider the ytem below, with D() k(+), w(t), G() +2, and H y() 2 ++2 2(+). Find the cloed loop tranfer function Y () R(), and range of k

More information

A Study on Simulating Convolutional Codes and Turbo Codes

A Study on Simulating Convolutional Codes and Turbo Codes A Study on Simulating Convolutional Code and Turbo Code Final Report By Daniel Chang July 27, 2001 Advior: Dr. P. Kinman Executive Summary Thi project include the deign of imulation of everal convolutional

More information

List coloring hypergraphs

List coloring hypergraphs Lit coloring hypergraph Penny Haxell Jacque Vertraete Department of Combinatoric and Optimization Univerity of Waterloo Waterloo, Ontario, Canada pehaxell@uwaterloo.ca Department of Mathematic Univerity

More information

New bounds for Morse clusters

New bounds for Morse clusters New bound for More cluter Tamá Vinkó Advanced Concept Team, European Space Agency, ESTEC Keplerlaan 1, 2201 AZ Noordwijk, The Netherland Tama.Vinko@ea.int and Arnold Neumaier Fakultät für Mathematik, Univerität

More information

k-means clustering mark = which(md == min(md)) nearest[i] = ifelse(mark <= 5, "blue", "orange")}

k-means clustering mark = which(md == min(md)) nearest[i] = ifelse(mark <= 5, blue, orange)} 1 / 16 k-means clustering km15 = kmeans(x[g==0,],5) km25 = kmeans(x[g==1,],5) for(i in 1:6831){ md = c(mydist(xnew[i,],km15$center[1,]),mydist(xnew[i,],km15$center[2, mydist(xnew[i,],km15$center[3,]),mydist(xnew[i,],km15$center[4,]),

More information

Theoretical Population Biology

Theoretical Population Biology Theoretical Population Biology 8 () 3 Content lit available at ScienceDirect Theoretical Population Biology journal homepage: www.elevier.com/locate/tpb On the ize ditribution of private microatellite

More information

Chapter 4. The Laplace Transform Method

Chapter 4. The Laplace Transform Method Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination

More information

2/7/2018. Strata. Strata

2/7/2018. Strata. Strata The strata option allows you to control how permutations are done. Specifically, to constrain permutations. Why would you want to do this? In this dataset, there are clear differences in area (A vs. B),

More information

Nonlinear Single-Particle Dynamics in High Energy Accelerators

Nonlinear Single-Particle Dynamics in High Energy Accelerators Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction

More information

11.2 Stability. A gain element is an active device. One potential problem with every active circuit is its stability

11.2 Stability. A gain element is an active device. One potential problem with every active circuit is its stability 5/7/2007 11_2 tability 1/2 112 tability eading Aignment: pp 542-548 A gain element i an active device One potential problem with every active circuit i it tability HO: TABIITY Jim tile The Univ of Kana

More information

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is EE 4G Note: Chapter 6 Intructor: Cheung More about ZSR and ZIR. Finding unknown initial condition: Given the following circuit with unknown initial capacitor voltage v0: F v0/ / Input xt 0Ω Output yt -

More information

A statistical analysis of the Evolutionary Trace method. Heidi Spratt

A statistical analysis of the Evolutionary Trace method. Heidi Spratt 1 A tatitical analyi of the Evolutionary Trace method Heidi Spratt 2 Content Lit of Illutration iii 1 Introduction 1 Evolution of Protein 1.1 Protein 3 1.2 Multiple Sequence Alignment 4 1.3 Model of Protein

More information

Annex-A: RTTOV9 Cloud validation

Annex-A: RTTOV9 Cloud validation RTTOV-91 Science and Validation Plan Annex-A: RTTOV9 Cloud validation Author O Embury C J Merchant The Univerity of Edinburgh Intitute for Atmo. & Environ. Science Crew Building King Building Edinburgh

More information

Chapter 12 Simple Linear Regression

Chapter 12 Simple Linear Regression Chapter 1 Simple Linear Regreion Introduction Exam Score v. Hour Studied Scenario Regreion Analyi ued to quantify the relation between (or more) variable o you can predict the value of one variable baed

More information

Çankaya University ECE Department ECE 376 (MT)

Çankaya University ECE Department ECE 376 (MT) Çankaya Univerity ECE Department ECE 376 (M) Student Name : Date : 13.4.15 Student Number : Open Source Exam Quetion 1. (7 Point) he time waveform of the ignal et, and t t are given in Fig. 1.1. a. Identify

More information

March 18, 2014 Academic Year 2013/14

March 18, 2014 Academic Year 2013/14 POLITONG - SHANGHAI BASIC AUTOMATIC CONTROL Exam grade March 8, 4 Academic Year 3/4 NAME (Pinyin/Italian)... STUDENT ID Ue only thee page (including the back) for anwer. Do not ue additional heet. Ue of

More information

Practice Problems - Week #7 Laplace - Step Functions, DE Solutions Solutions

Practice Problems - Week #7 Laplace - Step Functions, DE Solutions Solutions For Quetion -6, rewrite the piecewie function uing tep function, ketch their graph, and find F () = Lf(t). 0 0 < t < 2. f(t) = (t 2 4) 2 < t In tep-function form, f(t) = u 2 (t 2 4) The graph i the olid

More information

Fair Game Review. Chapter 7 A B C D E Name Date. Complete the number sentence with <, >, or =

Fair Game Review. Chapter 7 A B C D E Name Date. Complete the number sentence with <, >, or = Name Date Chapter 7 Fair Game Review Complete the number entence with , or =. 1. 3.4 3.45 2. 6.01 6.1 3. 3.50 3.5 4. 0.84 0.91 Find three decimal that make the number entence true. 5. 5.2 6. 2.65 >

More information

Codes Correcting Two Deletions

Codes Correcting Two Deletions 1 Code Correcting Two Deletion Ryan Gabry and Frederic Sala Spawar Sytem Center Univerity of California, Lo Angele ryan.gabry@navy.mil fredala@ucla.edu Abtract In thi work, we invetigate the problem of

More information

Channel Coding 2.

Channel Coding 2. Channel Coding Dr.-Ing. Dirk Wübben Intitute for Telecommunication and High-Frequency Technique Department of Communication Engineering Room: N3, Phone: /8-6385 wuebben@ant.uni-bremen.de Lecture Tueday,

More information

ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang

ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang

More information

R. W. Erickson. Department of Electrical, Computer, and Energy Engineering University of Colorado, Boulder

R. W. Erickson. Department of Electrical, Computer, and Energy Engineering University of Colorado, Boulder R. W. Erickon Department of Electrical, Computer, and Energy Engineering Univerity of Colorado, Boulder ZOH: Sampled Data Sytem Example v T Sampler v* H Zero-order hold H v o e = 1 T 1 v *( ) = v( jkω

More information

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou

More information

Problem Set 8 Solutions

Problem Set 8 Solutions Deign and Analyi of Algorithm April 29, 2015 Maachuett Intitute of Technology 6.046J/18.410J Prof. Erik Demaine, Srini Devada, and Nancy Lynch Problem Set 8 Solution Problem Set 8 Solution Thi problem

More information

EE 508 Lecture 16. Filter Transformations. Lowpass to Bandpass Lowpass to Highpass Lowpass to Band-reject

EE 508 Lecture 16. Filter Transformations. Lowpass to Bandpass Lowpass to Highpass Lowpass to Band-reject EE 508 Lecture 6 Filter Tranformation Lowpa to Bandpa Lowpa to Highpa Lowpa to Band-reject Review from Lat Time Theorem: If the perimeter variation and contact reitance are neglected, the tandard deviation

More information

Department of Mechanical Engineering Massachusetts Institute of Technology Modeling, Dynamics and Control III Spring 2002

Department of Mechanical Engineering Massachusetts Institute of Technology Modeling, Dynamics and Control III Spring 2002 Department of Mechanical Engineering Maachuett Intitute of Technology 2.010 Modeling, Dynamic and Control III Spring 2002 SOLUTIONS: Problem Set # 10 Problem 1 Etimating tranfer function from Bode Plot.

More information

Regression. What is regression? Linear Regression. Cal State Northridge Ψ320 Andrew Ainsworth PhD

Regression. What is regression? Linear Regression. Cal State Northridge Ψ320 Andrew Ainsworth PhD Regreion Cal State Northridge Ψ30 Andrew Ainworth PhD What i regreion? How do we predict one variable from another? How doe one variable change a the other change? Caue and effect Linear Regreion A technique

More information

Estimation of Peaked Densities Over the Interval [0,1] Using Two-Sided Power Distribution: Application to Lottery Experiments

Estimation of Peaked Densities Over the Interval [0,1] Using Two-Sided Power Distribution: Application to Lottery Experiments MPRA Munich Peronal RePEc Archive Etimation of Peaed Denitie Over the Interval [0] Uing Two-Sided Power Ditribution: Application to Lottery Experiment Krzyztof Konte Artal Invetment 8. April 00 Online

More information

STATISTICAL SIGNIFICANCE

STATISTICAL SIGNIFICANCE STATISTICAL SIGNIFICANCE EFFECT SIZE More to life than tatitical ignificance Reporting effect ize Turn out a lot of reearcher do not know what preciely p

More information

Numerical algorithm for the analysis of linear and nonlinear microstructure fibres

Numerical algorithm for the analysis of linear and nonlinear microstructure fibres Numerical algorithm for the anali of linear and nonlinear microtructure fibre Mariuz Zdanowicz *, Marian Marciniak, Marek Jaworki, Igor A. Goncharenko National Intitute of Telecommunication, Department

More information

Physics 218: Exam 1. Class of 2:20pm. February 14th, You have the full class period to complete the exam.

Physics 218: Exam 1. Class of 2:20pm. February 14th, You have the full class period to complete the exam. Phyic 218: Exam 1 Cla of 2:20pm February 14th, 2012. Rule of the exam: 1. You have the full cla period to complete the exam. 2. Formulae are provided on the lat page. You may NOT ue any other formula heet.

More information

MAE140 Linear Circuits Fall 2012 Final, December 13th

MAE140 Linear Circuits Fall 2012 Final, December 13th MAE40 Linear Circuit Fall 202 Final, December 3th Intruction. Thi exam i open book. You may ue whatever written material you chooe, including your cla note and textbook. You may ue a hand calculator with

More information

The Influence of the Load Condition upon the Radial Distribution of Electromagnetic Vibration and Noise in a Three-Phase Squirrel-Cage Induction Motor

The Influence of the Load Condition upon the Radial Distribution of Electromagnetic Vibration and Noise in a Three-Phase Squirrel-Cage Induction Motor The Influence of the Load Condition upon the Radial Ditribution of Electromagnetic Vibration and Noie in a Three-Phae Squirrel-Cage Induction Motor Yuta Sato 1, Iao Hirotuka 1, Kazuo Tuboi 1, Maanori Nakamura

More information

The machines in the exercise work as follows:

The machines in the exercise work as follows: Tik-79.148 Spring 2001 Introduction to Theoretical Computer Science Tutorial 9 Solution to Demontration Exercie 4. Contructing a complex Turing machine can be very laboriou. With the help of machine chema

More information

Inference for Two Stage Cluster Sampling: Equal SSU per PSU. Projections of SSU Random Variables on Each SSU selection.

Inference for Two Stage Cluster Sampling: Equal SSU per PSU. Projections of SSU Random Variables on Each SSU selection. Inference for Two Stage Cluter Sampling: Equal SSU per PSU Projection of SSU andom Variable on Eac SSU election By Ed Stanek Introduction We review etimating equation for PSU mean in a two tage cluter

More information

Correction for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002

Correction for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002 Correction for Simple Sytem Example and Note on Laplace Tranform / Deviation Variable ECHE 55 Fall 22 Conider a tank draining from an initial height of h o at time t =. With no flow into the tank (F in

More information

Reporting Checklist for Nature Neuroscience

Reporting Checklist for Nature Neuroscience Correponding Author: Manucript Number: 1 Manucript Type: Maanori Matuzaki Article Reporting Checklit for Nature Neurocience # Main Figure: # Supplementary Figure: # Supplementary Table: # Supplementary

More information

A Simple Example Binary Hypothesis Testing Optimal Receiver Frontend M-ary Signal Sets Message Sequences. = 4 for QPSK) E b Q 2. 2E b.

A Simple Example Binary Hypothesis Testing Optimal Receiver Frontend M-ary Signal Sets Message Sequences. = 4 for QPSK) E b Q 2. 2E b. Exercie: QPSK I Find the error rate for the ignal et n (t) = p 2E /T co(2pf c t + n p/2 + p/4), for n = 0,, 3 I Anwer: (Recall h P = d min 2 E b = 4 for QPSK) E Pr{e} = 2Q Q 2 = 2Q = 2Q 2E b h P E b 2

More information

CE3502. EMMA COMPARISON OF NUMBERS. Earth Week: Film, lecture International World Water Day: Monday Posters on water research at MTU

CE3502. EMMA COMPARISON OF NUMBERS. Earth Week: Film, lecture International World Water Day: Monday Posters on water research at MTU CE350. EMMA COMPARISON OF NUMBERS Earth Week: Film, lecture International World Water Da: Monda Poter on water reearch at MTU Majora Carter Greening the Ghetto Saturda, March 0 7:30pm Roza Center Majora

More information

EE 508 Lecture 16. Filter Transformations. Lowpass to Bandpass Lowpass to Highpass Lowpass to Band-reject

EE 508 Lecture 16. Filter Transformations. Lowpass to Bandpass Lowpass to Highpass Lowpass to Band-reject EE 508 Lecture 6 Filter Tranformation Lowpa to Bandpa Lowpa to Highpa Lowpa to Band-reject Review from Lat Time Theorem: If the perimeter variation and contact reitance are neglected, the tandard deviation

More information

Testing the Equality of Two Pareto Distributions

Testing the Equality of Two Pareto Distributions Proceeding of the World Congre on Engineering 07 Vol II WCE 07, July 5-7, 07, London, U.K. Teting the Equality of Two Pareto Ditribution Huam A. Bayoud, Member, IAENG Abtract Thi paper propoe an overlapping-baed

More information

RELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS

RELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS www.arpapre.com/volume/vol29iue1/ijrras_29_1_01.pdf RELIABILITY OF REPAIRABLE k out of n: F SYSTEM HAVING DISCRETE REPAIR AND FAILURE TIMES DISTRIBUTIONS Sevcan Demir Atalay 1,* & Özge Elmataş Gültekin

More information

Spring 2014 EE 445S Real-Time Digital Signal Processing Laboratory. Homework #0 Solutions on Review of Signals and Systems Material

Spring 2014 EE 445S Real-Time Digital Signal Processing Laboratory. Homework #0 Solutions on Review of Signals and Systems Material Spring 4 EE 445S Real-Time Digital Signal Proceing Laboratory Prof. Evan Homework # Solution on Review of Signal and Sytem Material Problem.. Continuou-Time Sinuoidal Generation. In practice, we cannot

More information

Assessing the Discriminatory Power of Credit Scores under Censoring

Assessing the Discriminatory Power of Credit Scores under Censoring Aeing the Dicriminatory Power of Credit Score under Cenoring Holger Kraft, Gerald Kroiandt, Marlene Müller Fraunhofer Intitut für Techno- und Wirtchaftmathematik (ITWM) Thi verion: Augut 27, 2003 Abtract:

More information

THE Z-TRANSFORM APPLIED TO BIRTH-DEATH MARKOV PROCESSES

THE Z-TRANSFORM APPLIED TO BIRTH-DEATH MARKOV PROCESSES he Zranform applied to birthdeath arov procee HE ZRASOR ALIED O BIRHDEAH ARKO ROCESSES H Department of Electronic and Communication niverity of DareSalaam O Box 3594 DareSalaam anania ABSRAC Birthdeath

More information

Common Weights in DEA Models with Reduced Singular Value Decomposition

Common Weights in DEA Models with Reduced Singular Value Decomposition Journal of Data Science 419-430, DOI: 10.6339/JDS.201804_16(2).0010 Common Weight in DEA Model with Reduced Singular alue Decompoition Haan Naeri 1, S.Emaeil Najafi 2,Abba Saghaei 3 1,2,3 Department of

More information

COPRIME ARRAYS AND SAMPLERS FOR SPACE-TIME ADAPTIVE PROCESSING

COPRIME ARRAYS AND SAMPLERS FOR SPACE-TIME ADAPTIVE PROCESSING COPRIE ARRAYS AND SAPLERS FOR SPACE-IE ADAPIVE PROCESSING Chun-Lin Liu 1 and P. P. Vaidyanathan 2 Dept. of Electrical Engineering, 136-93 California Intitute of echnology, Paadena, CA 91125, USA E-mail:

More information

diversity(datamatrix, index= shannon, base=exp(1))

diversity(datamatrix, index= shannon, base=exp(1)) Tutorial 11: Diversity, Indicator Species Analysis, Cluster Analysis Calculating Diversity Indices The vegan package contains the command diversity() for calculating Shannon and Simpson diversity indices.

More information

Dissimilarity and transformations. Pierre Legendre Département de sciences biologiques Université de Montréal

Dissimilarity and transformations. Pierre Legendre Département de sciences biologiques Université de Montréal and transformations Pierre Legendre Département de sciences biologiques Université de Montréal http://www.numericalecology.com/ Pierre Legendre 2017 Definitions An association coefficient is a function

More information

Chapter 8.2. Interval Estimation

Chapter 8.2. Interval Estimation Chapter 8.2. Iterval Etimatio Baic of Cofidece Iterval ad Large Sample Cofidece Iterval 1 Baic Propertie of Cofidece Iterval Aumptio: X 1, X 2,, X are from Normal ditributio with a mea of µ ad tadard deviatio.

More information

Combining allele frequency uncertainty and population substructure corrections in forensic DNA calculations

Combining allele frequency uncertainty and population substructure corrections in forensic DNA calculations Combining allele frequency uncertainty and population ubtructure correction in forenic DNA calculation arxiv:1509.08361v2 [tat.ap] 6 Oct 2015 Robert Cowell Faculty of Actuarial Science and Inurance Ca

More information

Lecture 9: Shor s Algorithm

Lecture 9: Shor s Algorithm Quantum Computation (CMU 8-859BB, Fall 05) Lecture 9: Shor Algorithm October 7, 05 Lecturer: Ryan O Donnell Scribe: Sidhanth Mohanty Overview Let u recall the period finding problem that wa et up a a function

More information

Unified Correlation between SPT-N and Shear Wave Velocity for all Soil Types

Unified Correlation between SPT-N and Shear Wave Velocity for all Soil Types 6 th International Conference on Earthquake Geotechnical Engineering 1-4 ovember 15 Chritchurch, ew Zealand Unified Correlation between SPT- and Shear Wave Velocity for all Soil Type C.-C. Tai 1 and T.

More information

CHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL

CHAPTER 4 DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL 98 CHAPTER DESIGN OF STATE FEEDBACK CONTROLLERS AND STATE OBSERVERS USING REDUCED ORDER MODEL INTRODUCTION The deign of ytem uing tate pace model for the deign i called a modern control deign and it i

More information

NCAAPMT Calculus Challenge Challenge #3 Due: October 26, 2011

NCAAPMT Calculus Challenge Challenge #3 Due: October 26, 2011 NCAAPMT Calculu Challenge 011 01 Challenge #3 Due: October 6, 011 A Model of Traffic Flow Everyone ha at ome time been on a multi-lane highway and encountered road contruction that required the traffic

More information

( ) y = Properties of Gaussian curves: Can also be written as: where

( ) y = Properties of Gaussian curves: Can also be written as: where Propertie of Gauian curve: Can alo be written a: e ( x μ ) σ σ π e z σ π where z ( ) x μ σ z repreent the deviation of a reult from the population mean relative to the population tandard deviation. z i

More information

Midterm Review - Part 1

Midterm Review - Part 1 Honor Phyic Fall, 2016 Midterm Review - Part 1 Name: Mr. Leonard Intruction: Complete the following workheet. SHOW ALL OF YOUR WORK. 1. Determine whether each tatement i True or Fale. If the tatement i

More information

Asymptotics of ABC. Paul Fearnhead 1, Correspondence: Abstract

Asymptotics of ABC. Paul Fearnhead 1, Correspondence: Abstract Aymptotic of ABC Paul Fearnhead 1, 1 Department of Mathematic and Statitic, Lancater Univerity Correpondence: p.fearnhead@lancater.ac.uk arxiv:1706.07712v1 [tat.me] 23 Jun 2017 Abtract Thi document i due

More information

Bogoliubov Transformation in Classical Mechanics

Bogoliubov Transformation in Classical Mechanics Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How

More information

Learning Multiplicative Interactions

Learning Multiplicative Interactions CSC2535 2011 Lecture 6a Learning Multiplicative Interaction Geoffrey Hinton Two different meaning of multiplicative If we take two denity model and multiply together their probability ditribution at each

More information

Jump condition at the boundary between a porous catalyst and a homogeneous fluid

Jump condition at the boundary between a porous catalyst and a homogeneous fluid From the SelectedWork of Francico J. Valde-Parada 2005 Jump condition at the boundary between a porou catalyt and a homogeneou fluid Francico J. Valde-Parada J. Alberto Ochoa-Tapia Available at: http://work.bepre.com/francico_j_valde_parada/12/

More information

Efficient Methods of Doppler Processing for Coexisting Land and Weather Clutter

Efficient Methods of Doppler Processing for Coexisting Land and Weather Clutter Efficient Method of Doppler Proceing for Coexiting Land and Weather Clutter Ça gatay Candan and A Özgür Yılmaz Middle Eat Technical Univerity METU) Ankara, Turkey ccandan@metuedutr, aoyilmaz@metuedutr

More information

Multivariate Analysis of Ecological Data

Multivariate Analysis of Ecological Data Multivariate Analysis of Ecological Data MICHAEL GREENACRE Professor of Statistics at the Pompeu Fabra University in Barcelona, Spain RAUL PRIMICERIO Associate Professor of Ecology, Evolutionary Biology

More information

The use of connected masks for reconstructing the single particle image from X-ray diffraction data

The use of connected masks for reconstructing the single particle image from X-ray diffraction data Mathematical Biology and Bioinformatic 2015. V. 10. Suppl. P. t1 t19. doi: 10.17537/2015.10.t1 The tranlation of original article Lunin V.Y., Lunina N.L., Petrova T.E. Matematichekaya biologiya i bioinformatika.

More information

c n b n 0. c k 0 x b n < 1 b k b n = 0. } of integers between 0 and b 1 such that x = b k. b k c k c k

c n b n 0. c k 0 x b n < 1 b k b n = 0. } of integers between 0 and b 1 such that x = b k. b k c k c k 1. Exitence Let x (0, 1). Define c k inductively. Suppoe c 1,..., c k 1 are already defined. We let c k be the leat integer uch that x k An eay proof by induction give that and for all k. Therefore c n

More information

UNITS FOR THERMOMECHANICS

UNITS FOR THERMOMECHANICS UNITS FOR THERMOMECHANICS 1. Conitent Unit. Every calculation require a conitent et of unit. Hitorically, one et of unit wa ued for mechanic and an apparently unrelated et of unit wa ued for heat. For

More information

General Relativity (sort of)

General Relativity (sort of) The Okefenokee Swamp Some book about relativity: Taylor & Wheeler Spacetime Phyic = TW1 (imple preentation, but deep inight) Taylor & Wheeler Introduction to Black Hole = TW (CO-level math, deep inight)

More information

LECTURE 17. Linear coupling (continued) Linear coupling (continued) Pretzel Orbits

LECTURE 17. Linear coupling (continued) Linear coupling (continued) Pretzel Orbits LETURE 7 Linear coupling (continued) oupling coefficient for ditribution of kew quadrupole and olenoid Pretzel Orbit Motivation and application Implication Linear coupling (continued) oupling coefficient

More information

ARTICLE Overcoming the Winner s Curse: Estimating Penetrance Parameters from Case-Control Data

ARTICLE Overcoming the Winner s Curse: Estimating Penetrance Parameters from Case-Control Data ARTICLE Overcoming the Winner Cure: Etimating Penetrance Parameter from Cae-Control Data Sebatian Zöllner and Jonathan K. Pritchard Genomewide aociation tudie are now a widely ued approach in the earch

More information

Correlation. January 11, 2018

Correlation. January 11, 2018 Correlation January 11, 2018 Contents Correlations The Scattterplot The Pearson correlation The computational raw-score formula Survey data Fun facts about r Sensitivity to outliers Spearman rank-order

More information

The Dynamics of Learning Vector Quantization

The Dynamics of Learning Vector Quantization The Dynamic of Learning Vector Quantization Barbara Hammer TU Clauthal-Zellerfeld Intitute of Computing Science Michael Biehl, Anarta Ghoh Rijkuniveriteit Groningen Mathematic and Computing Science Introduction

More information

Pulsed Magnet Crimping

Pulsed Magnet Crimping Puled Magnet Crimping Fred Niell 4/5/00 1 Magnetic Crimping Magnetoforming i a metal fabrication technique that ha been in ue for everal decade. A large capacitor bank i ued to tore energy that i ued to

More information

arxiv: v1 [math.ca] 23 Sep 2017

arxiv: v1 [math.ca] 23 Sep 2017 arxiv:709.08048v [math.ca] 3 Sep 07 On the unit ditance problem A. Ioevich Abtract. The Erdő unit ditance conjecture in the plane ay that the number of pair of point from a point et of ize n eparated by

More information

Control Systems Engineering ( Chapter 7. Steady-State Errors ) Prof. Kwang-Chun Ho Tel: Fax:

Control Systems Engineering ( Chapter 7. Steady-State Errors ) Prof. Kwang-Chun Ho Tel: Fax: Control Sytem Engineering ( Chapter 7. Steady-State Error Prof. Kwang-Chun Ho kwangho@hanung.ac.kr Tel: 0-760-453 Fax:0-760-4435 Introduction In thi leon, you will learn the following : How to find the

More information

Alternate Dispersion Measures in Replicated Factorial Experiments

Alternate Dispersion Measures in Replicated Factorial Experiments Alternate Diperion Meaure in Replicated Factorial Experiment Neal A. Mackertich The Raytheon Company, Sudbury MA 02421 Jame C. Benneyan Northeatern Univerity, Boton MA 02115 Peter D. Krau The Raytheon

More information

(b) 99%; n = 15; σ is unknown; population appears to be normally distributed.

(b) 99%; n = 15; σ is unknown; population appears to be normally distributed. MTH 345 Exam 3 Fall 2013 Jutify all anwer with neat and organized work. Clearly indicate your anwer. 100 point poible. 1. (12 pt.) Women height are normally ditributed with mean 63.6 in. and tandard deviation

More information

Advanced Method for Small-Signal Stability Assessment based on Neuronal Networks

Advanced Method for Small-Signal Stability Assessment based on Neuronal Networks 1 Advanced Method for Small-Signal Stability Aement baed on Neuronal Networ S. P. Teeuwen, I. Erlich, Member, IEEE, and M. A. El-Sharawi, Fellow, IEEE Abtract-- Thi paper deal with a new method for eigenvalue

More information

MULTI-LAYERED LOSSY FINITE LENGTH DIELECTRIC CYLINDIRICAL MODEL OF MAN AT OBLIQUE INCIDENCE

MULTI-LAYERED LOSSY FINITE LENGTH DIELECTRIC CYLINDIRICAL MODEL OF MAN AT OBLIQUE INCIDENCE Proceeding 3rd Annual Conference IEEE/EMBS Oct.5-8, 1, Itanbul, TURKEY MULTI-LAYERED LOSSY FINITE LENGTH DIELECTRIC CYLINDIRICAL MODEL OF MAN AT OBLIQUE INCIDENCE S.S. Şeker, B. Yeldiren Boğaziçi Univerity,

More information