similar to: plotting residuals

Is there an easy way of converting an aov.summary into a matrix in which the rows are the factor names and the columns are Df, Sum Sq, Mean Sq, F value and Pr. For example, convert Df Sum Sq Mean Sq F value Pr(>F) block 5 343.29 68.66 4.4467 0.015939 * N 1 189.28 189.28 12.2587 0.004372 ** P 1 8.40 8.40 0.5441 0.474904 K

Hi, I am working on table suing R and Latex. I am writing .Rnw file first in which i m reading input file and storing into dataframe. After filtering certain values from this dataframe. I am planning to display it. I don't want to use xtable since i need to change column names. *Sample .Rnw file* @ <<echo=FALSE>>= cat("\begin{table}[ht]")

I thought more about this last night after my email about putting an end to "Read the Manual" replies. It seems that there are a handful of R Super-Programmers. If I were one, I would get tired of 100's of questions like "How can I get my Excel file into R!!!". Concurrently, R is experiencing an exponential increase in new users to the extent that there are now

Hi all, I hope someone can help me with this. Suppose I import a text file and one of the columns looks like this: New York New York England Spain Spain Orlando New York England France I want to add a variable that is based on the previous one US US Europe Europe Europe US US Europe Europe How do that? Also, I would like to be able to export the data as

Dear R user, I am trying to reproduce the results in Montgomery D.C (2001, chap 13, example 13-1). Briefly, there are three suppliers, four batches nested within suppliers and three determinations of purity (response variable) on each batch. It is a two stage nested design, where suppliers are fixed and batches are random. y_ijk=mu+tau_i+beta_j(nested in tau_i)+epsilon_ijk Here are the

Dear Madame, Dear Sir, I am able to obtain the coefficients from a 'summary' of 'lm', but NOT from a 'summary' of 'aov'. The following example shows my steps. ## Initialize rm(list = ls()) # remove (almost) everything in the working environment utils::data(npk, package="MASS") # get data model <- yield ~ block + N*P*K ## Using lm npk.lm <-

I would appreciate it if someone could explain the results of the example from the aov() help file. The output given below is different from book Venables and Ripley - MASS The results In R1.5.1 (Under windows) is as follows: > N <- c(0,1,0,1,1,1,0,0,0,1,1,0,1,1,0,0,1,0,1,0,1,1,0,0) > P <- c(1,1,0,0,0,1,0,1,1,1,0,0,0,1,0,1,1,0,0,1,0,1,1,0) > K <-

Dear Prof. Ripley and Christoph, thank you very much for your comments. You have helped me a lot. Thanks, Tomas Goicoa >Dear Prof. Ripley > >Thank you for your email. Yes, this is of course the correct >syntax to save us the extra calculation. And I forgot the >"lower.tail = FALSE" for pf() in my example to obtain the >p-value. > >Thank you for the

Hi all, I'm sure this is simple but I don't get it. I have a table mytable<-c(rep("Disagree",37),rep("Agree",64)) table(mytable) this gives me Agree Disagree 64 37 but I didn't ask for it to be in alphabetic order. How can I get it in original order? Disagree Agree 37 64 Thanks, Jeff Jeffrey. M. Miller, PhD

Hi, Is it possible to configure NTLM HTTP proxy authentication using the winbind/squid "ntlm_auth" helper, to authenticate users against a Samba 3 server? I already have the NTLM auth working against a Windows 2003 Active Directory, but I also have a completely separate Samba 3 server that I would also like to configure NTLM proxy authentication against. Please advise, as I can't

Hi, I would like to be able to access the calculated results from the nls package. Using the example in R, fm3DNase1 we can reurn certain parts of the calculations: > coef(fm3DNase1) Asym xmid scal 2.345179 1.483089 1.041454 > resid(fm3DNase1) [1] -0.0136806237 -0.0126806237 0.0089488569 0.0119488569 -0.0025803222 [6] 0.0064196778 0.0026723396 -0.0003276604

Hi all, I would like to know a way of calculating the distance between 2 points when I already have the longitude and latitude of the points. For example : Point 1 : 52? 9'54.00"N 4?25'8.40"E Point 2 : 52?27'46.80"N 4?33'18.00"E Distance between point 1 and point in km .... Is there any functions already available for this ? Any help will be much

Hi! When checking validity of a model for a large number of experimental data I thought it to be interesting to check the information provided by the summary method programmatically. Still I could not find out which method to use to get to those data. Example (not my real world data, but to show the point): [BEGIN] > DNase1 <- subset(DNase, Run == 1) > fm1DNase1 <- nls(density ~

Well I have installed Ubuntu 8.40, Wine 9.59, and World of Warcraft, and it took me some time to get it working, but I finally did. Now, however, I'm having several problems with the way it runs. Namely: -the mouse is very slow and skippy at moving when in WoW -lots of objects in the game fail to render (the ground usually doesn't render either) -my FPS in the middle of nowhere, with

Hello, I'm writing from the otherside of the world from where my systems are, so details are coming in slow. We have a 6TB OCFS2 volume across 20 or so nodes all running OEL5.4 running ocfs2-1.4.4. The system has worked fairly well for the last 6-8 months. Something has happened over the last few weeks which has driven write performance nearly to a halt. I'm not sure how to proceed, and

Hi, (R version 2.15.0) I am running a pgm with 1 response (earlier standardized Y) and 44 independent vars (Xi) from the same data =a2: When I run the 'lm' function on single Xi at a time, the beta coefficient for let's say X1 is = -0.08 (se=0.03256) But when I run the same Y with 44 Xi-s with the 'step' function (because I left direction parameter empty, I assume a backward

Dear list: I am interested in the following sort of problem, as is found frequently in the field of QSAR. I have biological activity as a function of chemical structure, with structure defined in a categorical manner in that the SUBSTITUENT is the levels of the POSITION factor. For example, data from Kubinyi (http://www.kubinyi.de/dd-12.pdf) for this type of analysis is presented as follows:

All, I have just installed Samba on a Solaris 8 server and am trying to share two unix directories to an NT domain. The smb.conf file is set up as follows :- [spain]/usr/local/samba/lib $ more smb.conf [global] workgroup = coda username map = /usr/local/samba/lib/users.map hosts equiv = /etc/hosts.equiv [build] comment = autobuild path = /users/osproc/autobuild public = yes

Hi all, Does anyone know if an R function for the Mantel-Haenszel chi-square for a 2x2 table exists? I’ve also seen it called the randomization Q statistic. Note that I’m not looking for the Cochran-Mantel-Haenszel…I did see that out there as cmh.test. Thanks in advance, Jeff Internal Virus Database is out-of-date. Checked by AVG Free Edition. 11:27 AM [[alternative HTML version

Dear list, after fitting an lrm with the Design package (stored as "mymodel") I try running a summary, but I get the following error: dim(mydata) [1] 235 9 names(mydata) [1] "id" "VAR1" "VAR2" "VAR3" "VAR4" "VAR5" "VAR6" "VAR7" "VAR8" summary(mymodel) Error in `contrasts<-`(`*tmp*`,