Re: [R] Data Frame Transpose

2010-01-06 Thread Uwe Ligges



On 06.01.2010 03:14, Noli Sicad wrote:

Hi John

Thanks for your reply. I think I was posting properly the problem.

Here are the error, R script and console errors below.

Thanks. Noli

~~~
The error:
~~
Error in data.frame(CROP_ID = x[1, 1], CROPTYPE = x[1, 2], name =
colnames(x)[4:5],  :
  subscript out of bounds



Probably x does not have at least 5 columns ...
Since we do not have the full data, we cannot see what exactly happens.

Best,
Uwe Ligges



~~~

I have a dynamic subscripts for the Period, as result of linear
programming (LP) model run. How I generalise this line. Right now it
has 3 index only.

x01=y[,1], x02=y[,2], x03=y[,3])

This is sample the data.

PERIOD
1
1
1
1
2
2
2
2
3
3
3
4
4
5
5
5
5
5
6
6
6
6
6
6
7
8
9
10
10
10
10



R  script:

harvest.dat- read.dbf('C:\\Down2\\R_forestmgt\\Carbon\\forest_cut_m.dbf')

names(harvest.dat) = c(CROP_ID, CROPTYPE, PERIOD,CUT_AGE, AREA_CUT)

# Transpose 5 columns

fn- function(x) {
  y- t(x[,4:5])
  data.frame( CROP_ID=x[1,1], CROPTYPE=x[1,2], name=colnames(x)[4:5],
x01=y[,1], x02=y[,2], x03=y[,3])
  }

harvest.dat- do.call( rbind,
lapply(split(harvest.dat,list(harvest.dat$CROP_ID,harvest.dat$CROPTYPE)),fn)
)

write.csv(harvest.dat, forest_cut3.csv)
  ~

Scite console with r package
~

Rscript --vanilla --slave 
C:\Down2\R_forestmgt\Carbon\ForestCarbon_1_F_Clean7_transpose.R

[1] C:/Down2/R_forestmgt/Carbon
Loading required package: foreign
Loading required package: sp
Loading required package: methods
Loading required package: lattice
Warning messages:
1: package 'maptools' was built under R version 2.10.1
2: package 'foreign' was built under R version 2.10.1
3: package 'sp' was built under R version 2.10.1
Error in data.frame(CROP_ID = x[1, 1], CROPTYPE = x[1, 2], name =
colnames(x)[4:5],  :
  subscript out of bounds
Calls: do.call -  lapply -  FUN -  data.frame
Execution halted

Exit code: 1Time: 2.128

~

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Data Frame Transpose

2010-01-05 Thread Noli Sicad
Hi,

forests - read.csv(C:\\Down2\\R_forestmgt\\forest_cut-Age.csv)

m - forests

fn - function(x) {
  y - t(x[,2])
  data.frame( Croptype=x[1,1], Period =x[1,2], name=colnames(x)[2],
x01=y[,1])x01=y[,1], x02=y[,2], x03=y[,3] } ---Problem
here

m - do.call( rbind, lapply(split(m,list(m$Period,m$Croptype)),fn) )

m - m[order(m$Period,m$Croptype),]


I think I having a problem in here: x01=y[,1])x01=y[,1], x02=y[,2],
x03=y[,3]. how to address with my data. I have variable Period.

based on this http://www.mail-archive.com/r-h...@stat.math.ethz.ch/msg09264.html

P_ID Croptype  Period  Ini_Age  Area_Cut
83  SORI1   31  528.2465512
84  SORI1   32  74.55179899
85  SORI1   33  72.45778618
86  SORI1   34  139.5272947
82  SORI2   28  1.711642933
83  SORI2   29  2.50071
84  SORI2   30  432.5139327
93  ORM2   35  316.8422545
62  OTRM3   30  64.60526438
82  SORI3   27  26.93674606
3   SORM3   35  223.3658345
82  SORI4   26  2.50071
4   SORM4   34  1008.643
5   OTRI5   25  32.42603214
5   OTRM5   29  65.9031344
5   SORM5   32  223.1489321
5   SORM5   33  72.59203041
5   SORM5   35  222.8402746
6   OTRI6   22  2.49851
6   OTRI6   23  3.374626509
6   OTRI6   24  96.13462257
6   OTRM6   26  830.7463641
6   OTRM6   27  731.6228643
6   OTRM6   28  16.3519762
7   OTRM7   26  1636.5693
8   OTRM8   26  553.0050146
9   OTRM9   26  894.414033
10  OTRM10  24  38.72597099
10  OTRM10  25  308.6452707
10  OTRM10  26  786.1761969
10  SORM10  31  235.8360136

To this.

P_ID Croptype P1P2P3P4   P5P6P7
  P8P9P10
83  SORI31
84  SORI32
85  SORI33
86  SORI34
82  SORI28
83  SORI29
84  SORI30
93  SORM35
62  OTRM30
82  SORI27
3   SORM35
82  SORI26
4   SORM34
5   OTRI25
5   OTRM29
5   SORM32
5   SORM33
5   SORM35
6   OTRI22
6   OTRI23
6   OTRI24
6   OTRM26
6   OTRM27
6   OTRM28
7   OTRM26
8   OTRM26
9   OTRM

Thanks in advance. Noli

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] Data Frame Transpose

2010-01-05 Thread John Kane
Well, if nothing else, you have missing comma.  :)

