Displaying 9 results from an estimated 9 matches for "site_no".
2009 Apr 21
3
create objects in a loop and adding sqlQuery content to them
...something as I'm not able to create these
sitex_data objects with
the database values, neither using list nor assign...
Here some code snipplets:
>library (RODBC)
>channel <- odbcConnectExcel2007 ("biomass_data.xlsx")
>site_data <- sqlQuery(channel, "select site_no from [biomass_data
$] group by site_no")
#Here the values I want to loop through
>str(site_data)
'data.frame': 44 obs. of 1 variable:
$ site_no: num 4 7 9 10 15 16 17 18 19 20 ...
#Here my first try [error message on the line within the loop, saying
something like:
# ...
2009 Jun 04
2
wrong labels and colors of points in graph/plot
...ient2$intercept ~ coefficient2$average_height,
main="intercepts ::: height", ylab="intercepts", xlab="average height
per site [cm]", xlim=c(20,3020), ylim=c(-2,5), col=coefficient2$color)
highlight(coefficient2$intercept ~ coefficient2$average_height,
lbls=coefficient2$site_no,col="Red", cex = .6)
plot(coefficient2$intercept ~ coefficient2$average_dbh,
main="intercepts ::: dbh", ylab="intercepts", xlab="average dbh per
site [mm]", xlim=c(-10,360), ylim=c(-2,5),col=coefficient2$color )
highlight(coefficient2$intercept ~ coefficient...
2009 May 19
2
create string of comma-separated content of vector
...comma-separated content of a vector?
I've got the vector i with several numeric values as content:
>str(i)
num 99
and want to create a SQL statement to look like the following where
the part '(2, 4, 6, 7)' should be
the content of the vector i:
select * from [biomass_data$] where site_no in (2, 4, 6, 7)
Here my approach (which doesn't work):
site_all_data <= sqlQuery(channel, "select * from [biomass_data$]
where site_no in (",paste(i,sep=","),") )
sorry for spaming so much today to the mailing list...
-Katharina
--
Time flies like an arrow, f...
2010 May 18
2
Function that is giving me a headache- any help appreciated (automatic read )
...s the Upatoi Creek Near Columbus (ftbn)
#02342500 is the Uchee Creek Near Fort Mitchell (ftbn)
#02203000 is the Canoochee River Near Claxton (ftst)
#02196690 is the Horse Creek Gauge at Clearwater, S.C.
a <- "http://waterdata.usgs.gov/nwis/uv?format=rdb&period="
b <- "&site_no=021973269,02102908,02133500,02341800,02342500,02203000,02196690"
z <- paste(a, days, b, sep="")
L <- readLines(z)
#look for the data with USGS in front of it (this take advantage of
#the agency column)
L.USGS <- grep("^USGS", L, value = TRUE)
DF <- read.table(t...
2009 Jun 01
1
write values of points beside points in plot
...in a plot?
I'm plotting the following:
plot(coefficient$intercept ~ coefficient$average_BM_leaves_needles,
main="intercepts ::: BM_leaves_needles",
ylab="intercepts", xlab="average BM_leaves_needles per site [kg]")
and I want to have the value of coefficient$site_no written beside
each point within the plot (to recognize them)
sorry for troubling you with this question, but I somehow can't find
any information on it
(maybe because I'm not quite sure what to search for)
Thanks a lot,
Katharina
2011 Aug 26
1
Time Series data with data every half hour
I am working with data from the USGS with data every 30 minutes from
4/27/2011 to 8/25/2011.
I am having trouble with setting the frequency.
My R script is below:
> shavers=read.csv("shavers.csv")
> names(shavers)
[1] "agency_cd" "site_no" "datetime" "tz_cd" "Temp"
[6] "X04_00010_cd" "EC25" "X05_00095_cd" "DO"
"X06_00300_cd"
[11] "pH" "X07_00400_cd"
> as.POSIXct(shavers[1,3])
[1] "2...
2009 Oct 06
1
ggplot2 applying a function based on facet
...#02133500 is the Drowning Creek (ftbrbmcm)
#02341800 is the Upatoi Creek Near Columbus (ftbn)
#02342500 is the Uchee Creek Near Fort Mitchell (ftbn)
#02203000 is the Canoochee River Near Claxton (ftst)
a <- "http://waterdata.usgs.gov/nwis/uv?format=rdb&period="
b <- "&site_no=021973269,02102908,02133500,02341800,02342500,02203000"
z <- paste(a, days, b, sep="")
L <- readLines(z)
#look for the data with USGS in front of it (this take advantage of
the agency column)
L.USGS <- grep("^USGS", L, value = TRUE)
DF <- read.table(textConnect...
2012 Jan 11
2
Vegan(ordistep) error: Error in if (aod[1, 5] <= Pin) { : missing value where TRUE/FALSE needed
...dAgeConnectandGEOGraphy2.R")
print(myrda1)
Call: rda(formula = mygenind@tab ~ pcnmTRE_25_100_CS25 + pcnmTRE_25_10_CS25
+ pcnmTRE_25_2_CS25 +
pcnmTRE_25_5_CS25 + mydata$TreeCov + mydata$Hab_Config +
pcnmEYR_EO_100_CS25 +
pcnmEYR_EO_5000_CS25 + pcnmEYR_TH_10_CS25 + pcnmEYR_TH_2_CS25 +
mydata$Site_No + mydata$Landscape
+ Condition(pcnmCS_NULL + mydata$LAT.x + mydata$LONG.x), na.action =
"na.omit")
Inertia Proportion Rank
Total 1.8110 1.0000
Conditional 0.8681 0.4793 32
Constrained 0.0000 0.0000 0
Unconstrained 0.9429 0.5207 29
Inert...
2012 Jan 10
0
Error message in vegan ordistep
...dGEOGraphy2.R")
>
> print(myrda1)
Call: rda(formula = mygenind at tab ~ pcnmTRE_25_100_CS25 + pcnmTRE_25_10_CS25 + pcnmTRE_25_2_CS25 +
pcnmTRE_25_5_CS25 + mydata$TreeCov + mydata$Hab_Config + pcnmEYR_EO_100_CS25 +
pcnmEYR_EO_5000_CS25 + pcnmEYR_TH_10_CS25 + pcnmEYR_TH_2_CS25 + mydata$Site_No + mydata$Landscape
+ Condition(pcnmCS_NULL + mydata$LAT.x + mydata$LONG.x), na.action = "na.omit")
Inertia Proportion Rank
Total 1.8110 1.0000
Conditional 0.8681 0.4793 32
Constrained 0.0000 0.0000 0
Unconstrained 0.9429 0.5207 29...