similar to: Problem whit a piece of program

Hello, I'm working on the multiple categorical data (5-points scale) using linear mixed effect model and wondering if anyone knows about or works on the linear mixed effect model with ordered logit or probit link. I found that the "lmer" function in R is very flexible and supports various models, but not ordered logit/probit models. I may conduct my analysis by turning my DVs

Dear All, I'm having trouble on 2 laptops Lenovo B580 since upgrading to Centos6.5. ( Because it's a Lenovo I cannot switch the network card for a better supported network card. ) There on the latest kernel : root at jac network-scripts]# uname -a Linux jac.cawdekempen 2.6.32-431.1.2.0.1.el6.x86_64 #1 SMP Fri Dec 13 13:06:13 UTC 2013 x86_64 x86_64 x86_64 GNU/Linux With the help of

To enhance my understanding, and that of my students, I have a question about cox.zph in the survival package. If I have correctly gleaned the high-level point from the 1994 Biometrika paper of Grambsch and Therneau, it looks to me like cox.zph provides a mechanism to test for a simple trend in plots of a function of time, g(t) versus the scaled schoenfeld residuals and it also provides some

I downloaded R-Patched today (to see if another problem I want to ask about is still present or if its just me - as per the posting guide). I ran "tools/rsync-recommended" successfully. I then ran "configure --enable-R-shlib" successfully. Then make stops with the following error. gcc -shared -L/usr/local/lib -o tools.so text.o init.o Rmd5.o md5.o -L../../../../lib -lR

Hello comunity, I'm trying to find a similar function as decile of SPSS, NTILES (10) some of you know about that, I will appreciate your help in advance. In SPSS VARIABLES=Sales (A) /NTILES (10) /PRINT=NO for example if I have Next data input: case IdCust Sales Profit 1 265 140.81 314.31 2 266 1778.96 408.32 3 267 2663.66 820.85 4 268 2994.14 913.73 5 269 2185.58 686.64 6 271 105.21

Dear R useRs, I have a problem with nls.lm function of minpackl.lm package. I need to fit the Van Genuchten Model to a set of data of Theta and hydraulic conductivity with nls.lm function of minpack.lm package. For the first fit, the parameter estimates keep changing even after 1000 iterations (Th) and I have a following error message for fit of hydraulic conductivity (k); Reason for

Hi everyone (again), before you all start screaming that the reordering of factors has been discusse on several threads and is not particular to ggplot2, hear me out. I can easily reorder my x-axis factor in facet.grid() in ggplot2. What I cannot reorder are the factors represented on the strips. I can see that the graphs are changing, so I am afraid of what it is I am doing. Why is ggplot2

Hi, # For publications, I am not allowed to repeat the axes. I have tried to remove the axes using: # yaxt="n", but it did not work. I have not understood how to do this in ggplot2. Can you help me? # I also do not want loads of space between the graphs (see below script with Dummy Data). # If I could make it look like the examples on the (nice) examples page: #

Is names(model1$coef) what you're looking for? -----Original Message----- From: Kenneth Cabrera [mailto:krcabrer at epm.net.co] Sent: 29 November 2002 10:36 Cc: R-help at stat.math.ethz.ch Subject: [R] Obtaining the variable names of a glm object Hi, R users! Suppose I make a model like this:

I seem to remember reading somewhere about some style guide regarding R the number of comment characters (#) prior to the comment meaning something. Anyone know to what I'm referring? Where? ---------------------------------------------------------- SIGSIG -- signature too long (core dumped)

> From: Duncan Murdoch <murdoch at stats.uwo.ca> > To: ramasamy at cancer.org.uk > Date: Fri, 10 Mar 2006 08:33:09 -0500 > Subject: Re: [R] To improve my understanding of workspaces > Other than Emacs, I use the same work habits as Adai. An advantage of > this workflow is that almost everything is stored in text format, so it > is easy to compare different versions to

Dear all, I have created the following data (that you can run) in order to explain my problem: y <- rep(c(1,2), 8) id <- rep(1:8,each=2) x1 <- rep(c("-a","+a"), each = 8) x2 <- rep(c("-b","+b"), each = 2, times = 4) x3 <- rep(c("-c", "+c"), each = 4,2) df <- data.frame(cbind(id,y,x1,x2,x3)) If I do:

Hello there, I have a situation where I would like to select the first row of a particular factor for a data frame (data example below). So that is, I would like to select the first entry when the factor1 =A and then the first row when factor1=B etc. I have thousands of entries so I need some general way of doing this. I have a minimal example that should illustrate what I am trying to do. I am

Dear all, I apologize if this is a FAQ (seems a bit like one, but I didn't find anything). I'm looking for an easy way to cut one value out of a vector and shorten the vector accordingly. Something like: x <- c(1, 1, 0, 6, 2) throwaway(x[3]) which will return x = 1 1 6 2, with length(x) = 4. I know one could do this by hand, but then one would have to create a second vector y in a

hello all happy new year and hope you r having a good holiday. i would like to calculate the expectation of a particular random variable and would like to approximate it using a number of the functions contained in R. decided to do some experimentation on a trivial example. example ======== suppose x(i)~N(0,s2) where s2 = the variance the prior for s2 = p(s2)~IG(a,b) so the posterior is

Full_Name: Tanya Logvinenko Version: 1.7.0 OS: Windows 2000 Submission from: (NULL) (132.183.156.125) For unbalanced design, I ran into problem with ANOVA (aov function). The sum of squares for only for the second factor and total are computed correctly, but sum of squares for the first factor is computed incorreclty. Changing order of factors in the formula changes the ANOVA table. For the

Hi All, I have a 4-dimensional data. I'm using barchart() function from lattice package. The R code and data are below - code includes one for stack=TRUE and other for stack=FALSE. I would like to present the data in another form which would be plotting Factor3 levels (P, Q, R, S) as two stacked bars (side by side). Like, for each level of Factor1 there should be two bars: first bar showing

Dear R people, I have a couple of questions about post-doc analyses for 2 by 2 within subjects ANOVA. I conducted a psycholinguistic study that combined a 2 by 2 design and a latin square design. Specifically, I had 32 items each of which generated 4 conditions. Participants saw each of the 32 items only once: 8 in Condition A, 8 in B, 8 in C, and 8 in D. The table below serves as an example.

I can think of many brute-force ways to do this outside of R, but was wondering if there was a simple/elegant solution within R instead. I have a table that looks something like the following: Factor1 Factor2 Value A 11/11/2009 5 A 11/12/2009 4 B 11/11/2009 7 B 11/13/2009 8 >From that I need to generate all permutations of Factor1 and Factor2 and force a 0 for any combination that doesn?t

Hi All, I have googled and tried finding if someone has ever tried producing (Grouped + Stacked) Barplot. I couldn't find one. My data needs to be reshaped, but once it is done it would be something like this: Factor1 Factor2 Factor3 Value A X P 10 A X Q 20 A Y P 20 A Y Q 5 A Z P 20 A Z Q 10 B X P 20 B X Q 10 B