projects / valse.git / blame
| Commit | Line | Data |
|---|---|---|
| 39046da6 BA |
1 | vec_bin = function(X,r) |
| 2 | { | |
| 3 | Z = c() | |
| 4 | indice = c() | |
| 5 | ||
| 6 | j = 1 | |
| 7 | for(i in 1:length(X)) | |
| 8 | { | |
| 9 | if(X[i] == r) | |
| 10 | { | |
| 11 | Z[i] = 1 | |
| 12 | indice[j] = i | |
| 13 | j=j+1 | |
| 14 | } else | |
| 15 | Z[i] = 0 | |
| 16 | } | |
| 17 | ||
| 18 | return (list(Z=Z,indice=indice)) | |
| 35b42a4b B |
19 | } |
| 20 | ||
| 39046da6 BA |
21 | selectionmodele = function(vraisemblance) |
| 22 | { | |
| 23 | D = vraimsemblance[,2] | |
| 24 | D1 = unique(D) | |
| 25 | ||
| 26 | indice = rep(1, length(D1)) | |
| 27 | #select argmax MLE | |
| 28 | if (length(D1)>2) | |
| 29 | { | |
| 30 | for (i in 1:length(D1)) | |
| 31 | { | |
| 32 | A = c() | |
| 33 | for (j in 1:length(D)) | |
| 34 | { | |
| 35 | if(D[[j]]==D1[[i]]) | |
| 36 | a = c(a, vraimsemblance[j,1]) | |
| 37 | } | |
| 38 | b = max(a) | |
| 39 |