search for: nh4

...sure if this is an R-users question, but since most of you here are statisticians, I decided to give it a shot. I am using the lm() function in R to fit a dependent variable to a set of 3 to 5 independent variables. For this, I used the following commands: >model1<-lm(function=PBW~SO4+NO3+NH4) Coefficients: (Intercept) SO4 NO3 NH4 0.01323 0.01968 0.01856 NA and >model2<-lm(function=PBW~SO4+NO3+NH4+Na+Cl) Coefficients: (Intercept) SO4 NO3 NH4 Na Cl -0.0006987 -0.0119750 -0.0295042 0.0842989...

...day, I'm trying to get more time series in one plot. As there are bigger differences in values of variables I need logaritmic y axis. The code I use is the following: nvz_3_data <- read.csv('/home/tomas/R_outputs/nvz_3.csv') date <- (nvz_3_data$date) NO3 <- (nvz_3_data$NO3) NH4 <- (nvz_3_data$NH4) date_p <- as.POSIXct(date, "CET") par(mfrow=c(2,1), ylog = TRUE, yaxp = c(0.01, 100, 3)) plot(date_p, NO3, log = "y", type = "l", col = "darkred", main = "NVZ-1", xlab = "time", ylab = "NO3-" ) lines(dat...

...gSQL) # conn becomes an object which contains the DB connection: conn <- dbConnect(PgSQL(), host="localhost", dbname="BVS", user="postgres", password = "*******") query_duj_kal <- dbSendQuery(conn, "select zdroj as well, cas as date, fe2, fe3, mn, nh4, no2, no3, o2, teplota as temperature from analystats where zdroj like 'Dunaj Kalinkovo'") watchem_duj_kal <- dbGetResult(query_duj_kal) My intention with the data is to create a time series by a following script: date <- (watchem_duj_kal$date) NO3 <- (watchem_duj_kal$no3) N...

