Displaying 20 results from an estimated 500 matches similar to: "2SLS / TSLS / SEM non-linear"
2009 Aug 19
2
lmer with random slopes for 2 or more first-level factors?
I have data from a design in which items are completely nested within
subjects. Subject is the only second-level factor, but I have
multiple first-level factors (IVs). Say there are 2 such independent
variables that I am interested in. What is the proper syntax to fit a
mixed-effects model with a by-subject random intercept, and by-subject
random slopes for both the 2 IVs?
I can
2009 Feb 10
2
Mixed ANCOVA with between-Ss covariate?
Hi all,
I have data from an experiment with 3 independent variables, 2 are
within and 1 is between. In addition to the dependent variable, I have
a covariate that is a single measure per subject. Below I provide an
example generated data set and my approach to implementing the ANCOVA.
However the output confuses me; why does the covariate only appear in
the first strata? Presumably it should
2011 Feb 03
2
how to read the "Sum Sq" - column from summary.aov()
Dear R-Users,
I have a trivial problem, but extensive googling and ??'ing did not solve it: I want to obtain the sums of squares from a summary.aov() object for later use. Example:
> DV <- rnorm(100)
> IV1 <- as.factor(rep(c("male", "female"), times = 50))
> IV2 <- as.factor(rep(c("young", "old"), times = 50))
>
>
2009 Nov 06
2
Adjusting Yaxis (ylim) limits on a plotMeans(DV, IV1, IV2, error.bars="se")
Hello everyone,
I have tried to look for this everywhere and so far have no luck. I have a
plotMeans(DV, IV1, IV2, error.bars="se") graph that plots my data
(DV-continuous, IVs are factors, IV1 - two levels, IV2-four levels). I am
trying to increase a scale of my y-axis (to be consistent with my other
graphs), but unfortunately nothing works with "plotMeans" function, which
2007 Feb 19
1
Urgent: How to obtain the Consistent Standard Errors after apply 2SLS through tsls() from sem or systemfit("2SLS") without this error message !!!!!!!!!!!!!
Hi,
I am trying to obtain the heteroskedasticity consitent standard errors
(HCSE) after apply 2SLS. I obtain 2SLS through tsls from package sem or
systemfit:
#### tsls ####
library (sem)
Reg2SLS <-tsls(LnP~Sc+Ag+Ag2+Var+R+D,~I2+Ag+Ag2+Var+R+D)
summary (Reg2SLS)
#### systemfit ####
library (systemfit)
RS <- LnP~Sc+Ag+Ag2+Var+R+D
Inst <- ~I2+Ag+Ag2+Var+R+D
labels
2012 May 19
3
anovas ss typeI vs typeIII
Hi all,
I have been struggling with ANOVAs on R. I am new to R, so I created a simple data frame, and I do some analyses on R just to learn R and then check them on SPSS to make sure that I am doing fine. Here is the problem that I've run into:
when we use the aov function, it uses SS Type I as default (on SPSS it is Type III). Then I used the Anova function under cars package using the
2005 Jul 07
1
multivariate regression using R
Does anyone know if there is a way to run multivariate linear regression
in R? I tried using the lm function (e.g., lm(dv1, dv2~iv1+iv2+iv3),
but got error messages. Is my syntax wrong, or do I need a particular
package?
Thanks,
Jeff--
________________________________________________________
Jeffrey J. Lusk, Ph.D.
Postdoctoral Research Associate
Department of Forestry &
2008 Jun 22
1
two newbie questions
# I've tried to make this easy to paste into R, though it's probably
so simple you won't need to.
# I have some data (there are many more variables, but this is a
reasonable approximation of it)
# here's a fabricated data frame that is similar in form to mine:
my.df <- data.frame(replicate(10, round(rnorm(100, mean=3.5, sd=1))))
var.list <- c("dv1",
2008 Sep 25
2
levelplot/heatmap question
Hello!
I have data containing a large number of probabilities (about 60) of nonzero
coefficients to predict 10 different independent variables (in 10 different
BMA models). i've arranged these probabilities in a matrix like so:
(IV1) (IV2) (IV3) ...
p(b0) p(b0) p(b0)
p(b1) p(b1) p(b1)
p(b2) p(b2) p(b2)
...
where p(b1) for independent variable 1 is p(b1 !=
2001 Aug 22
1
limited formula length in tsls
Dear all,
Using the tsls package, I noticed that regression lists longer
than 64 character
are getting truncated. Looking at the original source,
tsls.formula <- function(model, instruments, data, subset, weights,
na.action, contrasts=NULL){
if (missing(na.action))
na.action <- options()$na.action
m <- match.call(expand.dots = FALSE)
if (is.matrix(eval(m$data,
2012 Nov 29
1
instrumental variables regression using ivreg (AER) or tsls (sem)
Dear friends,
I am trying to understand and implement instrumental variables
regression using R.
I found a small (simple) example here which purportedly illustrates the
mechanics (using 2-stage least-squares):
http://www.r-bloggers.com/a-simple-instrumental-variables-problem/
Basically, here are the R commands (reproducible example) from that
site:
# ------ begin R
library(AER)
2008 Nov 11
1
simulate data with binary outcome and correlated predictors
Hi,
I would like to simulate data with a binary outcome and a set of predictors that are correlated. I want to be able to fix the number of event (Y=1) vs. non-event (Y=0). Thus, I fix this and then simulate the predictors. I have 2 questions:
1. When the predictors are continuous, I can use mvrnorm(). However, if I have continuous, ordinal and binary predictors, I'm not sure how to simulate
2002 Jan 21
6
OpenSSH and OpenSSL snapshots
Hello,
In order to experiment with OpenCA, I have built an RPM for redhat 7.2
of a recent OpenSSL snapshot (the binary rpm on the OpenCA was built
with the different target directories and libraries. Unfortunately these
recent OpenSSL snapshots seems to break all OpenSSH tarballs and RPMs
that I have been able to find. None of them seem to compile
successfully, even the snapshots at
2008 Aug 22
1
filtering out data
Greetings,
Apologies for such a simple question, but I have been trying to figure this out for a few days now (I'm quite new to R), have read manuals, list posts, etc., and not finding precisely what I need. I am accustomed to working in Stata, but I am currently working through the multilevel package right now. In Stata, I would include the statement "if model1 == 1" at the end
2002 Jan 04
1
plotting missing data patterns
I have one other problem and then I'll stop and get back to writing. I
want to plot a missing data matrix (the R matrix a la Little and Rubin)
to graphically depict where the missing data lie in the matrix. Some
statistical packages produce a graphical depiction of the missing data
patterns by plotting the matrix, color-coding the plot with contrasting
colors for either missing or
2011 Apr 18
1
regression and lmer
Dear all,
I hope this is the right place to ask this question.
I am reviewing a research where the analyst(s) are using a linear
regression model. The dependent variable (DV) is a continuous measure.
The independent variables (IVs) are a mixture of linear and categorical
variables.
The
author investigates whether performance (DV - continuous linear) is a
function of age (continuous IV1 -
2012 Jun 14
3
p-values from lm()
Ein eingebundener Text mit undefiniertem Zeichensatz wurde abgetrennt.
Name: nicht verf?gbar
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20120614/e05ae92c/attachment.pl>
2012 Mar 21
1
How to do 2SLS in R
Hi List
I want to carry out structural mode. Following Example l have taken from
Basic Econometrics- Damodar Gujarati :
Advertising intensity function:
Ad/S = a0 + a1M + a2 (CD/S) + a3C + a4C2 + a5Gr + a6Dur – (1)
Concentration function:
C = b0 + b1 (Ad/S) + b2 (MES/S) -(2)
Price-cost margin function:
M = c0 + c1(K/S) + c2Gr + c3C + c4GD + c5(Ad/S) + c6 (MES/S)
2007 Apr 09
1
Dealing with large nominal predictor in sem package
Hi,
I am using tsls function from sem package to estimate a model which includes large number of data. Among its predictors, it has a nominal data which has about 10 possible values. So I expand this parameter into 9-binary-value predictors with the coefficient of base value equals 0. I also have another continuous predictor.
The problem is that, whenever I run the tsls, I will get 'System
2010 May 02
1
question about 2SLS
Hi All,
I am using R 2.11.0 on a Ubuntu machine. I estimated a model using "tsls"
from the package "sem". Is there a way to get Newey West standard errors for
the parameter estimates?
When estimating the model by OLS, I used "NeweyWest" from the package
"sandwich" to get HAC standard errors. But, I am not able to use the same
method with the results of the