Re: [R] [External] Re: Forestplot, grid graphics Viewplot grid.arange

Thank you for the suggestion Jim. Will try this the next time I do forest plots. Mary Putt, PhD, ScD Professor of Biostatistics Department of Biostatistics, Epidemiology & Informatics Pereleman School of Medicine University of Pennsylvania 215-573-7020 From: Michael Dewey Date: Fr

Re: [R] [External] Re: Forestplot, grid graphics Viewplot grid.arange

turn out really well using the forestplot package. Mary Mary Putt, PhD, ScD Professor of Biostatistics Department of Biostatistics, Epidemiology & Informatics Pereleman School of Medicine University of Pennsylvania 215-573-7020 From: Jim Lemon Date: Thursday, October 13, 2022 at 11:34 P

[R] Forestplot, grid graphics Viewplot grid.arange

rt() pushViewport(viewport(layout.pos.row = 1, layout.pos.col = 3)) fp_denmark upViewport(2) ### #Attempt 3 converting to grobs and use patchwork ### library(ggplotify) library(patchwork) fpd_grob <- grid2grob(print(fp_denmark)) p1 <- grid2grob(print(fp_denmark)) p

[R] help with stacked ggplot

order of my data is graphed alphabetical (and not what I am looking for); Asymmetric on the left, Common in the middle, then Private on the right. Is there an easy way to reorder this? Thanks for any help! Best, Mary [[alternative HTML version deleted]] ___

[R] How to create 10 minute time series from hourly data

but haven't had luck. Any input is greatly appreciated. ~Mary Ann Middleton, PhD Here is a sample of the data I have: Date Time date.time ms LEVEL TEMPERATURE Level_m 1 2016-05-31 15:25:00 2016-05-31 15:25:00 0 92.1767 25.171 9.401814 2 2016-05-31 16:25:00 2016-05-31 16:25:00 0 92

[R] Help with glm and glht for analysing count data

e right lines… Any help would be fantastic! Many thanks, Mary [[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

[R] unable to call certain functions

he FAQs related to being sure that the package was actually loaded/installed. Has anyone else had this problem and, if so, what was the fix? Appreciatively, Mary Pfauth Mary Pfauth Alaska Department of Environmental Conservation Air Quality Division 619 E. Ship Creek Ave. #249 Anchorage, AK 9950

Re: [R] Why 'gbm' is not giving me error when I change the response from numeric to categorical?

, 2013 at 3:50 PM, Marc Schwartz wrote: > > On Oct 4, 2013, at 2:35 PM, peter dalgaard wrote: > > > > > On Oct 4, 2013, at 21:16 , Mary Kindall wrote: > > > >> Y[Y < mean(Y)] = 0 #My edit > >> Y[Y >= mean(Y)] = 1 #My edit > > > > I

Re: [R] Not getting any result from 'gbm'?

Sorry David, The formula that I use here is fmla = as.formula(Y ~ x1 + x2 + x3 + x4 + x5 + x6 + x7 + x8) Thanks On Sat, Oct 5, 2013 at 2:02 AM, David Winsemius wrote: > > On Oct 3, 2013, at 3:07 PM, Mary Kindall wrote: > > > In the reproducible example given below, why I am

[R] Why 'gbm' is not giving me error when I change the response from numeric to categorical?

n.trees : argument of length 0 My question is, Is binarizing the response will have so much effect that it does not find anythin useful in the predictors? Thanks -- - Mary Kindall Yorktown Heights, NY USA [[alternative HTML version deleted]] ___

[R] How to create an ROC with three possible classes?

,-1],type="prob", aggregation= "average"); pred For two class case, I use to do the following but it is no longer valid for three classes. yhat <- pred[,2] y = test[, -1] plot.roc(y, yhat) Any help will be appreciated. Thanks -

[R] Which regression tree algorithm to use for large data?

