[R] Residual Plots

HI All, I Am creating a residual plot for my linear model. the code I created is : plot(eval$bty_avg,residuals,ylab=residuals, xlab=Score, main = Residual Analysis)Here data set is eval. eval$bty_avg is my response variable and residual is the var I have created using resid function to store the

Re: [R] geom_text in ggplot (position)

This problem is discussed in Winston Chang's R Graphics Handbook, section 3.9. Adapting his code to this example: test - data.frame(variables = c(PE_35, PE_49), value1=c(13,3), value2=c(75,31), value3=c(7,17),

Re: [R] ifelse and vs

Hi nusrat, The ifelse function returns the number of values that result from the logical expression in the first argument. If you use , you get one logical value. If you use you get logical values for the number of conditionals that you specify. For example: 1:10 0 1:10 11:20 [1] TRUE 1:10

Re: [R] Help with pedigree() function in kinship2

Your problem is that PatientID, FatherID, MotherID are factors. The authors of kinship2 (myself and Jason) simply never thought of someone doing this. Yes, that is an oversight. We will correct it by adding some more checks and balances. For now, turn your id variables into character or

Re: [R] Residual Plots

Hi Shivi82, The error message suggests that eval$bty_avg is a function. What does: str(eval) say about its components? Jim On Wed, May 13, 2015 at 7:33 PM, Shivi82 shivibha...@ymail.com wrote: HI All, I Am creating a residual plot for my linear model. the code I created is :

Re: [R] Mean rain per rain day

Hi Frederic, It looks like you have a number of daily rainfall values and you want to get the mean rainfall on days where there was some rain. What you probably want is: daily_rainfall-rpois(50,2) mean(daily_rainfall[daily_rainfall0]) You may have to add na.rm=TRUE if there are NAs in your data.

[R] rminer package

Hello, I googled and looked on r-bloggers website for how to plot several REC curves (rminer package) on the same plotting area but couldn't find anything useful. I want to compare several regression (not logic regression and not classification) predictive models analyzing the REC curves and the

[R] Help regarding multi-page web application

Hello, I am using shiny package to develop web application but got stuck while creating multi-page web application. I googled my problem but did not get proper help. I just found R is not supporting multi-page web application but it was to old question. *Is R supports multi-page web application

Re: [R] Count

Hi, On Wed, May 13, 2015 at 5:58 AM, venkadesan venky venkyno...@gmail.com wrote: Hello Team, I have data like this have Small doubts on the following calculation For example, Employee sizeCamp 1camp 2 Camp 3 11 0 0

[R] lmList Error in !unlist(lapply(coefs, is.null))

I have a dataframe with the following structure: 'data.frame': 13095 obs. of 1433 variables: $ my : Factor w/ 624 levels 19631,19632,..: 1 1 1 1 1 1 1 1 1 1 ... $ s1 : num NA NA NA NA NA NA NA NA NA NA ... Where my is a factor with the number of the month, s1,..,,s1426 vectors that

[R] Count

Hello Team, I have data like this have Small doubts on the following calculation For example, Employee sizeCamp 1camp 2 Camp 3 11 0 0 20 0 1 3

[R] Cross correlation between two time series over nested time periods?

Dear R users, I have two time series Calculate and plot cross correlation between two time series over nested time periods. Each point in either time series is for a week (not exactly a calendar week, but the first week in a calendar year always starts from Jan 1, and the other weeks in the

Re: [R] Fw: Downloading R

Hmmm. Follow the instructions at http://www.r-project.org/ After choosing a mirror, you will reach a page that says, in part: partial quote Source Code for all Platforms Windows and Mac users most likely want to download the precompiled binaries listed in the upper box, not the

Re: [R] Residual Plots

This example might help: tmp - data.frame(x=1:10, y=rnorm(10)) foo - lm(y~x, data=tmp) plot(tmp$x, residuals(foo)) It appears that eval$bty_avg is not what you think it is. -Don -- Don MacQueen Lawrence Livermore National Laboratory 7000 East Ave., L-627 Livermore, CA 94550 925-423-1062

Re: [R] binding two lists of lists of dataframes together

Hi David, I tried both solutions you provided(lapply(dlist, function(x) rbind(x,x) ) Map( rbind, smaller, smaller)) and they seem to transpose and reshape the data into a different structure. Whenever I added the str - I only get a NULL. You basically want a list of the same general

Re: [R] binding two lists of lists of dataframes together

On May 13, 2015, at 2:01 PM, Vin Cheng wrote: Hi David, I tried both solutions you provided(lapply(dlist, function(x) rbind(x,x) ) Map( rbind, smaller, smaller)) and they seem to transpose and reshape the data into a different structure. Whenever I added the str - I only get a

Re: [R] Code works on Mac but not Windows

Yes bar - function(y = 12345){ + message(y) + } foo - function(x = a, y = b){ + bar(y = y) + } bar() 12345 bar(y = abc) abc foo() b foo(y = xyz) xyz ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek / Research Institute for Nature and Forest team Biometrie Kwaliteitszorg /

Re: [R] Plotting times at night and getting plot limits correct

Hi Bob, Given the other answers I may be off target, but perhaps this will help: # create two nights worth of data Times-strptime( paste(c(2015-05-13,2015-05-14),paste(rep(c(18:23,0:6),2),:30:00,sep=)), %Y-%m-%d %H:%M:%S) # telescope the two nights into repeated hours

[R] error using the huge R package

Dear List I am trying to do an association network using some expression data I have, the data is really huge: 300 samples and ~30,000 genes. I would like to apply a gaussian graphical model to my data using the huge R package. Here is the code I am using dim(data) #[1] 317 32200 huge.out

[R] Post-hoc test on split-plot design

Hi Richard, Thank you very much for your cooperation. I installed the packages HH and run the following model: Bud.aov-aov(terms(Budburst~CO2*SoilTemp*Photoperiod+Greenhouse/(SoilTemp*Photoperiod),keep.order=TRUE),data=data) summary(Bud.aov) Df

Re: [R] Post-hoc tests on Split-plot design

Hi Richard, Thank you very much for your cooperation. I installed the packages HH and run the following model: Bud.aov-aov(terms(Budburst~CO2*SoilTemp*Photoperiod+Greenhouse/(SoilTemp*Photoperiod),keep.order=TRUE),data=data) summary(Bud.aov) Df

[R] matrix inf and zero's value replacement

Dear Group, I have the following matrix m [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,]021 Inf Inf Inf Inf Inf [2,]102 Inf Inf Inf Inf Inf [3,]210 Inf Inf Inf Inf Inf [4,]3210 Inf Inf Inf Inf [5,] Inf Inf Inf Inf0

Re: [R] matrix inf and zero's value replacement

On 13.05.2015 08:04, Ragia Ibrahim wrote: Dear Group, I have the following matrix m [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,]021 Inf Inf Inf Inf Inf [2,]102 Inf Inf Inf Inf Inf [3,]210 Inf Inf Inf Inf Inf [4,]3210 Inf

[R] Mean rain per rain day

Hi All, I want to compute Mean rain per rain day from rainfall data but i don't know how to go about that. Anyone who understand the approach I can use can help me. In addition, I would like to have RScript which can help me to compute it. thanks. Regards, Frederic. Frederic Ntirenganya Maseno

Re: [R] Code works on Mac but not Windows

Dear Glenn, Please keep the mailing list in cc. In this case I would drop the signature like values and use correct defaults. TermStructure - function(rates.data, method = ns) you can drop if(missing(method)) method = ns in that case. Best regards, ir. Thierry Onkelinx Instituut voor natuur-