similar to: How to draw a graph after model selection?

Hello r helpers! Below is the whole coding for my programme. Before proceed more further, let me explain for you. First of all, I need to compute trimmed mean. Till that step is ok. Then I need to compute ssdw which is sum of square deviation. If I do equal trimming at both tail of distribution that I chose, I will use the first ssd formulae which is "a". But if I am doing unequal

I wrote with C and created a file that I can read into S+, but R chokes because I've got multiple spaces between some variables and there is a space at the start of each data line. This happens becuase C was right-justifying the numbers as it printed them out. It looks like this: T Run Seed IntAct Change f0 f0T0 f0T1 f0T2 average0 variance0 relentropy0 f1 f1T0 f1T1 f1T2 average1 variance1

hello... I am trying to get this code to work, but as I get to the predict command, it displays an error due to the length of the data sets from the removal of the NA's. Here is the data, and the code that I am using so far, if you run it, you'll see the error pop up.... please help me to get around this problem. Thanks in advance. Location Dist. Size low1 .5 10.5 low2 .5 23 low3 .5 NA

I am a novice user of R and am stumbling on how to order a dataset produced during my session. I have a 1863 row X 14 column dataset that I want to put out to a file. I want the output sorted by the first column and then by the second column both in ascending order. The first column is character and the second is numeric (I hope). I used an "as.numeric" function to assign that

Hi R user,I am struggling to figure out on how I can calculate the average from the two tables in R. Any one can help me? really your help would be grateful- I am spending so much time to figure it out. It should not be so hard, I think. I have very big data but I have created a hypothetical data for simplification. for example I have : table 1 table 1: species occurance data

Hi everybody hacking in the 200-pre4 i found everything i needed, so i gave it a try ... and it works! But dont blame for the code, its just quick and dirty! --- ../syslinux-2.00-pre4/pxelinux.asm Sat Jun 15 07:25:51 2002 +++ ./pxelinux.asm Wed Aug 14 10:05:00 2002 @@ -17,6 +17,11 @@ ; version; incorporated herein by reference. ; ;

> Dear expeRts, > > Currently, my colors are as follows: > mycol <- > c("blue1","blue2","blue3","blue4","black","yellow4","yellow3","yellow2","y > ellow1") > heatmap(snp, Rowv=NA, Colv=NA, col=mycol) > > However, I would like to have the following colors: > bright blue ->

Hi, #working directory data1 #changed name data to data1.? Added some files in each of sub directories a1, a2, etc. ?indx1<- indx[indx!=""] lapply(indx1,function(x) list.files(x)) #[[1]] #[1] "a1.txt"??????? "mmmmm11kk.txt" #[[2]] #[1] "a2.txt"??????? "mmmmm11kk.txt" #[[3]] #[1] "a3.txt"??????? "mmmmm11kk.txt" #[[4]] #[1]

On Mon, 2007-05-14 at 23:41 +0200, vax9000 at gmail.com wrote: > Full_Name: vax, 9000 > Version: 2.4.0, 2.2.1 > OS: 2.4.0: Mac OS X; 2.2.1: Linux > Submission from: (NULL) (192.35.79.70) > > > To reproduce this bug, first go to the website "http://llmpp.nih.gov/DLBCL/" and > download the 14.8M data set "Web Figure 1 Data file". The direct link is >

I am trying to use t.test on the following data: date type INTERVAL nCASES MTF SDF MTO SDO nFST MF nOBS MO MB BIASCV BIASEV ME MAE RMSE CRCF 2001-06-15 avn GE1.00 4385 0.246 0.300 1.502 0.556 1367 1.373 4385 1.502 1.471 0.285 0.164 -1.256 1.266 1.399 0.056 2001-06-15 avn

Hallo I am teaching a very elementary intro course about R, and want to show the students how to find the distribution of the maximum of rolling two balanced die. Is there perhaps a more elegant way to do this, other than the way I am using below? (I would like to show them two ways - the one shown here, and perhaps and improved/elegant approach (showing off the "beauty" of R.

Dear R-people ... I'm a new user. I can't get predict.lm() to produce predictions for new independent data. There are some messages in archived help about this problem, but I still don't see my error after reviewing those. I understand that the new independent data must have the same name(s) as used when the model was made. In the example below, predict.lm produces the

How do I post text-string dates along the x-axis instead of the excel_date position values? I would like to plot a time-series of chemical values changing with time The first column is the internal date from excel. This should be the x-axis position of a graph of the plotted data. The second column is a text-string of date. This should be what is posted (instead of the excel_date). The third

Hi R users, I need some help in the followings: I'm doing factor analysis and I need to extract the loading values and the Proportion Var and Cumulative Var values one by one. Here is what I am doing: > fact <- factanal(na.omit(gnome_freq_r2),factors=5); > fact$loadings Loadings: Factor1 Factor2 Factor3 Factor4 Factor5 b1freqr2 0.246 0.486 0.145

Dear R help I am having difficulty doing a biplot of the first two factors of a factor analysis. I presume it is because the values in factor 2 for Milk and NUTS are not displayed in the component loadings. Loadings: Factor1 Factor2 RedMeat 0.561 -0.112 WhiteMeat 0.593 -0.432 Eggs 0.839 -0.195 Milk 0.679 Fish 0.300 0.951 Cereals -0.902 -0.267

Hi, I've come across a strange error when using the lrm.fit function and the subsequent predict function. The model is created very quickly and can be verified by printing it on the console. Everything looks good. (In fact, the performance measures are rather nice.) Then, I want to use the model to predict some values. I get the following error: "fit was not created by a Design

Dear R-helpers, I'm stuck with a little problem that surely has an easy solution but I can't think of a way to solve it. I'd really appreciate any help you can offer me! I'll provide a small example. Given a dataframe data.txt that looks like this: ID freq Var Var_mean Ratio_mean Var_median Ratio_median Var_sum Ratio_min Var_max Ratio_max Var_min

Dear all, The 'make check' step fails for the pacakge mva on IBM AIX. The tail of the Rout log file looks like: > for(factors in 2:4) print(update(Harman23.FA, factors = factors)) Call: factanal(factors = factors, covmat = Harman23.cor) Uniquenesses: height arm.span forearm lower.leg weight 0.170 0.107 0.166

Hi, In my analysis of impacts of insecticide-treated bednets on malaria, I look at the relationship between malaria incidence and mosquito behaviors. The condensed data set is copied here. Ordinary regression (lm) shows that Incidence was negatively related to Mortality. This makes sense because the latter reflected the strength of killing mosquitoes by insecticide-treated nets. Since the

Suharto, If you're interested in performance with subscripting, you might want to look at pqR (pqR-project.org). It has some substantial performance improvements for subscripting over R Core versions. This is especially true for the current development version of pqR (probably leading to a new release in about a month). You can look at a somewhat-stable snapshot of recent pqR development