similar to: R-alpha: some thoughts about factor()

I have a data frame something like: name wrist nLevel emot 1 4094 3.34 1 frustrated 2 4094 3.94 1 frustrated 3 4094 NA 1 frustrated 4 4094 3.51

# I originally intended to post this question here. # But it took some time for subscription authorization process. # So I have posted the same question to the forum when I was waiting # for the subscription confirmation notice. # It seems to me the forum is not so active, let me drop here too. Hello, The following WinSNMP trap receive program ''win_snmp.rb'' causes ruby

Hi all, I use kriging to interpolate the precipitation from stations, but the map of this results show lots of stripes. (please see the attachment)I think there's something wrong with the setting of the dimension of this matrix, however, I have no idea how to know or test to see if this setting is correct or not.I've tried to switch the latitude and longitude, but still got the same

HI to you all! I have problem with quilt.plot....this is the code: quilt.plot(long_gridded,lat_gridded,d18O_gridded,nrow=n,ncol=m,xaxt='n',yaxt='n',xlab=NA, ylab=NA,breaks=seq(d18O_gridded_1,d18O_gridded_2,length.out=number_colors),col=col,horizontal=FALSE,nlevel=number_colors - 1,legend.args=list(quote(delta^18*O~("‰")),col="black",cex=0.8,side=3,line=1))

Hi , I am new to R . I am facing difficulty how to make pairwise comparisions. For example. I have a file which looks like below a b c d x 3 6 7 6 y 7 8 6 5 z 5 4 7 8 Here I need to look for the each pairwise comparisions (ab,ac,ad,bc,bd,cd for each row) For instance ,looking at first row, for x i need to look for ab values and take the min(3,6) >5 ,if its satistfies the count should be

Hello, Is it possible to draw level lines on the surface of a bivariate function? In the following example, to draw surface and levels lines for a multivariate normal law, I use persp, trans3d, contourLines and lines, but if the lines are correctly drawn, some parts of them are, of course, visible even if they are drawn on a non visible "face". Any suggestion to avoid this problem

System info: Linux slackware R Version 1.5.1 (2002-06-17) ESS 5.1.21 Emacs 21.2.1 ------------------- Colleagues Thanks very much for advice. I was not aware of the Fields package (although I've been using Geostat and Spatial). The image.plot() function in Fields does a lovely job, with an added bonus of fixing the problem I was having of balancing the number of colour levels with the

I am using following codes to plot map library(sp) library(rgdal) library(maps) library(gplots) library(clim.pact) library(fields) source("/R/PlotGridded2DMap.R") source("/R/image.plot.fix.R") source("/R/image.plot.plt.fix.r") seasonal_plot<-function(input,lonll=-180,latll=-90,lonres=5.,latres=3.75,write_file=TRUE,The_title=NULL){ if(is.null(The_title)){

Hola: Estoy intentando generar un mapa de calor que muestre, dadas unas coordenadas geográficas determinadas, el ?calor? (cuanto mas alto el valor de esta variable mas rojo y cuanto mas bajo mas azul) de una determinada variable. Es decir con unos datos como estos: x y z -3,716575 40,41743 3.169 -3,710735

Dear R-helpers, I'm working on mass spectra in randomForest/R, and following the recommendations for the case of noisy variables, I don't want to use the default mtry (sqrt of nvariables), but I'm not sure up to which proportion mtry/nvariables it makes sense to increase mtry without "overtuning" RF. Let me tell my example: I have 106 spectra belonging to 4 classes, the

Hi, I'm newbie in R package, but i'm extremely gratefull for these list! is there possible to make kruskal.test summarize the sum of ranks?? e.g ##### z<-c(1.2,3.4,0,0,0,0,0) ranks(z) # like in octave package ans= 6 7 3 3 3 3 3 sum(ranks(z)) ans= 28 ### then apply multiple comparisions for many other data SE<-sqrt(n*(n*k)*(n*k+1)/12) # according Zar,JH. Bioestatistical

Christopher Lamb wrote: > A strategy that we have had luck with in the past is the concept of a > "vector boolean" type, which is the result of a comparison between two > vectors. It's then necessary to perform a reduction (and, or) on the > elements in order to produce a scalar boolean value that's suitable for > a branch. This also allows stuff like letting

Dear list, I try to compare the mean of a variable given a value of a factor with the mean of the same variable for all K-1 other non-missing values of this factor. This procedure I want to repeat for each level of the factor. Having read the recommendations of this list I want to avoid creating K-1 dummy variables and searched for options of the pairwise.t.test. But couldn't find a

I just tried rc4-cvs and I can't believe how good it sounds! I can hear very little difference between the wav and the q0-ogg on my above-average stereo... It really makes MP3 look (sound) *pathetic* in comparision! I did some comparisions on a few songs, from noisy 77-punk to Jean-Michel Jarre to EBM/Industrial, and q0/~64kbit RC4-vorbis sounds as good or better than V9/~128kbit LAME-mp3

Dear R help forum, I am using the function 'coxph' to obtain hazard ratios for the comparison of a standard treatment to new treatments. This is easily obtained by fitting the relevant model and then calling exp(coef(fit1)) say. I now want to obtain the hazard ratio for the comparison of two non-standard treatments. >From a statistical point of view, this can be achieved by dividing

Hi there, I was hoping to post this message. Im hoping I;ve got the right spot?! Hi there, sorry for basic question but im very new to R. Im trying to run a lme model with two categorical variables, each having 6 (for the explanatory variable C.f) and 5 levels (for expl. variable D.f) respectively. The lme runs fine, summary and anova commands are ok. But now i've tried running posthoc

Before I got into my question for the day I'd like to applaud all the helpful folks and time spent behind the asterisk project to get it where it is today. Great work and between this list, the doc list and the irc channel it's been a pleasure to deal with people willing to help others when and if they can. I will be moving in a few months and I'm concerned as to what kind of

Hi, I have taken one microarray experiment and trying to implement same statistical measures what they have done.I have taken datasets from GEO platform with accession number GSE1557. In the experiment,about half of double transgenic rats (dTGR) over-expressing the human renin and angiotensinogen genes die by age 7 weeks of terminal heart failure (THF); the other (preterminal) half develops

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*`,

You can use the graphics parameter "srt" to rotate displayed text by a specified number of degrees, e.g. srt=45 to put it on an angle, srt=90 to put it vertical. If you do this, may need to modify the call to text to increase ylim and change the plot location to give you more room. I'm working to update the 'balloonplot' function in the gregmisc package now to handle this