x01=y[,1]), x01=y[,1], x02=y[,2], x03=y[,3]
  --
 
 fn - function(x) {
   y - t(x[,2])
   data.frame( Croptype=x[1,1], Period =x[1,2],
 name=colnames(x)[2],
 x01=y[,1])x01=y[,1], x02=y[,2], x03=y[,3] }
 ---Problem
 here
 
 m - do.call( rbind,
 lapply(split(m,list(m$Period,m$Croptype)),fn) )
 
 m - m[order(m$Period,m$Croptype),]
 
 
 I think I having a problem in here: x01=y[,1])x01=y[,1],
 x02=y[,2],
 x03=y[,3]. how to address with my data. I have variable
 Period.
 
 based on this 
 http://www.mail-archive.com/r-h...@stat.math.ethz.ch/msg09264.html
 
 P_ID Croptype  Period  Ini_Age  Area_Cut
 83      SORI    1   
    31      528.2465512
 84      SORI    1   
    32      74.55179899
 85      SORI    1   
    33      72.45778618
 86      SORI    1   
    34      139.5272947
 82      SORI    2   
    28      1.711642933
 83      SORI    2   
    29      2.50071
 84      SORI    2   
    30      432.5139327
 93      ORM    2   
    35      316.8422545
 62      OTRM    3   
    30      64.60526438
 82      SORI    3   
    27      26.93674606
 3       SORM    3 
      35      223.3658345
 82      SORI    4   
    26      2.50071
 4       SORM    4 
      34      1008.643
 5       OTRI    5 
      25      32.42603214
 5       OTRM    5 
      29      65.9031344
 5       SORM    5 
      32      223.1489321
 5       SORM    5 
      33      72.59203041
 5       SORM    5 
      35      222.8402746
 6       OTRI    6 
      22      2.49851
 6       OTRI    6 
      23      3.374626509
 6       OTRI    6 
      24      96.13462257
 6       OTRM    6 
      26      830.7463641
 6       OTRM    6 
      27      731.6228643
 6       OTRM    6 
      28      16.3519762
 7       OTRM    7 
      26      1636.5693
 8       OTRM    8 
      26      553.0050146
 9       OTRM    9 
      26      894.414033
 10      OTRM    10   
   24      38.72597099
 10      OTRM    10   
   25      308.6452707
 10      OTRM    10   
   26      786.1761969
 10      SORM    10   
   31      235.8360136
 
 To this.
 
 P_ID Croptype P1        P2   
     P3        P4   
    P5        P6 
   P7
   P8        P9     
   P10
 83      SORI    31
 84      SORI    32
 85      SORI    33
 86      SORI    34
 82      SORI       
     28
 83      SORI       
     29
 84      SORI       
     30
 93      SORM       
     35
 62      OTRM       
             30
 82      SORI       
             27
 3       SORM     
               35
 82      SORI       
                
     26
 4       SORM     
                
       34
 5       OTRI     
                
               25
 5       OTRM     
                
               29
 5       SORM     
                
               32
 5       SORM     
                
               33
 5       SORM     
                
               35
 6       OTRI     
                
                
       22
 6       OTRI     
                
                
       23
 6       OTRI     
                
                
       24
 6       OTRM     
                
                
       26
 6       OTRM     
                
                
       27
 6       OTRM     
                
                
       28
 7       OTRM     
                
                
               26
 8       OTRM     
                
                
                
       26
 9       OTRM
 
 Thanks in advance. Noli
 
 __
 R-help@r-project.org
 mailing list
 https://stat.ethz.ch/mailman/listinfo/r-help
 PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
 and provide commented, minimal, self-contained,
 reproducible code.
 


  __
Yahoo! Canada Toolbar: Search from anywhere on the web, and bookmark your 
favourite sites. Download it now
http://ca.toolbar.yahoo.com.

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] Data Frame Transpose

2010-01-05 Thread Noli Sicad
Hi John

Thanks for your reply. I think I was posting properly the problem.

Here are the error, R script and console errors below.

Thanks. Noli

~~~
The error:
~~
Error in data.frame(CROP_ID = x[1, 1], CROPTYPE = x[1, 2], name =
colnames(x)[4:5],  :
 subscript out of bounds
~~~

I have a dynamic subscripts for the Period, as result of linear
programming (LP) model run. How I generalise this line. Right now it
has 3 index only.

x01=y[,1], x02=y[,2], x03=y[,3])

This is sample the data.

PERIOD
1
1
1
1
2
2
2
2
3
3
3
4
4
5
5
5
5
5
6
6
6
6
6
6
7
8
9
10
10
10
10



R  script:

harvest.dat - read.dbf('C:\\Down2\\R_forestmgt\\Carbon\\forest_cut_m.dbf')

names(harvest.dat) = c(CROP_ID, CROPTYPE, PERIOD,CUT_AGE, AREA_CUT)

# Transpose 5 columns

fn - function(x) {
 y - t(x[,4:5])
 data.frame( CROP_ID=x[1,1], CROPTYPE=x[1,2], name=colnames(x)[4:5],
x01=y[,1], x02=y[,2], x03=y[,3])
 }

harvest.dat - do.call( rbind,
lapply(split(harvest.dat,list(harvest.dat$CROP_ID,harvest.dat$CROPTYPE)),fn)
)

write.csv(harvest.dat, forest_cut3.csv)
 ~

Scite console with r package
~
Rscript --vanilla --slave 
C:\Down2\R_forestmgt\Carbon\ForestCarbon_1_F_Clean7_transpose.R
[1] C:/Down2/R_forestmgt/Carbon
Loading required package: foreign
Loading required package: sp
Loading required package: methods
Loading required package: lattice
Warning messages:
1: package 'maptools' was built under R version 2.10.1
2: package 'foreign' was built under R version 2.10.1
3: package 'sp' was built under R version 2.10.1
Error in data.frame(CROP_ID = x[1, 1], CROPTYPE = x[1, 2], name =
colnames(x)[4:5],  :
 subscript out of bounds
Calls: do.call - lapply - FUN - data.frame
Execution halted
Exit code: 1Time: 2.128
~

__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.