f <- (structure(list(X = structure(96:97, .Label = c("119DAmm",
"119DN",
"119DNN", "119DO", "119DOC", "119Flow",
"119Nit", "119ON", "119OPhos",
"119OrgP", "119Phos", "119TKN",
"119TOC", "148DAmm", "148DN",
"148DNN", "148DO", "148DOC", "148Flow",
"148Nit", "148ON", "148OPhos",
"148OrgP", "148Phos", "148TKN",
"148TOC", "179DAmm", "179DN",
"179DNN", "179DO", "179DOC", "179Flow",
"179Nit", "179ON", "179OPhos",
"179OrgP", "179Phos", "179TKN",
"179TOC", "185DAmm", "185DN",
"185DNN", "185DO", "185DOC", "185Flow",
"185Nit", "185ON", "185OPhos",
"185OrgP", "185Phos", "185TKN",
"185TOC", "190DAmm", "190DN",
"190DNN", "190DO", "190DOC", "190Flow",
"190Nit", "190ON", "190OPhos",
"190OrgP", "190Phos", "190TKN",
"190TOC", "198DAmm", "198DN",
"198DNN", "198DO", "198DOC", "198Flow",
"198Nit", "198ON", "198OPhos",
"198OrgP", "198Phos", "198TKN",
"198TOC", "202DAmm", "202DN",
"202DNN", "202DO", "202DOC", "202Flow",
"202Nit", "202ON", "202OPhos",
"202OrgP", "202Phos", "202TKN",
"202TOC", "215DAmm", "215DN",
"215DNN", "215DO", "215DOC", "215Flow",
"215Nit", "215ON", "215OPhos",
"215OrgP", "215Phos", "215TKN",
"215TOC", "61DAmm", "61DN", "61DNN",
"61DO", "61DOC", "61Flow", "61Nit",
"61ON", "61OPhos", "61OrgP",
"61Phos", "61TKN", "61TOC", "BCOPhos",
"BCOrgP", "BCPhos", "BCTKN",
"BCTOC", "HCDAmm", "HCDN", "HCDNN",
"HCDO", "HCDOC", "HCFlow",
"HCNit", "HCON", "HCOPhos", "HCOrgP",
"HCPhos", "HCTKN", "HCTOC",
"SCDAmm", "SCDN", "SCDNN", "SCDO",
"SCDOC", "SCFlow", "SCNit",
"SCON", "SCOPhos", "SCOrgP", "SCPhos",
"SCTKN", "SCTOC"), class = "factor"),
RiverMile = structure(c(8L, 8L), .Label = c("119",
"148",
"179", "185", "190", "198",
"202", "215", "61", "BC",
"HC",
"SC"), class = "factor"), Constituent = structure(5:6,
.Label c("DAmm",
"DN", "DNN", "DO", "DOC",
"Flow", "Nit", "ON", "OPhos",
"OrgP",
"Phos", "TKN", "TOC"), class =
"factor"), X2.1.06 = c(2.5,
2494), X4.1.06 = c(2.4, NA), X5.1.06 = c(2.3, NA), X6.1.06 = c(2.1,
NA), X7.1.06 = c(2, NA), X8.1.06 = c(1.9, NA), X9.1.06 = c(2,
NA), X10.1.06 = c(2.1, NA), X11.1.06 = c(2.9, NA), X12.1.06 = c(2.6,
NA), X1.1.07 = c(2.1, 4229.55), X2.1.07 = c(2.1, 1823.5),
X3.1.07 = c(2.8, 3617.5), X4.1.07 = c(2.3, NA), X5.1.07 = c(3.5,
NA), X6.1.07 = c(2.8, 10974.5), X7.1.07 = c(2.5, 9652), X8.1.07 = c(2.4,
6700.5), X9.1.07 = c(2.2, 11438.5), X10.1.07 = structure(c(34L,
33L), .Label = c("-0.0214", "-0.3660",
"#VALUE!", "0.0000",
"0.0010", "0.0040", "0.0060",
"0.0080", "0.0086", "0.0090",
"0.0100", "0.0110", "0.0120",
"0.0140", "0.0220", "0.0230",
"0.0240", "0.0280", "0.0300",
"0.0720", "0.0890", "0.1000",
"0.1200", "0.1280", "0.1400",
"0.1500", "0.1880", "0.1910",
"0.2000", "0.2080", "0.2200",
"0.3900", "1141.0000", "2.2000",
"2.3000", "2.5000", "2.6000",
"2606.0000", "3.0000", "3.4000",
"4.5000", "4.6000", "4710.0000",
"4720.0000", "4729.0000",
"6.9768", "7.7065", "7.7338",
"8.0097"), class = "factor"),
X11.1.07 = c(2.2, 10968), X12.1.07 = c(2.2, 7361.5), X1.1.08 = c(2,
6797)), .Names = c("X", "RiverMile",
"Constituent", "X2.1.06",
"X4.1.06", "X5.1.06", "X6.1.06",
"X7.1.06", "X8.1.06", "X9.1.06",
"X10.1.06", "X11.1.06", "X12.1.06",
"X1.1.07", "X2.1.07", "X3.1.07",
"X4.1.07", "X5.1.07", "X6.1.07",
"X7.1.07", "X8.1.07", "X9.1.07",
"X10.1.07", "X11.1.07", "X12.1.07",
"X1.1.08"), row.names = c(NA,
-2L), class = "data.frame"))
plot(f[2,4:26], f[1,4:26])
#Error in plot.window(...) : need finite 'xlim' values
#In addition: Warning messages:
#1: In min(x) : no non-missing arguments to min; returning Inf
#2: In max(x) : no non-missing arguments to max; returning -Inf
#I want a scatterplot like what I thought would happen with the plot command
and the rows and columns that contain the data. If there is a NA in either
case there is no point . plotted. This is how plot works with the data in
columns, I think. This is a small example of a much larger data set are
there any suggestions. What am I missing?
thanks
stephen
--
Let's not spend our time and resources thinking about things that are so
little or so large that all they really do for us is puff us up and make us
feel like gods. We are mammals, and have not exhausted the annoying little
problems of being mammals.
-K. Mullis
[[alternative HTML version deleted]]
Try: x <- unlist(f[2, 4:26]) y <- unlist(f[1, 4:26]) plot(x, y) On Tue, May 6, 2008 at 12:43 PM, stephen sefick <ssefick at gmail.com> wrote:> f <- (structure(list(X = structure(96:97, .Label = c("119DAmm", "119DN", > "119DNN", "119DO", "119DOC", "119Flow", "119Nit", "119ON", "119OPhos", > "119OrgP", "119Phos", "119TKN", "119TOC", "148DAmm", "148DN", > "148DNN", "148DO", "148DOC", "148Flow", "148Nit", "148ON", "148OPhos", > "148OrgP", "148Phos", "148TKN", "148TOC", "179DAmm", "179DN", > "179DNN", "179DO", "179DOC", "179Flow", "179Nit", "179ON", "179OPhos", > "179OrgP", "179Phos", "179TKN", "179TOC", "185DAmm", "185DN", > "185DNN", "185DO", "185DOC", "185Flow", "185Nit", "185ON", "185OPhos", > "185OrgP", "185Phos", "185TKN", "185TOC", "190DAmm", "190DN", > "190DNN", "190DO", "190DOC", "190Flow", "190Nit", "190ON", "190OPhos", > "190OrgP", "190Phos", "190TKN", "190TOC", "198DAmm", "198DN", > "198DNN", "198DO", "198DOC", "198Flow", "198Nit", "198ON", "198OPhos", > "198OrgP", "198Phos", "198TKN", "198TOC", "202DAmm", "202DN", > "202DNN", "202DO", "202DOC", "202Flow", "202Nit", "202ON", "202OPhos", > "202OrgP", "202Phos", "202TKN", "202TOC", "215DAmm", "215DN", > "215DNN", "215DO", "215DOC", "215Flow", "215Nit", "215ON", "215OPhos", > "215OrgP", "215Phos", "215TKN", "215TOC", "61DAmm", "61DN", "61DNN", > "61DO", "61DOC", "61Flow", "61Nit", "61ON", "61OPhos", "61OrgP", > "61Phos", "61TKN", "61TOC", "BCOPhos", "BCOrgP", "BCPhos", "BCTKN", > "BCTOC", "HCDAmm", "HCDN", "HCDNN", "HCDO", "HCDOC", "HCFlow", > "HCNit", "HCON", "HCOPhos", "HCOrgP", "HCPhos", "HCTKN", "HCTOC", > "SCDAmm", "SCDN", "SCDNN", "SCDO", "SCDOC", "SCFlow", "SCNit", > "SCON", "SCOPhos", "SCOrgP", "SCPhos", "SCTKN", "SCTOC"), class = "factor"), > > RiverMile = structure(c(8L, 8L), .Label = c("119", "148", > "179", "185", "190", "198", "202", "215", "61", "BC", "HC", > "SC"), class = "factor"), Constituent = structure(5:6, .Label > c("DAmm", > "DN", "DNN", "DO", "DOC", "Flow", "Nit", "ON", "OPhos", "OrgP", > "Phos", "TKN", "TOC"), class = "factor"), X2.1.06 = c(2.5, > 2494), X4.1.06 = c(2.4, NA), X5.1.06 = c(2.3, NA), X6.1.06 = c(2.1, > NA), X7.1.06 = c(2, NA), X8.1.06 = c(1.9, NA), X9.1.06 = c(2, > NA), X10.1.06 = c(2.1, NA), X11.1.06 = c(2.9, NA), X12.1.06 = c(2.6, > NA), X1.1.07 = c(2.1, 4229.55), X2.1.07 = c(2.1, 1823.5), > X3.1.07 = c(2.8, 3617.5), X4.1.07 = c(2.3, NA), X5.1.07 = c(3.5, > NA), X6.1.07 = c(2.8, 10974.5), X7.1.07 = c(2.5, 9652), X8.1.07 = c(2.4, > > 6700.5), X9.1.07 = c(2.2, 11438.5), X10.1.07 = structure(c(34L, > 33L), .Label = c("-0.0214", "-0.3660", "#VALUE!", "0.0000", > "0.0010", "0.0040", "0.0060", "0.0080", "0.0086", "0.0090", > "0.0100", "0.0110", "0.0120", "0.0140", "0.0220", "0.0230", > "0.0240", "0.0280", "0.0300", "0.0720", "0.0890", "0.1000", > "0.1200", "0.1280", "0.1400", "0.1500", "0.1880", "0.1910", > "0.2000", "0.2080", "0.2200", "0.3900", "1141.0000", "2.2000", > "2.3000", "2.5000", "2.6000", "2606.0000", "3.0000", "3.4000", > "4.5000", "4.6000", "4710.0000", "4720.0000", "4729.0000", > "6.9768", "7.7065", "7.7338", "8.0097"), class = "factor"), > X11.1.07 = c(2.2, 10968), X12.1.07 = c(2.2, 7361.5), X1.1.08 = c(2, > 6797)), .Names = c("X", "RiverMile", "Constituent", "X2.1.06", > "X4.1.06", "X5.1.06", "X6.1.06", "X7.1.06", "X8.1.06", "X9.1.06", > "X10.1.06", "X11.1.06", "X12.1.06", "X1.1.07", "X2.1.07", "X3.1.07", > "X4.1.07", "X5.1.07", "X6.1.07", "X7.1.07", "X8.1.07", "X9.1.07", > "X10.1.07", "X11.1.07", "X12.1.07", "X1.1.08"), row.names = c(NA, > -2L), class = "data.frame")) > > plot(f[2,4:26], f[1,4:26]) > > #Error in plot.window(...) : need finite 'xlim' values > #In addition: Warning messages: > #1: In min(x) : no non-missing arguments to min; returning Inf > #2: In max(x) : no non-missing arguments to max; returning -Inf > > #I want a scatterplot like what I thought would happen with the plot command > and the rows and columns that contain the data. If there is a NA in either > case there is no point . plotted. This is how plot works with the data in > columns, I think. This is a small example of a much larger data set are > there any suggestions. What am I missing? > > thanks > > stephen > -- > Let's not spend our time and resources thinking about things that are so > little or so large that all they really do for us is puff us up and make us > feel like gods. We are mammals, and have not exhausted the annoying little > problems of being mammals. > > -K. Mullis > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at 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. >
cast the two vectors as.matrix-- see here: plot(as.matrix(f[2,4:26]), as.matrix(f[1,4:26])) y stephen sefick wrote:> > f <- (structure(list(X = structure(96:97, .Label = c("119DAmm", "119DN", > "119DNN", "119DO", "119DOC", "119Flow", "119Nit", "119ON", "119OPhos", > "119OrgP", "119Phos", "119TKN", "119TOC", "148DAmm", "148DN", > "148DNN", "148DO", "148DOC", "148Flow", "148Nit", "148ON", "148OPhos", > "148OrgP", "148Phos", "148TKN", "148TOC", "179DAmm", "179DN", > "179DNN", "179DO", "179DOC", "179Flow", "179Nit", "179ON", "179OPhos", > "179OrgP", "179Phos", "179TKN", "179TOC", "185DAmm", "185DN", > "185DNN", "185DO", "185DOC", "185Flow", "185Nit", "185ON", "185OPhos", > "185OrgP", "185Phos", "185TKN", "185TOC", "190DAmm", "190DN", > "190DNN", "190DO", "190DOC", "190Flow", "190Nit", "190ON", "190OPhos", > "190OrgP", "190Phos", "190TKN", "190TOC", "198DAmm", "198DN", > "198DNN", "198DO", "198DOC", "198Flow", "198Nit", "198ON", "198OPhos", > "198OrgP", "198Phos", "198TKN", "198TOC", "202DAmm", "202DN", > "202DNN", "202DO", "202DOC", "202Flow", "202Nit", "202ON", "202OPhos", > "202OrgP", "202Phos", "202TKN", "202TOC", "215DAmm", "215DN", > "215DNN", "215DO", "215DOC", "215Flow", "215Nit", "215ON", "215OPhos", > "215OrgP", "215Phos", "215TKN", "215TOC", "61DAmm", "61DN", "61DNN", > "61DO", "61DOC", "61Flow", "61Nit", "61ON", "61OPhos", "61OrgP", > "61Phos", "61TKN", "61TOC", "BCOPhos", "BCOrgP", "BCPhos", "BCTKN", > "BCTOC", "HCDAmm", "HCDN", "HCDNN", "HCDO", "HCDOC", "HCFlow", > "HCNit", "HCON", "HCOPhos", "HCOrgP", "HCPhos", "HCTKN", "HCTOC", > "SCDAmm", "SCDN", "SCDNN", "SCDO", "SCDOC", "SCFlow", "SCNit", > "SCON", "SCOPhos", "SCOrgP", "SCPhos", "SCTKN", "SCTOC"), class > "factor"), > > RiverMile = structure(c(8L, 8L), .Label = c("119", "148", > "179", "185", "190", "198", "202", "215", "61", "BC", "HC", > "SC"), class = "factor"), Constituent = structure(5:6, .Label > c("DAmm", > "DN", "DNN", "DO", "DOC", "Flow", "Nit", "ON", "OPhos", "OrgP", > "Phos", "TKN", "TOC"), class = "factor"), X2.1.06 = c(2.5, > 2494), X4.1.06 = c(2.4, NA), X5.1.06 = c(2.3, NA), X6.1.06 = c(2.1, > NA), X7.1.06 = c(2, NA), X8.1.06 = c(1.9, NA), X9.1.06 = c(2, > NA), X10.1.06 = c(2.1, NA), X11.1.06 = c(2.9, NA), X12.1.06 = c(2.6, > NA), X1.1.07 = c(2.1, 4229.55), X2.1.07 = c(2.1, 1823.5), > X3.1.07 = c(2.8, 3617.5), X4.1.07 = c(2.3, NA), X5.1.07 = c(3.5, > NA), X6.1.07 = c(2.8, 10974.5), X7.1.07 = c(2.5, 9652), X8.1.07 > c(2.4, > > 6700.5), X9.1.07 = c(2.2, 11438.5), X10.1.07 = structure(c(34L, > 33L), .Label = c("-0.0214", "-0.3660", "#VALUE!", "0.0000", > "0.0010", "0.0040", "0.0060", "0.0080", "0.0086", "0.0090", > "0.0100", "0.0110", "0.0120", "0.0140", "0.0220", "0.0230", > "0.0240", "0.0280", "0.0300", "0.0720", "0.0890", "0.1000", > "0.1200", "0.1280", "0.1400", "0.1500", "0.1880", "0.1910", > "0.2000", "0.2080", "0.2200", "0.3900", "1141.0000", "2.2000", > "2.3000", "2.5000", "2.6000", "2606.0000", "3.0000", "3.4000", > "4.5000", "4.6000", "4710.0000", "4720.0000", "4729.0000", > "6.9768", "7.7065", "7.7338", "8.0097"), class = "factor"), > X11.1.07 = c(2.2, 10968), X12.1.07 = c(2.2, 7361.5), X1.1.08 = c(2, > 6797)), .Names = c("X", "RiverMile", "Constituent", "X2.1.06", > "X4.1.06", "X5.1.06", "X6.1.06", "X7.1.06", "X8.1.06", "X9.1.06", > "X10.1.06", "X11.1.06", "X12.1.06", "X1.1.07", "X2.1.07", "X3.1.07", > "X4.1.07", "X5.1.07", "X6.1.07", "X7.1.07", "X8.1.07", "X9.1.07", > "X10.1.07", "X11.1.07", "X12.1.07", "X1.1.08"), row.names = c(NA, > -2L), class = "data.frame")) > > plot(f[2,4:26], f[1,4:26]) > > #Error in plot.window(...) : need finite 'xlim' values > #In addition: Warning messages: > #1: In min(x) : no non-missing arguments to min; returning Inf > #2: In max(x) : no non-missing arguments to max; returning -Inf > > #I want a scatterplot like what I thought would happen with the plot > command > and the rows and columns that contain the data. If there is a NA in > either > case there is no point . plotted. This is how plot works with the data > in > columns, I think. This is a small example of a much larger data set are > there any suggestions. What am I missing? > > thanks > > stephen > -- > Let's not spend our time and resources thinking about things that are so > little or so large that all they really do for us is puff us up and make > us > feel like gods. We are mammals, and have not exhausted the annoying little > problems of being mammals. > > -K. Mullis > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at 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. > >----- Yasir H. Kaheil, Ph.D. Catchment Research Facility The University of Western Ontario -- View this message in context: http://www.nabble.com/General-Plotting-Question-tp17087346p17087502.html Sent from the R help mailing list archive at Nabble.com.
Try also: x1=as.numeric(f[2,4:26]) y1=as.numeric(f[1,4:26]) plot(x1,y1) HTH, Jorge On Tue, May 6, 2008 at 12:43 PM, stephen sefick <ssefick@gmail.com> wrote:> f <- (structure(list(X = structure(96:97, .Label = c("119DAmm", "119DN", > "119DNN", "119DO", "119DOC", "119Flow", "119Nit", "119ON", "119OPhos", > "119OrgP", "119Phos", "119TKN", "119TOC", "148DAmm", "148DN", > "148DNN", "148DO", "148DOC", "148Flow", "148Nit", "148ON", "148OPhos", > "148OrgP", "148Phos", "148TKN", "148TOC", "179DAmm", "179DN", > "179DNN", "179DO", "179DOC", "179Flow", "179Nit", "179ON", "179OPhos", > "179OrgP", "179Phos", "179TKN", "179TOC", "185DAmm", "185DN", > "185DNN", "185DO", "185DOC", "185Flow", "185Nit", "185ON", "185OPhos", > "185OrgP", "185Phos", "185TKN", "185TOC", "190DAmm", "190DN", > "190DNN", "190DO", "190DOC", "190Flow", "190Nit", "190ON", "190OPhos", > "190OrgP", "190Phos", "190TKN", "190TOC", "198DAmm", "198DN", > "198DNN", "198DO", "198DOC", "198Flow", "198Nit", "198ON", "198OPhos", > "198OrgP", "198Phos", "198TKN", "198TOC", "202DAmm", "202DN", > "202DNN", "202DO", "202DOC", "202Flow", "202Nit", "202ON", "202OPhos", > "202OrgP", "202Phos", "202TKN", "202TOC", "215DAmm", "215DN", > "215DNN", "215DO", "215DOC", "215Flow", "215Nit", "215ON", "215OPhos", > "215OrgP", "215Phos", "215TKN", "215TOC", "61DAmm", "61DN", "61DNN", > "61DO", "61DOC", "61Flow", "61Nit", "61ON", "61OPhos", "61OrgP", > "61Phos", "61TKN", "61TOC", "BCOPhos", "BCOrgP", "BCPhos", "BCTKN", > "BCTOC", "HCDAmm", "HCDN", "HCDNN", "HCDO", "HCDOC", "HCFlow", > "HCNit", "HCON", "HCOPhos", "HCOrgP", "HCPhos", "HCTKN", "HCTOC", > "SCDAmm", "SCDN", "SCDNN", "SCDO", "SCDOC", "SCFlow", "SCNit", > "SCON", "SCOPhos", "SCOrgP", "SCPhos", "SCTKN", "SCTOC"), class > "factor"), > > RiverMile = structure(c(8L, 8L), .Label = c("119", "148", > "179", "185", "190", "198", "202", "215", "61", "BC", "HC", > "SC"), class = "factor"), Constituent = structure(5:6, .Label > c("DAmm", > "DN", "DNN", "DO", "DOC", "Flow", "Nit", "ON", "OPhos", "OrgP", > "Phos", "TKN", "TOC"), class = "factor"), X2.1.06 = c(2.5, > 2494), X4.1.06 = c(2.4, NA), X5.1.06 = c(2.3, NA), X6.1.06 = c(2.1, > NA), X7.1.06 = c(2, NA), X8.1.06 = c(1.9, NA), X9.1.06 = c(2, > NA), X10.1.06 = c(2.1, NA), X11.1.06 = c(2.9, NA), X12.1.06 = c(2.6, > NA), X1.1.07 = c(2.1, 4229.55), X2.1.07 = c(2.1, 1823.5), > X3.1.07 = c(2.8, 3617.5), X4.1.07 = c(2.3, NA), X5.1.07 = c(3.5, > NA), X6.1.07 = c(2.8, 10974.5), X7.1.07 = c(2.5, 9652), X8.1.07 > c(2.4, > > 6700.5), X9.1.07 = c(2.2, 11438.5), X10.1.07 = structure(c(34L, > 33L), .Label = c("-0.0214", "-0.3660", "#VALUE!", "0.0000", > "0.0010", "0.0040", "0.0060", "0.0080", "0.0086", "0.0090", > "0.0100", "0.0110", "0.0120", "0.0140", "0.0220", "0.0230", > "0.0240", "0.0280", "0.0300", "0.0720", "0.0890", "0.1000", > "0.1200", "0.1280", "0.1400", "0.1500", "0.1880", "0.1910", > "0.2000", "0.2080", "0.2200", "0.3900", "1141.0000", "2.2000", > "2.3000", "2.5000", "2.6000", "2606.0000", "3.0000", "3.4000", > "4.5000", "4.6000", "4710.0000", "4720.0000", "4729.0000", > "6.9768", "7.7065", "7.7338", "8.0097"), class = "factor"), > X11.1.07 = c(2.2, 10968), X12.1.07 = c(2.2, 7361.5), X1.1.08 = c(2, > 6797)), .Names = c("X", "RiverMile", "Constituent", "X2.1.06", > "X4.1.06", "X5.1.06", "X6.1.06", "X7.1.06", "X8.1.06", "X9.1.06", > "X10.1.06", "X11.1.06", "X12.1.06", "X1.1.07", "X2.1.07", "X3.1.07", > "X4.1.07", "X5.1.07", "X6.1.07", "X7.1.07", "X8.1.07", "X9.1.07", > "X10.1.07", "X11.1.07", "X12.1.07", "X1.1.08"), row.names = c(NA, > -2L), class = "data.frame")) > > plot(f[2,4:26], f[1,4:26]) > > #Error in plot.window(...) : need finite 'xlim' values > #In addition: Warning messages: > #1: In min(x) : no non-missing arguments to min; returning Inf > #2: In max(x) : no non-missing arguments to max; returning -Inf > > #I want a scatterplot like what I thought would happen with the plot > command > and the rows and columns that contain the data. If there is a NA in > either > case there is no point . plotted. This is how plot works with the data > in > columns, I think. This is a small example of a much larger data set are > there any suggestions. What am I missing? > > thanks > > stephen > -- > Let's not spend our time and resources thinking about things that are so > little or so large that all they really do for us is puff us up and make > us > feel like gods. We are mammals, and have not exhausted the annoying little > problems of being mammals. > > -K. Mullis > > [[alternative HTML version deleted]] > > ______________________________________________ > 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. >[[alternative HTML version deleted]]
On 5/6/2008 12:43 PM, stephen sefick wrote:> f <- (structure(list(X = structure(96:97, .Label = c("119DAmm", "119DN", > "119DNN", "119DO", "119DOC", "119Flow", "119Nit", "119ON", "119OPhos", > "119OrgP", "119Phos", "119TKN", "119TOC", "148DAmm", "148DN", > "148DNN", "148DO", "148DOC", "148Flow", "148Nit", "148ON", "148OPhos", > "148OrgP", "148Phos", "148TKN", "148TOC", "179DAmm", "179DN", > "179DNN", "179DO", "179DOC", "179Flow", "179Nit", "179ON", "179OPhos", > "179OrgP", "179Phos", "179TKN", "179TOC", "185DAmm", "185DN", > "185DNN", "185DO", "185DOC", "185Flow", "185Nit", "185ON", "185OPhos", > "185OrgP", "185Phos", "185TKN", "185TOC", "190DAmm", "190DN", > "190DNN", "190DO", "190DOC", "190Flow", "190Nit", "190ON", "190OPhos", > "190OrgP", "190Phos", "190TKN", "190TOC", "198DAmm", "198DN", > "198DNN", "198DO", "198DOC", "198Flow", "198Nit", "198ON", "198OPhos", > "198OrgP", "198Phos", "198TKN", "198TOC", "202DAmm", "202DN", > "202DNN", "202DO", "202DOC", "202Flow", "202Nit", "202ON", "202OPhos", > "202OrgP", "202Phos", "202TKN", "202TOC", "215DAmm", "215DN", > "215DNN", "215DO", "215DOC", "215Flow", "215Nit", "215ON", "215OPhos", > "215OrgP", "215Phos", "215TKN", "215TOC", "61DAmm", "61DN", "61DNN", > "61DO", "61DOC", "61Flow", "61Nit", "61ON", "61OPhos", "61OrgP", > "61Phos", "61TKN", "61TOC", "BCOPhos", "BCOrgP", "BCPhos", "BCTKN", > "BCTOC", "HCDAmm", "HCDN", "HCDNN", "HCDO", "HCDOC", "HCFlow", > "HCNit", "HCON", "HCOPhos", "HCOrgP", "HCPhos", "HCTKN", "HCTOC", > "SCDAmm", "SCDN", "SCDNN", "SCDO", "SCDOC", "SCFlow", "SCNit", > "SCON", "SCOPhos", "SCOrgP", "SCPhos", "SCTKN", "SCTOC"), class = "factor"), > > RiverMile = structure(c(8L, 8L), .Label = c("119", "148", > "179", "185", "190", "198", "202", "215", "61", "BC", "HC", > "SC"), class = "factor"), Constituent = structure(5:6, .Label > c("DAmm", > "DN", "DNN", "DO", "DOC", "Flow", "Nit", "ON", "OPhos", "OrgP", > "Phos", "TKN", "TOC"), class = "factor"), X2.1.06 = c(2.5, > 2494), X4.1.06 = c(2.4, NA), X5.1.06 = c(2.3, NA), X6.1.06 = c(2.1, > NA), X7.1.06 = c(2, NA), X8.1.06 = c(1.9, NA), X9.1.06 = c(2, > NA), X10.1.06 = c(2.1, NA), X11.1.06 = c(2.9, NA), X12.1.06 = c(2.6, > NA), X1.1.07 = c(2.1, 4229.55), X2.1.07 = c(2.1, 1823.5), > X3.1.07 = c(2.8, 3617.5), X4.1.07 = c(2.3, NA), X5.1.07 = c(3.5, > NA), X6.1.07 = c(2.8, 10974.5), X7.1.07 = c(2.5, 9652), X8.1.07 = c(2.4, > > 6700.5), X9.1.07 = c(2.2, 11438.5), X10.1.07 = structure(c(34L, > 33L), .Label = c("-0.0214", "-0.3660", "#VALUE!", "0.0000", > "0.0010", "0.0040", "0.0060", "0.0080", "0.0086", "0.0090", > "0.0100", "0.0110", "0.0120", "0.0140", "0.0220", "0.0230", > "0.0240", "0.0280", "0.0300", "0.0720", "0.0890", "0.1000", > "0.1200", "0.1280", "0.1400", "0.1500", "0.1880", "0.1910", > "0.2000", "0.2080", "0.2200", "0.3900", "1141.0000", "2.2000", > "2.3000", "2.5000", "2.6000", "2606.0000", "3.0000", "3.4000", > "4.5000", "4.6000", "4710.0000", "4720.0000", "4729.0000", > "6.9768", "7.7065", "7.7338", "8.0097"), class = "factor"), > X11.1.07 = c(2.2, 10968), X12.1.07 = c(2.2, 7361.5), X1.1.08 = c(2, > 6797)), .Names = c("X", "RiverMile", "Constituent", "X2.1.06", > "X4.1.06", "X5.1.06", "X6.1.06", "X7.1.06", "X8.1.06", "X9.1.06", > "X10.1.06", "X11.1.06", "X12.1.06", "X1.1.07", "X2.1.07", "X3.1.07", > "X4.1.07", "X5.1.07", "X6.1.07", "X7.1.07", "X8.1.07", "X9.1.07", > "X10.1.07", "X11.1.07", "X12.1.07", "X1.1.08"), row.names = c(NA, > -2L), class = "data.frame")) > > plot(f[2,4:26], f[1,4:26]) > > #Error in plot.window(...) : need finite 'xlim' values > #In addition: Warning messages: > #1: In min(x) : no non-missing arguments to min; returning Inf > #2: In max(x) : no non-missing arguments to max; returning -Inf > > #I want a scatterplot like what I thought would happen with the plot command > and the rows and columns that contain the data. If there is a NA in either > case there is no point . plotted. This is how plot works with the data in > columns, I think. This is a small example of a much larger data set are > there any suggestions. What am I missing?You need to look closely at your data. It's not all numeric; column 23 is a factor, because of that "#VALUE!" entry. Duncan Murdoch