Displaying 20 results from an estimated 2000 matches similar to: "lmer function and Inverse mills ratio"
2009 Jul 11
2
Heckman Selection Model/Inverse Mills Ratio
I have so far used the following command
glm(formula = s ~ age + gender + gemedu + gemhinc + es_gdppc +
imf_pop + estbbo_m, family = binomial(link = "probit"))
My question is
1. How do i discard the non significant selection variables (one out of the
seven variables above is non-significant) and calculate the Inverse Mills
Ratio of the significant variables
2. I need the inverse
2009 Sep 08
0
Inverse Mills in clustered (multilevel) cross-sectional panel data
Dear R saviors,
kindly address to this problem, I would really appreciate any takers. I am
trying to resolve this issue of IMR in clustered (multilevel)
cross-sectional panel data for more than two months now,.
The characteristics of my dataset are as follows:
- some 900 000 individuals
- total of 60 countries
- cross-sectional time series at the country level max 10 years, not all
2009 Jul 12
0
ERROR message while using <-invMillsRatio()
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
2009 Jul 10
0
GLM for Probit for Panel Data
Hello
I am working on a panel data, my panel variable is the variable "yearctry",
let me explain what I mean, yearctry is calculated based on the year and the
ISD phone code of a country, eg, for the year 2000 say and for country USA
say (code = 001), my yearctry variable will then be 2000001, there are 2000
observations (ie 2000 individual responses with yearctry = 2000001), I have
65
2009 Jul 12
2
Heckman Selection MOdel Help in R
Hi Saurav!
On Sun, Jul 12, 2009 at 6:06 PM, Pathak,
Saurav<s.pathak08 at imperial.ac.uk> wrote:
> I am new to R, I have to do a 2 step Heckman model, my selection equation is
> below which I was successful in running but I am unable to proceed further,
>
>
>
> I have so far used the following command
>
> glm(formula = s ~ age + gender + gemedu + gemhinc + es_gdppc +
2009 Oct 18
2
How to create MULTILEVELS in a dataset??
Dear R users
I have a data set which has five variables. One depenedent variable y, and 4
Independent variables (education-level, householdincome, countrygdp and
countrygdpsquare). The first two are data corresponding to the individual
and the next two coorespond to the country to which the individual belongs
to. My data set does not make this distinction between individual level and
country
2009 Oct 17
0
how to cluster data for use with lmer
Dear R users
My data set is e
> names(e)
[1] "yearctry" "discent" "age" "gender"
"gemeduc" "gemhhinc" "ref_group" "fearfail_ref" "knowent_ref"
"nbgoodc_ref"
[11] "nbstatus_ref" "estbbuso_ref" "lngdp" "lngdpsq"
2009 Jul 15
0
DECLARING A PANEL VARIABLE???
Hi
I am working on a panel data, my data are clustered/grouped by the variable
"yearctry", I am running the regression below, but I cant make the
regression recognise "yearctry" as the panel variable in the regression
myProbit<- glm(s ~ age + gender + gemedu + gemhinc + es_gdppc +
imf_pop + estbbo_m, family = binomial(link = "probit"), data =
adpopdata)
Can
2009 Jul 16
1
PROBIT REGRESSION FOR GROUPED/CLUSTERED DATA
Hello all
I have been working to fix this for weeks now, It should be simple to fix.
Please help
Let me explain what I am doing, I have a data set for 65 countries over a
period of 9 years (2000-2008). Each country has on an average say 2000
interviews, so that the total set has roughly 65*9*2000 data
points/observations (of course there are missing vales as well). Now let me
explain how are the
2009 Jul 01
2
getOptions("max.print") in R
I am typing the following on the command prompt:
>variab = read.csv(file.choose(), header=T)
>variab
It lists 900,000 ( this is the total number of observations in "variab" )
minus 797124 observations and prompts the following message
[ reached getOption("max.print") -- omitted 797124 entries ]]
Is there a way to see the entire set of data, ie all of 900,000 obs, and
2009 Jun 29
1
Stata file Import and Analysis in R
Hi
I have a stata data set (.dta file) of size 600 MB, I need to import it in
R and do a 2SLS multilevel analysis on the data set. I would eb grateful if
help is provided for the first part of how to import this big file from
Stata to R and then how to open the imported Stata file in R?
Kindly help
Thanks in advance
--
Dr.Saurav Pathak
PhD, Univ.of.Florida
Mechanical Engineering
Doctoral
2009 Jun 29
2
Large Stata file Import in R
Hi
I am using Stata 10 and I need to import a data set in stata 10 to R, I have
saved the dataset in lower versions of Stata as well by using saveold
command in Stata.
My RAM is 4gb and the stata file is 600MB, I am getting an error message
which says :
"Error: cannot allocate vector of size 3.4 Mb
In addition: There were 50 or more warnings (use warnings() to see the first
50)"
Thus
2009 Jun 30
1
Stata file and R Interaction :File Size Problem in Import
Hi
I am using Stata 10 and I need to import a data set in stata 10 to R, I
have saved the dataset in lower versions of Stata as well by using saveold
command in Stata.
My RAM is 4gb and the stata file is 600MB, I am getting an error message
which says :
"Error: cannot allocate vector of size 3.4 Mb
In addition: There were 50 or more warnings (use warnings() to see the
first
2004 Aug 19
1
sample selection problem, inverse mills ratio (Heckman, Lewbel, ...)
-----Ursprüngliche Nachricht-----
Von: Wildi Marc, wia
Gesendet: Mittwoch, 18. August 2004 10:11
An: r-help@lists.R-project.org
Betreff:
Hi
Does anybody know from an R-package devoted to sample selection problems (Heckman's lambda, Lewbel, ...)?
Thanks and best regards
Marc Wildi
[[alternative HTML version deleted]]
2011 Nov 30
1
How can I pick a matrix from a function? (Out Product of Gradient)
Hi all,
I would like to use optim() to estimate the equation by the log-likelihood
function and gradient function which I had written. I try to use OPG(Out
Product of Gradient) to calculate the Hessian matrix since sometime Hessian
matrix is difficult to calculate. Thus I want to pick the Gradient matrix
from the gradient function.
Moreover, could R show the process of calculation on gradient
2003 Jan 14
4
density plot - beginner's question
Hi,
I am trying to plot densities given on a two dimensional grid. My
data is in the an external file, and is arranged in three columns:
x, y, density
how may i get a plot of this? i would like to get (1) a three
dimensional plot and (2) a color coded two dimensional plot.
I have tried using
image(x, y, density)
but i am asked to put the data in ascending order. i am not sure
how i may
2007 Jan 24
2
modify rectangle color from image
Hi,
I need some suggestion on how I could modify the color on some
rectangle that I have created using "image".
In other words, I have a 5x5 matrix, say, m.
m <- matrix(rnorm(25), nrow=5)
I create a grid of rectangles by:
image(m)
Now I want to change the color of rectangle (3,3) to blue.
I don't know how this could be done, and searching the web has
2005 Feb 28
2
3d scatterplots of more than 1 data set
hi,
i am need to plot two or more sets of data in a 3d scatterplot,
each set with different color.
i tried Rcmdr, and the 3d scatterplot facility, based on rgl. that
is what i need. but i cannot seem to code different sets of data
differently. any help will be very helpful.
i tried scatterplot3d, but it is difficult to get the right angle in
it. i need to be able to rotate the axes, and
2016 Apr 15
0
Odds Ratio and OR CI
Please don't post in HTML. Your post is almost unreadable
Also,lease have a look at
http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example and/or http://adv-r.had.co.nz/Reproducibility.html for some suggestions about what you should include in your question.
Welcome to R.
John Kane
Kingston ON Canada
> -----Original Message-----
> From: rosita21 at
2011 Oct 07
1
loess question
Hi All,
I am trying to use loess to smooth a 2D image, and also obtain the
standard error for every pixel. I see that the standard error does not
make sense. For example, running the following:
library(stats)
x <- array(c(1:100), dim=c(100,100))
y <- t(x)
v <- exp(-((x-50)^2+(y-50)^2)/30^2)
s <- v*0.02
g_noise <- rnorm(10000, mean = 0, sd = s)
f <- v + g_noise
f.loess <-