[R] Changing mu in svyttest for one sample test

Hi All, I am doing a one-sample t-test for survey data with replicates (svrepdesign). I was wondering if there is a way to change the mu in svyttest for one sample test? I can't find a way to do this online for svyttest but I was able to do it for ttest. for example; in the following, I can

Re: [R] Survey Data - comparing multiple groups

. At this point, I am still exploring R and if I am able to provide additional feedback on this issue, I will post. I really appreciate your help. Thank you. From: Abby Spurdle Sent: July 27, 2019 9:39:49 PM To: Mavra Ahmed Cc: r-help@r-project.org Subject: Re: [R

Re: [R] Survey Data - comparing multiple groups

rom the survey package, the svyranktest is shown with svydesign only and not with svrepdesign. I could consider using other tests such as svyglm but am also having issues with it. Thank you, Mav From: Abby Spurdle Sent: July 25, 2019 9:48:05 PM To: Ma

[R] Survey Data - comparing multiple groups

Hi Everyone, I am trying to assess whether the mean of carbohydrate intake (dependent variable = continuous) is significantly different between five groups (groups = people categorized into five different groups) (independent variable= categorical). This is survey data with replicates and

Re: [R] Summary Table for Quantile Regression for Complex Survey

ace it here? Thanks for your help, Mavra ____ From: Koenker, Roger W Sent: July 16, 2019 11:37:54 AM To: Mavra Ahmed Cc: r-help@r-project.org Subject: Re: [R] Summary Table for Quantile Regression for Complex Survey try summary(newer, se = �nid�) Roger Koenker

[R] Summary Table for Quantile Regression for Complex Survey

Hello All, I am running quantile regression for a complex survey design and am running the following code: newdf<-withReplicates(df, quote(coef(rq(api00~api99, tau=0.5, weights=.weights This provides me with the theta and SE but I would really like the p-values and residual degrees of

[R] Graphic functions in the ade4TKGUI graphic user interface in R

, plabels.boxes.draw = TRUE, : non convenient selection for dfxy (can not be converted to dataframe) Any suggestions. Thanks Ahmed Attia, Ph.D. Agronomist & Crop Modeler __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see h

[R] Maximum number of iterations (maxit) does not work in hydroPSO-modFit-optimr in r

ist(numiter=3)) res= optimr(par=xopt,my.function,method="L-BFGS-B", lower=xmin,upper=xmax,control=list(maxit=3,trace=T)) #Only the hjkb form dfoptim that has argument maxfeval and it works. res=hjkb(par=xopt,my.function, lower=xmin,upper=xmax,control=list(max

Re: [R] r-data partitioning considering two variables (character and numeric)

Thanks Bert, worked nicely. Yes, genotypes with only one ID will be eliminated before partitioning the data. Best regards Ahmed Attia On Mon, Aug 27, 2018 at 8:09 PM, Bert Gunter wrote: > Just partition the unique stand_ID's and select on them using %in% , say: > > id

[R] r-data partitioning considering two variables (character and numeric)

did not produce the desired output, the data are partitioned within the stand_ID. For example, one row of stand_ID 7 goes to training and two rows of stand_ID 7 go to testing. How can I partition the data by Genotype and stand_ID together?. Ahmed Attia __ R-h

Re: [R] remove rows of a matrix by part of its row name

Thanks a lot. The code works great. Regards Ahmed ** Ahmed Serag Analytical Chemistry Department Faculty of Pharmacy Al-Azhar University Cairo Egypt From: Rui Barradas <ruipbarra...@sapo.pt> Sent: Tuesday, May 22, 2018 2

[R] remove rows of a matrix by part of its row name

73.1 Any ideas or suggestion please ? Regards ** Ahmed Serag Analytical Chemistry Department Faculty of Pharmacy Al-Azhar University Cairo Egypt [[alternative HTML version deleted]] __ R-help@r-project.org mailing

Re: [R] netCDF to GeoTIFF by layer in r

Hi Ben, Adding 'band' to raster provided what I want; importnetcdf <- raster(file.nc,band=4) writeRaster(importnetcdf,file.tiff,"GTiff",overwrite=FALSE) Thanks for your help. Ahmed Ahmed Attia, Ph.D. Agronomist & Soil Scientist On Wed, Mar 28, 2018 at 11:53 AM,

[R] Fwd: netCDF to GeoTIFF by layer in r

6800, 43200, 72576 (nrow, ncol, ncell) resolution : 0.0083, 0.0083 (x, y) extent : -180, 179.9998, -56, 83.5 (xmin, xmax, ymin, ymax) coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0 data source : C:\Users\Ahmed\Documents\file.tif names :

Re: [R] LSD-test

The model should be class aov or lm and my model class is aovlist. tried tidy from broom library but did not work. To make it class aov, I had to remove the error term; model <- aov(Rotationdata_R$`GY(Mg/ha)`~Rep+code*as.factor(Nitrogen),data=Rotationdata_R) Ahmed Attia, Ph.D. Agronomist &am

[R] LSD-test

ass "c("aovlist", "listof")" to a data.frame What is wrong with this. Ahmed Attia, Ph.D. Agronomist & Soil Scientist __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-hel

[R] Using nls.lm to fit a non-continuous dates range

odel a_start<-8 b_start<-2*log(2)/a_start m<-nls(df$Y~a*exp(-b*df$days),start=list(a=a_start,b=b_start)) Error in nls(df$Y ~ a * exp(-b * df$days), start = list(a = a_start, b = b_start)) : parameters without starting value in 'data': Y, days Could you please help Thank you Ahm

Re: [R] Convert date to continuous variable in R

ew=TRUE) plot(LAI~Date, data=LAI_simulation,xlab="",ylab="", xaxt="n",type="p",cex=0.2,col=1,xaxs="i",yaxs="i",ylim=c(0,6)) axis.Date(side=3, at=seq(min(LAI_simulation$Date), max(LAI_simulation$Date),by="year"), format="

[R] Convert date to continuous variable in R

LAI~Date, data=LAI_simulation) LAI_simulation$Date<- as.numeric(LAI_simulation$Date,origin="2009-10-07") #R does not see the origin axis.Date(side=3,at=LAI_simulation$Date) Ahmed Attia, Ph.D. Agronomist & Soil Scientist [[alternative

Re: [R] Predict

= 2016, frequency =365) plot(forecast(ets(Stand_Height),10)) a=seq(as.Date("2009-12-01"),by="weeks",length=11) axis(1, at = a, labels = format(a, "%Y %b %d"), cex.axis=0.6) #Error :$ operator is invalid for atomic vectors Thanks Ahmed Attia, Ph.

[R] Predict

<- predict(lm(height~Date, data=Stand_Height)); dl=leafbiom$Date[i-1]; de=leafbiom$Date[i]; SH[i]=sum(SH[leafbiom$Date==de$Date==dl]) } SH The SH output is the prediction of Stand height in 25 observations only and provides NA for the remaining 2616 iterations. Thanks for your help. Ah

[R] For loop

quot;i"])/2) } #Error in charToDate(x) : Thanks for your help Ahmed Attia, Ph.D. Agronomist & Soil Scientist __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the postin

Re: [R] Help/ Mathematics

uot;2012-02-12"] - leafbiom97$LeafBiog[leafbiom97$Date == "2010-03-15"]+ Litterfall_Ahmed97$littperiod[Litterfall_Ahmed97$Date =="2011-04-08"])/ (sum(GPP_Ahmed13$GPP[GPP_Ahmed13$Date >= "2010-03-12" & GPP_Ahmed13$Date <= "2012-04-12"])/2)

[R] Help/ Mathematics

format. Which library or function Should I use to tell R do these calculations of these variables at different dates. I appreciate your help. Ahmed Attia, Ph.D. Agronomist & Soil Scientist __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more,

[R] Variable 'A' is not a factor Error message

I am running a DOE with the following code library(Rcmdr) library(RcmdrMisc) library(RcmdrPlugin.DoE) # Define plackett burman experiment PB.DOE <- pb(nruns= 12 ,n12.taguchi= FALSE ,nfactors= 12 -1, ncenter= 0 , replications= 1 ,repeat.only= FALSE ,randomize= TRUE

[R] PSM With Panel Data

that how I can do this. Any help in this regard will be much appreciated With Best Regards,  Ahmed Arif  [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman

[R] R help - "UCL at different levels"

Hello Experts, I need your help to understand the how can we calculate the UCL values at different levels dynamically? Customer No Month Sales UCL 100 201501 500 550.75 100 201502 100 550.75 101 201501 400 425.25 101 201502 50 425.25 In the above scenario i need calculate

[R] R packages(maxLik ())‏

Hi all.. I am currently finish recent research to study master's degree in statistics. And in fact, I faced two problems in the practical side using r packages . In the first, generation of the new distributional data( weibull-lomax dist ) and I've successfully overcame them by the following

[R] pch size in a legend

Hi R users, I would like to adjust the pch size in a legend without changing the text size, pt.cex does not do the job. R 2.15.2 32 bit. legend(0,2100, legend=c(2009,2010,2012,2013,2014), col = 1,cex=1,lty=NA,pch=c(1,2,6,7,8),lwd=2,bty=n) Thanks Ahmed Attia, Ph.D. Agronomist Soil Scientist

Re: [R] Specifying strip names in xyplot (strip.custom)

Hi I attach the data and a xyplot I made using strip.default; # This is my code; trellis.device(col=F) xyplot(GrainM~Nrate|S*Loc, data = Ahmed, layout=c(5,3,1), xlab=expression(bold(paste(N rate (kg , ha^-1,, #xlab=expression(paste(Herbicide rate (g isoxaflutole , ha^-1, ))), #ylab

[R] Specifying strip names in xyplot (strip.custom)

. In the code below, factor.levels represent three Loc Thanks for your help. Ahmed trellis.device(col=F) xyplot(GrainM~Nrate|S*Loc, data = Ahmed, layout=c(5,3,1), xlab=expression(bold(paste(N rate (kg , ha^-1,, ylab=expression(bold(paste(Grain yield (Mg ,ha^-1,, ylim=c(10,20

[R] Prooblem using RMpi on OSX moutain lion

I have posted the following question on stackoverflow http://stackoverflow.com/questions/24940918/error-using-rmpi-on-osx-moutain-lion (but i didn't get any feedback). I am posting it here, hoping that a more R-focued user group might be able to help me out.. Thank you

[R] Fwd: Plot the means with simultaneous horizontal and vertical error bars

to draw the horizontal lines. Thank you so much. AA Ahmed M. Attia Research Assistant Dept. of SoilCrop Sciences Texas AM University ahmed.at...@ag.tamu.edu Cell phone: 001-979-248-5215 FAX: 001-308-455-4024 __ R-help@r-project.org mailing list

[R] Fwd: Plot the means with simultaneous horizontal and vertical error bars

would be the average response point (y Lint) with vertical and horizontal error bars. Each value of the x-axis (X) represents a range of values that in x Water. The x Water column should be used to draw the horizontal lines. Thank you so much. AA Ahmed M. Attia Research Assistant Dept

Re: [R] Fwd: Plot the means with simultaneous horizontal and vertical error bars

Ahmed M. Attia Research Assistant Dept. of SoilCrop Sciences Texas AM University ahmed.at...@ag.tamu.edu Cell phone: 001-979-248-5215 FAX: 001-308-455-4024 On Fri, Jul 18, 2014 at 3:30 PM, William Dunlap wdun...@tibco.com wrote: Is the original poster having trouble translating the synopsis

[R] Tick marks at the 0:0 point in a plot ()

Dear R usres, Sorry to bother you with my basic question. I have a quick question about tick mark in R. A reviewer wants the tick marks of the graph to start at 0:0 point. My current code below produces the tick mark a little far from the corner; plot(Simu2~Irrig, data=Ahmed,pch=17, col=1,subset

Re: [R] Tick marks at the 0:0 point in a plot ()

Dear William, Worded beautifully. Yes, I wanted to get the tick marks at the left corner. I had to just add xaxs=i and yaxs=i Thanks Ahmed Ahmed M. Attia Research Assistant Dept. of SoilCrop Sciences Texas AM University ahmed.at...@ag.tamu.edu Cell phone: 001-979-248-5215 FAX: 001-308-455

[R] Tick marks at the 0:0 point in a plot ()

Dear R usres, Sorry to bother you with my basic question. I have a quick question about tick mark in R. A reviewer wants the tick marks of the graph to start at 0:0 point. My current code below produces the tick mark a little far from the corner; plot(Simu2~Irrig, data=Ahmed,pch=17, col=1,subset

[R] xyplot y scale label help

Hi R usres, I would like to specify y label units in xyplot. I tried the yaxp, but it does not work. Following are my codes; trellis.device(col=F) xyplot(GY~N|S,data = Ahmed, layout=c(5,1,1), xlab=expression(paste(N rate (kg , ha^-1,))), ylab=expression(paste(Biomass yield (Mg ,ha^-1

Re: [R] regsubsets

..etc.) if Best subset regression do not give satisfactory results. I think it's interesting to   compare between all of them.   On Sunday, March 30, 2014 2:54 PM, Mohammad Ahmed mwh...@yahoo.com wrote: I am not able to find the criterion under which the best subset model for a certain size

[R] regsubsets

I am not able to find the criterion under which the best subset model for a certain size is chosen for the function regsubsets (leaps package in R) . For leaps function, we can determine the method ( cp,r2 or adj r2). For example, x - predictors y - response g - regsubsets(x,y,nvmax=1,nbest=1)

[R] Help: Split-Split model nested in a nested structure

and year 2 separately. We are using the following ANOVA model. We do not know whether we are missing something here. Help will be highly arreciated. Regards Zia Ahmed, CIMMYT model-aov(YIELD~ENV*SOWING*N_RATE*GEN+Error(FARMERS/N_RATE/GEN), data=mydata) summary(model) # District: two DIST

[R] sourcing scripts

: memory.size(max=TRUE) memory.limit() gc() memory.limit(16000) memory.size() gc() rm(list=ls(all=TRUE)) source(H:\\HIRAM\\Script_1.R) cheers Ahmed Kenawy __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] sourcing scripts

Hi Barry, I removed the lines related to memory size from other scripts and kept it only in the master script. Now it works. Cheers Ahmed Kenawy From: b.rowling...@gmail.com [b.rowling...@gmail.com] On Behalf Of Barry Rowlingson [b.rowling

[R] FW: sourcing scripts

Hi Barry, Thanks for your reply. Actually, I tried both max=TRUE and min=TRUE. Neither of them works. Ahmed Kenawy From: b.rowling...@gmail.com [b.rowling...@gmail.com] On Behalf Of Barry Rowlingson [b.rowling...@lancaster.ac.uk] Sent: Wednesday

[R] linear line in x, y plot

, data=Ahmed) NWI - seq(-0.08, -0.02, len = -0.02) lines(NWI, predict(reg1, list(NWI = NWI)),xlim=c(-0.08,-0.02),ylim=c(0,1),pch=1,col=1,lwd=2, lty=1) When I wrote the following code abline(reg1,pch=2,col=2,lwd=2, lty=2,xlim=c(-0.06,-0.02),ylim=c(0.3,0.8)) the line shows up, but I can not control

[R] Text in a plot without x axis

Dear R users, Can I add a text to a plot without x axis? I saw a function that enables placing the text with mouse text(location, text to place, pos, ...) locator(1) but it did not work out. -- Ahmed M. Attia Research Assistant Dept. Of SoilCrop Sciences Texas AM University ahmed.at

[R] Custom x axis in a boxplot

Dear R users, How can I customize x axis in a boxplot? xaxt=n to suppress x axis, but I would like to customize x axis for different groups. Thanks -- Ahmed M. Attia Research Assistant Dept. Of SoilCrop Sciences Texas AM University ahmed.at...@ag.tamu.edu Cell phone: 001-979-248-5215

[R] Nash-Sutcliffe coefficient

Hi Could you please inform how can I get the Nash-Sutcliffe coefficient in r? Comparing observed vs. simulated values. Appreciated -- Ahmed M. Attia Research Assistant Dept. Of SoilCrop Sciences Texas AM University ahmed.at...@ag.tamu.edu Cell phone: 001-979-248-5215

[R] Fwd: axis lab font in r

-- Forwarded message -- From: Ahmed Attia ahmedati...@gmail.com Date: Tue, Sep 17, 2013 at 11:17 AM Subject: axis lab font in r To: r-help-requ...@r-project.org I have an a question about the axis lab font in r. I use font.lab, but r changes the font of x axis lab only. How I

[R] replacing Na's with values on different records

I'm sure I'm missing something really obvious in the for loop... Here is simplified data for 3 patients, we need filling in Na's with same WT for each patient, AUC halved for visit 3, doubled for visit 5 for the same patient, based on visit 4 for(i in unique(u3s$ID)){

[R] Negative page rank values

Hello, I'm using igraph for some simple network analysis of a directed graph. When I calculate the page-rank using page.rank(g), the returned vector contains negative values. The vector sums to 1, which tells me that the algorithm working fine. How do I interpret the negative values and is

Re: [R] selection based on dates

I'm still a novice at R, so this may be a bit convoluted but it works: colnames(x) = c(date, id, value) do.call(rbind, (dlply(as.data.frame(x), .(id), function (y) y[-c(which(as.Date(y$date, %m/%d/%y) == min(as.Date(y$date, %m/%d/%y)), arr.ind=T)),]))) - FA -- Original Message --

Re: [R] Select only rows that don't contain one number

x[!apply(x, 1, function (y) any(y==-1)),] -- Original Message -- From: Dimitri Liakhovitski dimitri.liakhovit...@gmail.com To: r-help r-help@r-project.org Sent: 7/30/2013 10:06:02 AM Subject: [R] Select only rows that don't contain one number Hello! I have a data frame:

[R] simulating help

Dear all Could you please give me a link to a paper or article that explains how can I calculate type I error and power of a test using my data set by using R ? Also how can I do the same thing on a simulated data? I want this to compare between two tests I will be so grateful if you could help

[R] Importing data by odbcConnectExcel in 64 bit

, as.is = TRUE) Error in odbcConnectExcel(F:\\Cotton2012\\validation.xlsx) : odbcConnectExcel is only usable with 32-bit Windows -- Ahmed M. Attia Research Assistant Dept. Of SoilCrop Sciences Texas AM University ahmed ahmeda...@zu.edu.eg.at...@ag.tamu.edu Cell phone: 001-979-248-5215

Re: [R] kernlab::kkmeans initial centers

the number of clusters and lets the method assigns random cluster memberships. I am also wondering about the output centers. Are they the centers in the eigenvectors space? If so, what kind of initial centers I am supposed to pass to the method? --ahmed On Wed, Apr 3, 2013 at 12:54 AM, Pascal Oettli kri

Re: [R] kernlab::kkmeans initial centers

' not found --ahmed On Wed, Apr 3, 2013 at 10:30 AM, Pascal Oettli kri...@ymail.com wrote: Hi, The output of the provided example shows you which kind of matrix you need for cluster centers. In the example, 3 centers per variable (4) are defined. Thus, the initial cluster centers are defined

[R] kernlab::kkmeans initial centers

Hi, I am trying to pass initial cluster assignments to the kkmeans methodhttp://rss.acs.unt.edu/Rdoc/library/kernlab/html/kkmeans.htmlof kernlab. It is not clear to me how I can set the parameter *centers* with initial cluster centers as stated in the documentation? thanks, --ahmed

[R] Fwd: Help with nonlinear regression

-- Forwarded message -- From: r-help-ow...@r-project.org Date: Mon, Jan 14, 2013 at 12:31 AM Subject: Help with nonlinear regression To: ahmedati...@gmail.com The message's content type was not explicitly allowed -- Forwarded message -- From: Ahmed Attia

[R] Fwd: Help with nonlinear regression

-- Forwarded message -- From: Ahmed Attia ahmedati...@gmail.com Date: Mon, Jan 14, 2013 at 2:31 AM Subject: Help with nonlinear regression To: r-help@r-project.org Dears- I have a question about nonlinear regression. Asymptotic Regression Model below is the equation

[R] Getting the R squared value in asymptotic regression model

Please help getting the R squared value in asymptotic regression model I use the code below model1-nls(GN1~SSasymp (nrate,a,b,c), data = data.1 ) and R produced the modell coefficients without the R squared value? -- Ahmed M. Attia Research Assistant Dept. Of SoilCrop Sciences Texas AM

[R] Help

Help Please Hello, I want to find the whole hat matrix not only the hat values. Is there anyway that could be extracted from lm function ?. If not, please suggest something else. Direct computations using chol2inv or solve are not stable if t(X)%*%X has   high determinant. In this case lm is

Re: [R] Apply same linear model to subset of dataframe

, main=paste(DV, DV[i], sep==)) lapply(fit$model[, -1], function(x) plot(x, fit$resid)) fits[[i]] - fit } Many thanks Ross From: Jean V Adams jvad...@usgs.gov Date: Tuesday, 6 November 2012 19:20 To: Ross Ahmed rossah...@googlemail.com Cc: r-help@r-project.org Subject: Re: [R] Apply same

Re: [R] Apply same linear model to subset of dataframe

Thanks Jean This works for the plots, but it only stores the last lm() computed Ross From: Jean V Adams jvad...@usgs.gov Date: Tuesday, 6 November 2012 14:12 To: Ross Ahmed rossah...@googlemail.com Cc: r-help@r-project.org Subject: Re: [R] Apply same linear model to subset of dataframe

Re: [R] Apply same linear model to subset of dataframe

Superb Jean, many thanks Ross From: Jean V Adams jvad...@usgs.gov Date: Tuesday, 6 November 2012 19:20 To: Ross Ahmed rossah...@googlemail.com Cc: r-help@r-project.org Subject: Re: [R] Apply same linear model to subset of dataframe DV - c(mpg, drat, gear) IV - list(c(cyl, disp, hp), c(wt

[R] Apply same linear model to subset of dataframe

],lm.gear$residuals) Many thanks Ross Ahmed [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

Re: [R] ROC plot for KNN

Qian Liu littleduck24 at gmail.com writes: Hi I need some help with ploting the ROC for K-nearest neighbors. Since KNN is a non-parametric classification methods, the predicted value will be either 0 or 1. It will not be able to test for different cutoff to plot ROC. What is the package

[R] Problem indexing a factor variable

Hi, I have a dataframe gexp_1 with 115 rows of samples and 27000 columns of gene expression measurements with each column corresponding to one gene. I now have a smaller vector of genes tp for which I need to pull out the data form the data frame. # first gene from the list with 16 genes x

Re: [R] predictive accuracy

in the area of using R to solve challenging issues. Best Regards Ahmed -Original Message- From: Mike Marchywka [mailto:marchy...@hotmail.com] Sent: Thursday, May 26, 2011 8:55 PM To: gunter.ber...@gene.com; El-Tahtawy, Ahmed Cc: r-help@r-project.org Subject: RE: [R] predictive accuracy

[R] predictive accuracy

I am trying to develop a prognostic model using logistic regression. I built a full , approximate models with the use of penalization - design package. Also, I tried Chi-square criteria, step-down techniques. Used BS for model validation. The main purpose is to develop a predictive model

Re: [R] predictive accuracy

...@me.com] Sent: Thursday, May 26, 2011 10:54 AM To: El-Tahtawy, Ahmed Cc: r-help@r-project.org Subject: Re: [R] predictive accuracy On May 26, 2011, at 7:42 AM, El-Tahtawy, Ahmed wrote: I am trying to develop a prognostic model using logistic regression. I built a full , approximate models

[R] Fwd: Fwd: R help

On 11/21/10 11:24, Ahmed Attia wrote: -- Forwarded message -- From: Ahmed Attia ahmedati...@gmail.com Date: Sun, Nov 21, 2010 at 11:06 AM Subject: R help To: r-help@r-project.org Dear All, I'm a beginner user in R and I would like to make a quadratic plus plateau

[R] Help with linear plus plateau fit model

Dear R users, If you could please provide an example with the data set for the linear plus plateau model in R. I've 5 N rates and I need to make this model between N rates and grain Yield. Thanks so much Ahmed -- Ahmed M. Attia Assistant Lecturer El-Khattara farm Station Agronomy Dept

[R] R help

Dear All, I'm a beginner user in R and I would like to make a quadratic and plateau model in R. Can you help please with an example? Thanks so much -- Ahmed M. Attia Assistant Lecturer El-Khattara farm Station Agronomy Dept., Zgazig Univ., Egypt Visiting Scientist Haskell Agricultural

Re: [R] Ranked Set Sampling

it provided the author agree to that. Thanks again Ahmed On Sun, Oct 3, 2010 at 3:29 PM, David Winsemius dwinsem...@comcast.netwrote: On Oct 3, 2010, at 3:01 PM, David Winsemius wrote: On Oct 3, 2010, at 2:13 PM, Jeffrey Spies wrote: This is certainly not my area of expertise, but like Peter

[R] Ranked Set Sampling

Dear R Users; Are you aware of any package that calculates Ranked Set Sample? If you have a code that you are willing to share, I will acknowledge that in my work. Thanks much Ahmed [[alternative HTML version deleted]] __ R-help@r-project.org

[R] Ranked Set Sampling

Ahmed [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained

[R] stratified Wilcoxon's rank sum test with the locally best weights

ANOVA. I have been directed to use SAS with PROC FREQ with the SCORES=modridit option and CMH2 (SAS only has the locally best weight option??) I will never go back to SAS, and believe a code for the above test already exists?? Thanks in advance Ahmed

[R] Error in affypdnn package

parameters, (may take some time)... || |Error in object$call : $ operator not defined for this S4 class Please help me to get rid of this error(Error in object$call : $ operator not defined for this S4 class). Best regards, Ahmed Ryadh Hasan PhD Student, Knowledge-Based

[R] ecdf error

Dear friends, Ecdf used to work nicely before. e.g. ecdf(df1v[,c(1:6)],group=df1v$trt,q=.5,col=1:2,label.curves=list(keys=l ines), datadensity=rug) but now it does not work at all, even for a simpler code??? By the way, I used the same libraries. Thanks for your help. Ahmed

[R] duplicating records

could not find much help ( R code) to accomplish this duplication. Thanks for your help. Ahmed [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

[R] How do I calculate r^2 of SSasymp non-linear regression model in R?

Is it possible to get r^2 value of SSasymp non-linear regression model in R. Does anyone have any suggestions ? Thanks Zia __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] colorkey: Placement of labels and ticks in center positions of color key of a levelplot ?

Is it possible to place labels and ticks in center positions of each cut of color key of levelplot ? Here is my code that I used to produce a map (see attachment). The two labels and ticks were placed in left sides each segment of the color key. I want to show them center of the

[R] Resent: colorkey: Placement of labels and ticks in center positions of color key of a levelplot ?

Is it possible to place labels and ticks in center positions of each cut of color key of levelplot ? Here is my code that I used to produce a map (see attachment). The two labels and ticks were placed in left sides each segment of the color key. I want to show them center of the

[R] putting circles/buffers arround points on a 2D graph/plot

Hi, I'll try to keep my question brief and to the point, to anyone who helps 'THANK YOU!' I want to get a circle/ring/buffer or some other form of drawn line around certain points on a x,y plot, the points usually don't form a circle, often they form a 'wobbly' shape, so ideally I'd like

Re: [R] plot data with a colour scale - more details!

values ranging from 0.3-0.9). At this point, the discontinous shading would be a polishing step, and I would be just content with the color scale corresponding to the range of hydropahty values. Thanks. -Ray On Mon, Jan 19, 2009 at 8:55 AM, Ahmed, Rayhan raah...@utmb.edu wrote

Re: [R] plot data with a colour scale

Could you please provide further details - I am still still unable to produce this image. Thanks! -ray [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] GAM vs. MGCV packages

Hi all, I am a new R- user and I am going through the R-manuals, but I could not find an answer for my question. I am confused about when to use the GAM package and when to use the MGCV package?? My Model is a GAM model of continuous outcome and many non-linear continuous predictors (using