Bab 1 Mewakili fungsi matematika Seperti judulnya, dalam buku ini kita akan menggunakan bahasa komputer R untuk mengimplementasikan operasi kalkulus bersama dengan operasi terkait seperti grafik dan penyelesaian.

Topik kalkulus pada dasarnya adalah tentang fungsi matematika dan operasi yang dilakukan pada mereka. Konsep “fungsi matematika” adalah sebuah ide. Jika kita akan menggunakan bahasa komputer untuk bekerja dengan fungsi matematika, kita perlu menerjemahkannya ke dalam beberapa entitas dalam bahasa komputer. Artinya, kita memerlukan konstruksi bahasa untuk merepresentasikan fungsi dan kuantitas yang diambil fungsi sebagai input dan dihasilkan sebagai output.

Jika Anda kebetulan memiliki latar belakang pemrograman, Anda mungkin berpikir bahwa pilihan representasi sudah jelas. Untuk jumlah, gunakan angka. Untuk fungsi, gunakan fungsi bahasa R. Memang itulah yang akan kita lakukan, tetapi seperti yang akan Anda lihat, situasinya sedikit lebih rumit dari itu. Sedikit, tetapi kerumitan kecil itu perlu ditangani sejak awal.

1.1 Angka, besaran, dan nama Kompleksitas yang disebutkan di bagian sebelumnya berasal dari penggunaan dan situasi dunia nyata di mana kita ingin dapat menerapkan ide matematis fungsi. Input yang diambil oleh fungsi dan output yang dihasilkannya belum tentu berupa angka. Seringkali, mereka adalah kuantitas .

Perhatikan contoh besaran berikut: uang, kecepatan, tekanan darah, tinggi, volume. Bukankah masing-masing besaran ini adalah bilangan? Tidak terlalu. Pertimbangkan uang. Kita dapat menghitung, menambah, dan mengurangi uang – sama seperti kita menghitung, menambah, dan mengurangi dengan angka. Tetapi uang memiliki properti tambahan: jenis mata uang. Ada Dolar 1 , Renminbi, Euro, Rand, Krona dan sebagainya. Uang sebagai ide adalah abstraksi dari jenis yang kadang-kadang disebut dimensi . Dimensi lainnya adalah panjang, waktu, volume, suhu, luas, sudut, luminositas, suku bunga bank, dan sebagainya. Ketika kita berhadapan dengan suatu besaran yang memiliki dimensi, tidak cukup memberikan angka untuk mengatakan “berapa” dari dimensi yang kita miliki. Kita juga perlu memberikan unit. Banyak dari ini akrab dalam kehidupan sehari-hari. Misalnya, dimensi panjang diukur dengan satuan meter, atau inci, atau mil, atau parsec, dan seterusnya. Dimensi uang diukur dengan euro, dolar, renminbi, dan sebagainya. Area sebagai dimensi diukur dengan meter persegi, atau kaki persegi, atau hektar, atau ekar, dll.

Sayangnya, sebagian besar teks matematika memperlakukan dimensi dan satuan dengan mengabaikannya. Hal ini menyebabkan kebingungan dan kesalahan ketika upaya dilakukan untuk menerapkan ide-ide matematika ke kuantitas di dunia nyata.

Dalam buku ini kita akan menggunakan fungsi dan kalkulus untuk bekerja dengan kuantitas dunia nyata. Kita tidak bisa mengabaikan dimensi dan satuan. Sayangnya, bahasa komputer utama seperti R dan Python dan JavaScript tidak menyediakan cara sistematis untuk menangani dimensi dan unit secara otomatis. Di R, misalnya, kita dapat dengan mudah menulis

x <- 7 yang menyimpan kuantitas dengan nama x. Tapi bahasanya mengacungkan sesuatu seperti

y <- 12 meters ## Error: :1:9: unexpected symbol ## 1: y <- 12 meters ## ^ Kurangnya notasi komputer yang tepat untuk mewakili kuantitas berdimensi dengan satuannya, kita memerlukan beberapa cara lain untuk melacak berbagai hal.

Berikut adalah apa yang akan kita lakukan. Saat kita merepresentasikan besaran di komputer, kita akan menggunakan nama besaran untuk mengingatkan kita apa dimensi dan satuannya. Misalkan kita ingin, misalnya, untuk mewakili pendapatan tahunan keluarga. Pilihan yang masuk akal adalah menggunakan nama seperti incomeatau income_per_yearatau family_income. Kami bahkan mungkin memberikan unit dalam nama, misalnya family_income_euros_per_year. Tetapi biasanya yang kami lakukan adalah mendokumentasikan unit secara terpisah untuk pembaca manusia.

Gaya penggunaan nama besaran sebagai pengingat dimensi dan satuan besaran ini berada di luar tradisi notasi matematika. Dalam matematika sekolah menengah, Anda menemukan x dan kamu dan t dan θ . Anda bahkan mungkin pernah melihat subskrip dulu lebih spesifik, misalnya x 0 . Seseorang yang terlatih dalam fisika akan mengetahui penggunaan huruf secara tradisional dan tipikal. x kemungkinan adalah posisi, t kemungkinan adalah waktu, dan x 0 adalah posisi pada waktu nol.

Di banyak bidang aplikasi, ada banyak kuantitas berbeda yang akan direpresentasikan. Sebagai ilustrasi, pertimbangkan buku teks tentang aliran turbulen fluida: Aliran Turbulen oleh Stephen B. Pope (ISBN 978-0521598866). Dalam materi pendahuluan buku, ada bagian sepanjang 14 halaman berjudul “Nomenklatur.” Gambar 1.1 menunjukkan sebagian kecil dari bagian tersebut, hanya simbol yang dimulai dengan huruf kapital R dan S.

Bagian dari satu halaman dari bagian nomenklatur sepanjang 14 halaman dari Turbulent Flows oleh Stephen B. Pope Gambar 1.1: Bagian dari satu halaman dari bagian nomenklatur sepanjang 14 halaman dari Turbulent Flows oleh Stephen B. Pope

Setiap halaman di bagian nomenklatur padat dengan simbol. Inti dari kebanyakan dari mereka adalah satu atau dua huruf, misalnya
R atau R e . Tetapi hanya ada begitu banyak huruf yang tersedia, bahkan jika huruf Yunani dan lainnya dimasukkan ke dalam berbagai hal. Dan lebih banyak kekhususan ditambahkan dengan menggunakan subskrip, superskrip, topi (mis
^ S saya j ), batang (mis
¯ S saya j ), bintang, dan sebagainya. Favorit saya di halaman ini adalah R e δ ⋆ di mana subskrip δ ⋆ sendiri memiliki superskrip ⋆ .

Bagi para ahli yang bekerja secara ekstensif dengan manipulasi aljabar, sistem tata nama yang rumit ini mungkin tidak dapat dihindari dan bahkan optimal. Tetapi ketika para ahli menerapkan ide-ide mereka sebagai program komputer, mereka harus menggunakan pendekatan yang jauh lebih biasa untuk menamai sesuatu: string karakter seperti Reynolds.

Mungkin ada keuntungan dari pendekatan duniawi. Dalam buku ini, kita akan bekerja dengan kalkulus, yang digunakan di sejumlah bidang. Kita akan menggambar contoh dari banyak bidang. Jadi kami ingin kuantitas kami memiliki nama yang mudah dikenali dan diingat. Hal-hal seperti incomedan blood_pressuredan farm_area. Jika kita mencoba untuk tetap pada xdan ydan z, notasi komputer kita akan menjadi tidak dapat dipahami oleh pembaca manusia dan dengan demikian mudah mengalami kesalahan.

1.2 Fungsi bahasa R R, seperti kebanyakan bahasa komputer, memiliki konstruksi pemrograman untuk mewakili operasi yang mengambil satu atau lebih input dan menghasilkan output. Di R, ini disebut “fungsi.” Di R, semua yang Anda lakukan melibatkan fungsi, baik secara eksplisit maupun implisit.

Mari kita lihat versi R dari fungsi matematika, eksponensial. Fungsi diberi nama expdan kita dapat melihat pemrograman yang dikandungnya

exp ## function (x) .Primitive(“exp”) Tidak banyak di sana, kecuali notasi komputer. Dalam R, fungsi dapat dibuat dengan kata kunci function. Misalnya, untuk membuat fungsi yang menerjemahkan pendapatan tahunan menjadi pendapatan harian, kita dapat menulis:

as_daily_income <- function(yearly_income) { yearly_income / 365 } Nama yang dipilih untuk fungsi, as_daily_income, adalah arbitrer. Kita bisa saja menamai fungsi itu apa saja. (Ini adalah praktik yang baik untuk memberikan nama fungsi yang mudah dibaca dan ditulis dan mengingatkan Anda tentang fungsi tersebut.)

Setelah kata kunci function, ada sepasang tanda kurung. Di dalam tanda kurung adalah nama yang diberikan ke input ke fungsi. Hanya ada satu masukan as_daily_incomeyang kami panggil yearly_incomehanya untuk membantu mengingatkan kami apa fungsi yang dimaksudkan untuk dilakukan. Tapi kita bisa memanggil input apa saja.

