similar to: saving object function

Hi all, I appreciate the help this list has given me before. I have a question which has been perplexing me. I have been working on doing a Bayesian calculating inserting studies sequentially after using a non-informative prior to get a meta-analysis type result. I created a function using three iterations of this, my code is below. I insert prior mean and precision (I add precision manually

My question is about the "Signif. codes" and the p-value, specifically, the output when I run summary(nameofregression.lm) So you get this little key: Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 And on a regression I ran, next to the intercept data, I get '***' Coefficients: > > Estimate Std. Error t value Pr(>|t|) > >

Is there anyway to make the lines in the dispersion command come forward in a plot and allow the fill in the dispersion parameter be transparent so all of the lines I am using to note confidence intervals are shown? ------------------------------------------- Joe King, M.A. Ph.D. Student University of Washington - Seattle 206-913-2912 jp@joepking.com

I wasn't sure about the subject so I am sorry about the vagueness, but if I have a vector of values, how do I get a certain type of vectors to be bold or italics? So let x be a vector from -10 to 10, how would I tell it to bold the x values that are greater than 0? Or like absolute value of 5? Is it possible to do this in like an if else statement? Could I use x<-seq(-10,10,1)

In the "filled bands" part of xYplot of the Hmisc package, is there a way to have multiple bands with multiple lines? or does it just allow one for now? So I had an example bit ago had a made up line and CI, now if I wanted to make a second line with a CI filled in can I put them on the same plot? x<-seq(1,10,1) y<-seq(1,10,1) ci<-y*.10 ciupper<-y+ci

Dears, I'm a college student and In doing my statistics homework. I use R with SEM package as my tool for sem analysis, but my teacher told me AMOS is more suitable for such analysis. Could someone help tell me whether it is true that some commercial software is better accepted in academic fields? Sorry if I should not post such topics here. -- Best Regards, Reeyarn T. Lee Accounting

I am wanting to plot a 95% confidence band using segplot, yet I am wanting to have groups. For example if I have males and females, and then I have them in different races, I want the racial groups in different panels. I have this minor code, completely made up but gets at what I am wanting, 4 random samples and 4 samples of confidence, I know how to get A & B into one panel and C&D in to

Ok so I am trying to add a subtitle to my xYplot using Hmisc and I have some dummy code of x<-seq(1,10,1) y<-rev(seq(1,10,1)) ci<-y*.10 ciupper<-y+ci cilower<-y-ci this code works fine: xYplot(Cbind(y,ciupper,cilower)~x, type=c("b"),plot.points = TRUE, method='filled bands', main='main title') but when I add sub=. at the end and use this

Hello, I want to model my data with the following model: Y1=X1*coef1+X2*coef2 Y2=X1*coef2+X2*coef3 Note: coef2 appears in both lines Xi, Yi is input versus output data respectively How can I do this in R? I got this far: lm(Y1~X1+X2,mydata) now how do I add the second line of the model including the cross dependency? Your help is greatly appreciated! Cheers, Bastiaan

Dear all, I''d like to ask help on R code to get the same results as the following Splus code: >indata<-importData("/home/data_new.csv") >indata[1:5,4] [1] 0930 1601 1006 1032 1020 I tried the following R code: > indata<-read.csv("/home/data_new.csv") > indata[1:5,4] [1] 930 1601 1006 1032 1020 I''d like the first

BAYESIAN INFERENCES FOR MILKING TEMPERAMENT IN CANADIAN HOLSTEINS Hi All, I have a data set "x" with several variables. Sample of the data is shown below V1 v2 v3 v4 5 6 9 10 3 4 7 10 4 6 10 18 I want the frequency of each data point sorted by their occurrence. Below is the output that I want 10 =3 6=2 4=2 9=1 5=1 7=1 3=1 How do

Hello all, if I have multiple imputed data sets, is there a command or function in R in any package you know of to combine those, I know one common MI approach is rubins rules, is there a way to do this using his rules or others? I know theres ways, like using Amelia from Gary King's website to create the imputed data sets, but how to make them into one or combine them for analysis.

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20091122/430e3297/attachment-0001.pl>

So many matrices are square symmetrical (i.e. variance-covariance matrices), is there any way to get R to split the matrix on its diagonal and just return one diagonal? So if I have mat<-matrix(c(1,4,3,4,1,2,3,2,1), nrow = 3, ncol=3, byrow=TRUE) is there anyway to get the lower right diagonal instead of the entire symmetric matrix? -------------------------------------------

Is there a way to get confidence intervals around an ordered regression like polr() in the MASS package? ------------------------------------------- Joe King, M.A. Ph.D. Student University of Washington - Seattle 206-913-2912 <mailto:jp@joepking.com> jp@joepking.com ------------------------------------------- Ad astra per aspera - "Through hardships to the stars"

I have a loop that I would like to use to extract the "stddev" for each itteration so I can average the "stddev" for all the runs. It would be helpful to know how to extract the "stddev" for each run from the VarCorr. Thanks MCruns<-1000 sighatlvec<-rep(NA,MCruns) sighatbvec<-rep(NA,MCruns) sighatevec<-rep(NA,MCruns) for(mc in 1:MCruns) {

Hi all, I am trying to run a two factor anova, but one of the factors is a random factor, now I am also running in SPSS and it seems its dividing by the wrong term to get the appropriate F term. here is my data. In SPSS the F scores about double the ones in R, how can I specify one of my factors as a random factor or change it to where it does the right model fitting? I am using the lm command

Hi all, i am still using the Hmisc package and I like the filled bands part, is there a way though to have different groups have different color of bands, maybe a lighter version of the color of the line that is used? ------------------------------------------- Joe King, M.A. Ph.D. Student University of Washington - Seattle 206-913-2912 jp@joepking.com

Hi all, in R I have?Sys.timezone() function to get the current working Time zone. However I want to have a vector to get the list of all available time zones, like say, LETTERS gives me all letters. Is there any function like that? Thanks in advance,

Hi everbody this is my first post, I have two test server .(Both of them is idle) db1 : RHEL4 OCFS2 db2 : RHEL3 OCFS I test the IO both of them The result is below. db1(Time Spend)db2(Time Spend)OS Test Command dd (1GB) (Yazma)0m0.796s0m18.420stime dd if=/dev/zero of=./sill.t bs=1M count=1000 dd (1GB) (Okuma)0m0.241s8m16.406stime dd of=/dev/zero if=./sill.t bs=1M count=1000 cp