Email : natalieekaren@gmail.com
RPubs : https://rpubs.com/karennatalie/
Jurusan : Statistika
Address : ARA Center, Matana University Tower
Jl. CBD Barat Kav, RT.1, Curug Sangereng, Kelapa Dua, Tangerang, Banten 15810.
contoh :
integrand <- function(x) {2/((x+4)*sqrt(x))}
integrate(integrand, lower = 0, upper = Inf)## 3.141593 with absolute error < 0.00012
library(mosaicCalc)## Loading required package: mosaicCore
## Registered S3 method overwritten by 'mosaic':
## method from
## fortify.SpatialPolygonsDataFrame ggplot2
##
## Attaching package: 'mosaicCalc'
## The following object is masked from 'package:stats':
##
## D
K = antiD(9*x^2+4 ~ x)
K## function (x, C = 0)
## 3 * x^3 + 4 * x + C
phi = 22/7
Luaslingkaran = function(p,r)
{Luas = phi*r^2
return(cat("Luas Lingkaran :", Luas))
}
Luaslingkaran(22/7,16)## Luas Lingkaran : 804.5714
KelilingLingkaran = function(phi,r)
{keliling = 2*phi*r
return(cat("Keliling Lingkaran :", keliling))
}
KelilingLingkaran(22/7,16)## Keliling Lingkaran : 100.5714
VolumeBola = function(phi,r)
{Volume =round(4/3*phi*r^3, digits = 1)
return(cat("Volume Bola :", Volume))
}
VolumeBola(22/7,16)## Volume Bola : 17164.2
Data frekuensi
IPK = c(3.7,2.6,4,3.5,3,3.1)
Frek = c(3,5,1,7,4,2)
Frekuensi = round(Frek, digit=0)
N = data.frame(IPK,
Frekuensi)
N| IPK | Frekuensi |
|---|---|
| 3.7 | 3 |
| 2.6 | 5 |
| 4.0 | 1 |
| 3.5 | 7 |
| 3.0 | 4 |
| 3.1 | 2 |
ekstrak.data = function (x,frek)
{
min = min(x)
max = max(x)
rata= round(sum(x*frek)/sum(frek), digits=2)
return(cat(c("minimum =", min,
"maksimum =", max,
"rata-rata =", rata
)))
}
ekstrak.data(N$IPK,N$Frekuensi)## minimum = 2.6 maksimum = 4 rata-rata = 3.22
ekstrak.data = function(x, Fr)
{
median = median(x)
return(cat(c("median=", median)))
}
ekstrak.data(N$IPK, N$Frekuensi)## median= 3.3
Mode <- function(IPK)
{
mode <- unique(IPK)
mode[which.max(tabulate(match(IPK, mode)))]
}
Mode(IPK)## [1] 3.7
Variansi <- function(IPK)
{variansi = var(IPK)
return(cat(c("Variansi=",variansi)))
}
Variansi(IPK)## Variansi= 0.261666666666667
StandDev <- function(IPK)
{standev = sd(IPK)
return(cat(c("Standar Deviasi=",standev)))
}
StandDev(IPK)## Standar Deviasi= 0.511533641774094