Bagian tersebut function(yearly_income)menentukan bahwa hal yang sedang dibuat akan menjadi a functiondan kita memanggil input yearly_income. Setelah bagian ini muncul badan fungsi. Tubuh berisi ekspresi R yang menentukan perhitungan yang akan dilakukan. Dalam hal ini, ekspresinya sangat sederhana: bagi yearly_incomedengan 365 – jumlah (perkiraan) hari dalam setahun.

Sangat membantu untuk membedakan antara nilai input dan peran yang akan dimainkan input dalam fungsi. Nilai pendapatan tahunan mungkin 61362(dalam, katakanlah, dolar). Untuk berbicara tentang peran yang dimainkan input, kami menggunakan kata argument . Misalnya, kita mungkin mengatakan sesuatu seperti, ” as_yearly_incomeadalah fungsi yang membutuhkan satu argumen.”

Setelah kata kunci functiondan tanda kurung di mana argumen didefinisikan datang sepasang kurung kurawal {dan }berisi beberapa pernyataan R. Pernyataan-pernyataan ini adalah badan fungsi dan berisi instruksi untuk perhitungan yang akan mengubah input menjadi output.

Inilah fitur mengejutkan dari bahasa komputer seperti R … Nama yang diberikan untuk argumen tidak masalah sama sekali, asalkan digunakan secara konsisten di badan fungsi. Jadi programmer mungkin telah menulis fungsi R dengan cara ini:

as_daily_income <- function(x) { x / 365 } atau bahkan

as_daily_income <- function(ghskelw) { ghskelw / 365 } Semua versi as_daily_income()wasiat yang berbeda ini melakukan hal yang persis sama dan digunakan dengan cara yang persis sama, terlepas dari nama yang diberikan untuk argumen tersebut. 2 Seperti ini:

as_daily_income(61362) ## [1] 168.1151 Seringkali, fungsi memiliki lebih dari satu argumen. Nama-nama argumen dicantumkan di antara tanda kurung setelah kata kunci function, seperti ini:

as_daily_income <- function(yearly_income, duration) { yearly_income / duration } Dalam kasus seperti itu, untuk menggunakan fungsi kita harus menyediakan semua argumen. Jadi, dengan definisi dua argumen terbaru dari as_daily_income(), penggunaan berikut menghasilkan pesan kesalahan:

as_daily_income(61362) ## Error in as_daily_income(61362): argument “duration” is missing, with no default Sebagai gantinya, tentukan kedua argumen:

as_daily_income(61362, 365) ## [1] 168.1151 Satu lagi aspek argumen fungsi di R … Argumen apa pun dapat diberi nilai default . Sangat mudah untuk melihat cara kerjanya dengan sebuah contoh:

as_daily_income <- function(yearly_income, duration = 365) { yearly_income / duration } Dengan nilai default untuk durationfungsi dapat digunakan dengan satu atau dua argumen:

as_daily_income(61362) ## [1] 168.1151 as_daily_income(61362, duration = 366) ## [1] 167.6557 Baris kedua adalah perhitungan yang tepat untuk tahun kabisat.

Untuk menutup, mari kembali ke expfungsi, yang merupakan bawaan dari R. Argumen tunggal untuk expdiberi nama xdan isi fungsi agak samar: .Primitive(“exp”).

Akan sering terjadi bahwa fungsi yang kita buat akan memiliki tubuh yang tidak melibatkan ekspresi matematika tradisional seperti x / d . Seperti yang akan Anda lihat di bab-bab selanjutnya, di dunia modern banyak fungsi matematika yang terlalu rumit untuk diwakili oleh notasi aljabar.

Jangan ditunda oleh badan fungsi non-aljabar seperti itu. Pada akhirnya, apa yang perlu Anda ketahui tentang suatu fungsi untuk menggunakannya hanyalah tiga hal:

Apa argumen untuk fungsi dan apa yang mereka perjuangkan. Apa jenis hal yang sedang diproduksi oleh fungsi. Bahwa fungsinya berfungsi seperti yang diiklankan, misalnya menghitung apa yang akan kita tulis secara aljabar sebagai e x . 1.3 Literasi penggunaan argumen Ingatlah bahwa nama-nama yang dipilih oleh pemrogram suatu fungsi bersifat arbitrer. Anda akan menggunakan fungsi dengan cara yang persis sama meskipun namanya berbeda. Demikian pula, saat menggunakan fungsi, Anda dapat memilih sendiri ekspresi apa yang akan menjadi nilai argumen.

Misalnya, Anda ingin menghitung 100 e − 2.5 . Mudah:

