search for: bliese

...tion(x,y){ panel.xyplot(x,y) lines(x,fitted(lm(y~poly(x,1),na.action=na.omit))) lines(x,fitted(lm(y~poly(x,2),na.action=na.omit)),lty=3)}, xlab="Time", ylab="Average Reaction Time") Happy to send the very small TEMP file for the curious... Any suggestions? Paul MAJ Paul Bliese, Ph.D. Walter Reed Army Institute of Research Phone: (301) 319-9873 Fax: (301) 319-9484 paul.bliese at na.amedd.army.mil

...ot;c:\\temp\\file1b.csv",header=T)) What is tedious about this solution is that we have to change the file name in step 3 every time. Is there a way to have R identify all the files in a directory and create one big data frame? I'm working in Windows with R 1.6.2. Thanks Paul MAJ Paul Bliese, Ph.D. Walter Reed Army Institute of Research Phone: (301) 319-9873 Fax: (301) 319-9484 paul.bliese at na.amedd.army.mil

R Version 2.2.0 Platform: Windows When I use barplot but select a ylim value greater than zero, the graph is distorted. The bars extend below the bottom of the graph. For instance the command produces a problematic graph. barplot(c(200,300,250,350),ylim=c(150,400)) Any help would be appreciated. Paul [[alternative HTML version deleted]]

....6708204 0.5 0.2236068 Where, of course, column 1 is linear, column 2 Quadratic and 3 cubic. My experience is that the coding scheme used in R works "better" than the integer scheme discussed in Kirk for many regression type analyses. Can anyone enlighten me as to why? Thanks, Paul Bliese Walter Reed Army Institute of Research -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the...

We have been using Least Angle Regression (lars) to help identify predictors in models where the outcome is continuous. To do so we have been relying on the lars package. Theoretically, it should be possible to use the lars procedure within a general linear model (glm) framework - we are particular interested in a logistic regression model. Does anyone have examples of using lars with logistic

...ference. The reason why the contrasts look different is that R is rescaling the variables to have unit length. The unit length of the first row, for instance is 20, so R divides -3, -1, 1, 3 by sqrt(20) and returns: > c(-3,-1,1,3)/sqrt(20) [1] -0.6708204 -0.2236068 0.2236068 0.6708204 Paul Bliese -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat....

...ion in MASS to create a matrix containing random normal variables and dichotomize these variables into 0,1; however, this is a less than ideal solution as my observed correlation matrix is downwardly biased and the amount of the bias is related to the prevalence of each variable. Thanks, Paul D. Bliese Heidelberg, Germany COMM: +49-6221-172626