...ur script into mine like this: - library(Rdbi) library(RdbiPgSQL) conn <- dbConnect(PgSQL(), host="localhost", dbname="BVS", user="postgres", password = "******") query_duj_kal <- dbSendQuery(conn, "select zdroj as well, cas as date, fe2, fe3, mn, nh4, no2, no3, o2, teplota as temperature from analystats where zdroj like 'Dunaj Kalinkovo' order by date") watchem_duj_kal <- dbGetResult(query_duj_kal) date <- (watchem_duj_kal$date) date_p <- strptime(date, "%m-%d-%Y") no3 <- (watchem_duj_kal$no3) nh4 <- (watc...

...uppose that each column in pred has NAs in them. Is there some way to fix this? Thanks > str(pred) 'data.frame':   200 obs. of  13 variables:  $ mnO2: num  9.8 8 11.4 4.8 9 13.1 10.3 10.6 3.4 9.9 ...  $ Cl  : num  60.8 57.8 40 77.4 55.4 ...  $ NO3 : num  6.24 1.29 5.33 2.3 10.42 ...  $ NH4 : num  578 370 346.7 98.2 233.7 ...  $ oPO4: num  105 428.8 125.7 61.2 58.2 ...  $ PO4 : num  170 558.8 187.1 138.7 97.6 ...  $ Chla: num  50 1.3 15.6 1.4 10.5 ...  $ a1  : num  0 1.4 3.3 3.1 9.2 15.1 2.4 18.2 25.4 17 ...  $ a2  : num  0 7.6 53.6 41 2.9 14.6 1.2 1.6 5.4 0 ...  $ a3  : num  0 4.8 1....

How can I read a space-delimited file, where the data values for each case are folded before column 80, and so appear on two lines for each case? The first few cases look like this loc type bio H2S sal Eh7 pH buf P K Ca Mg Na Mn Zn Cu NH4 OI DVEG 676 -610 33 -290 5.00 2.34 20.238 1441.67 2150.00 5169.05 35184.5 14.2857 16.4524 5.02381 59.524 OI DVEG 516 -570 35 -268 4.75 2.66 15.591 1299.19 1844.76 4358.03 28170.4 7.7285 13.9852 4.19019 51.378 OI DVEG 1052 -610 32 -282 4.20 4.18 18.716 1154.27 1750.36 4041.27 26455.0 17.8066 15.327...

...;|'Al'|'CO3'|'HCO3'|'Alk-Tot'|'As'|'Ba'|'Be'|'Bi'|'Ca'|'Cd'|'Cl'|'Co'|'Cr'|'Cu'|'DO'|'Fe'|'Hg'|'K'|'Mg'|'Mn'|'Mo'|'Na'|'NH4'|'NO3-NO2'|'Oil-grease'|'Pb'|'pH'|'Sb'|'SC'|'Se'|'SO4'|'Sr'|'TDS'|'Tl'|'V'|'Zn' 'D-1'|'2007-12-12'|0.000|0.106|1.000|231.000|231.000|0.011|0.000|0.002|0.000|100.000|0.000|1.43...

Hi, my name is Luis, and I have a problem with a dataset. Its name is algae and count the collection of data in a lake and respective proliferation of algae. The parameters that it has are: "mxPH", "mnO2", "Cl", "NO3" "NH4", "oPO4", "PO4", "Chla" and "a1" all numerics. a1 - algae1 If I try to do SVM with him, it doesn't generate a square matrix, why? If I try to do same thing in another dataset as Glass of the package "mlbench" with the data(Glass) that has...

...fferent direction to which I am now facing. Earlier Posting:--- I am sorry if this is more of a stats question than an R-question, but I have found it difficult to get a clear answer by other means. Q. Would it be "wrong" to specify a nested model and retain a common intercept, e.g. lm(NH4 ~ Site/TideCode + 1) I am aware (?) that my Site-coefficients are now calculated relative to my reference Site (treatment.contrasts), *but* that my TideCode levels now relate to their reference level within Site. Is that correct? Thank you in advance for help. Regards, Mark Difford. Mark D...

...are defined from the same values, -9999. I used the code below. Below the code are parts of the results. Cheers, Ulrich Data import: chemicS <- read.table("ChemieUlli_4_Quellen.csv", header = TRUE, sep = ",",na.strings = "-9999") Count EC NO3 NO2 NH4 3504 630.0000 33.00 0.001 0.01 3505 NA 26.66 <NA> <NA> 3506 NA 0.72 <NA> <NA> 3507 NA NA <NA> <NA> 3508 NA NA <NA> <NA> 3509 NA NA <NA> <NA> 3510 1210.0000 14.00...

Hi everyone! I tried to look up in the previous topics whether there was something similar to my question but I think there's nothing corresponding... Anyway! Here is my problem: I have a plot (qplot()) in which I want to insert axis labels with the following text: kg NH4+ N-equiv. * ha-1 yr-1 (well... 4 smaller and down and + as a power (the same for the -1) but I am sure you understood what I meant). So, in my code, I typed : qplot(...., xlab="Amount of N-fertilizer" expression(kg NH[4]^+ N-equiv. * ha^-1 yr^-1)) but I don't get what I expect......

...h R and it's the first time I'm using the R-help list... I hope I'm in the right place, if not: Sorry!! I need to do non linear regressions on a data set which columns are: "river.name" "Portata" "PTG.P" "PO4.P" "NT.N" "NH4.N" "NO3.N" "BOD5" "SiO2" I need to predict every variable (PTG, PO4, NT, ..., which are concentration of substances in water) starting from the "Portata" variable (which is the water flow) The functions that I'm using are: > fz1<-fun...

Dear R-help, I am sorry if this is more of a stats question than an R-question, but I have found it difficult to get a clear answer by other means. Q. Would it be "wrong" to specify a nested model and retain a common intercept, e.g. lm(NH4 ~ Site/TideCode + 1) I am aware (?) that my Site-coefficients are now calculated relative to my reference Site (treatment.contrasts), *but* that my TideCode levels now relate to their reference level within Site. Is that correct? Thank you in advance for help. Regards, Mark Difford. Mark D...

Dear List-Members, Is the application of multiple comparison procedures (using the multcomp package) to the output from a rank-based ANOVA straightforward, or do I need to take heed ? That is, is it as simple as: glht( aov(rank(NH4) ~ Site, data=mydat), linfct=mcp(Site="Tukey") ) Thanks in advance for your help. Regards, Mark Difford. -- View this message in context: http://www.nabble.com/Multiple-comparisons--rank-based-anova-tf4516025.html#a12881037 Sent from the R help mailing list archive at Nabble.com.

..."CrTot", "CuDis", "CuTot", "DO", "FeDis", "FeTot", "HgDis", "HgTot", "MgDis", "MgTot", "MnDis", "MnTot", "MoDis", "MoTot", "N", "NH4", "NaDis", "NaTot", "NiDis", "NiTot", "PDis", "PTot", "PbDis", "PbTot", "SC", "SO4", "SbDis", "SbTot", "SeDis", "SeTot", "SrDis", &quo...

..."CrTot", "CuDis", "CuTot", "DO", "FeDis", "FeTot", "HgDis", "HgTot", "MgDis", "MgTot", "MnDis", "MnTot", "MoDis", "MoTot", "N", "NH4", "NaDis", "NaTot", "NiDis", "NiTot", "PDis", "PTot", "PbDis", "PbTot", "SC", "SO4", "SbDis", "SbTot", "SeDis", "SeTot", "SrDis", &quo...