Pinnwand für den Kurs Wahrscheinlichkeit und Statistik install.packages("languageR", repos ="http://cran.r-project.org") library(languageR) Versuch: vgrenze=function(n,c){ int=dbinom(0:n,n,0.5) mittel= n/2 + 1 interval=1:(n/2) for (i in 1:(n/2)) { indices<- seq(mittel-i, mittel+i) ; interval[i]<- sum(int[indices]) } for(x in 1:(n/2)){ if(interval[x] >= c){ return(x) } } } a) f1 = function(k){dbinom(k,50,1/37)} plot(f1,0,50) b) dbinom(0,20,1/37) c) dhyper(3,4,28,5) d) qgeom(0.9,1/37) #Teil 1 varianz <- function(x){ E <- sum(x)/length(x) E2 <- (x - E)^2 out <- sum(E2)/length(E2) return(out) } #Bzw Stichprobenvarianz (was auch R berechnet): sVarianz <- function(x){ E <- sum(x)/length(x) E2 <- (x - E)^2 out <- sum(E2)/(length(E2)-1) return(out) } standard <- function(x){ return(varianz(x)^(1/2)) } coVarianz <- function(x,y){ E1 <- sum(x)/length(x) E2 <- sum(y)/length(y) F1 <- (x-E1) F2 <- (y-E2) G <- F1*F2 out <- sum(G)/length(G) return(out)q } correl <- function(x,y){ return(coVarianz(x,y)/(sd(x)*sd(y))) } #Teil 2 # P(A|B)=P(A,B)/P(B) library(languageR) #Der Zähler: z = verbs[verbs$RealizationOfRec=="NP"& verbs$AnimacyOfTheme == "inanimate" & verbs$AnimacyOfRec == "animate",] nrow(z) #Der Nenner: n= verbs[ verbs$AnimacyOfTheme == "inanimate" & verbs$AnimacyOfRec == "animate",] nrow(n) condProb =nrow(z)/nrow(n) condProb verbs5 = verbs[exp(verbs$LengthOfTheme)>=5,] nchar(as.character(verbs[2,2])) abstand <- function(n){ return(sum(1:n-1)+1) } verbsVerb5 = verbs[nchar(as.character(verbs$Verb))>=6,] x=6 f <- function(x){ return(x+4) } fak <- function(x){ y=1 for(i in 1:x){ y=y*i } if(x==0){ y=1 } return(y) } ueber <- function(x,y){ return(fa Neue Spalten "Variablen" definieren mit $: verbs$cross = cross