similar to: Problem with dummy variables controlling for years

Hello everyone, It seems that the answer to my question should be obvious, but I''ve not yet found a way to do what I want. Any advice on the best method would be greatly appreciated. Have the following database table: CREATE TABLE datasets ( id int not null auto increment, name varchar not null, value varchar not null, year varchar not null, state

Hallo, I have this dataset with repeated measures. There are two within-subject factors, "formant" (2 levels: 1 and 2) and "f2 Ref" (25 levels: 670, 729, 788, 846, 905, 1080, 1100, 1120, 1140, 1170, 1480, 1470, 1450, 1440, 1430, 1890, 1840, 1790, 1740, 1690, 2290, 2210, 2120, 2040, 1950), and one between-subject factor, lang (2 levels:1 and 2). The response variable

I'm attempting to do model selection with AIC, using a glm and a lognormal distribution, but: fit1<-glm(BA~Year,data=pdat.sp1.65.04, family=gaussian(link="log")) ## gives the same result as either of the following: fit1<-glm(BA~Year,data=pdat.sp1.65.04, family=gaussian) fit1<-lm(BA~Year,data=pdat.sp1.65.04) fit1 #Coefficients: #(Intercept) Year2004 # -1.6341

I'm trying to plot a time series (quarterly observations, seasonal effect, T=56),regression line and predicted point on the same plot. I'm using following commands: > plot(qdts,xlim=c(1982,1997)) > lines(fitted(m2)~time,data=qd,col="red") > points(predict(m2,newdata=nqd),col="blue") where: qdts<-ts(qd$durables,start=c(1982,1),freq=4)

Rodolfo, I haven't looked up how to do this with sip phones, but the zap channels can be configured in groups that will hunt through the group until a non-busy line is found. http://www.voip-info.org/tiki-index.php?page=Asterisk%20ZAP%20channels#comments Here is a link to PBX hunting with the dial plan. http://www.voip-info.org/tiki-index.php?page=PBX+Hunt+Groups I haven't tried

Hi R Gurus! I conducted regression and anova followings : Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 6.07e-01 5.95e-02 10.19 < 2e-16 *** nemp 2.87e-06 1.04e-07 27.63 < 2e-16 *** as.factor(corridor1)A -8.81e-02 2.13e-02 -4.14 3.6e-05 *** as.factor(corridor1)B

Hi, I am having a problem figuring out how to extract a subset of rows. I have a list with 68 similar data.frames. Each data.frame is 500 rows by 5 columns. I want to take one row from each data.frame based upon the data in a particular column (i.e. it matches a symbol). For example: > str(database) List of 68 $ X2011.01.11:'data.frame': 500 obs. of 5 variables: ..$ Symbol

Hi! I have some questions about MARS model's coefficient of determination. I use the MARS method in my master's thesis and I have noticed some problems with the MARS model's R^2. You can see the following example that the MARS model's R^2 is too big when i have used mars() -function for MARS model building, and when I have made MARS-model using a linear regression, it

here ,I want to plot two lines in one figure.But I have two problems 1) how to move one of the y-axis to be the right ? I tried to the command"axis(2)",But I failed. 2) how to add the axis information correctly.Since I have use the cmommand "axis(1,at=1:6,labels=gradeinfo$gradenam)" but it seems that the correct information that I want is superposition with the old axis

Howdy I apologize for duplicated posting. But I decided to correct my previous posting. I had the regression results using r <- lm(Y ~ nemp + as.factor(devt), data=d). First, there is the result of anova(r). Here I could not find regression coefficients. Response: Y Df Sum Sq Mean Sq F value Pr(>F) nemp 1 58.2 58.2 1233.23 < 2e-16 ***

I would be very grateful to anyone who could point to the error of my ways in the following. I have a dataframe called net1, as such: > str(net1) `data.frame': 114192 obs. of 9 variables: $ server : Factor w/ 122 levels "AB93-99","AMP93-1",..: 1 1 1 1 1 1 1 1 1 1 ... $ ts :'POSIXct', format: chr "2006-06-30 12:31:44"

Dear list members, I am trying to estimate parameters of the AR(1)-GARCH(1,1) model. I have one additional dummy variable for the AR(1) part. First I wanted to do it using garchFit function (everything would be then estimated in one step) however in the fGarch library I didn't find a way to include an additional variable. That would be the formula but, as said, I think it is impossible to add

Am I trying to perform multiple linear regressions on each 'VARIABLE2'. I figured out that there are different ways, using the following code: (data is given at the end of this message) reg <- lapply(split(TRY, VARIABLE2), function(X){lm(X2 ~ X3, data=X)}) lapply(reg, summary) Which produces the following: $`1` Call: lm(formula = X2 ~ X3, data = X) Residuals: Min

Dear list, I'm trying to mimic the analysis of Wedderburn (1974) as cited by McCullagh and Nelder (1989) on p.328-332. This is the leaf-blotch on barley example, and the data is available in the `faraway' package. Wedderburn suggested using the variance function mu^2(1-mu)^2. This variance function isn't readily available in R's `quasi' family object, but it seems to me

>>>>> Constantin Ahlmann-Eltze via R-devel >>>>> on Mon, 10 Aug 2020 10:05:36 +0200 writes: > Thanks Ben for verifying the issue. It is always reassuring to hear > when others can reproduce the problem. > I wrote a small patch that fixes the issue > (https://github.com/r-devel/r-svn/pull/11): > diff --git

I am a novice R user, and I am having difficulty understanding R's influence plots. I am trying to remove outliers from a particular variable, "sib." I am able to generate influence plots and further outlier information such as below (which is a shortened example). For my analyses, I end up excluding the points R refers to, 7, 18, 26, and 105. However, my question is, how can I

Thanks Ben for verifying the issue. It is always reassuring to hear when others can reproduce the problem. I wrote a small patch that fixes the issue (https://github.com/r-devel/r-svn/pull/11): diff --git a/src/nmath/qnbinom.c b/src/nmath/qnbinom.c index b313ce56b2..d2e8d98759 100644 --- a/src/nmath/qnbinom.c +++ b/src/nmath/qnbinom.c @@ -104,6 +104,7 @@ double qnbinom(double p, double size,

Hi I have been trying so many different things to get my Inverse Mills Ratio going for a Two stage Heckman Model, I have tried the following so far (the commands are listed below till teh point where I get an error), I get an error in the last sentence (marked in bold below), if this were successful then I could have used the IMR as a control in my OLS (which would be the OLS for the outcome

A couple of suggestions: 1. This screams out for do.call. Try jj <- do.call("rbind", t1). 2. Use rowSums() instead of apply(..., 1, sum). Andy _____ From: r-help-bounces at stat.math.ethz.ch on behalf of Mike Nielsen Sent: Sat 7/8/2006 7:20 PM To: r-help at stat.math.ethz.ch Subject: Re: [R] Combining a list of similar dataframes into a single dataframe [Broadcast] Well,

Hi all, as I told earlier, I am installing an email server with postfix/dovecot/dspam (for now that's all, though some programs may come later). For now, what works good enough is that I receive mails and spams are tagged. Dovecot version 1.1.7-r1 compiled from the gentoo depository with debug. The path of a mail is 1/ to be passed to dspam by postfix as a virtual transport with this in