Displaying 18 results from an estimated 18 matches for "mu_s".
Did you mean:
mm_s
2012 Apr 16
1
eval a SYMSXP from C
Can someone offer some advice on how to properly evaluate a SYMSXP
from a .Call ?
I have the following in R:
variable xn, with an attribute "mu" which references the variable mu
in the global environment.
I know "references" is a loose term; mu was defined in this fashion as
a way to implement deferred binding:
foo <- function(x,mu) {
attr(x,"mu") <-
2009 Oct 17
2
Recommendation on a probability textbook (conditional probability)
I need to refresh my memory on Probability Theory, especially on
conditional probability. In particular, I want to solve the following
two problems. Can somebody point me some good books on Probability
Theory? Thank you!
1. Z=X+Y, where X and Y are independent random variables and their
distributions are known.
Now, I want to compute E(X | Z = z).
2.Suppose that I have $I \times J$ random number
2010 Jun 08
0
GMM: "The covariance matrix of the coefficients is singular"
...ng this message? Any ideas on how to get
rid of it?
My code:
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
Function g, where E[g]=0:
g <- function(theta,x){
N <- length(x)
#Transform parameters in order to make them x>0, x<0 or 0<x<1
a <- exp(theta[1])
b <- -exp(theta[2])
mu_s <- exp(theta[3])
prob <- exp(theta[4])/(1+exp(theta[4]))
sigma <- exp(theta[5])
sigma_s <- exp(theta[6])
nt1 <- diff(x) - a - b*x[-N]
m1 <- nt1
m2 <- nt1^2 - (sigma^2-2*prob*(mu_s^2+sigma_s^2))/(1-2*prob)
m3 <- nt1^3
m4 <- nt1^4 -
(-2*prob*mu_s^4-12...
2006 Oct 05
2
Writing Text into Plots
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1
Hello,
I have a simple question on the text() method in plots, e.g.:
text(3,4,adj=1,cex=1.0, expression(alpha == beta))
I know there exists a lot more like frac(), etc which could be used for
expression. But a help(frac) doesn't return any results - where do I
have to look for a documentation of possible text commands?
More in detail I'm
2009 Nov 15
2
lme model specification
Dear all
this is a question of model specification in lme which I'd for which I'd greatly appreciate some guidance.
Suppose I have data in long format
gene treatment rep Y
1 1 1 4.32
1 1 2 4.67
1 1 3 5.09
. . . .
. . . .
. . . .
1 4 1 3.67
1 4 2 4.64
1 4 3 4.87
.
2011 Jan 03
1
Greetings. I have a question with mixed beta regression model in nlme.
*Dear R-help:
My name is Rodrigo and I have a question with nlme package
in R to fit a mixed beta regression model. The details of the model are:
Suppose that:*
*j in {1, ..., J}* *(level 1)*
*i in {1, ..., n_j}* *(level 2)*
*y_{ij} ~ Beta(mu_{ij} * phi_{ij}; (1 - mu_{ij}) * phi_{ij})
y_{ij} = mu_{ij} + w_{ij}
*
*with*
*logit(mu_{ij}) = Beta_{0i} + Beta_{1i} * x1_{ij} + b2 * x2_{ij}
2011 Jan 03
0
Greetings. I have a question with mixed beta regression model in nlme (corrected version).
*Dear R-help:
My name is Rodrigo and I have a question with nlme package
in R to fit a mixed beta regression model. I'm so sorry. In the last
email, I forgot to say that W is also a unknown parameter in the mixed
beta regression model. In any case, here I send you the correct formulation.
**
Suppose that:*
*j in {1, ..., J}* *(level 1)*
*i in {1, ..., n_j}* *(level 2)*
*y_{ij} ~
2010 May 18
1
Maximization of quadratic forms
Dear R Help,
I am trying to fit a nonlinear model for a mean function $\mu(Data_i,
\beta)$ for a fixed covariance matrix where $\beta$ and $\mu$ are low-
dimensional. More specifically, for fixed variance-covariance matrices
$\Sigma_{z=0}$ and $\Sigma_{z=1}$ (according to a binary covariate $Z
$), I am trying to minimize:
$\sum_{i=1^n} (Y_i-\mu_(Data_i,\beta))' \Sigma_{z=z_i}^{-1} (Y_i-
2003 Mar 06
6
type III Sum Sq in ANOVA table - Howto?
Hello,
as far as I see, R reports type I sums of squares. I'd like to get R to
print out type III sums of squares.
e.g. I have the following model:
vardep~factor1*factor2
to get the type III sum of squares for factor1 I've tried
anova(lm(vardep~factor2+factor1:factor2),lm(vardep~factor1*factor2))
but that didn't yield the desired result.
Could anyone give me a hint how to proceed?
2002 Dec 10
1
autoregressive poisson process
Dear R users,
I am trying to find a package that can estimate
an autoregressive model for discrete data. I am
imagining a Poisson or Gamma process in which the
mean (say mu) follows a process such as
mu_t = a + b*x + c*mu_{t-1}
Suppose I have data on the time-series Poisson
outcomes and x and would like to obtain ML estimates
for b and c.
Does anyone know of a package that can do this
2002 May 16
1
glm(y ~ -1 + c, "binomial") question
This is a question about removing the intercept in a binomial
glm() model with categorical predictors. V&R (3rd Ed. Ch7) and
Chambers & Hastie (1993) were very helpful but I wasn't sure I
got all the answers.
In a simplistic example suppose I want to explore how disability
(3 levels, profound, severe, and mild) affects the dichotomized
outcome. The glm1 model (see below) is
2004 Oct 08
1
nlme vs gls
Dear List:
My question is more statistical than R oriented (although it originates
from my work with nlme). I know statistical questions are occasionally
posted, so I hope my question is relevant to the list as I cannot turn
up a solution anywhere else. I will frame it in the context of an R
related issue.
To illustrate the problem, consider student achievement test score data
with multiple
2009 Apr 26
1
simulate arima model
I am new in R.
I can simulate Arma, using Arima.sim
However, I want to simulate an Arima Model. Say (1-B)Zt=5+(1-B)at. I do not
know how to deal with 5 in this model.
Can any one could help me?
Thank you very much!
Regards,
--
View this message in context: http://www.nabble.com/simulate-arima-model-tp23239027p23239027.html
Sent from the R help mailing list archive at Nabble.com.
2011 Sep 25
1
Question about syntax in lm function
I encounters some codes in ggplot2 manual and confused with one of its lm syntax. The code is here:
library(ggplot2) d <- subset(diamonds, carat < 2.5 & rbinom(nrow(diamonds), 1, 0.2) == 1) d$lcarat <- log10(d$carat) d$lprice <- log10(d$price) detrend <- lm(lprice ~ lcarat, data = d) d$lprice2 <- resid(detrend) mod <- lm(lprice2 ~ lcarat * color, data = d) # ***
what
2011 Nov 08
2
why NA coefficients
Hi, I am trying to run ANOVA with an interaction term on 2 factors (treat has 7 levels, group has 2 levels). I found the coefficient for the last interaction term is always 0, see attached dataset and the code below:
> test<-read.table("test.txt",sep='\t',header=T,row.names=NULL)
> lm(y~factor(treat)*factor(group),test)
Call:
lm(formula = y ~ factor(treat) *
2001 Nov 14
0
OPEN ssh pkg
Hi,
I have a few corrections for OPENssh-2.9p1 solaris packaging.
I use to configure:
setenv CFLAGS -I/opt/local/include
setenv LDFLAGS "-R/opt/JGMAzlib/lib -L/opt/JGMAzlib/lib"
./configure --prefix=/opt/OPENssh --with-tcp-wrappers --sysconfdir=/etc/opt/OPENssh --localstatedir=/var/opt/OPENssh --with-pid-dir=/var/opt/OPENssh/run --with-ssl-dir=/opt/OPENssl --with-ipv4-default
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
2003 Dec 01
0
No subject
<----------------------------------------------------------------------->
Changes to user passwords are captured by a special DLL, which traps and
then stores the password changes in encrypted form in a private area.
On each synchronization schedule, the synchronization service first examines
the SAM file for changes, and then checks this private area for passwords
to be synchronized. Once