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:3

[R] Forestplot, grid graphics Viewplot grid.arange

rt(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)) p2 <- grid2grob(prin

[R] help with stacked ggplot

the 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.1498

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

specified seems to only work with one variable? Apologies if my explanation is poor, I would be more than happy to give you more information if it would help. I’m not sure if anything I’ve done if right or even if I’m on the right lines… Any help would be fantastic! Many thanks, Mary

[R] unable to call certain functions

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 99501 (907) 269-6879 mary.pfa

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 dwinsem...@comcast.netwrote: On Oct 3, 2013, at 3:07 PM, Mary Kindall wrote: In the reproducible example given below, why I

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 marc_schwa...@me.com wrote: On Oct 4, 2013, at 2:35 PM, peter dalgaard pda...@gmail.com 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 have no clue about gbm, but I don't think

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

in the predictors? 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

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

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 -- - Mary Kindall Yorktown Heights, NY USA [[alternative HTML version deleted]] __ R-help@r

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

? Or it is 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 version deleted

[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: CApath:

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.

[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,]

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

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

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:

[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:

[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){

Re: [R] Ellipse in PCA with parameters a and bdefined.

Ok... so, in my model my a is built using the standard deviation of the first principal component and b with the second, so my x and Y should be : PCA $ scores [, 1], PCA $ scores [, 2] but in this way I do not get out a confidence interval set on my parameters but many ellipses. Thanks Mary

Re: [R] Ellipse in PCA with parameters a and bdefined.

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.

[R] Ellipse in PCA with parameters a and bdefined.

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 from the R help mailing list archive

[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

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-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide

Re: [R] Iterative loop using repeat

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-loop-using-repeat-tp4654625p4654777.html Sent from the R help mailing

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:

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

[R] Iterative loop using repeat

) 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-using-repeat-tp4654625.html Sent from the R

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:

[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

[R] POSIXct to ts

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 21:00:00 2008-07-11 22:00:00

Re: [R] POSIXct to ts

() 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. Cheers, Mary Ann - Original

[R] fractal package

correct? Please HELP me.. Thank you so much. Sincerely yours, Mary Joy Regidor Student MSU-IIT, Iligan City Philippines [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

[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

[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.frame Any suggestions appreciated. Mary Ann  F1 F1 [[alternative

[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

, 2012 at 6:05 PM, Mary Kindall mary.kind...@gmail.com 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 redundant legend. Your basic problem is that you're using the fill

[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 var1var2

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

+ 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-help@r-project.org mailing list https

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

Custom Legend) g - On Mon, Jan 16, 2012 at 6:55 PM, Felipe Carrillo mazatlanmex...@yahoo.comwrote: Mary: Here's one way. ## change the variable name to whatever title you want on your legend data = melt(data, id=inputs,variable_name=customName) data g - ggplot(data,aes(x

[R] relative frequency plot using ggplot or other function

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 version deleted

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 jto...@gmail.com 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 groups

[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

Re: [R] mode of frequency distribution table

#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 mary.kind...@gmail.com wrote: In a frequency distribution table (bell shaped), how can we find the most

[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/listinfo/r

[R] (Edited) cbind alternate for data frames

] 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/listinfo/r

[R] aggregate function

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/listinfo/r-help

Re: [R] aggregate function

at 11:26 AM, Mary Kindall mary.kind...@gmail.com wrote: Hi I have a data frame with values in following format. param case1 1 a 2 b 2 c 2 d 3 e 4 f how to use aggregate so that it I only one row

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

, 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 the posting guide http

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

, 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 read the posting guide

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

(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 HTML version deleted

[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.names = FALSE)[1

[R] download.file

) != 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/mailman/listinfo/r-help PLEASE do read

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

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 do read the posting guide http

[R] download files using ftp: avoid error

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

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

2011 15:07:15 Mary Kindall wrote: I am planning to download a large number of files from some website. I am using the following script. files2down = c('aaa', 'bbb', ) for (i in 1: len) { print(paste('downloading file', i, ' of total ', len)); url = paste

[R] installing packages systemwide

I installed some downloaded packages in R. I always do $sudo R CMD INSTALL anRpackage.tar.gz 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] save and load in R

can I do that. Regards - M -- - 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

Re: [R] save and load in R

, Mary Kindall wrote: I have a list of txt files that I want to convert into .rdata R data object. filenames 1. ./file1.txt 2. ./file2.txt 3. ./file3.txt 4. ./file4.txt 5. ./file5.txt 6. ./file6.txt 7. ./file7.txt 8. ./file8.txt 9. ./file9.txt 10. ./file10.txt I saved these files

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

) Reduce(merge.all, lapply(L, function(x) aggregate(x[2], x[1], toString))) Cheers !!! - M On Tue, Jun 14, 2011 at 11:27 AM, Gabor Grothendieck ggrothendi...@gmail.com wrote: On Tue, Jun 14, 2011 at 11:21 AM, Mary Kindall mary.kind...@gmail.com wrote: I resolved it. There was a problem in type

[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] combine the data frames into comma separated list.

NA 11 4NANA Thanks -- M On Mon, Jun 13, 2011 at 6:35 PM, Dr. D. P. Kreil (Boku) david.kr...@boku.ac.at wrote: Hi, try ?merge Best, David. On 13 June 2011 23:48, Mary Kindall mary.kind...@gmail.com wrote: Hi R users, I am new to R and am trying to merge

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

:14 PM, Dr. D. P. Kreil (Boku) david.kr...@boku.ac.at wrote: ?write.csv Cheers, David. On 14 June 2011 01:07, Mary Kindall mary.kind...@gmail.com wrote: Thanks for reply. The following code is working but only patially. How to get the condensed values separated by comma. dataframe1

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 ggrothendi...@gmail.com wrote: On Mon, Jun 13, 2011 at 5:17 PM, Mary Kindall mary.kind...@gmail.com wrote: Hi R users, I am new to R

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:

[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

[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.character() to change this? Thanks Mary

[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 -

[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(clin_b)[c(5,7,10:36 )], collapse

[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-project.org mailing

[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

, 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://www.R-project.org/posting-guide.html

[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, Mary

[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 tL

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

(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. Sincerely, Mary A. Marion [[alternative HTML version deleted

[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] 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

[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

[R] (no subject)

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/listinfo/r

[R] Frequency table to histogram

$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 HTML version deleted

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 provide a copy of the first few lines of your

[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

[R] R code -column headings

) 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 __ R-help@r-project.org mailing list https

[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-999#number of data pointsn

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

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-999#number of data pointsn-dim(manpower)[1]#alpha level to use

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

, 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 option, which, if set to TRUE, returns logarithms of probabilities

[R] rmultinomial() function

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-project.org/posting-guide.html and provide commented

[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, Mary Helen Black

[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