Deutsches R-Forum

Gehen wir es zeilenweise durch

s_c_j <- function(j){
Ich vermute, dass s_c_y eine Funktion von j und von X sein soll, da X ja in der Funktion verwendet wird. Daher würde ich hier allgemein s_c_j <- function(j, X) schreiben.

Hier sollte in a schon zehn Ergebnisse zu den zehn Werten in X stehen. Ich vermute, dass die for-Schleife dazu gedacht war, zehnmal was mit den einzelnen Werten von X zu machen. Wenn ja, dann ist sie wahrscheinlich unnötig. Schwer zu sagen, weil der Rest nicht soviel Sinn macht.

Alarm! Hier wird X neu definiert und das ist bestimmt nicht gut. Bis hierhin war X eine Serie von zehn Zufallszahlen, jetzt wird X für den Rest der Funktion die Zahlen von eins bis zehn und wenn wir die Funktion verlassen ist X schon wieder die Serie mit Zufallszahlen. Wenn das so gewollt ist, halte ich das für einen irreführenden Programmierstil, wahrscheinlich ist es aber nicht so gemeint?

Also eine for-Schleife ist dazu da, einen Anweisungsblock mehrfach zu wiederholen. Eine Funktion endet aber beim ersten Aufruf eines return, sodass ein mehrfaches Ausführen von return von vorneherein unmöglich ist.

MathHam hat geschrieben: Di Jan 25, 2022 10:36 amGehen wir es zeilenweise durch

s_c_j <- function(j){
Ich vermute, dass s_c_y eine Funktion von j und von X sein soll, da X ja in der Funktion verwendet wird. Daher würde ich hier allgemein s_c_j <- function(j, X) schreiben.

Ich dachte es wäre hier sinnvoll, nur das j in der Funktion anzugeben da sich das X dann sowieso anpasst, sobald ich den Vektor X ändere.

Hier sollte in a schon zehn Ergebnisse zu den zehn Werten in X stehen. Ich vermute, dass die for-Schleife dazu gedacht war, zehnmal was mit den einzelnen Werten von X zu machen. Wenn ja, dann ist sie wahrscheinlich unnötig. Schwer zu sagen, weil der Rest nicht soviel Sinn macht.


Nein, in a ist das Ergebnis nur eine Zahl. Hier werden die Realisationen von X in die Funktion c_j eingesetzt und aufsummiert. Das Ergebnis wird dann durch die Anzahl von Realisationen geteilt. Also als Antwort erhalte ich hier eine Zahl (mein Code gibt auch die richtige Zahl aus)

X <- rnorm(10)
j <- 1
       
a <- sqrt(2)*cos(j*pi*X)
print(a)

Alarm! Hier wird X neu definiert und das ist bestimmt nicht gut. Bis hierhin war X eine Serie von zehn Zufallszahlen, jetzt wird X für den Rest der Funktion die Zahlen von eins bis zehn und wenn wir die Funktion verlassen ist X schon wieder die Serie mit Zufallszahlen. Wenn das so gewollt ist, halte ich das für einen irreführenden Programmierstil, wahrscheinlich ist es aber nicht so gemeint?


Ich glaube hier liegt mein Problem. Die verschiedenen s_c_j die nun für jedes j separat berechnet wurden sind nun sozusagen ein konstanter Faktor und in das phi_j [...9]

Haben Sie mit den zusätzlichen Informationen eine Ahnung wie ich meine Funktion plotten kann? Es ist nicht so leicht das Thema hier schriftlich zu erläutern...

MathHam hat geschrieben: Di Jan 25, 2022 11:57 am

Code: Alles auswählen

#Beispiel
X<-rnorm(10)

s_c_1<-sqrt(2)*cos(1*pi*X) #alle 10 Realisierungen von s_c_1
s_c_1done<-(sum(s_c_1)/10) #fertiger Wert für s_c_1

s_c_2<-sqrt(2)*cos(2*pi*X) #alle 10 Realisierungen von s_c_2
s_c_2done<-(sum(s_c_2)/10) #fertiger Wert für s_c_2

s_c_3<-sqrt(2)*cos(3*pi*X) #alle 10 Realisierungen von s_c_2
s_c_3done<-(sum(s_c_3)/10) #fertiger Wert für s_c_2

X <- rnorm(10)

s_c <- function(j, v = X) mean(sqrt(2)*cos(j*pi*v))

print(s_c(1))
print(s_c(2))
print(s_c(3))

Code: Alles auswählen

#phi_j<-sqrt(2)*cos(j*pi*x) mit x = x-Achse (keine Zufallsvariable!)

xs <- seq(0,10,.01)
ys <- sapply(xs, \(x) s_c(1, x))
plot(xs, ys, type = "l", main = "j = 1")