search for: nieuwenhuis

...ustered according to `Country', but not so for the predictors. The error message I received: Error in tapply(dat, dat$COUNTRY, function(x) glm(participate ~ FEMALE + : arguments must have same length > Could you please help me solving this little problem? Thanks in advance: Rense Nieuwenhuis [[alternative HTML version deleted]]

...;, no links seems to be present. All checks are OK, but no vignette or whatever appears in the help-files of the package. Clearly, I'm missing something. I would really appreciate some pointers on how to integrate an 'additional' pdf file to my package. Thanks in advance, Rense Nieuwenhuis

...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 ranef() and fixef(), but it would be a waste of time if such a function would already exist. Does any of you such a function? Regards, Rense Nieuwenhuis

...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? Thanks in advance, Rense Nieuwenhuis > version _ platform powerpc-apple-darwin8.6.0 arch powerpc os darwin8.6.0 system powerpc, darwin8.6.0 status major 2 minor 3.0 year 2006 month 04 day 24 svn rev 37909 language R ve...

...How is this possible? Is this a rounding-problem or is it the way R creates the linear model? Using SPSS results in b-values of zero. I'm using R 2.3 on MacOS 10.4.6, on a G4 processor. The syntax I used is printed below. Any explanations of this behaviour are really appreciated. Rense Nieuwenhuis #-------------------------- x <- c(0,1,1) y <- c(2,1,3) plot(x,y) b <- 1 length <- 100 for(i in 1:length) { x[1] <- i/length model <- lm(y~x) b[i] <- model$coefficients[2] } plot(b) #---------------- [[alternative HTML version deleted]]

...spond to the order in which the polygons are drawn. Since some of the full world-map isn't drawn this time, the color-vector doesn't correspond anymore. This results in the coloring of the wrong countries. Does anybody know of a way to solve this? Thanks very much in advance, Rense Nieuwenhuis [[alternative HTML version deleted]]

...er between men and women (so, on the first level). In order to model this, I want the 'Sex'-variable to be random on the first level, as described in Snijders & Bosker, page 110. Does anybody know if this is possible and how this can be done using R? Many thanks in advance. Rense Nieuwenhuis PS. Please excuse me for not providing a self-contained example. I couldn't find a data-set in the lme4-package that fitted my question. [[alternative HTML version deleted]]

...t to extract from the summary (or the full model) the part stating the random parameters. More specifically, I want to extract from the summary the following: (Intercept) 0.2869401 (Intr) typeSngl 0.2791040 -0.617 Residual 0.7302233 How could this be done? Thanks for the effort, Rense Nieuwenhuis Linear mixed-effects model fit by REML Data: Exam AIC BIC logLik 9158.56 9234.241 -4567.28 Random effects: Formula: ~type | school Structure: General positive-definite, Log-Cholesky parametrization StdDev Corr (Intercept) 0.2869401 (Intr) typeSngl...

...eral: VarFunc) of lme returns a single value, representing a coëfficiënt along which variance is increasing / decreasing? - In general: how can a variance-component / heteroscedasticity be made dependent on some level-2 variable (school level in my examples) ? Many thanks in advance, Rense Nieuwenhuis [[alternative HTML version deleted]]

...parameters. More specifically, I want to extract from the summary the following: (Intercept) 0.2869401 (Intr) typeSngl 0.2791040 -0.617 Residual 0.7302233 How could this be done, and how can the same be done using the lmer- function from the lme4-package? Thanks for the effort, Rense Nieuwenhuis Linear mixed-effects model fit by REML Data: Exam AIC BIC logLik 9158.56 9234.241 -4567.28 Random effects: Formula: ~type | school Structure: General positive-definite, Log-Cholesky parametrization StdDev Corr (Intercept) 0.2869401 (Intr) typeSngl...

...also looked at NIS, but the information I found was all several years old so I have concerns that it may not be the correct choice. I'd really appreciate thoughts and comments. I'm not looking for detailed instructions, but rather some advice on the overall plan. Sincerely, -- Cliff Nieuwenhuis Foresite Software LLC 608-393-1019

...t to extract from the summary (or the full model) the part stating the random parameters. More specifically, I want to extract from the summary the following: (Intercept) 0.2869401 (Intr) typeSngl 0.2791040 -0.617 Residual 0.7302233 How could this be done? Thanks for the effort, Rense Nieuwenhuis Linear mixed-effects model fit by REML Data: Exam AIC BIC logLik 9158.56 9234.241 -4567.28 Random effects: Formula: ~type | school Structure: General positive-definite, Log-Cholesky parametrization StdDev Corr (Intercept) 0.2869401 (Intr) typeSngl...

...nts moved at the same time, so weight-factors need to be applied to contexts Does any of you have ideas on how to estimate such models? Preferably using the NLME / LMER packages. In the literature these kinds of models are referred to as multiple membership models. Thanks in advance, Rense Nieuwenhuis

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

https://bugs.freedesktop.org/show_bug.cgi?id=100228 Bug ID: 100228 Summary: [NV137] unknown chipset (137000a1) Product: xorg Version: unspecified Hardware: x86-64 (AMD64) OS: Linux (All) Status: NEW Severity: normal Priority: medium Component: Driver/nouveau Assignee: nouveau

...fects (standardized and unstandardized) in the output. Can I create these effects using R? I know how to calculate them 'by hand', but I'd like to be able to do it automatically. If one of you is able to give some information, I'd be very gratefull, with high regards, Rense Nieuwenhuis

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

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