Displaying 18 results from an estimated 18 matches for "fx1".
Did you mean:
0x1
2006 Jul 03
1
gamm
...g.4 %in% g.3 %in% g.2 %in% g.1
(Intercept) Residual
StdDev: 0.6621173 1.007227
Variance function:
Structure: fixed weights
Formula: ~invwt
Fixed effects: y.0 ~ X - 1
Value Std.Error DF t-value p-value
X(Intercept) 2.0870364 0.3337787 392 6.252755 0.0000
Xs(x0)Fx1 -0.0000325 0.1028794 392 -0.000316 0.9997
Xs(x1)Fx1 0.3831694 0.0957323 392 4.002509 0.0001
Xs(x2)Fx1 1.4584330 0.3909237 392 3.730736 0.0002
Xs(x3)Fx1 -0.0123951 0.0143162 392 -0.865809 0.3871
Correlation:
----------------------------...
2012 Jul 27
2
How can I access an element of a string?
....org]
> > On Behalf Of jpm miao
> > Sent: Thursday, July 26, 2012 9:12 PM
> > To: r-help
> > Subject: [R] How can I access the title of a table read via read.csv?
> >
> > Hi,
> >
> > I have a table which I can read via read.csv:
> >
> > fx1<-read.csv(file="A_FX_M.csv", header=TRUE)
> >
> > TIME REER NTD JPY GBP HKD
> > 1 198001 124.26 36.030 237.96 2.263980 4.8366
> > 2 198002 126.59 36.030 244.05 2.290426 4.8765
> > 3 198003 128.33 36.026 248.62 2.206045 4.9960
> >...
2002 Aug 10
0
?subexpressions, D, deriv
...ormula for a higher derivative) possible? and if so, could
someone give me some clues as to how to go about it?
Thank you
----------------------------------------- the original code
> invroot2pi
[1] 0.3989423
> fx
expression((invroot2pi/sigma) * exp(-0.5 * (((x - mu)/sigma)^2)))
> fx1<-D(fx,"x") # first derivative
> fx1
-(invroot2pi/sigma) * (exp(-0.5 * (((x - mu)/sigma)^2)) * (0.5 *
(2 * (1/sigma * ((x - mu)/sigma)))))
> fx2<-D(fx1,"x") # second derivative
> fx2
-(invroot2pi/sigma) * (exp(-0.5 * (((x - mu)/sigma)^2)) * (0.5 *
(2 *...
2012 May 03
1
conducting GAM-GEE within gamm4?
...oups Name Variance Std.Dev.
Xr s(dist_slag) 1611344 1269.39
Xr.0 s(Depth) 98622 314.04
Number of obs: 10792, groups: Xr, 8; Xr.0, 8
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
X(Intercept) -13.968 5.145 -2.715 0.00663 **
Xs(dist_slag)Fx1 -35.871 33.944 -1.057 0.29063
Xs(Depth)Fx1 3.971 3.740 1.062 0.28823
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
Correlation of Fixed Effects:
X(Int) X(_)F1
Xs(dst_s)F1 0.654
Xs(Dpth)Fx1 -0.030 0.000
>
How do I ensure that autocorrel...
2012 Jul 27
2
How can I access the title of a table read via read.csv?
Hi,
I have a table which I can read via read.csv:
fx1<-read.csv(file="A_FX_M.csv", header=TRUE)
TIME REER NTD JPY GBP HKD
1 198001 124.26 36.030 237.96 2.263980 4.8366
2 198002 126.59 36.030 244.05 2.290426 4.8765
3 198003 128.33 36.026 248.62 2.206045 4.9960
4 198004 127.85 36.063 251.67 2.215330 4.9760
5 198005 124....
2010 May 19
1
Displaying smooth bases - mgcv package
...basis functions used by a
thin plate regression spline in a gamm model. I've been searching the help
files, but I can't really figure out how to get the plots of the basis
functions. Anybody an idea?
Some toy code :
require(mgcv)
require(nlme)
x1 <- 1:1000
x2 <- runif(1000,10,500)
fx1 <- -4*sin(x1/50)
fx2 <- -10*(x2)^(1/4)
y <- 60+ fx1 + fx2 + rnorm(1000,0,5)
test <- gamm(y~s(x1)+s(x2))
plot(test$gam)
This gives me the conditional estimates, but not the basis functions for the
different splines. I'd like to have the basis functions for s(x1) and s(x2).
Cheers...
2007 Sep 21
1
problem with 'integrate'
...essage - even with stop.on.error = FALSE. Is there any way I can
work around this? Perhaps there is something similar to
"suppressWarnings()" that I might be able to use?
questionableFunction <- function(B = 10)
{
fx4 <- function(B, y, z){
1 + exp(B*y - B*z);
}
fx1 <- function(x, r0 = 0.12, k=-0.10, a=0.5, tshift=8.36){
(r0*exp(r0*t0 + (k/B)*log(fx4(B, t0, tshift)) - r0*x - (k/B)*log(fx4
(B, x, tshift))))
}
integrate(fx1, 8.36, 10, stop.on.error = FALSE);
}
Thanks in advance for your help...
Dan Rabosky
Department of Ecology and Evolutionary Biol...
2006 Jun 06
1
gamm error message
...function:
Structure: Different standard deviations per stratum
Formula: ~1 | Z
Parameter estimates:
TRUE FALSE
1.0000000 0.9721987
Fixed effects: y ~ X.0 - 1
Value Std.Error DF t-value p-value
X.0(Intercept) 0.04426102 0.1041540 98 0.4249573 0.6718
X.0s(X)Fx1 -0.01053900 0.1039558 98 -0.1013796 0.9195
Correlation:
X.0(I)
X.0s(X)Fx1 0.002
Standardized Within-Group Residuals:
Min Q1 Med Q3 Max
-1.96955692 -0.81101014 0.07004034 0.72528662 2.44244302
Number of Observations: 100
Number of Gr...
2011 Jul 19
2
Incorrect degrees of freedom for splines using GAMM4?
...pwhtc +pwhtsq +edu2+ edu3 +husbgs
+husbcol+ husbmiss
+ +currpmh +pastpmh +neverpmh, random= ~(1|id)
+(1|cogtest), data=global)
Using > summary(global.gamm4$mer), I get the following output for my spline
term, indicating that I use the expected 4 degrees of freedom.
Xs(adjpatx)Fx1 0.1018943 0.1073225 0.949
Xs(adjpatx)Fx2 -0.0708114 0.1123845 -0.630
Xs(adjpatx)Fx3 0.7459511 0.6836413 1.091
Xs(adjpatx)Fx4 -0.2062321 0.0923569 -2.233
However, when I use > summary(global.gamm4$gam). I get an estimate of
degrees of freedom that is not 4:
Approximate significance...
2013 Jun 07
1
gamm in mgcv random effect significance
...1, int=~1),correlation=corAR1())
Any help would be appreciated.
Thanks.
Will Shadish
********************************************
g1
$lme
Linear mixed-effects model fit by maximum likelihood
Data: data
Log-likelihood: -437.696
Fixed: fixed
X(Intercept) Xz Xint Xs(xc)Fx1
0.6738466 -2.5688317 0.0137415 -0.1801294
Random effects:
Formula: ~Xr - 1 | g
Structure: pdIdnot
Xr1 Xr2 Xr3 Xr4
Xr5 Xr6 Xr7 Xr8 Residual
StdDev: 0.0004377781 0.0004377781 0.0004377781 0.0004377781 0.000437778...
2002 Sep 23
0
arima() in package ts.
...e of the innovations).
Can it be computed directly from the residuals? How?
(Since the whole MLE structure is very complicated, I'm
pretty sure the answer is ``No it can't'', but I'd like
a little more insight into what's going on here.)
(3) I tried
x <- rnorm(100)
fx1 <- arima(x,c(1,0,0),c(1,0,0))
fx2 <- arima(x,c(2,0,0))
Note that seasonality is not specified and so defaults to 1.
I would have thought that the coeficients of the second
model, say phi.21 and phi.22, would be related to the coefficients
of the first (``multiplicative'') model, sa...
2012 Jul 19
3
write list to ascii
Dear all,
apologies for this (perhaps recurrent) question but I did not found a question when searching mailing lists.
How to write a list of a simple kind, e.g.:
abc <- list(one=(1:2), two=(1:5))
# to a file? I understand that write() & co. cannot work but when I try
sink("aa.txt", append=T, split=T)
abc
sink()
# the output is indeed "split by rows" in the
2016 Apr 04
1
PA-RISC (hppa) video cards init failure loading the device driver kernel module
Dear Ilia,
this page summarize all the possible options of the "nouveau" kernel module:
https://nouveau.freedesktop.org/wiki/KernelModuleParameters/
I will test the options that you have suggested as soon as possible.
Try a PCI video card is a very good advice because the
AGP/FX1(chipset) is on the top of the "usual suspects" list;
unfortunately the c8000 has 3.3 PCI slots and all the (very old) PCI
video cards in my hands have 5V PCI connector
Many thanks for your suggestions and advice.
Simone Mannori - Italy
On 4 April 2016 at 04:18, Ilia Mirkin <imir...
2016 Apr 04
0
[Radeon] PA-RISC (hppa) video cards init failure loading the device driver kernel module
...ummarize all the possible options of the "nouveau" kernel module:
>
> https://nouveau.freedesktop.org/wiki/KernelModuleParameters/
>
> I will test the options that you have suggested as soon as possible.
>
> Try a PCI video card is a very good advice because the
> AGP/FX1(chipset) is on the top of the "usual suspects" list;
> unfortunately the c8000 has 3.3 PCI slots and all the (very old) PCI
> video cards in my hands have 5V PCI connector
>
> Many thanks for your suggestions and advice.
>
> Simone Mannori - Italy
>
>
> On 4...
2004 Aug 06
0
toner cartridges
...em # 14</font></font></center>
</td>
<td>
<center><font color="#000099"><font size=+1>Leserfax 500, 700</font></font></center>
</td>
<td>
<center><font color="#000099"><font size=+1>FX1</font></font></center>
</td>
<td>
<center><font color="#000099"><font size=+1>$59</font></font></center>
</td>
</tr>
<tr>
<td>
<center><font color="#000099"><font size=+1&g...
2002 Jul 12
1
toner cartridges
...em # 14</font></font></center>
</td>
<td>
<center><font color="#000099"><font size=+1>Leserfax 500, 700</font></font></center>
</td>
<td>
<center><font color="#000099"><font size=+1>FX1</font></font></center>
</td>
<td>
<center><font color="#000099"><font size=+1>$59</font></font></center>
</td>
</tr>
<tr>
<td>
<center><font color="#000099"><font size=+1&g...
2016 Apr 03
3
PA-RISC (hppa) video cards init failure loading the device driver kernel module
Dear "nouveau" developers,
I know that many very competent guys have already spent a lot of time
and efforts on this issue without success. I have started to play with
"hppa" two weeks ago and, with the support of the linux-parisc mailing
list people, now I have a - almost fully - working workstation (hp
c8000).
Everythings work perfectly BUT the video card. No matter the
2009 Jan 28
2
t.test in a loop
Hi All,
I've been having a little trouble with creating a loop that will run a a
series of t.tests for inspection,
Below is the code i've tried, and some checks i've looked at.
I've used the get(paste()) idea as i was told previously that the use of the
eval should try and be avoided.
I've run a single syntax to check that my systax is correct and works
without any problems