similar to: PairedData: a new package for analysing paired (numerical) data

Dear R-Users, ade4TkGUI is a new package available on CRAN. It implements a Tcl/Tk graphical user interface (GUI) for the ade4 package. Only the most basic functions of ade4 have a GUI in this first version : classical one-table data analysis methods (PCA, COA, MCA, PCO, etc.), one table with groups of rows (BGA, WGA, DA), and two-tables analysis methods (Coinertia analysis, CCA, PCAIV).

Dear R-Users, ade4TkGUI is a new package available on CRAN. It implements a Tcl/Tk graphical user interface (GUI) for the ade4 package. Only the most basic functions of ade4 have a GUI in this first version : classical one-table data analysis methods (PCA, COA, MCA, PCO, etc.), one table with groups of rows (BGA, WGA, DA), and two-tables analysis methods (Coinertia analysis, CCA, PCAIV).

Hello All I am new to this list. I have a problem where for a single sample drawn from normal population, null hypothesis is that variance = k (say). Alternative hypothesis is variance > k. Now if we know the true variance, then I would like to calculate the sample size required to produce certain power (for some significance). How do I do this ? I thought of using pwr.chisq.test, and I

Hello, i' m looking for a function to compute geary statistic for spatial correlation. Thanks. -- St?phane DRAY --------------------------------------------------------------- Biom?trie et Biologie ?volutive - Equipe "?cologie Statistique" Universite Lyon 1 - Bat 711 - 69622 Villeurbanne CEDEX - France Tel : 04 72 43 27 56 Fax : 04 78 89 27 19 04 72 43 27 57 E-mail

Hi Guys, I am having a real hard time trying to figure out for microarry. Here is my code One-Sample t-Test dim(data.sub) [1] 10000 140 ##there are 10000 probesets and 140 columns hist(data.sub) ## Histogram. Identify if the probesets are normal distributed q<-rnorm(10000) ##generate 10000 random, normal distributed values qqplot(data.sub,q)) ##Show the plot of the probeset

Dear all, I would like to compute power of a square non symmetric matrix. This is a part of a simulation study. Matrices are quite large (e.g., 900 by 900), and contains many 0 (more than 99 %). I have try the function mtx.exp of the Biodem package: library(Biodem) m <- matrix(0, 900, 900) i <- sample(1:900, 3000, replace = T) j <- sample(1:900, 3000, replace = T) for(x in 1:3000)

Hello, I am trying to play with RODBC library and Excel Files. In my file (doubs.xls) there are 2 spreadsheets: > library(RODBC) > connection<-odbcConnect("Excel Files") > sqlTables(connection) TABLE_CAT TABLE_SCHEM TABLE_NAME TABLE_TYPE REMARKS 1 F:\\Th?se\\R\\Doubs NA Faune$ SYSTEM TABLE NA 2 F:\\Th?se\\R\\Doubs NA Milieu$

Hello list, i have a glm with family=binomial, link=logit but there is over-dispersion. So, in order to take into account for this problem i choose to do a glm with family=quasibinomial(). I'm not an expert on this subject and i ask if someone could validate my approach (i'm not sure for the tests) : quasi_glm(myformula,quasibinomial(),start=mystart) summary(quasi) # test t for

Hello list, (PC Windows 2000, R 1.2.3) i want to compile C source file for R on windows 2000. I'm working with Metrowerks CodeWarrior for creating and compiling C files. I've see that R.lib file must be obtain with R.ddl by Visual C++. I don't have Visual C++ (and i don't want to buy it !!) but i want to obtain R.lib file. Is there a solution to my problem ??? Any help would

Hi there, I'm trying to fit GEE models on binomial data. My problem is how to "select" a model if ever it is possible and how I decide a factor has to be drop out or not. More specifically: 1) Are QIC (Pan 2001) reliable and can be used in a model selection framework (Burhnam and Anderson 1998)? If yes, how can the QIC be estimated from available information in a gee object? 2) If

