similar to: Help with normal distribution in random samples...

Hi all, I have data from an experiment with 3 independent variables, 2 are within and 1 is between. In addition to the dependent variable, I have a covariate that is a single measure per subject. Below I provide an example generated data set and my approach to implementing the ANCOVA. However the output confuses me; why does the covariate only appear in the first strata? Presumably it should

Hi All, I've been having a little trouble with creating a loop that will run a a series of t.tests for inspection, Below is the code i've tried, and some checks i've looked at. I've used the get(paste()) idea as i was told previously that the use of the eval should try and be avoided. I've run a single syntax to check that my systax is correct and works without any problems

On my machine (Mac OS 10.5.6, R 2.8.1) the following plot is drawn such that a substantial portion of the descender of the title is covered by the plot: library(ggplot2) qplot(x=1:10,y=1:10,main='p q j g') -- Mike Lawrence Graduate Student Department of Psychology Dalhousie University www.thatmike.com Looking to arrange a meeting? Do so at: http://www.timetomeet.info/with/mike/ ~

Hi, when performing a spearman_test stratified by a given factor in package "coin", how is it possible to obtain the value of Rho, the Spearman correlation coefficient ? Thanks in advance Gilles (F) [[alternative HTML version deleted]]

Hello, I have a very simple data set i imported from excel including 96 averages in a column along with 96 standard errors associated with those averages (calculated in excel). I plotted the 95 averages using r and I am wondering if it is possible to plot the second column of standard errors while applying error bars to each value so they represent the error corresponding to each average? thanks,

To understand some functions about spatial point patterns in "spatstat" ,I should know some background about it, and the best way is to read the monograph, and "Statistical Analysis of Spatial Point Patterns" (2nd edt.) is a better choise. But I can not find it anywhere I can. Who can help me? Thank you! ----- unangu at gmail.com -- View this message in context:

I'm trying to create a colored map that would show the number of students per state. My data frame consists of two columns - state and count. I'm using the following code library(maps) map("usa") library(plotrix) state.col<-color.scale(gre$count,0,0,c(0,1)) map("state",fill=TRUE,col=state.col) I'm getting a map, but the values are not being mapped to correct

Hopefully this is a pretty simple question:   Is there a function in R that calculates the "mode" of a sample?   That is, I would like to be able to determine the value that occurs the most frequently in a data set.   I tried the default R "mode" function, but it appears to provide a storage type or something else.    I tried the RSeek and some R documentation that I

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

Hello, here I have some code, which behaves quite strange, IMHO. There is a "res.txt" which will be collected before printing it. There is a paste-statement, which has a comment at the end of the line, which is this one: "# !!!!!HERE!!!!!". If you throw out the first hash-mark "#" on that line, the printout behaves like if it is done for each of the for-loop

Dear R users, I have a matrix "final" which looks like this: final oSO4 oNO3 mSO4 mNO3 [1,] 3.3728 0.2110 1.9517421 1.01883602 [2,] 0.8249 0.0697 1.5970292 0.11368781 [3,] 0.2636 0.1004 0.6012445 0.24356332 [4,] 8.0072 0.3443 6.1016998 3.63207149 [5,] 13.5079 0.6593 12.4011068 1.55323386 [6,] 6.1293 0.1989 5.7620926 0.12884845 [7,] 0.6004 0.0661

Hi guRus, I'm looking for advice on a good way to approach analysis of some multi-level data I've obtained. I had humans classify words and measured response time. Words were 10 positive words ("happy", "joy", etc) and 10 negative words ("sad","grumpy", etc). Words were also presented in either white, red or green color. All variables were

Hi all, I've noticed that many computational neuroscience research groups use MATLAB. While it's possible that MATLAB may have some features unavailable in R, I suspect that this may instead simply be a case of costly tradition, where researchers were taught MATLAB as students and pay for it as researchers because it's all they know. I'd like to attempt to break the cycle by

Dear List, One persistent feedback I am getting to people who are newly introduced to R ( especially in this cost cutting recession) is - 1) The website looks a bit old. While the current website does have a lot of hard work behind it, should n't a world class statistics package have a better website instead. You can check out www.knime.org which is an open source software , and free, and

Hello everyone I'm working with R 2.8.1 on a windows machine I have a question regarding time series analysis The first question is how does R expect the input file to be structured? I'm working with a *.txt file similar to the abbreviated one here: Date,stage 4/2/1953,7.56 4/3/1953,7.56 4/4/1953,7.54 4/5/1953,7.53 4/6/1953,7.5 4/7/1953,7.47 4/8/1953,7.44 4/9/1953,7.41 4/10/1953,7.37

the below... class User < ActiveRecord::Base has_and_belongs_to_many :teams class Team < ActiveRecord::Base has_and_belongs_to_many :users has_one :captain, :class_name => ''User'' produces the error... Mysql::Error: Unknown column ''users.team_id'' in ''where clause'': SELECT * FROM users WHERE (users.team_id = 1) LIMIT 1 i

This question might seem silly, because I felt that it MUST be in the mailing list archives or help files somewhere, but I simply couldn't find it. I want to make some simple CDF (cumulative distribution function) plots to check whether distributions are Gaussian / normal. But in order to check how "normal" the distribution is, I really need the y-axis to be Gaussian as well

Hi all, Where f(x) is a logistic function, I have data that follow: g(x) = f(x)*.5 + .5 How would you suggest I modify the standard glm(..., family='binomial') function to fit this? Here's an example of a clearly ill-advised attempt to simply use the standard glm(..., family='binomial') approach: ######## # First generate some data ######## #define the scale and location of

Hi, poma wrote: > - "unsigned char c;" does not solve the problem > > - "c >= 0 && c <= ' '" solves the problem for the current git What the ... ?!? My proposal of "(c >= 0 ... )" was only made to avoid the need for a type change of variable "c". Both proposals should be equivalent ... on the first morning glimpse they

Hi All - "i'm new to linux" My Specs: Suse 9.0 & Windows 2003 Server (PDC) My needs: authenticating with Windows 2003 Server My Problem: Samba 2.x.x not working with Windows 2003 Server (worked with Windows 2000 Server!?) My Solution: Try Samba 3.0.1... Unfortunately - I could noy locate a RPM for Samba 3.0.1 so far... Have any one got any idea where I would be able to