impossible to apply rpart on my data. I am using the following rpart command: varimp = rpart(fmla, dat=tmpData, method = "class")$variable.importance Thanks -- - Mary Kindall Yorktown Heights, NY USA [[alternative HTML versi

[R] "error setting certificate" issue in R-studio

Hi all, I'm a newbie to R and I can't figure out how to fix this error I keep getting in R-studio: > "error setting certificate verify locations:\n CAfile: \n CApath: none\n" Error in twInterfaceObj$doAPICall(cmd,params,"GET",...) : Error: error setting certificate verify locations: CAfile: CAp

Re: [R] Generate positive definite matrix with constraints

Thank you but I want to create a matrix (4,4) or (8,8) not just only a (2,2) matrix... -- View this message in context: http://r.789695.n4.nabble.com/Generate-positive-definite-matrix-with-constraints-tp4667449p4667509.html Sent from the R help mailing list archive at Nabble.com. ___

Re: [R] Generate positive definite matrix with constraints

Sorry in the previuos message I've inverted: I would like generate ... "same eigenvectors and different eigenvalues" -- View this message in context: http://r.789695.n4.nabble.com/Generate-positive-definite-matrix-with-constraints-tp4667449p4667450.html Sent from the R help mailing list ar

[R] Generate positive definite matrix with constraints

Hi, I have a question for my simulation problem: I would like to generate a positive (or semi def positive) covariance matrix, non singular, in wich the spectral decomposition returns me the same values for all dimensions but differs only in eigenvectors. Ex. sigma [,1] [,2] [1,] 5

Re: [R] How to menage decimals in order to obtain accurate estimates

I've checked my iteration loop with other dataset, i've not problems with integer numbers or other dataset with decimals but not all decimals dataset. I would like it works with all type of data. Thanks -- View this message in context: http://r.789695.n4.nabble.com/How-to-menage-decimals-in-or

Re: [R] How to menage decimals in order to obtain accurate estimates

I'm sorry, I have a data set contains numbers with more decimal places, the data are available in R "data (milk)". I made a loop to get a robust correlation matrix but in calculations, decimals cause me problems ( especially I do not know how R rounds numbers by default) so the diagonal should

[R] How to menage decimals in order to obtain accurate estimates

Hi, I'm trying to calculate some functions loop, unfortunately my data are decimal numbers that not allow me to get an accurate estimate. I do not want to use neither "round "nor "signif" but I would prefer to have all the information available, (the function loop is generic so you should adapt to

Re: [R] Using loop for a random vector (Montecarlo method)

Hi, I've understand how to do a permutation matrix but how I can utilize it? I would like to have a sample of 5 unit with rotation parameters! The same with 10 unit or 20 unit... Thanks -- View this message in context: http://r.789695.n4.nabble.com/Using-loop-for-a-random-vector-Montecarlo-m

[R] Using loop for a random vector (Montecarlo method)

Hi, I would like to replicate a sort of Montecarlo experiment: I have to generate a random variable N(0,1) with 100 observations but I have to contaminate this values at certain point in order to obtain different vectors with different samples: tab<-function(N,n,E,L){ for(i in 1:100){ X1<-rnorm(

Re: [R] Ellipse in PCA with parameters "a" and "b"defined.

thank you David, it was my first idea, but i don't know if it is right statistically speaking!!! -- View this message in context: http://r.789695.n4.nabble.com/Ellipse-in-PCA-with-parameters-a-and-b-defined-tp4656215p4656274.html Sent from the R help mailing list archive at Nabble.com. _

Re: [R] Ellipse in PCA with parameters "a" and "b"defined.

l set on my parameters but many ellipses. Thanks Mary -- View this message in context: http://r.789695.n4.nabble.com/Ellipse-in-PCA-with-parameters-a-and-b-defined-tp4656215p4656242.html Sent from the R help mailing list archive at Nabble.com. __

[R] Ellipse in PCA with parameters "a" and "b"defined.

er is what I have to choose as X and Y into the formula? if "a" and "b" are scalars or vectors should be x and y scalars or vectors too? Thanks Mary -- View this message in context: http://r.789695.n4.nabble.com/Ellipse-in-PCA-with-parameters-a-and-b-defined-tp4656215.html Sent

[R] import data (txt) separated by many blank spaces and lines with variable names that correspond to multiple observations.

Hello, I have read other posts and in part I tried to solve my problem. I have a txt file but the values ​​as well as being separated by many blank spaces and lines eg.var_a var_b var_c 1 2 4 5

Re: [R] Getting discrete colors on plot

adient of colors? Also how would I specify the colors that I want, such as the RGB equivalents of green, yellow, and red? My real data will have many more records. -Mary __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-h

Re: [R] Iterative loop using "repeat"

b<-b2 } return(b) } But my question now is what is the difference between the recursive process "while do" and suggested earlier with "if else"? Thank you for all precious tips suggested me Mary -- View this message in context: http://r.789695.n4.nabble.com/Iterative

Re: [R] Iterative loop using "repeat"

Thanks for the help:your reply has focused my problem in particular I'm trying to do a converge process; the robustm() function is: function (x,z) 1# {eigen<-eigen(x) 2# d<-madmatrix(z) ##I have created this function

Re: [R] Iterative loop using "repeat"

Thanks for the help, For Mario: you're absolutely right in fact the process never stopped!!! in particular I'm trying to do a converge process; the robustm() function is: function (x,z) 1# {eigen<-eigen(x) 2# d<-madmatrix(z)

Re: [R] Iterative loop using "repeat"

I'm sorry I read it again the answers and I realized that the condition "if(abs(b2[i,i]-b[i,i])<=0.001)" must be tested for the entire matrix and not only for the diagonal, I misspelled. thanks -- View this message in context: http://r.789695.n4.nabble.com/Iterative-loop-using-repeat-tp46546

[R] Iterative loop using "repeat"

ntified as output (in this case b2) and the remarket in the process creating b3, b4, etc. up to find the minimum difference between the two matrices and returns me the last.) Sorry if I am not so clear, thanks Mary -- View this message in context: http://r.789695.n4.nabble.com/Iterative-loop-us

Re: [R] PROBLEMS INSTALLING LIFEMETRICS

I have Statconn. Lifemetrics is a software that is written in R. I have a problem to install this add-ins in Excel. I don't know very good r language and so I don't know where is the error to install it. -- View this message in context: http://r.789695.n4.nabble.com/PROBLEMS-INSTALLING-LIFEMETR

[R] FASTA annot issue

them back and I need to use MEGA to keep my data uniform. Any help would be greatly appreciated Thank you Mary [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read

Re: [R] POSIXct to ts

ot; calculation from decompose() to generate an acf of the data after the seasonal patterns and trends are stripped ( ac f ( n a.omit(x.ts.decomp$random) ) , which is ultimately what I need, however, that isn't solid justification for choosing that calculation.  Any pointers appreciated.

[R] POSIXct to ts

  I am stuck trying to go from POSIXct to as.ts Any suggestions on how to tackle this would be greatly appreciated.   Sincerely, Mary Ann A sample of the data looks like this: 'data.frame': 26925 obs. of 2 variables: $ date : POSIXct, format: "2008-07-11

[R] fractal package

coercing argument of type 'double' to logical   I cannot proceed because of the warning messages appeared. Is my installation correct? Please HELP me.. Thank you so much. Sincerely yours, Mary Joy Regidor Student MSU-IIT, Iligan City Philippines [[alternative HTML version deleted]]

[R] pvalue calculate

I have a value a=300 observation (x) = sample(1:50) How to find a p-value from this. I need to show that "a" is different fom mean(x). Thanks -- - Mary Kindall Yorktown Heights, NY USA [[alternative HTML version deleted]] __

[R] output Shapiro-Wild results to a table

the following with errors. > write.table(x=F1.norm,file="I:/R_Work/F1/Shapiro.csv", sep=",") Error in > as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = > stringsAsFactors) : cannot coerce class '"htest"' into a data.fram

[R] which function/method to find agreement between two

Hi I have data in the following format: itemsperson1 person2 - - -- car honda,toyotahonda bikesuzuki suzuki pant Lee Levis, Lee shirt Van_housen Hollister house rented

Re: [R] ggplot- using geom_point and geom_line at the same time

0:07 AM, Hadley Wickham wrote: > On Mon, Jan 16, 2012 at 6:05 PM, Mary Kindall > wrote: > > Thanks for reply > > I wanted to have legend name with spaces. Right now I am using the > > following code but it produce two legends. I have to use Gimp to cut the > &

Re: [R] ggplot- using geom_point and geom_line at the same time

crete('my Custom Legend') g <- g + scale_shape_discrete("my Custom Legend") g - On Mon, Jan 16, 2012 at 6:55 PM, Felipe Carrillo wrote: > Mary: > Here's one way. > ## change the variable name to whatever title you want on your legend > data = melt(data, id=&qu

[R] ggplot- using geom_point and geom_line at the same time

lt;- g + geom_point() g <- g + scale_x_continuous(name='Number of inputs') g <- g + scale_y_continuous('Conversion time (sec.)') g -- - Mary Kindall Yorktown Heights, NY USA [[alternative HTML version deleted]] ___

[R] ggplot- using geom_point and geom_line at the same time

Hi I am plotting line chart using ggplot and want to use geom_line and geom_point simultaneously. I get the plot but now I have two legends. None of the legend is representing the true values. I need the legend with shape and color both. Thanks > con = textConnection("inputs var1var2va

Re: [R] relative frequency plot using ggplot or other function

/55078149a733dd1a0b42a57faf847036.png http://had.co.nz/ggplot2/graphics/90983232ced45a93d9fbbe40afffd69a.png Thanks On Thu, Jan 12, 2012 at 12:13 PM, Justin Haynes wrote: > On Thu 12 Jan 2012 09:02:27 AM PST, Mary Kindall wrote: > >> Hi >> I have a data frame in the following form. There are two grou

[R] relative frequency plot using ggplot or other function

28 2 200 0.0022677303 0.02948625 3 300 0.0061182673 0.01739936 4 400 0.0152237225 0.02569902 5 500 0.0300215262 0.03639880 6 600 0.0597610250 0.07717765 Thanks -- - Mary Kindall Yorktown Heights, NY USA [[alternative HTML versi

Re: [R] mode of frequency distribution table

100 300 > 8 5 4 2 1 1 1 1 1 1 > #extract any range from sorted table > > t[order(as.numeric(t),decreasing = TRUE)][1:3] > x > 5 6 4 > 8 5 4 > > I hope this helps. > > Andrija > > > On Sun, Jan 8, 2012 at 1:48 PM, Mary Kindall > wrot

[R] mode of frequency distribution table

. Thanks. -- - Mary Kindall Yorktown Heights, NY USA [[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

[R] (Edited) cbind alternate for data frames

,1:4] > y = table2[1:100,3:6] > > z = cbind(x,y) //hanging the machine > > write.table(z,'out.txt) > > > > -- > - > Mary Kindall > Yorktown Heights, NY > USA > > [[alternative HTML version deleted]]

[R] cbind alternate

:100,5] z = cbind(x,y) //hanging the machine write.table(z,'out.txt) -- - Mary Kindall Yorktown Heights, NY USA [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listi

Re: [R] aggregate function

either of these: > > aggregate(x$case1, x$param, FUN=paste, collapse=",") > aggregate(x$case1, x[["param"]], FUN=paste, collapse=",") > > ) > > require(data.table) >>> x <- data.table(x) >>> x[ >>> >> + ,

[R] aggregate function

should be param case1 1 a 2 b,c,d 3 e 4 f Thanks M -- - Mary Kindall Yorktown Heights, NY USA [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listi

[R] read.table : fill missing entry with "unAvailable" [edit]

read.table(DataFile, header=FALSE, fill_missing_entry = 'unAvailable') Or remove the the complete row in which the missing entry appear. thanks -- - Mary Kindall Yorktown Heights, NY USA -- - Mary Kindall Yorktown Heights, NY USA [[alternative

[R] read.table : fill missing entry with "unAvailable" [edit]

ader=FALSE, fill_missing_entry = 'unAvailable') 1 thanks -- - Mary Kindall Yorktown Heights, NY USA [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do re

[R] read.table : fill missing entry with "unAvailable"

ader=FALSE, fill_missing_entry = 'unAvailable') thanks -- - Mary Kindall Yorktown Heights, NY USA [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

[R] list.dir() function

Hi I have an organism directory that contains two folders galGal3 and hg19 and many other files. orgDir = '/home/mary/org' When I try to use list.dir() function, it gives me the same answer, no matter what is the value of full.names argument. > list.dirs(path = indexDir, full.nam

[R] download.file

File, 'wget' , quiet = TRUE) != 0) { file.remove(destFile) } } thanks -- - Mary Kindall Yorktown Heights, NY USA [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mai

[R] Question: Self selection bias and censoring in R

, R code etc. on how to update my simple Z test for proportions to accommodate this would be most welcome, or alternative approaches in R. Best regards, Mary __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] download files using ftp: avoid error

; Program execution continues, whether or not the download from url was > > successful. However, wget is, I believe, not available on Windows. > > > > Rgds, > > Rainer > > > > > > On Friday 16 September 2011 15:07:15 Mary Kindall wrote: > > > I am plannin

[R] download files using ftp: avoid error

2down[i], sep='') destfile = paste (dest, 'inDir', files2down[i], sep='/' ) download.file(url, destfile, quiet = FALSE) } It works fine as long as the file is present. When the file is not present, it exit from loop. Is there a way to continue looping if error occ

[R] installing packages systemwide

I installed some downloaded packages in R. I always do $sudo R CMD INSTALL By default it is storing these packages into my directory /home/mary/R/x86_64-pc-linux-gnu-library/2.13/. However I want them to be systemwide into /usr/local/lib/R/site-library/ folder. I tried $sudo R R

Re: [R] save and load in R

the output files } newnames <- paste("file", 1:numFiles, sep="") #output variables to load into for ( i in 1:numFiles) { load(file = outfilenames[i]); assign(newnames[i], dataFile) #assign into corresponding output variables; } # Regards - M On Sun,

[R] save and load in R

ch a way that the data is loaded into a variable that is same as the file name without extension. file1 = load (file = './file1.Rdata') file2 = load (file = './file2.Rdata') file3 = load (file = './file3.Rdata') file4 = load (file = './file4.Rdata') How can I

[R] [Resolved] combine the data frames into comma separated list.

7;,'z','z'), Target4 = c('xx','yy','zz','tt','uu'))); L <- list(dataframe1, dataframe2, dataframe3, dataframe4) merge.all <- function(...) merge(..., all = TRUE) Reduce(merge.all, lapply(L, function(x) aggregate(x[2], x[1], toStr

Re: [R] combine the data frames into comma separated list.

Superb Gabor, Though I dont know what is happening, but yes it is workin and without fail. Thanks - M On Mon, Jun 13, 2011 at 8:20 PM, Gabor Grothendieck wrote: > On Mon, Jun 13, 2011 at 5:17 PM, Mary Kindall > wrote: > > Hi R users, > > I am new to R and am trying to merge

Re: [R] combine the data frames into comma separated list.

=TRUE, incomparables=''), by = 'Src', all=TRUE, incomparables='') ddply(test, ~ Src , colwise(paste, .(Target1)), collapse ="," ); Thanks On Mon, Jun 13, 2011 at 7:14 PM, Dr. D. P. Kreil (Boku) < david.kr...@boku.ac.at> wrote: > ?write.csv &

Re: [R] combine the data frames into comma separated list.

x ccc 43 yy ddd 55 zz 66 tt 76 uu 82 aaa 94 10 4 11 4 Thanks -- M On Mon, Jun 13, 2011 at 6:35 PM, Dr. D. P. Kreil (Boku) < david.kr...@bok

[R] combine the data frames into comma separated list.

- 5 - -zz 6 - - tt,uu Basically I am trying to make a consolidated table. Help appreciated. Thanks M - Mary Kindall Yorktown Heights USA [[alternative HTML version

[R] combine the data frames into comma separated list.

- 5 - -zz 6 - - tt,uu Basically I am trying to make a consolidated table. Help appreciated. Thanks M - Mary Kindall Yorktown Heights USA [[alternative HTML version

Re: [R] Tests for the need of cluster analysis

to determine whether cluster analysis is valid for my data, I will learn R and use it on my data. My data are panel data consisting of 49 states and 25 years. Currently, I am estimating models with fixed state and time effects. Thanks for any help you can give me. Cheers, Mary On Mon, May 2

[R] Tests for the need of cluster analysis

Does R have the capability to perform tests for the need of clustering analysis (e.g., in prabclus)? I am using panel data with two-way fixed effects but am unsure about whether I should be using a cluster option as well to estimate my model.-- View this message in context: http://r.789695.n4.nab

[R] Strategy for maintaining R in student PC labs

Dear All At institutes where there are PC labs used by Masters and undergraduate students and software is maintained /managed by a central IT services, do you always update R every time a new version is released - around beginning of April and October? Do you also ask them to install all patch rel

[R] How to use escape characters in a string

Hi, I want to use a character as below in R, as.character("`X^`R\`S") [1] "`X^`R`S" Warning messages: 1: '\`' is an unrecognized escape in a character string 2: unrecognized escape removed from "`X^`R\`S" But I found errors. Can I use some option in as.c

[R] Fwd: how to write a long string to many lines in a fasta text file

Hi, I want to write a long string to many lines in R, but I do not know how to do it. Hope to get your help. Thanks. Mary [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[R] error using frailty term

Trying to fit a model using the Surv() with a frailty term but get the following cryptic error message: -- Error in frailty.brent(sqrt(x), y, lower = 0) : Ties for max(y), I surrender Calls: runplots ... coxpenal.fit -> e

[R] randomForest question--problem with ntree

identical. Have changed mtry from mtry=5 to mtry=6 successfully. Have seen same problem on both a Windows machine and our linux system running 2.8 and 2.9. Sample code follws. Thanks in advance for help. Mary > m1<-as.formula(paste("as.factor(EAD)~", paste(names(cli

[R] saving only part of a computation in integrate

Hello, I am using Beta <-integrate(dnorm,mean=0,sd=1,-Inf,2.3552) > Beta 0.9907436 with absolute error < 4.2e-06 1-Beta results in an error How can I store into Beta only .9907436? Thank you. Sincerely, Mary A. Marion __ R-help@r-pr

[R] formatting in r

etc. were added by hand. I need an example of how I can mix alphanumeric and numeric data on output. Can you assist? R is proving to be a really fine computational language that is easy to use. Thank you. Sincerely, Mary A. Marion u xbar alpha zcrit zcrit2 zstatpvalue 140

[R] graphs

t;black") curve(dnorm(x,mean=150, sd=15),from=75, to=225, col='orange', add=TRUE) Thank you. Sincerely, Mary A. Marion __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

[R] mathematical notation in R

Hello, I have an R function which includes the statement LBAuo<- -. Rather than printing out - I want it to print out - infinity. As of yet I have not been able to do that. I am not in a graphics window just outputting a set of variables. Can you help? Thanks. Sincerely, Mar

[R] IF STATEMENTS

Hello, I am working on using if statements. What is the error message telling me here and how do I correct for it? I have tried various combinations of quotes. Thank you. Sincerely, Mary A. Marion #Find critical values crit<-function(n,alpha,type) { if (type==twoSided) { alpha2=alpha/2

[R] math symbols in R

, Mary A. Marion __ 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, reproducible code.

[R] naming of columns in R dataframe consisting of mixed data (alphanumeric and numeric)

void the following error: Error in names(tresults2) <- c("Statistic ", "Value") : 'names' attribute [2] must be the same length as the vector [1] Am I to use c( , ) or list( , ) with a dataframe? Thank you for your help. Sincerel

[R] "taking the log then later exponentiate the result" query

Hi, I am trying to figure out the observed acceptance rate and M, using generalised rejection sampling to generate a sample from the posterior distribution for p. I have been told my code doesn't work because I need to "take the log of the expression for M, evaluate it and then exponentia

[R] Generalised Rejection Sampling

Hi, I am trying to figure out the observed acceptance rate and M, using generalised rejection sampling to generate a sample from the posterior distribution for p. I have been told my code doesn't work because I need to "take the log of the expression for M, evaluate it and then exponent

[R] (no subject)

along the x axis? Any help would be much appreciated! Mary _ [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listi

Re: [R] Frequency table to histogram

, obviously it's far greater than that.. Louisa > Date: Tue, 7 Apr 2009 10:00:55 -0400 > From: villa...@dms.umontreal.ca > To: statsstud...@hotmail.com > CC: r-help@r-project.org > Subject: Re: [R] Frequency table to histogram > > Hi Mary, > > Could you please provi

[R] Frequency table to histogram

temp1<-rep(temp$Chest,100) hist(temp1) But this creates a v. odd histogram so I know it must be wrong! Would be very grateful if any one could help me out! Mary _ [[alternative HT

[R] print of objects in R

Hello, I have been watching my output as I create functions and do other things in r. One thing I don't like is the [1,] type notation at the beginning of a line. I have been able to change that to a number such as 1 2 etc. using as.data.frame(object). How can I stop the printing of a lin

[R] R code -column headings

t;Midwest", "South", "West" ) colnames(x)=c("Public", "Private") as.data.frame(as.table(x)) 2] add a column heading for region when I use addmargins(x) addmargins(x) #margin.table(x,1) #margin.table(x,2) Thank you. Sincerely, Mary A. Marion

Re: [R] Residual resampling for non linear reg model

es, : singular gradientafter > R has only produced a few bootstraps.Could anyone suggest where I am going > wrong? Would greatly appreciate it!Mary > manpower<-read.table("http://www.creem.st-and.ac.uk/len/classes/mt3607/data/manhours_surgical.dat";, > header=TRUE)#attach

[R] Residual resampling for non linear reg model

gradientafter R has only produced a few bootstraps.Could anyone suggest where I am going wrong? Would greatly appreciate it!Mary manpower<-read.table("http://www.creem.st-and.ac.uk/len/classes/mt3607/data/manhours_surgical.dat";, header=TRUE)#attach dataattach(manpower)B<-99

Re: [R] [PS] rmultinomial() function

[EMAIL PROTECTED]> Date: Wednesday, March 19, 2008 7:36 am Subject: RE: [PS] [R] rmultinomial() function To: Mary Black <[EMAIL PROTECTED]>, r-help@r-project.org > Mary -- > > The dmultinomial function (try ?Multinomial, noting that it is an > uppercase M) has a "log" o

[R] rmultinomial() function

scribed what the difference between those generators is. Any clarification would be greatly appreciated! Thanks, Mary __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-p

[R] R version of SAS Proc Varclus

I am interested in finding an R version of SAS "Proc Varclus". SAS's Proc Varclus implements an oblique cluster analysis based on principal components. How can I find out if R has a package that runs the same algorithm implemented in SAS "Proc Varclus"? Tha

[R] R version of SAS Proc Varclus

I am interested in finding an R version of SAS "Proc Varclus". SAS's Proc Varclus implements an oblique cluster analysis based on principal components. How can I find out if R has a package that runs the same algorithm implemented in SAS "Proc Varclus"? Thank you,

[R] extract year or month from date

Hi, I am having trouble extracting just the year or the month or the day from a date such as 5/7/2007 which is May 7th 2007. Is there any particular function to extract just the year from this format? When I am reading this data from a text file it is reading it correctly in the same format but d

[R] read sas data into R

Hi I am having trouble using read.ssd. Can someone help? The code that I have written is *sashome<-"C:/Mary/Datasets"* *read.ssd(read.ssd(file.path(sashome, "core", "sashelp"), "surv_1v",* *sascmd = file.path(sashome, "sas.exe"))* Here the