similar to: Adding .PDF files to a package

Dear helpers, I've come up with what is probably a simple problem, but I cannot find the solution. I have a data-set containing survey-data from several countries. What I want to do is to perform some regression analyses, for each country separately. The question is, how to do this nicely (thus without repeating the same syntax with another `subset' argument). I thought of the

Dear expeRts, I am trying to plot a lme-object {package nlme) in such a way, that on a selected level the x-axis represents the value on a selected predictor and the y-axis represents the predicted-outcome variable. The graphs would than consist of several lines that each represent one group. I can't find such a plotting function. I could write such a function myself, based on

Dear R users, I am trying to get data from the clipboard into R on MacOSX. I tried the following, but got an error message: read.delim("clipboard") Error in file(file, "r") : unable to open connection In addition: Warning message: unable to contact X11 display Obviously, I'm not running R using X11. I'm wondering, can I import data from the clipboard on MacosX?

Dear R users, experimenting with the lm function in R, I've encountered some behaviour I don't understand with my limited knowledge of regression. I made a data-'set' of three measurements (see syntax below). Using lm (linear model) to fit the regression-line, I expected to find an intercept of 2.0 and a slope of 0, but in fact the slope is slightly below zero. Amazed by

Dear ListeRs, I am trying to fit a heteroscedastic multilevel model using lmer{lme4- package). Take, for instance, the (fictive) model below. lmer(test.result ~ homework + Sex -1 + (1 | School)) Suppose that I suspect the error terms in the predicted values to differ between men and women (so, on the first level). In order to model this, I want the 'Sex'-variable to be random on

Dear helpers, I'm estimating multilevel regression models, using the lme-function from the nlme-package. Let's say that I estimated a model and stored it inside the object named 'model'. The summary of that model is shown below: Using summary(model)$tTable , I receive the following output: > summary(model)$tTable Value Std.Error DF t-value

Dear helpeRs, I'm working with the map-package and came upon a problem which I couldn't solve. I hope onee of you can. If not, this can be seen as a suggestion for new versions of the package. I'm trying to create a map of some European countries, filled with colors corresponding to some values. Let's say I have the following countries and I assign the following colors

Dear All, I'm trying to model heteroscedasticity using a multilevel model. To do so, I make use of the nlme package and the weigths-parameter. Let's say that I hypothesize that the exam score of students (normexam) is influenced by their score on a standardized LR test (standLRT). Students are of course nested in "schools". These variables are contained in the

LS, I'm estimating multilevel regression models, using the lme-function from the nlme-package. Let's say that I estimated a model and stored it inside the object named 'model'. The summary of that model is shown below: Using summary(model)$tTable , I receive the following output: > summary(model)$tTable Value Std.Error DF t-value

Hi, I'm running some data simulations using (mixed effects)* regression models that show difficulty to converge. Therefore, I seek a way of capturing warnings (of false convergence) inside a loop. Inside that loop, I modify data and estimate a model. I do so many times with slightly different modifications of the data. Next, I extract some of the model parameters and store these in a matrix.

Hi, I'm trying to import categorical data from SPSS to R using the script: xxx <-spss.get("xxx.por", use.value.labels=TRUE) but unfortunately am getting an error message 'error reading portable-file dictionary'. I have successfully imported data in the past. What could be the problem with this data? Thanks Simo

Hi there, as a student sociology, I'm starting to learn about SEM. The course I follow is based on LISREL, but I want to use the SEM-package on R parallel to it. Using LISREL, I found it to be very usable to be able to see the total direct and total indirect effects (standardized and unstandardized) in the output. Can I create these effects using R? I know how to calculate them

Hi, I am trying to learn some basic statistics stuff but I cannot find any elementary statistics exercises using R language. Using RKward would be even better... I need that in analysing sociological data, obtained through questionnairres - findind corelations between variables, relations between different types of data, etc. Could anyone recommend simple tutorials/exercises, available on www

Dear all, I'm investigating on how to estimate a specific kind of cross- classified multilevel model. I think it is often referred to as a multiple-membership model. To problem is this: I want to study changing attitudes of people, for which I use (amongst other things) contextual data of the neighborhood they live in. I have two surveys (panel-design) and I know the neighborhood

LS, I'm estimating multilevel regression models, using the lme-function from the nlme-package. Let's say that I estimated a model and stored it inside the object named 'model'. The summary of that model is shown below: Using summary(model)$tTable , I receive the following output: > summary(model)$tTable Value Std.Error DF t-value

Hi R users, I need to draw a map of select European countries with country names shown on the map. Does anyone know how to do this in R? Also, is it possible to draw a historical map of European countries using R? Thanks a lot for your help. Maomao [[alternative HTML version deleted]]

Hi, I cannot figure out how to draw a certain plot: could someone help me out? I have this data.frame from a survey my.data that looks like something like this: col1 col2 col3 col4 1 5 5 4 5 2 3 5 3 1 3 2 3 4 5 4 3 1 1 2 5 5 5 4 5 6 4 2 5 5 .... Each row represents a single questionnaire

I have a small network with mixed OS's -- some Linux (mix of distros, Ubuntu, PCLinuxOS, CentOS) and some Windows (XP and Win7). I'm replacing my server and have installed CentOS 6.2 on it. There are no Windows servers on the network. I'd like to have the server manage user accounts and passwords. Ideally, I'd be able to log in to any computer on my network with the same

Hi, Fitting all possible models (GLM) with 10 predictors will result in loads of (2^10 - 1) models. I want to do that in order to get the importance of variables (having an unbalanced variable design) by summing the up the AIC-weights of models including the same variable, for every variable separately. It's time consuming and annoying to define all possible models by hand. Is there a

Hi Is it possible to subset an n-dimensional array by a vector of n dimensions? E.g. assume I have > x <- array(1:24, dim=2:4) > x[1,1,2] [1] 7 > dims <- c(1,1,2) I would like a function that I can supply x and dims as parameters to, and have it return 7. Also, I would like to do something like: > x[1,1,] [1] 1 7 13 19 > dims2<- c(1,1,NA) And have a function of x and