similar to: How to create MULTILEVELS in a dataset??

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

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"

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

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

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

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 +

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

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

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

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 <-

Hi, For some reason I have been unable to use the predict function when I desire the standard error to be calculated too. For example, when I try the following: l<- loess(d~x+y, span=span, se=TRUE) p<- predict(l, se=TRUE) I get the following error message: Error in vector("double", length) : vector size cannot be NA In addition: Warning message: In N * M1 : NAs produced by

hi all, i am trying to make multi-color plots. that is, i generally use, plot(x, y, type="n") text(x, y, labels=class) here, the vector class denotes the class of each point. there are usually 3-4 classes of points. how may i display the different classes in different colors? thanks for any help. -- saurav

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

Hi, I am building an application where i need to have an admin interface from where the admin can assign permission for different methods (create, edit etc.) in controller based on roles . It should be similar to django administration, where we can assign permission on methods based on the role. How can i do the implementaton ? Can anybody give me some pointers on how to start ? Thanks and

Hi all, I need a help. I have got a user registration form, where the user signs up and then an activation link is sent to the user email id for activating his/her account. Now I need to send an another email to the user, only after he logs in for the first time in my site, then the second email should be sent. Can anybody give me some suggestion on how to do it ?? NB: When the user clicks on

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

hi, i have just begun using R, so please bear with me. i am trying to use cmdscale and display the result. i read the data using read.table(), calculate the proximity matrix using dist() and the display the result using the cmdscale(). this is very fine. in addition, i want the display to distinguish between two classes of records in my data. i have my data records marked as "1" or

Hello, I am trying to create a lucid guest VM with a vmchannel between the guest and the hypervisor. Details below: KERNEL: # uname -a Linux saurav-desktop 2.6.38-2-generic #29~lucid1-Ubuntu SMP Mon Feb 7 13:35:14 UTC 2011 i686 GNU/Linux QEMU: # /usr/bin/qemu --version QEMU emulator version 0.13.0, Copyright (c) 2003-2008 Fabrice Bellard LIBVIRT: # /usr/sbin/libvirtd --version

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

hi all, I have an application that runs on a diskless card. The application stores its data on the flash and also as a backup on the filesystem of a workstation which is NFS mounted on the card. Can I achieve data integrity/relaiability if I use ext3 on NFS filesystem in case of power failure/card jackout? What about the data in the flash? What should be the ext3 filesystem mode? Any help is