Looking for help with the following problem. Given a sample of zeros and ones, for example: > VECTOR1<-rep(c(1,0),c(15,10)) > VECTOR1 [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 0 0 0 0 0 0 0 0 0 How would I create a new sample (VECTOR2) also containing zeros and ones, in which the phi-coefficient between the two sample vectors was drawn from a population with a known

...elect any individuals from group 2. My fundamental question is: "What are people's thoughts about using bootstaping in nested data? Does one have to sample with replacement taking into consideration the group structure in the data?" If so, any suggestions on how to do this? Paul Bliese US Army Medical Research Unit - Europe

...or possible solutions. Paul -----Original Message----- From: r-help-bounces at stat.math.ethz.ch [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of David Duffy Sent: Monday, September 12, 2005 1:34 AM To: r-help at stat.math.ethz.ch Subject: [R] Simulate phi-coefficient > From: "Bliese, Paul D LTC USAMH" <paul.bliese at us.army.mil> > > Given a sample of zeros and ones, for example: > > VECTOR1<-rep(c(1,0),c(15,10)) > How would I create a new sample (VECTOR2) also containing zeros and > ones, in which the phi-coefficient between the two sample vect...

...n the example, notice that the SE estimates are much closer to the SE values in the original model without weights. > model<-glm(y1~x1,weights=w) > varmat<-infjack.glm(model, groups=1:5) > se<-sqrt(diag(varmat)) > se (Intercept) x1 0.7827224 0.4969040 MAJ Paul Bliese, Ph.D. Walter Reed Army Institute of Research Phone: (301) 319-9873 Fax: (301) 319-9484 paul.bliese at na.amedd.army.mil -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info&quot...

...1 2 2 4 1 3 3 3 1 4 4 2 1 5 5 1 NA > TEMP[,c(1,3)][TEMP[,c(1,3)]==1&!is.na(TEMP[,c(1,3)])]<-10 Error in "[<-.data.frame"(*tmp*, TEMP[, c(1, 3)] == 1 & !is.na(TEMP[, : matrix subscripts not allowed in replacement Thanks! Paul Bliese Walter Reed Army Institute of Research

Recently, I was working with some lagged designs where a vector of observations at one time was used to predict a vector of observations at another time using a lag 1 design. In the work, I noticed a lot of negative correlations, so I ran a simple simulation with 2 matched points. The crude simulation example below shows that the correlation can be -1 or +1, but interestingly if you do this

Does anyone know of another version of the Mcnemar test that provides: 1. Odds Ratios 2. 95% Confidence intervals of the Odds Ratios 3. Sample probability 4. 95% Confidence intervals of the sample probability Obviously the Odds Ratios and Sample probabilities are easy to calculate from the contingency table, but I would appreciate any help on how to calculate the confidence

...m TEMP2<-make.univ(x=TEMP,dvs=TEMP[,c(10,13,16)]) #transforming it to univariate form again xyplot(MULTDV~TIME|as.factor(SUBNUM),data=TEMP2,type=c("p","r","g"),col="blue",col.line="black",xlab="Time",ylab="SAT") # taken from Bliese Paul ? Multilevel Modeling in R Now I want to be able to identify those SUBNUM that have a downwards trend like (for eg SUBNUM 7) without picking them by inspecting the plot. Is there some way how I could access these coefficients for this apparently regression line? Thank you very much and have...

...these basically the same?). I would have expected lme and gee to provide more similar answers. Any thoughts on why gee and lme are not more similar would be appreciated as would any clarification about when ignoring nonindepedence leads to too small versus too large SE values.... Paul MAJ Paul Bliese, Ph.D. Walter Reed Army Institute of Research Phone: (301) 319-9873 Fax: (301) 319-9484 paul.bliese at na.amedd.army.mil

I have a case where I would like to change multiple columns containing numbers to factors. I can change each column one at a time as in: TEMP.FACT$EXPOS01<-factor(TEMP.FACT$EXPOS01,levels=c(1,2,3),labels=c("No ne","Low Impact","MedHigh Imp")) TEMP.FACT$EXPOS02<-factor(TEMP.FACT$EXPOS02,levels=c(1,2,3),labels=c("No ne","Low

...nction to start with in the first place? I'm wondering whether there might be anyone knowing how to conduct a meta-analysis based on correlations in R? Thanks in advance Sebastian P.S.: Of course, I'm dreaming of such a step-by-step-script like the absolutely marvellous ones provided by Bliese for multilevel-analysis in R :-) --------- Dipl.-Psych. Sebastian Stegmann Managing Editor, British Journal of Management Goethe University Institute of Psychology Department of Social Psychology Kettenhofweg 128 60054 Frankfurt am Main Germany http://www.sozialpsychologie.uni-Frankfurt.de/ Phone:...

Recent changes to read.spss() in the foreign package return a dataframe containing additional attributes. For example, >TEMP<-read.spss(choose.files(), to.data.frame=T,use.value.labels=F) > str(TEMP) `data.frame': 780 obs. of 8 variables: $ EXPOS01: atomic 1 1 2 1 2 3 2 4 2 1 ... ..- attr(*, "value.labels")= Named num 5 4 3 2 1 .. ..- attr(*,

...ta(mydata)" command returns an error saying the file cannot be found. It appears to be looking in "c:\documents...\temp" instead of in the data subdirectory assocatied with the package. Functions and help files seem to work fine. Operating system: Windows 2000 Any thoughts? Paul Bliese Walter Reed Army Institute of Research -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the...