Hi there, I'm trying to fit GEE models on binomial data. My problem is how to "select" a model if ever it is possible and how I decide a factor has to be drop out or not. More specifically: 1) Are QIC (Pan 2001) reliable and can be used in a model selection framework (Burhnam and Anderson 1998)? If yes, how can the QIC be estimated from available information in a gee object? 2) If

A full set of course notes on 'Analysing spatial point patterns in R' is now available at <http://www.csiro.au/resources/pf16h.html> The course is based on the package 'spatstat'. It includes a brief introduction to R, a detailed introduction to the 'spatstat' package, and a discussion of statistical methodology. Adrian Baddeley

Dear all, I'm trying to produce two dotcharts side-by-side within a Sweave document. When I'm compiling this example: \documentclass{article} \begin{document} <<fig=T,width=8,height=4>>= par(mfrow = c(1, 2), cex = 0.7) for(i in 1:2) dotchart(1:10) @ <<fig=T,width=8,height=4>>= par(mfrow = c(1, 2), cex = 0.7) for(i in 1:2) hist(1:10) @

Dear Mag. Ferri Leberl, I'm using something like: ----------------------- tex.tex --------------------------- \documentclass{article} \usepackage{graphicx} \usepackage{fancyvrb} \newcommand{\Rlogo}{\protect\includegraphics[height=1.8ex,keepaspectratio]{Rlogo.pdf}} \newcommand{\myinput}[1] {\begin{scriptsize} \VerbatimInput[frame=single,label=#1]{#1} \end{scriptsize}} \title{The R logo,

Dear Mag. Ferri Leberl, I'm using something like: ----------------------- tex.tex --------------------------- \documentclass{article} \usepackage{graphicx} \usepackage{fancyvrb} \newcommand{\Rlogo}{\protect\includegraphics[height=1.8ex,keepaspectratio]{Rlogo.pdf}} \newcommand{\myinput}[1] {\begin{scriptsize} \VerbatimInput[frame=single,label=#1]{#1} \end{scriptsize}} \title{The R logo,

Hi all, In the next release of R (4.4) the option to obtain a paired t.test with the formula interface for "long" data has been removed: t.test(x ~ group, paired = TRUE) # now results in an error. Exploring how one might obtain a paired t.test, there seems to be some inconsistency between the formula and default interfaces as well as within the formula interface. ``` extra.1 <-

There's a funny inconsistency in how t.test handles paired=T or paired=F. If x and y parameters are lists, paired=F works, but paired=T doesn't. > lg=read.csv("my.csv") > a = subset(lg, condition=="a")["score"] > b = subset(lg, condition=="b")["score"] > t.test(a,b) > t.test(a,b, paired=TRUE) Error in `[.data.frame`(y, yok)

Hello All I am new to this list. I have a problem where for a single sample drawn from normal population, null hypothesis is that variance = k (say). Alternative hypothesis is variance > k. Now if we know the true variance, then I would like to calculate the sample size required to produce certain power (for some significance). How do I do this ? I thought of using pwr.chisq.test, and I

Hi, Assuming that your dataset is similar to the one below: set.seed(25) dat1<- data.frame(Algae.Mass=sample(40:50,10,replace=TRUE),Seagrass.Mass=sample(30:70,10,replace=TRUE),Terrestrial.Mass=sample(80:100,10,replace=TRUE),Other.Mass=sample(40:60,10,replace=TRUE),Site.X.Treatment=rep(c("ALA1A","ALA1U"),each=5),stringsAsFactors=FALSE) library(reshape2)

Dear R-help, After making a typo (reminiscent of FORTRAN 77, I guess) I found the following: > identical(all.equal(2^(-10:10), 2**(-10:10)), TRUE) [1] TRUE I have tried to find the documentation about the ** operator but I was unsuccesful this way: > sessionInfo() Version 2.3.1 (2006-06-01) powerpc-apple-darwin8.6.0 attached base packages: [1] "methods" "stats"