heute ist der Geburtstag von
Pierre-Simon Laplace (28.03.1749 - 05.03.1827)


Auf das Thema antworten  [ 2 Beiträge ] 
WS08: morgen is es soweit 
Autor Nachricht

Registriert: 10/2007
Beiträge: 24 + 13
Mit Zitat antworten
Beitrag WS08: morgen is es soweit
die Übung sollte fertig sein!

also hier mal der R-Code:

Code:
"ueb" <- function ()
{

   werner_bcd <<- read.table("werner_bcd.txt")
   werner_bcd <<- matrix(scan("werner_bcd.txt",na.strings="*"), ncol=9,byrow=T)
   dimnames(werner_bcd)  <<- list(werner_bcd[,1],
   c("Id","Age","Height","Weight","Brthpill","Cholstrl","Albumin","Calcium","Uricacid"))

   werner_bcd <- as.data.frame(werner_bcd)

   tfvec   <-   !is.na(werner_bcd[,8])
   tfvec   <-   !is.na(werner_bcd[,9]) & tfvec

   Ca <- werner_bcd[tfvec,8]
   Ur <- werner_bcd[tfvec,9]   

   par(mfrow=c(3,1))

# 1. Histogramm mit "eindimensionalen Streuungsdiagramm"

   hist(Ur, col="darkgreen",
      breaks=seq(min(Ur),max(Ur)+3,by=3),
      main="URICACID Verteilung", xlab = "URICACID", ylab="Anzahl"
      )
   rug(Ur, side=1)
   

# 2. Schätzungen



   print(paste("Lokation:"))
   print(summary(Ur))

   print(paste("Streuung:", diff(range(Ur))))
   
   print(paste("Schiefe:", skew(Ur)))

   print(paste("Kurtosis:", kurt(Ur)))

#3 Korrelationskoeffizient
   




   print(paste("Korrelationskoeffizient von Calcium und Uricacid: p=",(cov(Ur,Ca)/(var(Ur)*var(Ca)))))

   plot(Ca,Ur,      
      xlab="Calcium", ylab="Uricacid") # Streuungsdiagramm

   abline(lsfit(Ca,Ur),col="brown")
   legend(86,100,c("Calcium-Uracid", "Regression"))

#4 Boxplot

   boxplot(werner_bcd[,c(2:length(werner_bcd))],
      col="darkgreen",
      outline=T)
   title("Werner Blood Chemistry Data")


}


Habs leider noch nicht geschafft in Latex gscheit das File zu schreiben, bin aber am Arbeiten.

Für die, die in Latex schreiben wollen, die Anleitung:

1. Das angehängte File in \MikTex\tex\latex kopieren
2. Unter MiKTeX Options "Refresh FNDB" drücken
3. in latex \usepackage{sweave} dazuschreiben
4. code schreiben.

wenn R-Code implementiert werden soll geht das folgender maßen:

<<>>=
2+2
@

die Datei muss unter *.Rnw abgespeichert werden

In R dann mit >Sweave("*.Rnw") aufrufen, sweave erstellt automatisch ein *.tex file. Dieses dann wieder in LaTeX ausführen und viel Spaß!

wenn jemand ein fertiges pdf file hat bitte posten... auch wenn dies mit word ist! ich versuch derweil noch in latex was zu schreiben.

lg Peter

p.s.: kann das file leider nicht hochladen. ihr solltet es in eurem R-Ordner unter \R\share unter dem Namen "sweave.sty" finden.

viel spaß!

_________________
Die Mathematik ist die Mutter und Tochter aller Wissenschaften


Mo 27-10-2008 19:24:25
Diesen Beitrag melden
Profil

Registriert: 10/2007
Beiträge: 24 + 13
Mit Zitat antworten
Beitrag Re: morgen is es soweit
hab vergessen die schiefe und kurt funktion zu posten

skew <-function(x)
{
round(sum(((x-mean(x))/sd(x))^3)/length(x),digits=2)
}

kurt <- function(x)
{
round(sum(((x-mean(x))/sd(x))^4)/length(x)-3, digits=2)
}

_________________
Die Mathematik ist die Mutter und Tochter aller Wissenschaften


Di 28-10-2008 12:00:56
Diesen Beitrag melden
Profil
Beiträge der letzten Zeit anzeigen:  Sortiere nach  
Auf das Thema antworten   [ 2 Beiträge ] 


Wer ist online?

Mitglieder in diesem Forum: 0 Mitglieder und 3 Gäste


Du darfst neue Themen in diesem Forum erstellen.
Du darfst Antworten zu Themen in diesem Forum erstellen.
Du darfst deine Beiträge in diesem Forum nicht ändern.
Du darfst deine Beiträge in diesem Forum nicht löschen.
Du darfst keine Dateianhänge in diesem Forum erstellen.

Suche nach:
Gehe zu:  
cron
Powered by phpBB © phpBB Group.  |  Designed by STSoftware for PTF  |  © Czechnology 2007 - 2024  |  Deutsche Übersetzung durch phpBB.de