similar to: manova

Hi, I am using the sim1000G R package to simulate data for case/control study. I can not figure out how to manipulate this code to be able to generate 10% or 50% causal SNPs in R. This is whole code provided as example on GitHub: library(sim1000G) vcf_file = "region-chr4-357-ANK2.vcf.gz" #nvariants = 442, ss=1000 vcf = readVCF( vcf_file, maxNumberOfVariants = 442 ,min_maf =

Check out the proposed prices when this is approved. BELL CANADA REPORT ON THE ECONOMIC EVALUATION FOR THE TARIFF REVISION OF Bell Canada's Access Services Tariff Item 315 ? Zero-Dialed Emergency Call Routing Service (0-ECRS) *2 March 2006 TABLE OF CONTENTS Page 1.0 GENERAL 3 1.1 Purpose of the Study 3 2.0 SERVICE DESCRIPTION 3 2.1 Service Characteristics 3 2.2

Hi, I am trying to do multiple regression on a set of data using backward stepwise regression....however backward stepwise regression is critised for overfitting data. To actually observe the bias and to come up with a better method to use..Could you all stats experts kindly give me pointers to any alternative procedure (or references) to use over backward stepwise regression from your

Hi, I am trying to get samba 3.0.7 working with our win2k DC. I installed samba from the ports collection, so the kerberos library looks to be the heimdel version. I can use kinit to create a ticket and it authenticates against the DC just fine. However when I attempt to use "net ads join" it fails with the following response: [2004/11/09 16:32:30, 0] utils/net_ads.c:ads_startup(183)

Fully up-to-date CentOS 5.2 and Firefox browser. A previous thread had a problem with plugins disappearing. I've discovered, by accident, a causal relationship. Running a Thunderbird mail client (could happen with others? Haven't tested with Evolution yet.) I click a link that opens the browser. BANG! The plugins.dat is "emptied". Restored from my rsync backup, repeat, BANG!

Does anyone know of a package that includes the Modified Sims test [Gewerke, 1983, Sims, 1972]? This test is used in econometrics and is a kind of alternative to the Granger test [Granger, 1969], which is in the package lmtest. Thanks in advance, chris Refernces: Gewerke, J., R. Meese, and W. Dent (1983), "Comparing Alternative Tests of Causality in Temporal Systems: Analytic Results and

Hi Since a stepwise procedure for variable selection (as e.g. in SPSS) for a LDA is not implemented in R and anyway I cannot be sure, that all the required assumptions for e.g. a procedure using a statistic based on wilks' lambda, hold (such as normality and variance homogeneity) I would like to ask you, what you would recommend me: shall I e.g. define a criterion such as the error-rate

I added a bibliographic reference to one of my .Rnw documents, to wit: \begin{thebibliography}{9} \bibitem{Jose} Jos{\'{e}} C. Pinheiro and Douglas M. Bates, \emph{Mixed-Effects Models in S and S-PLUS}, Springer, 2000. \end{thebibliography} just before the \end{document}. When I run Sweave on the result these lines disappear without a trace. Is this intentional? Subquestion: is

Hi Berina, I'm not an expert on genetics. I haven't looked at the package. And I've only glanced at your question. So, this is probably not the best response. But as no one else has responded, here's some comments: (1) Have you checked if there's a function in the package to do what you want? The remainder of these questions assume that you have, and the answer is no. (2)

Dear R-list, i have 6 different sets of samples. Each sample has about 5000 observations, with each observation comprised of 150 baseline covariates (X), 125 of which are dichotomous. Roughly 20% of the observations in each sample are "treatment" and the rest are "control" units. i am doing propensity score matching, i have already estimated propensity scores(predicted

Hey all, I encountered a problem drawing a histogram. You can view the picture here: http://dl.dropbox.com/u/4836866/Bad_Histogramm.png What happens: the bars are drawn with different starting points, thus no straight zero-line is there. And bars are overlapping. (or sometimes apart from each other.) How it happens: hist(volcano, breaks=10) # and any other data This also happens with

Dear R-help, I'm using prop.test() to compute a confidence interval for a proportion under R version 2.4.1, as in: prop.test(x = 340, n = 400)$conf [1] 0.8103309 0.8827749 I have two questions: 1) from the source code my understanding is that the confidence interval is computed according to Wilson, E.B. (1927) Probable inference, the law of succession, and statistical inference. J. Am.

Hello, I've run into a problem which is both maddening and rather hard to replicate, therefore I wondered if someone might know of a plausible explanation for it. I couldn't find anything in the archives, though maybe I'm searching for the wrong thing. I'm calling some C code using .C, and get the vector I'm interested in back as the 7th location in the returned list.

Dear All, I am trying to perform MANOVA. I have table with 504 columns(species) and 36 rows) with two grouping (season and location) Zx <- Z[c(4:504)] Zxm <- as.matrix(Z) m<- manova(Zxm~Season*location, data=Z) when I do summary.aov, I get respond for each species but summary.manova summary.manova(m) :" residuals have rank" 24<501. What can it be the reason for this error

Hi all, I am experimenting the function "manova" in R. I tried it on a few data sets, but I did not understand the result: I used "summary(manova_result)" and "summary(manova_result, test='Wilks')" and they gave a bunch of numbers... But I need the Sum-of-Squares of BETWEEN and WITHIN matrices... How do I read off from the R's manova results? Any

Hi, I would like to conduct a MANOVA. I know that there 's the manova() funciton and the summary.manova() function to get the appropriate summary of test statistics. I just don't manage to specify my model in the manova() call. How to specify a model with multiple responses and one explanatory factor? If I type:

Hello, I had a couple questions about manova modeling in R. I have calculated a manova model, and generated a summary.manova output using both the Wilks test and Pillai test. The output is essentially the same, except that the Wilks lambda = 1 - Pillai. Is this normal? (The output from both is appended below.) My other question is about the use of MANOVA. If I have one variable which has a

I'm trying to learn to use manova(), and don't understand why none of the following work: > data(iris) > fit <- manova(~ Species, data=iris) Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : incompatible dimensions > fit <- manova(iris[,1:4] ~ Species, data=iris) Error in model.frame(formula, rownames, variables, varnames, extras,

Dear colleagues, I had a question wrt the car package. How do I evaluate whether a simpler multivariate regression model is adequate? For instance, I do the following: ami <- read.table(file = "http://www.public.iastate.edu/~maitra/stat501/datasets/amitriptyline.dat", col.names=c("TCAD", "drug", "gender", "antidepressant","PR",

Hello all, I have a question on manova in R: I'm using the function "manova()" from the stats package. Is there anything like a stepwise (backward or forward) manova in R (like there is for regression and anova). When I enter: step(Model1, data=Mydata) R returns the message: Error in drop1.mlm(fit, scope$drop, scale = scale, trace = trace, k = k, : no 'drop1'