kommazahlen in script berechnen

[manes]

hallo forum,

für eine kleine statistische übung wollte ich einen einzeiler schreiben, aber unter der bash funktioniert das nur für ganze zahlen und bringt ansonsten falsche ergebnisse oder fehlermeldungen.

Code: Alles auswählen

#!/bin/bash
echo "Wie groß sollte unter bestimmten Bedingungen die Stichprobe der Untersuchung sein? Gib bitte die folgenden Daten ein!"
read -p "Größe der Grundgesamtheit N: " N
read -p "Breite des Konfidenzintervalls e: " e
read -p "Größe des Stichprobenanteils p: " p
read -p "Standardnormalverteilung Z: " Z
echo "Die Stichprobe sollte etwa" $(((((Z*Z)*p*(1-p))/(e*e))/(1+(((((Z*Z)*p*(1-p))/(e*e))-1)/N)))) "Personen umfassen"

offenbar muß ich auf bc zurückgreifen, aber dessen syntax verstehe ich leider garnicht. hat jemand einen tip für mich, wie ich die bc für obenstehende formel benutzen kann?
merci im vorabbereich!
manes

edit: typo

[GoKi]

Auf die Schnelle:

Code: Alles auswählen

#!/bin/bash
echo "Wie groß sollte unter bestimmten Bedingungen die Stichprobe der Untersuchung sein? Gib bitte die folgenden Daten ein!"
read -p "Größe der Grundgesamtheit N: " N
read -p "Breite des Konfidenzintervalls e: " e
read -p "Größe des Stichprobenanteils p: " p
read -p "Standardnormalverteilung Z: " Z
n=$(echo "scale = 7; ((((($Z*$Z)*$p*(1-$p))/($e*$e))/(1+((((($Z*$Z)*$p*(1-$p))/($e*$e))-1)/$N))))" | bc)
echo "Die Stichprobe sollte etwa $n Personen umfassen"

[manes]

danke Goki, es funktioniert. auch wenn ich es auf den ersten blick noch nicht kapiere...
manes

[Meillo]

Bei bc müssen Variablennamen in Kleinbuchstaben sein. Wenn du das berücksichtigst kannst du die Formel 1:1 übernehmen. Allerdings darfst du nicht vergessen bc mit -l zu starten, sonst rechnet er nur in Ganzzahlen (oder gibt nur Ganzzahlen aus -- eines von beiden).

Ansonsten ist awk eine gute Wahl für Gleitkommaarithmetik.

[manes]

Bei bc müssen Variablennamen in Kleinbuchstaben sein. (...) bc mit -l zu starten, sonst rechnet er nur in Ganzzahlen (oder gibt nur Ganzzahlen aus -- eines von beiden).Ansonsten ist awk eine gute Wahl für Gleitkommaarithmetik.

hier irrt der meister! nach GoKis vorlage berechnet bc auch ohne option -l und ohne kleinbuchstaben korrekte ergebnisse und gibt auch nicht nur ganzzahlige ergebnisse aus. und ich hab mir die "softwareschätze" vom ulm.ccc.de angetan und war auch selber angetan, aber bei awk wars mit meiner konzentration vorbei. aber danke auch dafür!
manes

[GoKi]

Das mit den Variablen stimmt schon. Nur sieht bc bei meiner Variante die Variablen nicht mehr. Diese werden durch echo expandiert und als Zeichenkette in bc eingelesen.

Die Option -l hat den Effekt, dass direkt ein scale gesetzt wird (scale = 20), also mit Nachkommazahlen berechnet wird. Für die einfache Formel oben brauchen wir aber ansonsten keine Funktionen aus der Math-Library.

[Meillo]

Bei bc müssen Variablennamen in Kleinbuchstaben sein. (...) bc mit -l zu starten, sonst rechnet er nur in Ganzzahlen (oder gibt nur Ganzzahlen aus -- eines von beiden).Ansonsten ist awk eine gute Wahl für Gleitkommaarithmetik.

hier irrt der meister! nach GoKis vorlage berechnet bc auch ohne option -l und ohne kleinbuchstaben korrekte ergebnisse und gibt auch nicht nur ganzzahlige ergebnisse aus.

Das mit den Gleitkommazahlen liegt an `scale', das ich noch nicht kannte. Hatte bisher immer -l gewählt, dann ist scale automatisch gesetzt.

Bei den Kleinbuchstaben ändert GoKi nichts an der von mir beschriebenen Situation. Wenn du genau hinschaust sorgt er dafür, dass bc gar keine Variablen zu sehen bekommt. Stattdessen werden von der Shell alle Werte direkt eingesetzt und bc berechnet nur noch den Ausdruck (der keine Variablen mehr enthält).

[EDIT: Da war GoKi schon schneller.]

und ich hab mir die "softwareschätze" vom ulm.ccc.de angetan und war auch selber angetan, aber bei awk wars mit meiner konzentration vorbei. aber danke auch dafür!

Freut mich, dass dir der Vortrag gefallen hat.

In awk würde man es beispielsweise so machen:

Code: Alles auswählen

awk 'BEGIN { Z='$Z'; p='$p'; e='$e'; N='$N'; print ((((Z*Z)*p*(1-p))/(e*e))/(1+(((((Z*Z)*p*(1-p))/(e*e))-1)/N))) }'