100 * exp(-2.5) ## [1] 8.2085 Tapi kemungkinan besar itu − 2.5 dimaksudkan untuk mewakili sesuatu yang lebih umum. Misalnya, mungkin Anda menghitung berapa banyak obat yang masih ada di dalam tubuh sepuluh hari setelah dosis 100 mg diberikan. Akan ada tiga kuantitas yang terlibat bahkan dalam perhitungan sederhana ini: dosis, jumlah waktu sejak dosis diambil, dan apa yang disebut “konstanta waktu” untuk eliminasi obat melalui hati atau mekanisme lain. (Untuk mengikuti contoh ini, Anda tidak perlu mengetahui apa itu konstanta waktu. Tetapi jika Anda tertarik, berikut ini contohnya. Misalkan suatu obat memiliki konstanta waktu 4 hari. Ini berarti bahwa 63% obat akan dieliminasi selama periode 4 hari.)

Dalam menulis perhitungan, ada baiknya untuk memperjelas dan menjelaskan arti dari setiap besaran yang digunakan dalam perhitungan. Jadi, alih-alih 100 * exp(-2.5), Anda mungkin ingin menulis:

dose <- 100 # mg duration <- 10 # days time_constant <- 4 # days dose * exp(- duration / time_constant) ## [1] 8.2085 Lebih baik lagi, Anda bisa mendefinisikan fungsi yang melakukan perhitungan untuk Anda:

drug_remaining <- function(dose, duration, time_constant) { dose * exp(- duration / time_constant) } Kemudian, melakukan perhitungan untuk situasi tertentu yang dijelaskan di atas adalah masalah penggunaan fungsi:

drug_remaining(dose = 100, duration = 10, time_constant = 4) ## [1] 8.2085 Dengan menggunakan nama deskriptif yang baik dan secara eksplisit melabeli argumen yang mana, Anda menghasilkan dokumentasi yang jelas dan terpelajar tentang apa yang ingin Anda lakukan dan bagaimana orang lain, termasuk “Anda di masa depan”, harus mengubah hal-hal untuk mewakili situasi baru.

1.4 Sehubungan dengan … Kami telah menggunakan fungsi R untuk mewakili perhitungan jumlah dari input seperti dosis dan konstanta waktu obat.

Tetapi fungsi R memainkan peran yang jauh lebih besar dari itu. Fungsi digunakan untuk hampir semua hal, mulai dari membaca file data hingga menggambar grafik hingga mengetahui jenis komputer apa yang digunakan. Yang menarik bagi kami di sini adalah penggunaan fungsi untuk mewakili dan mengimplementasikan operasi kalkulus. Operasi ini memiliki nama yang mungkin atau mungkin belum Anda kenal: diferensiasi, integrasi, dll.

Ketika kalkulus atau operasi matematika serupa sedang dilakukan, Anda biasanya harus menentukan variabel atau variabel mana yang sedang dilakukan operasi “sehubungan dengan.” Untuk mengilustrasikannya, pertimbangkan operasi sederhana secara konseptual untuk menggambar grafik suatu fungsi. Lebih khusus, mari kita menggambar grafik berapa banyak obat yang tersisa dalam tubuh sebagai fungsi waktu sejak dosis diberikan. Farmakokinetik dasar dari proses ini dikemas dalam drug_remaining()fungsi. Jadi yang ingin kita lakukan adalah menggambar grafik drug_remaining().

Ingat bahwa drug_remaining()memiliki tiga argumen: dose, duration, dan time_constant. Grafik tertentu yang akan kita gambar menunjukkan sisa obat sebagai fungsi durasi. Artinya, operasi grafik akan sehubungan denganduration . Kami akan mempertimbangkan, katakanlah, dosis 100 mg obat dengan konstanta waktu 4 hari, mungkin dengan melihat interval durasi dari 0 hari hingga 20 hari.

Dalam buku ini, kita akan menggunakan operasi yang disediakan oleh paket mosaicand mosaicCalcuntuk R. Operasi dari paket ini memiliki notasi yang sangat spesifik untuk diekspresikan sehubungan dengan . Notasi tersebut menggunakan karakter tilde, ~. Berikut cara menggambar grafik yang kita inginkan, menggunakan slice_plot()operasi paket:

slice_plot( drug_remaining(dose = 100, time_constant = 4, duration = t) ~ t, domain(t = 0:20))

Grafik yang tepat akan diberi label dengan benar, misalnya sumbu horizontal dengan “Waktu (hari)” dan sumbu vertikal dengan “Sisa obat (mg)”. Anda akan melihat bagaimana melakukannya di bab berikutnya, yang mengeksplorasi operasi grafik fungsi secara lebih rinci.