Dear Everyone,
I 've just begun to use the library ncdf and I would like to compare meteorological observational data with forecast data, so to make verification. The netcdf files I'm using contain data of many different parameters in many different stations. I could read easily that I needed, but naturally I do not need the data of all the stations. On the other hand, the order of the stations is not the same in the observation files and in the forecast files. Let's take that I have a list of those stations (with station numbers) where I would like to make the verification. I read the observations in all the possible stations and I receive a data.frame (first column with the station numbers, second column with - let' say - the temperature data...and naturally we can have more columns with different parameters). I make the same with the forecast data, as I wrote the orders of the station numbers in the two dataset are different, and naturally there are some stations that you can find in one data.frame but not in the other. How could I make (or rewrite) my two data.frames (observation and forecast), where the first coulumn is totally the same as in the station list (even the order of the stations)?? For example, I have this data.frame as obsesrvation: 12866 14.4 12844 14.1 12843 16.5 12860 14.9 12851 9.8 12846 15.3 ...and have this data.frame as forecast: 12830 12.808611 12836 12.725081 12843 15.241580 12844 15.185887 12846 13.723515 12851 8.498717 12860 15.715260 12866 14.262023 12870 12.968392 ...and my list of the stations I needed 12836 12843 12846 12860 12870 Thank you very much for your help or suggestions in advance! Gabriella Csima [hidden email] [[alternative HTML version deleted]] ______________________________________________ [hidden email] 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. |
? merge
> obs <- read.table(textConnection("12866 14.4 + + 12844 14.1 + + 12843 16.5 + 12860 14.9 + 12851 9.8 + + 12846 15.3"), col.names=c('station', 'obs')) > fore <- read.table(textConnection("12830 12.808611 + 12836 12.725081 + 12843 15.241580 + 12844 15.185887 + 12846 13.723515 + 12851 8.498717 + 12860 15.715260 + 12866 14.262023 + 12870 12.968392"), col.names=c('station', 'fore')) > closeAllConnections() > > # use 'merge' to group them > > x <- merge(obs, fore, by="station", all=TRUE) > > x station obs fore 1 12830 NA 12.808611 2 12836 NA 12.725081 3 12843 16.5 15.241580 4 12844 14.1 15.185887 5 12846 15.3 13.723515 6 12851 9.8 8.498717 7 12860 14.9 15.715260 8 12866 14.4 14.262023 9 12870 NA 12.968392 > You can then order by the first column. On Thu, May 20, 2010 at 6:22 AM, Csima Gabriella <[hidden email]> wrote: > Dear Everyone, > > > > I 've just begun to use the library ncdf and I would like to compare meteorological observational data with forecast data, so to make verification. The netcdf files I'm using contain data of many different parameters in many different stations. I could read easily that I needed, but naturally I do not need the data of all the stations. On the other hand, the order of the stations is not the same in the observation files and in the forecast files. > > Let's take that I have a list of those stations (with station numbers) where I would like to make the verification. I read the observations in all the possible stations and I receive a data.frame (first column with the station numbers, second column with - let' say - the temperature data...and naturally we can have more columns with different parameters). I make the same with the forecast data, as I wrote the orders of the station numbers in the two dataset are different, and naturally there are some stations that you can find in one data.frame but not in the other. > > > > How could I make (or rewrite) my two data.frames (observation and forecast), where the first coulumn is totally the same as in the station list (even the order of the stations)?? > > > > For example, I have this data.frame as obsesrvation: > > > > 12866 14.4 > > 12844 14.1 > > 12843 16.5 > 12860 14.9 > 12851 9.8 > > 12846 15.3 > > > > > > ...and have this data.frame as forecast: > > > > 12830 12.808611 > 12836 12.725081 > 12843 15.241580 > 12844 15.185887 > 12846 13.723515 > 12851 8.498717 > 12860 15.715260 > 12866 14.262023 > 12870 12.968392 > > > ...and my list of the stations I needed > > > > 12836 > 12843 > 12846 > 12860 > 12870 > > > > > Thank you very much for your help or suggestions in advance! > > Gabriella Csima > > [hidden email] > > > > [[alternative HTML version deleted]] > > ______________________________________________ > [hidden email] 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. > -- Jim Holtman Cincinnati, OH +1 513 646 9390 What is the problem that you are trying to solve? ______________________________________________ [hidden email] 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. |
Free forum by Nabble | Edit this page |