3:36 p.m.
Dear all, happy new & greetings from Paris !
I experience some problem with the matrix functions. Returning to an older
script with new data it does not work properly. Well, here is some code:
###### lissage spatial des PPV
list datalist = 13 14 15 16 17 18 19 20 21 22 23 24
list ppv4 = 25 26 27 28
list ppv8 = 29 30 31 32 33 34 35 36
matrix V4 = {ppv4}
matrix V8 = {ppv8}
Everything is ok with the list commands, but when i want to create the
matrices V4 and V8, from lists ppv4 and ppv8 gretl does not do it and says:
Types de données non conformes pour l'opération
Data type non conformant for the operation
I dont really understand what happens because the script worked smoothly before.
The data are declared as cross-section, under 1.8.0 and running under win XP.
Note also that i tried to declare the matrices directly with variables number
with the same results.
Need helped or will have to reverse to an older version !
thanxs
F.
--
Franck Nadaud
Economiste
CIRED
UMR 8568 CNRS - EHESS
45 bis avenue de la Belle Gabrielle
94736 Nogent-sur-Marne Cedex
TEL 33-1-43-94-73-94
FAX: 33-1-43-94-73-70
MOB: 06-07-39-92-75
France
4:11 p.m.
On Mon, 2 Feb 2009, Franck Nadaud wrote:
I experience some problem with the matrix functions. Returning to an
older
script with new data it does not work properly. Well, here is some code:
###### lissage spatial des PPV
list datalist = 13 14 15 16 17 18 19 20 21 22 23 24
list ppv4 = 25 26 27 28
list ppv8 = 29 30 31 32 33 34 35 36
matrix V4 = {ppv4}
matrix V8 = {ppv8}
Everything is ok with the list commands, but when i want to
create the matrices V4 and V8, from lists ppv4 and ppv8 gretl
does not do it and says:
Types de données non conformes pour l'opération
Could you send a copy of your complete script, please? I tried a
simple example of this with current gretl:
<script>
open data4-10
list ppv4 = 1 2 3 4 5
list ppv8 = 4 5 6 7 8 9
matrix V4 = {ppv4}
matrix V8 = {ppv8}
</script>
and it works OK, so I suppose there must be something more
complicated going on. (You would get the "types" error message if
the identifier "V4" had been used for something other than a
matrix at an earlier point in the script.)
Allin Cottrell
2:53 p.m.
New subject: problem with matrix & loops
Dear Allin,
thank you for your response, I just read it today morning and so I could test
the small script you sent and compare with mine.
Well, here are the results:
1) you script works perfectly well;
2) the first part of mine also... now I can define my matrices by named list
assignement and it works ok.
I guess there is some arcane win XP problem around here... gretl does not seem
to be culprit !
But now I have problems with the second segment of my script :
############################################################################
matrix W1=zeros(257,258)
list contigus1 = 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52
matrix VW1={contigus1}
loop i = 1..257 --quiet
scalar row=$i
loop j = 1..15 --quiet
scalar col=VW1[$i,$j]
matrix W1[row,col]=1
end loop
end loop
matrix WL1=W1[,1:257]
delete W1
### calcul de la matrice de lissage
matrix DEGW1=1./sumr(WL1)
matrix INVDEG=I(257).*DEGW1
matrix W1W=WL1*INVDEG
### lissage d'ordre 1 des donnees a l'aide de la matrice de contiguite
matrix LISS1=W1W*X
loop m=1..12
series lw1$m=LISS1[,$m]
end loop
############################################################################
When i run it: gretl says :
gretl version 1.8.0
Current session: 2009/02/03 15:38
? matrix W1=zeros(257,258)
Replaced matrix W1
? list contigus1 = 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52
Replaced list contigus1
? matrix VW1={contigus1}
Replaced matrix VW1
? loop i = 1..257 --quiet
scalar row=$i
loop j = 1..15 --quiet
scalar col=VW1[$i,$j]
matrix W1[row,col]=1
end loop
end loop
Index value 0 is out of bounds
> matrix W1[row,col]=1
> loop j = 1..15
i dont see why he finds an index value 0 while the loops are indexed from 1 to
257 for i and 1 to 15 for j ?
I could continue with the other script but i need these computations for
comparison.
To be brief, the script computes from a data set (matrix X), the matrix of
spatially smoothed data (each observation is replaced by the mean of its
immediate neighbors). A contiguity (neighborhood) matrix is read from the list
of indexes of the neighbors of each information stored in the matrix VW1 =
{contigus1}.
But as i said, i cannot create the matrix. Ooops, yes should have stored it
once and for all, long ago ! But once i have fixed that, i will send the tests
to the community.
thanx for help and cheers.
Franck
3:13 p.m.
New subject: problem with matrix & loops
On Tue, 3 Feb 2009, Franck Nadaud wrote:
thank you for your response, I just read it today morning and so
I could test the small script you sent and compare with mine.
Well, here are the results:
1) you script works perfectly well;
That's good to hear.
2) the first part of mine also... now I can define my matrices
by named list assignement and it works ok.
Hmm, but what accounts for the change in this respect?
I guess there is some arcane win XP problem around here... gretl
does not seem to be culprit !
This is very unlikely to be a problem with the operating system.
I'd still be interested to see the entire script that failed.
But now I have problems with the second segment of my script :
############################################################################
matrix W1=zeros(257,258)
list contigus1 = 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52
matrix VW1={contigus1}
loop i = 1..257 --quiet
scalar row=$i
loop j = 1..15 --quiet
scalar col=VW1[$i,$j]
matrix W1[row,col]=1
end loop
end loop
############################################################################
When i run it: gretl says :
gretl version 1.8.0
Current session: 2009/02/03 15:38
? matrix W1=zeros(257,258)
Replaced matrix W1
? list contigus1 = 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52
Replaced list contigus1
? matrix VW1={contigus1}
Replaced matrix VW1
? loop i = 1..257 --quiet
> scalar row=$i
> loop j = 1..15 --quiet
> scalar col=VW1[$i,$j]
> matrix W1[row,col]=1
> end loop
> end loop
Index value 0 is out of bounds
>> matrix W1[row,col]=1
>> loop j = 1..15
i dont see why he finds an index value 0 while the loops are
indexed from 1 to 257 for i and 1 to 15 for j ?
You're trying to set W1[row,col]. The "row" value is obtained
from the loop index i, and therefore should run from 1 to 257, as
you say. But "col" is obtained as VW1[i,j]. From looking at your
script, we don't know what values this will take on, and it
appears that for some i and j VW1[i,j] = 0.
Allin.
5843
days inactive
5844
days old
3 comments
2 participants
participants (2)
-
Allin Cottrell -
Franck Nadaud