similar to: nlme Random Effects Specification

Hi all, I am running regressions with many covariates, most of which remain the same each time (control variables). Instead of writing 30 demographic variables every regression, is there a way I could call them all at once using a variable called, perhaps "demog"? I have tried: > demog <- list(age1, age2, age3) but I get an error when I try to call a list in a regression. I also

Dear All, I am using rhierMnlRwMixture in the bayesm package for the analysis of data from a choice experiment. I am trying to follow the margarine example set out in the bayesm manual (p.28). However, after several attempts I keep getting an error message with regards to my Z matrix as below. > Error in Z %*% t(matrix(olddelta, ncol = nz)) : >requires numeric/complex matrix/vector

In response to those who asked for a better explanation of what the Vilno software does, here's a simple example that gives some idea of what it does. LABRESULTS is a dataset with multiple rows per patient , with lab sodium measurements. It has columns: PATIENT_ID, VISIT_NUM, and SODIUM. DEMO is a dataset with one row per patient, with demographic data. It has columns: PATIENT_ID, GENDER.

On Windows, if you load a dll, this is locked. Hence, for package installations, close all R instances, start one without loading packages and then update packages. Best, Uwe Ligges On 26.08.2017 15:18, Bill Denney wrote: > Hi, > > > > When installing packages in Windows (currently using Windows 10 with all > service packs), occasionally, I get a warning similar to the

I think that this response should be added to R for Windows FAQ 3.5. -- Sent from my phone. Please excuse my brevity. On August 26, 2017 11:45:55 PM PDT, Uwe Ligges <ligges at statistik.tu-dortmund.de> wrote: >On Windows, if you load a dll, this is locked. >Hence, for package installations, close all R instances, start one >without loading packages and then update packages. >

Hi, When installing packages in Windows (currently using Windows 10 with all service packs), occasionally, I get a warning similar to the following: package 'Rcpp' successfully unpacked and MD5 sums checked Warning in install.packages : unable to move temporary installation 'C:\Users\William Denney\Documents\R\win-library\3.4\file32701900456\Rcpp' to

Below you will find the output from a failed multi-level model run. I am trying to estimate the following model: Pr(PLFP=1)= logistic regression -> B1_j * bm + B2_j * wm + B3_j bf + B4_j wf + B5 yrsed+ B6 age+ B7 age^2+e_ij B1_j = G01 + G11 bmxd + d1 B2_j = G02 + G12 wmxd + d2 B3_j = G03 + G13 bfxd + d3 B4_j = G04 + G14 wfxd + d4 d1-d4 freely correlated Note that there is no

Hello All, I am generating some plots where the title is generated with substitute and paste. An example: nval <- 20 plot(0,0) title(substitute(paste("n = ", n), list(n = nval))) But when compared to: plot(0,0) title("n = 20") the title in the first plot looks slightly different (it is not in bold). How can I get the two titles to look exactly the same? Also, how can I

Hello everybody, I have been recently using reshape to convert "long" data to "wide" data. Everything was going well until I reached some problematic datasets. It has taken me a couple of weeks to finally figure out what might be happening. The problem is reproducible with test cases, and on two versions of R (Windows 2.2.0 and x86-64 Fedora Core 3 R 2.2.0). The data

Dear all, I am using the nlme package to find the best model using non-linear mixed effects model. After the all the other things needed, I did > M1<- nlme(NShrubs~A*exp(-B*CED),data=DataSg,fixed=list (A~SITE,B~1), random=B~1,start=list(fixed=c (1.66,1.66,1.66,1.66,1.66,1.66,1.66,1.66,0.0038)),control=lmc) and it was ok. However, the data showed to need random parameters at both

Hi all, I'm running a multidimensional scaling model on data that are first transformed into a similarity matrix using the program daisy in the cluster package. If I already provide the similarity matrix, the plotting routine will give me the variable names on the plot itself. But, it seems that when I use the daisy program to create the similarity matrix, I can't get the program

Hello, I wuold like to propose my modifications of the original cor.test to you : I tried to calcolate the correct p-value for Spearman and Kendall's test with ties. Let me know what you think. Thanks you for your time. Antonietta di Salvatore test <- function(x, ...) UseMethod("test") test.default <- function(x, y, alternative = c("two.sided",

Dear List I have math test scores for male and female students where gender is a dummy code (female =1). I also have a variety of other demographic variables. However to begin, I want to create a very simple stripchart where female math scores are a blue circle and male scores are a red triangle. I am having difficulty using conditional statements to accomplish this. Thank you. ------

Hi everyone, I'm having some difficulty getting "simple effects" for the ols() function in the rms package. The example below illustrates my difficulty -- I'll be grateful for any help. #make up some data exD <- structure(list(Gender = structure(c(1L, 2L, 1L, 2L, 1L, 1L, 1L, 2L, 1L, 2L, 2L, 2L, 1L, 2L), .Label = c("F", "M"), class = "factor"),

Dear R people, I am having a problem with modeling the following SAS code in R: Class ID Gr Hemi Region Gender Model Y = Gr Region Hemi Gender Gr*Hemi Gr*Region Hemi*Region Gender*Region Gender*Hemi Gr*Hemi*Region Gender*Hemi*Region Gr*Gender*Hemi*Region Random Intercept Region Hemi /Subject = ID (Gr Gender) I.e., ID is a random effect nested in Gr and Gender, leading to ID-specific

Hi folks, I am tearing my hair out on this one. I am using an example from Pinheiro and Bates. ### this works data(Orange) mod.lis <- nlsList(circumference ~ SSlogis(age, Asymp, xmid, scal), data=Orange ) ### This works mod <- nlme(circumference ~ SSlogis(age, Asymp, xmid, scal), data=Orange, fixed = Asymp + xmid + scal ~ 1, start =

I've attached a perl script that I wrote to build smbpasswd from /etc/passwd. Please note that this writes a samba-2.0 format smbpasswd file. This script will: # - create accounts that exist in /etc/passwd and not in smbpasswd # - remove accounts that exist in smbpasswd and not in /etc/passwd # - preserve existing password, flags and LCT fields in smbpasswd # - sets new account passwords to

Pfizer is hiring a statistician supporting early drug discovery for Global Medicinal Chemistry and Pharmacokinetics, Dynamics and Metabolism (PDM). If interested, go to http://pfizercareers.com/ and search for job ID 985944. The location is Groton, Connecticut. Responsibilities The statistician will have a consulting-type role supporting a large pool of scientists in two large platform line

Full_Name: Dirk Koschuetzki Version: 2.1.1 OS: source code Submission from: (NULL) (194.94.136.34) Hello, >From the source code (R-2.1.1, file: .../R-2.1.1/src/library/stats/R/) ****************************** cor.test.default <- function(x, y, alternative = c("two.sided", "less", "greater"), method = c("pearson", "kendall",

Hi all, I have collected response time data from 178 participants ('sub') for each combination of 4 within-Ss factors ('con','int','tone','cue'). Additionally, I have recorded the gender of each participant, so this forms a between-Ss factor ('gender'). Normally this would be analyzed using aov: