J4Help Trovare combinazione di valori che porti ad una media e DS nota

[Tilion]

Ciao a tutti :D,
ho un problema che devo risolvere il prima possibile, ma non so da dove partire!
Chiedo qui sperando nell'aiuto di qualche anima pia.

Ho un set di 52 valori (sono età di pazienti ospedalizzati). Ora, io ho una data media e deviazione standard di 37 di questi 52 valori, MA NON SO QUALI SIANO. Come faccio a trovare quella combinazione di 37 valori che mi dà esattamente quella media e DS? Idee?

In b4 "no vai via" & "j4s is not your personal army" :lol:

[Recidivo]

No vai viah.

[Frigg]

j4s is not your personal army

[Bicio]

ciao, sei nuovo? :uhm:

[gnappinox1]

Presentati con il tuo vero nome :sisi:

[Skywolf]

Con un programmino fai la prova con tutti i possibili valori finchè non viene fuori il numero giusto.

Il numero di combinazioni è il coefficiente binomiale C(52,37) che vale 52! / (37!15!) (ovvero 4481381406320), auguri vivissimi ^____^

[Raist²]

non ho il pc abbastanza potente

[Moloch]

puoi provare a

[gnappinox1]

42

[Skywolf]

Con un programmino fai la prova con tutti i possibili valori finchè non viene fuori il numero giusto.

Il numero di combinazioni è il coefficiente binomiale C(52,37) che vale 52! / (37!15!) (ovvero 4481381406320), auguri vivissimi ^____^

volevo però aggiungere, per dare una nota d'ottimismo, che in realtà statisticamente non occorre provare tutte le combinazioni ma solo la metà di 4481381406320 !
no fattoriale intended :jfs:

[Fruttolo]

Cazzo sei?
È gheisvillage qui?

[Raist²]

FAI UN SISTEMA!

[Tilion]

Con un programmino fai la prova con tutti i possibili valori finchè non viene fuori il numero giusto.

Il numero di combinazioni è il coefficiente binomiale C(52,37) che vale 52! / (37!15!) (ovvero 4481381406320), auguri vivissimi ^____^

:lol: :lol: :lol: sono giusto un paio di tentativi

[GenghisKhan]

QUI la soluzione

[fulviuz]

grazie a dio non sonasega di matematica

[Sinex/]

Ahhaha

[Ciome]

edit: no, va be, non c'ho cazzi di incazzarmi anche per te

[Mdk]

Rete neurale is the way

[gnappinox1]

Algoritmo genetico :sisi:

[tigerwoods]

Problema che si risolve in due secondi di matlab :caffe:

[GenghisKhan]

+1 :sisi:

[koba44]

Problema che si risolve in due secondi di matlab :caffe:

Ok ma non puoi usare una distribuzione normale.

[tigerwoods]

azzo centra? :uhm:

[koba44]

azzo centra? :uhm:

Dovrai dargli una distribuzione o no, se hai a disposizione solo m e sigma?

Se la distribuzione è normale, potrebbero saltare fuori valori negativi, per esempio...

EDIT: okok, normale troncata. https://en.wikipedia.org/wiki/Truncated_distribution

Con Matlab: https://it.mathworks.com/help/stats/prob.truncatabledistribution.truncate.html?s_tid=g n_loc_drop
alla voce Generate Random Numbers from a Truncated Distribution

[tigerwoods]

clc
close all
clear all

n = 52;
k = 37;

data = zeros(n,2);

%sim data
for i = 1:n
data(i,1) = randi(100,1); %age 1-100 uniform dist
data(i,2) = i;
end
data_index = sort(randperm(n,k)); %randomise k samples from data
avg = mean(data(data_index,1));
sigma2 = var(data(data_index,1));

%compute combinations
n_comb = nchoosek(n,k);
combs = nchoosek(data(:,2),k);
solution = zeros(size(combs,1),k);

%get solutions
tol = 1e-10;
j = 1;
for i = 1:size(combs,1)
test_avg = mean(data(combs(i,:),1));
test_sigma2 = var(data(combs(i,:),1));
if abs(test_avg - avg) < tol && abs(test_sigma2 - sigma2) < tol
solution(j,:) = combs(i,:);
j = j+1;
end
end


provo a farlo partire vediamo se ci mette molto :asd:

[koba44]

Aaah cerchi esattamente 52 valori in brute-force!

Io avevo interpretato che l'OP volesse un generatore di campioni... che poi sarebbe la cosa più logica....

[koba44]

Brute force, eh!... zzzz... asd

[Skywolf]

sbaglio o manca la parte in cui gli si fornisce i 52 campioni di cui l'OP evidentemente già dispone e da cui deve estrarne 37 con media e varianza note? :uhm:

[koba44]

sbaglio o manca la parte in cui gli si fornisce i 52 campioni di cui l'OP evidentemente già dispone e da cui deve estrarne 37 con media e varianza note? :uhm:
Nessuno ha capito cosa precisamente l'OP voglia, nemmeno lui.

[Fruttolo]

Io ho capito che vuole il cazzo.

[Zodd00]

Io ho capito che vuole il cazzo.

Fruttolo the best

[Fruttolo]

52 campioni di grigio.

[koba44]

Tigerwoods, è ancora lì che gira il coso...

[Mad_One]

ma siamo sicuri che ci sia una sola soluzione? :uhm:

[GenghisKhan]

Io sono sicuro che non me ne fotte un cazzo :sisi:

[tigerwoods]

Tigerwoods, è ancora lì che gira il coso...
l'ho stoppato dopo un po' :asd:

ma siamo sicuri che ci sia una sola soluzione? :uhm:

le soluzioni possono essere anche 4481381406320

[koba44]

ma siamo sicuri che ci sia una sola soluzione? :uhm:
Le soluzioni sono infinite.

[tigerwoods]

Le soluzioni sono infinite.
No