search for: blomberg

...iables. The design is a 2*2 factorial with 10 replicates per treatment. Looking at the code for summary.manova, the error involves a problem with qr(). Does anyone have a suggestion as to how to deal with this error? The analysis works fine when there are only 32 variables. Thanks, Simon. Simon Blomberg, PhD Depression & Anxiety Consumer Research Unit Centre for Mental Health Research Australian National University http://www.anu.edu.au/cmhr/ Simon.Blomberg at anu.edu.au +61 (2) 6125 3379

...here the corSymm correlation structure is modelled, rather than known. I have also searched the R archives, but no luck. I think it should be of the form gls(response~A*B*C, data=dat, correlation=corSymm(...?)) but I don't understand the arguments to corSymm. Thanks in advance, Simon. Simon Blomberg, PhD Depression & Anxiety Consumer Research Unit Centre for Mental Health Research Australian National University http://www.anu.edu.au/cmhr/ Simon.Blomberg at anu.edu.au +61 (2) 6125 3379

Could anyone please explain how to set up contrasts between means in R. I want to know if "before I conduct an experiment and believe the mean for 1 and 2 will be different from means 3 and 4, Is this true?" That is what I have to prove or disprove, I thought that contrasts would be the way to go. Thanks for the help. Igor [[alternative HTML version deleted]]

Dear Lister, I understand merge() can be used to join 2 data frames based on 1 variable. But how about merge based on more than 2 variables? Thank you so much! -- WenSui Liu (http://spaces.msn.com/statcompute/blog) Senior Decision Support Analyst Health Policy and Clinical Effectiveness Cincinnati Children Hospital Medical Center [[alternative HTML version deleted]]

Hello, can anyome tell me how to access the full script of the panel.smooth function so that I can change the thickness of the smoothing line or its colour ? All I could access is : > panel.smooth function (x, y, col = par("col"), bg = NA, pch = par("pch"), cex = 1, col.smooth = "red", span = 2/3, iter = 3, ...) { points(x, y, pch = pch, col = col, bg =

Let me modify my question about user-defined covariance structures for LME models: Can somebody tell me how I can see the code for the definition of the correlation structures that come with the NLME package. Specifically I like to see the code for the functions coef, corMatrix, and intialize for any of the pre-defined correlation structures, and use this as a template to define a new correlation

...ces to "R Development Core Team", and I received no hits. So I tentatively conclude that our paper is the first Nature paper to cite R. A great many thanks to the R Development Core Team for R, and Prof. Bates for lmer. Cheers, Simon. (I'm off to the pub to celebrate.) -- Simon Blomberg, B.Sc.(Hons.), Ph.D, M.App.Stat. Centre for Resource and Environmental Studies The Australian National University Canberra ACT 0200 Australia T: +61 2 6125 7800 email: Simon.Blomberg_at_anu.edu.au F: +61 2 6125 075...

Forgive me please if this has been answered previously, but the r-help archive doesn't have a "search" capability. It's a pity, as I'm sure that there's a wealth of informative help in there. I am unable to get the "cloud" function to produce a fundamental 3D scatterplot. I have a 3-column numeric matrix: >is.matrix(sanity.MIF); is.numeric(sanity.MIF);

>From time to time people request symbolic computations beyond what D() and deriv() etc can provide. A brief look at the internet shows that there are many more or less developed computer algebra packages freely available. Therefore, I wondered if it would be an idea to try to 'integrate' one of these packages in R, which I guess can be done in more or less elegant ways... I do not know

I am looping over many data files and reading in the data with F <- read.table(filename) to read in a 22000 by 15 matrix. Works fine on the first matrix F, but I get the following error when the second file is read into F: Error in "row.names<-.data.frame"(*tmp*, value = row.names) : duplicate row.names are not allowed I have tried picking a column of the matrix and

Hello 1/ 'priors' is a table looking like: "W123" "T678" "S789" 23 42 11 12 35 9 etc 2/ WBS <- labels(priors) gives me a result of class list and length 1 looking like: "W123" "T678" "S789" I want to read W123 into X[1] as W, T687 into X[2] as T and S789 into X[3] as S using

Hi R-people, I am wanting to run Factorial ANOVA followed by Scheffe tests on some spatial subjective data. I'm comparing X-Y independent coordinates against x-y dependent coordinates. There are only four independent spatial coordinates that form a square. I am wondering whether I am doing the right thing, because there doesn't seem to be a simple way of doing this. I have attempted to

Dear there, I am wondering how to test whether a simple linear regression model (e.g. y=1.05x) is significantly different from a 1 to 1 line (i.e. y=x). Thanks. Regards, Jin [[alternative HTML version deleted]]

...a Windows XP Pro host. Thanks, Simon. > sessionInfo() Version 2.3.1 (2006-06-01) i386-pc-linux-gnu attached base packages: [1] "methods" "stats" "graphics" "grDevices" "utils" "datasets" [7] "base" -- Simon Blomberg, B.Sc.(Hons.), Ph.D, M.App.Stat. Centre for Resource and Environmental Studies The Australian National University Canberra ACT 0200 Australia T: +61 2 6125 7800 email: Simon.Blomberg_at_anu.edu.au F: +61 2 6125 0757 CRICOS Provider # 00120C

Hello, My name is Rodrigo, I am using R program and I have a trouble. I am trying to do a dendrogram with genetics information. Let me explain... The Similarity Matrix was already did, and with this matrix I want to construct a dendrogram. So, the distance is done. I need to transform this matrix (that I have) in a dendrogram, I woud be very grateful if someone could help me. PS: I am sending

...basic understanding of statistics. Jason -- Assistant Professor of Statistics Mathematics and Computer Science Department University of Richmond, Virginia 23173 (804) 289-8081 fax:(804) 287-6664 http://www.mathcs.richmond.edu/~wowen "This is R. There is no if. Only how." Simon (Yoda) Blomberg

...S=en_AU.UTF-8;LC_PAPER=en_AU.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_AU.UTF-8;LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base loaded via a namespace (and not attached): [1] tools_2.8.0 > Simon. -- Simon Blomberg, BSc (Hons), PhD, MAppStat. Lecturer and Consultant Statistician Faculty of Biological and Chemical Sciences The University of Queensland St. Lucia Queensland 4072 Australia Room 320 Goddard Building (8) T: +61 7 3365 2506 http://www.uq.edu.au/~uqsblomb email: S.Blomberg1_at_uq.edu.au Policie...

Thanks for telling me that you could not get my message, I hope this work better... so my question was: I built a population matrix to which I applied the fonction eigen in order to find the main parameters about my population. I know that the first eigen value correspond to lambda or exponential growth rate of my population. My problem is that I want to have the 95% confidence interval of the

...s the rate I get from the root class of my queue. In both cases, the rate is calculated by rrdtool, I submit the number of bytes that have been transmitted. Shouldn''t the rate of packets leaving the root class be the same as the rate given by snmp? Thank you for the help, regards, David Blomberg

...and 1:20 where 1:20 is the default. Then I use iptables to mark all packets I send out on eth1. I then filter marked packets into class 1:10. I expected all packets sent on eth1 to end up in class 1:10, but some packets still go to 1:20. Did I do it wrong? Thank you for any help. regards, David Blomberg My script: --------------------------------------------------------------------------- /sbin/tc qdisc add dev eth1 root handle 1: htb default 20 r2q 1 /sbin/tc class add dev eth1 parent 1: classid 1:1 htb rate 200kbit burst 15k /sbin/tc class add dev eth1 parent 1:1 classid 1:10 htb rate 100kbit...