Displaying 17 results from an estimated 17 matches for "po4".
Did you mean:
p4
2013 Apr 12
2
Creating "%d/%m/%Y %H:%M:%S" format from separate date and time columns
...ween
contacts, I need to calculate a contact end time. The contact starts
at the date and time shown in V4 and V5, and lasts for the duration
shown IN SECONDS in V6:
> data2<- read.csv(file=file.choose(), header=F, sep=" ")
> head(data2)
V1 V2 V3 V4 V5 V6
1 3 PO4 CO1 2011-04-29 07:27:21 28
2 3 PO4 CO1 2011-04-24 05:57:39 20
3 3 PO4 CO1 2011-04-14 10:29:49 4
4 3 PO4 CO1 2011-04-16 07:27:31 63
5 3 PO4 CO1 2011-04-18 15:46:20 1
6 3 PO4 CO1 2011-04-18 15:45:57 1
To start with I have tried to make the start data and time into one new column:
startt<...
2011 Mar 18
1
general question about dropping terms of glm model fits
...ar list!
as I am currently helping someone with their statistical analysis of a
count survey, I stumbled upon a very basic question upon model optimization:
when fitting a model like:
model<-lmer(abundance~(x+y+z)^3,family=poisson,...)
in which x,y,z are continuous abiotic parameters such as po4
concentration, no2-concentration, which terms / interaction terms would
you recommend removing FIRST?
the ones of lowest significance (i.e. the ones with highest p-value) OR
the ones with the most complex interaction structure (even though
p-values may be low-ish)?
another question just poppe...
2007 Jun 15
1
interpretation of F-statistics in GAMs
...alues among tests (like to
standardize them by teh sum of F within each test?) I append two summaries
below.
Thanks in advance,
Robert
### example 1 ###
Family: gaussian
Link function: identity
Formula:
dep[sel, i] ~ s(date, k = 3) + s(depth, k = kn) + s(temp, k = kn) +
s(light, k = kn) + s(PO4, k = kn) + s(DIN, k = kn) + s(prop.agpla,
k = kn)
Parametric coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 5.1048 0.0384 132.9 <2e-16 ***
---
Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
App...
2011 Oct 09
1
sapply(pred,cor,y=resp)
...hem. 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.9 18.9 7.5 1.4 3.2 0 2.5 0 ...
$ a4 : num...
2011 Mar 07
4
png inside loop
...(at
least I am not able to explain the behaviour, therefore I assume it to
be a strange behaviour):
attach(water) # I know, this is not recommended
names(water[3:10])
[1] "temp" "pH" "DO" "BOD" "COD" "no3" "no2" "po4"
for (i in names(water)[3:10]){
fname<-paste("Henni/GFX/fem",i,".png",sep="")
mname<-paste("Henni/GFX/mal",i,".png",sep="")
png(fname,1000,1000)
xyplot(N_female~eval(parse(text=i)) |group,xlab=i,ylab="A...
2009 Feb 03
3
non linear regression with nls
Hello,
I'm a beginner with 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...
2009 Oct 26
1
issue with levels of a factor after subsetting
...","Sin?","Atrato","San Juan","Mira","Micay",
"Patia","Canal del Dique","Iscuand?","Guapi"),]
And I plot the two graphs using:
xyplot(Valor~A?o|Var,groups=Rio,data=car4[car4$Var%in%c("NT","PO4","HDD","CTE",
"SST","OCT"),],layout=c(2,3),subscripts=T,scale=list(y=list(relation="free")),type="b")
xyplot(Valor~A?o|Var,groups=Rio,data=car5[car5$Var%in%c("NT","PO4","HDD","CTE",
"SST&q...
2006 May 17
1
Fix for augPred/gsummary problem (nlme library)
...he colon.] Remove these, or change them, using, e.g.:
class(obj$DateTime) <- NULL
class(obj$AnyVariable) <- 'numeric' ## leaves the actual labels/units intact so that you can later restore them.
2) Execute your lme() statement &c on the object, e.g.:
test.1 <- lme(Chla ~ PO4, random=~1|Site, data=obj) ## or simply: lme(obj)
augPred(test.1)
plot(augPred(test.1))
(Note that if you are using a data.frame() as your data object you will need to supply a 'primary' statement to augPred(), e.g. augPred(test.1, primary=~PO4).
Regards,
Mark Difford.
--------------...
2011 Jan 10
3
Help with Data Transformation
...LHR020GW-01E2 N WG INO BRLS NONE E300 CL = 23590.9
2269 LHR020GW-01E2 N WG INO BRLS NONE E300 PO4 ND 50
2270 LHR020GW-01E2 N WG INO BRLS NONE E300 SO4 = 22460
2272 LHR020GW-01E2 N WG MET BRLS...
2011 Feb 27
3
accessing variables inside a function, inside a loop
...v",sep="\t",header=T)
source("functions.R")
attach(water)
names(water):
[1] "site" "sample" "temp" "pH" "DO" "BOD"
[7] "COD" "no3" "no2" "po4" "N_male" "N_female"
vars<-names(water)[3:10]
par(mfrow=c(4,2))
for (i in vars)
{
comp.plot(i,N_male,DATA=water)
}
comp.plot<-function(parameter,y,DATA){ #this is actually in functions.R
x<-DATA[[parameter]]
y<-DATA[[y]]
plot(y~para...
2007 Oct 03
0
datasets
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 the parameters:
&...
2011 Jan 13
0
Help with Data Transformation - RESOLVED
...;HG","HGACIDLAB","HGEXTINO","HGEXTORG","HGNONMOB""HGSEMIMOB","K",
"MEHG","MG","MN","MO","NH3N","NI","NO2N","NO3","NO3N","PB""PO4","S","SB","SE","SO4",
"SOLID","SSC","TL","TOC","V","Zn","ZN"
- These are currently stored in the "PARLABEL" column.
* For each sample ID I would like to create a si...
2009 Feb 10
3
summary of a list
...tList), uppervalues (upperList) and lower values (lowerList).
A sample of the list of function I use in the loop is the following:
FormList <- list(PTG.P ~ fz1(Portata, a, b), PTG.P ~ fz2(Portata, a, b), PTG.P ~ fz3(Portata,a, b, d, e),
PTG.P ~ fz4(Portata, a, b), PTG.P ~ fz5(Portata, a, b, d), PO4.P ~ fz1(Portata, a, b),
PO4.P ~ fz2(Portata, a, b), ...
And the loop I use is:
resultList <- list()
for (i in 1:length(formList))
{
resultList[[i]] <- nls(formList[[i]], data=subset(dati, Fiume=="Laveggio"), start=startList[[i]],
nls.control(maxiter=1000, warnOnly=TRUE), a...
2013 Mar 29
1
Create values based on a table of conditions
..."CO8", "CO9", "PI1", "PI2", "PI3", "PI4", "PI5", "PI6",
"PI7", "PI8", "PO1", "PO10", "PO11", "PO12", "PO13", "PO2", "PO3",
"PO4", "PO5", "PO6", "PO7", "PO8", "PO9", "RD1", "RD2", "WB1"
), class = "factor"), collar = c(36L, 64L, 25L, 28L, 11L, 46L,
58L, 12L, 4L, 52L, 2L, 63L, 14L, 12L, 12L, 57L, 51L, 51L, 24L,
42L, 44L, 25L, 23L, 7...
2013 Apr 26
2
Remove reciprocal data from a grouped animal social contact dataset
...y animals that could have
contacted each other:
> animals
animals
1 CO1
2 CO2
3 CO3
4 CO4
5 CO5
6 CO6
7 CO7
8 CO8
9 CO9
10 CO10
11 CO11
12 CO12
13 CO13
14 CO14
15 CO15
16 CO16
17 CO17
18 PO1
19 PO2
20 PO3
21 PO4
22 PO5
23 PO6
24 PO7
25 PO8
26 PO9
27 PO10
28 PO11
29 PO12
30 PO13
31 PI1
32 PI2
33 PI3
34 PI4
35 PI5
36 PI6
37 PI7
38 PI8
39 RD1
40 RD2
41 WB1
42 WB2
Because both collars may have recorded the single contact, I need t...
2007 Jun 24
2
matlab/gauss code in R
...ence).
>
> best,
> Simon
>
> >
> >
> > ### example 1 ###
> >
> > Family: gaussian
> > Link function: identity
> >
> > Formula:
> > dep[sel, i] ~ s(date, k = 3) + s(depth, k = kn) + s(temp, k = kn) +
> > s(light, k = kn) + s(PO4, k = kn) + s(DIN, k = kn) + s(prop.agpla,
> > k = kn)
> >
> > Parametric coefficients:
> > Estimate Std. Error t value Pr(>|t|)
> > (Intercept) 5.1048 0.0384 132.9 <2e-16 ***
> > ---
> > Signif. codes: 0 '***' 0.001...
2008 Jun 30
4
Rebuild of kernel 2.6.9-67.0.20.EL failure
Hello list.
I'm trying to rebuild the 2.6.9.67.0.20.EL kernel, but it fails even without
modifications.
How did I try it?
Created a (non-root) build environment (not a mock )
Installed the kernel.scr.rpm and did a
rpmbuild -ba --target=`uname -m` kernel-2.6.spec 2> prep-err.log | tee
prep-out.log
The build failed at the end:
Processing files: kernel-xenU-devel-2.6.9-67.0.20.EL
Checking