[R] Multiplication of regression coefficient by factor type variable

2018-11-17 Thread Julian Righ Sampedro
s are treated as 'numeric' variables. Though the ols estimates are obviously very different. Thank you in advance for your views on the problem. Best Regards, julian [[alternative HTML version deleted]] __ R-help@r-project.org mailing li

Re: [R] testing whether clusters in a PCA plot are significantly different from one another

2017-01-09 Thread Marchesi, Julian
Dear Micheal So I would be much better off just reporting the PCA as is and conclude what i can from plot cheers Julian Julian R. Marchesi Deputy Director and Professor of Clinical Microbiome Research at the Centre for Digestive and Gut Health, Imperial College London, London W2 1NY Tel

Re: [R] testing whether clusters in a PCA plot are significantly different from one another

2017-01-06 Thread Marchesi, Julian
many thanks david for such a swift response, really appreciate your help cheers Julian Julian R. Marchesi Deputy Director and Professor of Clinical Microbiome Research at the Centre for Digestive and Gut Health, Imperial College London, London W2 1NY Tel: +44 (0)20 331 26197 and Professor

[R] testing whether clusters in a PCA plot are significantly different from one another

2017-01-06 Thread Marchesi, Julian
Rplot_PCA.pdf Description: Rplot_PCA.pdf __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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] Random effects in GAMs

2016-02-13 Thread Julian Chen
I am now trying to use random effects in GAMs developed by Professor Simon Wood. Prof Wood uses s(...,bs="re") to account for the random effects. Random intercepts models or random slopes models are two different types of mixed linear models or general random effects model (Cameron and Trivedi, 20

Re: [R] Random effects in GAMs

2016-02-13 Thread Julian Chen
uggestions? # gam(crash~s(main.traffic.volume, ID, bs=“re")+s(minor.traffic.volume, ID, bs="re"), family=nb(),data=mydata3sg) # Many thanks! On Fri, Feb 12, 2016 at 1:57 PM, Julian

[R] Call R from Matlab

2015-02-27 Thread Julian Keil
happens if I call R from Matlab, not when I start R by itself or use the terminal to start R. I already installed R new, but it did not help. Can anyone enlighten my on this matter? Thanks a lot! Julian Dr. Julian Keil AG Multisensorische Integration Psychiatrische Univer

[R] dbHasCompleted() always returns TRUE for POSTGRES Database with RJDBC?

2014-10-20 Thread julian . bothe
occur (and a Error: "java.lang.OutOfMemoryError: GC overhead limit exceeded" when trying to set a larger size). Does anybody else have or had this problem? All the best Julian adapted Example from help (?dbHasCompleted ), Works ### require("RSQLite"

[R] Error with RJDBC when trying dbUnloadDriver()

2014-09-12 Thread julian . bothe
Postgres_drv) #Error in .valueClassTest(standardGeneric("dbUnloadDriver"), "logical", : # invalid value from generic function 'dbUnloadDriver', class "NULL", expected "logical" -- All the best Julian [[alternativ

[R] Multiple Imputation of longitudinal data in MICE and statistical analyses of object type mids

2014-07-22 Thread Julian Schulze
NULL # save as.mids as.mids(imp_long,.imp=1, .id=2) # doesnt work as.mids(imp_long) # doesnt work Best, Julian __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the postin

Re: [R] rfImpute

2014-06-24 Thread Julian Karch
I am having the same problem trainImp <- rfImpute(train[,!colnames(train)=="score"],train$score) yields Error in randomForest.default(xf, y, ntree = ntree, ..., do.trace = ntree, : NA not permitted in predictors any ideas? -- View this message in context: http://r.789695.n4.nabble.com/r

Re: [R] problem quitting

2013-04-10 Thread Julian Wells
the issue seemed to crop up after I installed RStudio. Best wishes, Julian Wells Begin forwarded message: > From: r-help-requ...@r-project.org > Subject: R-help Digest, Vol 122, Issue 5 > Date: 5 April 2013 11:00:08 GMT+01:00 > To: r-help@r-project.org > Reply-To: r-help@r-proje

[R] Problem quitting

2013-04-03 Thread Julian Wells
27;Permission denied' I'm running R 2.15.3 GUI 1.53 Leopard build 64-bit on a MacBook Pro Processor 2.3 GHz Intel Core i7 Software Mac OS X Lion 10.7.5 (11G63) Julian Wells __ R-help@r-project.org mailing list https://stat.ethz.ch/mailma

Re: [R] Hclust tree to Figtree w/ branch lengths

2013-01-19 Thread Julian Banchier
Hi, This worked: fit <- hclust(d, method="ward") library(ape) p <- (as.phylo(fit)) write.tree(p, file="MyNewick.tre") Cheers, J -- View this message in context: http://r.789695.n4.nabble.com/Hclust-tree-to-Figtree-w-branch-lengths-tp4655990p4656057.html Sent from the R help mailing list arch

[R] Hclust tree to Figtree w/ branch lengths

2013-01-18 Thread Julian Banchier
uot;) d <- dist(data, method = "euclidean") library(ape) fit <- as.phylo(hclust(d, method="ward") library(ctc) write.table(as.phylo(hc2Newick(fit), file="C:/path/output.txt",row.names=FALSE,col.names=FALSE)) Would anyone have an idea about how to make this work?

[R] Indexing in summaryBy

2012-05-15 Thread Julian Wucherpfennig
I'm trying to use a self-written function with the summaryBy function (doBy package). I have lots of data from Monte Carlo experiments comparing different estimators across different (combinations of) parameter values, similar to the following form: colnames(mydata) <- c("X", "b0", "b1", # parame

[R] p-value for hazard ratio in Cox proportional hazards regression?

2011-12-10 Thread Thierry Julian Panje
Hi, I'm new to R and using it for Cox survival analysis. Thanks to this great forum I learned how to compute the HR with its confidence interval. My question would be: Is there any way to get the p-value for a hazard ratio in addition to the confidence interval? Thanks, Thierry -- Thier

[R] transforming a matrix of logical to 0 and 1 while keeping the dim of matrix

2011-07-20 Thread Julian TszKin Chan
eed to copy the dim information to y too. attributes(y)<-attributes(x) Is this a correct way to do it in R? Is there any single step function which can do the something? Thanks Regards, TszKin Julian __ R-help@r-project.org mailing list https:

[R] txtProgressBar and parallelized foreach

2011-06-14 Thread Julian TszKin Chan
? How can i create a progress bar in parallelized foreach? Thanks cl<-makeCluster(core) registerDoSNOW(cl) pb<-txtProgressBar(min = 0, max = n, style = 3) foreach( i = 1:n ,.combine=c) %dopar%{ setTxtProgressBar(pb, i) Sys.sleep(0.1) } close(pb) stopCluster(cl) Regards, TszKin

[R] problems with geom_vline, histograms, scale=free and facets in ggplot

2011-06-10 Thread Julian Burgos
facet_grid(code2 ~ code1) p + geom_vline(aes(xintercept = n.mean), mean.data,col="red") # END R CODE ## I can´t figure out why the second example does not work when having the "scales=free" argument, and why I get the extra face

[R] Question about ggplot2

2011-05-26 Thread Julian TszKin Chan
Hi all, Is there any way for me to to string in the argument of qplot or ggplot? for example qplot(x='carat',y='price',data=diamonds,geom=c('point','smooth')) instead of qplot(x=carat,y=price,data=diamonds,geom=c('point','smooth'))

Re: [R] Multi-dimensional non-linear fitting - advice on best method?

2011-04-26 Thread Julian Gilbey
On Sun, Apr 24, 2011 at 07:02:48PM -0400, Ravi Varadhan wrote: > Julian, > > You have not specified your problem fully. What is the nature of f? Is f a > scalar function or is it a vector function (2-dim)? It's something like this (only a bit worse): given x, work out alpha

Re: [R] Multi-dimensional non-linear fitting - advice on best method?

2011-04-26 Thread Julian Gilbey
her > it is a y1 or a y2 observation. You can also specify a within-pair > correlation. Will take a look, thanks! Julian __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting gu

[R] Multi-dimensional non-linear fitting - advice on best method?

2011-04-23 Thread Julian Gilbey
estimated, and I'll assume the errors to be normal for sake of simplicity. If there were only one response per input, I would use the nls() function, but what can I do in this case? Many thanks, Julian __ R-help@r-project.org mailing list

[R] Conflict between gam::gam and mgcv::gam

2011-03-07 Thread Julian Shaw
) library(mgcv) g3 = mgcv::gam( y ~ s(x) ) plot(g3) detach(package:mgcv,unload=TRUE) library(gam) g4= gam::gam( y ~ s(x,5) ) plot(g4,se=TRUE) Julian Shaw Notice: This e-mail message is intended only for the named recipient(s) above. It may contain

[R] Question about foreach (with doSNOW), is that a bug?

2011-02-24 Thread Julian TszKin Chan
luster(cl) # However if i put require(lmtest) inside the foreach loop, it works again! cl<-makeCluster(2) registerDoSNOW(cl) foreach(i=1:2) %dopar% { require(lmtest) coeftest(lm(length ~ age, data=Mandible, subset=(age <= 28))) } stopCluster(cl)

Re: [R] R-/Text-editor for Windows?

2011-01-31 Thread Julian Burgos
ed to, and it is not for everyone. But, it is totally customizable, and combined with org mode and org babel it becomes a tool as no other. I´ve tried about every editor, including Tinn-R, Notepad++, and Eclipse with StatEt, but after switching to Emacs/ESS there is no turning back. Julian -

[R] post-hoc comparisons in GAMs (mgcv) with parametric terms

2011-01-25 Thread Julian Burgos
hoc test on the categorical predictor in the models where it is significant. Any suggestions? Thanks, Julian -- Julian Mariano Burgos Hafrannsóknastofnunin/Marine Research Institute Skúlagata 4, 121 Reykjavík, Iceland Sími/Telephone : +354-5752037 Bréfsími/Telefax:  +354-5752001 Netfang/Email

[R] rprofile.site

2011-01-07 Thread Julian Pritsch
ource(.trPaths[5], echo=TRUE, max.deparse.length=150) It doesn't do that after I configure (R>configure>permanent(rprofile.site)) but after i closed R/Tinn-R I have to configure R again. Is that necessary? Thank you, Julian __ R-help@r-pro

[R] prediction intervals for (mcgv) gam objects

2010-12-29 Thread Julian
As I understand it,  predict.lm(l ,newdata=nd ,interval="confidence") yields confidence bands for the predicted mean of new observations and lm.predict(l ,newdata=nd ,interval="prediction") yields confidence bands for new observations themselves, given an lm object l.   However with regard to {m

[R] Resource for learning C/R interface

2010-12-18 Thread Julian TszKin Chan
Hi all, Is there any tutorial for learning C/R interface ? Thanks Regards, Julian [[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] How can i select a set of element in a list ?

2010-12-08 Thread Julian TszKin Chan
ength(x) - length(y) ]] > y [[1]] [1] 2 [[2]] [1] 3 [[3]] [1] 4 [[4]] [1] 5 Regards, Julian __ 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] Is there any function can perform "outer lapply" ?

2010-12-06 Thread Julian TszKin Chan
x,y, MoreArgs = MoreArgs, SIMPLIFY = SIMPLIFY, USE.NAMES = USE.NAMES) } Regards, Julian __ 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.htm

Re: [R] help me understand how things work.

2010-09-16 Thread Julian Burgos
r some reason R prints the attributes of the object when you do the inverse but not when you do the square root (I´m curious about why...If anyone has an answer please pitch in). If you only want to get a number, do this: > x=as.numeric(dist(temp)) > class(x) [1] "numeric" > 1/

[R] Multichannel singular spectrum analysis (M-SSA) for R?

2010-07-26 Thread Julian Smith
Does anyone know of any tools to do multichannel(multivariate) singular spectrum analysis (M-SSA) in R? Seems like this would be a useful thing to have. -Hilbert [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://

Re: [R] Strange behaviour when using diff with POSIXt and POSIXlt objects

2010-05-21 Thread Julian Burgos
ocale: [1] LC_COLLATE=Icelandic_Iceland.1252 LC_CTYPE=Icelandic_Iceland.1252 [3] LC_MONETARY=Icelandic_Iceland.1252 LC_NUMERIC=C [5] LC_TIME=Icelandic_Iceland.1252 attached base packages: [1] grDevices datasets splines graphics stats tcltk utils [8] methods base other attache

Re: [R] Strange behaviour when using diff with POSIXt and POSIXlt objects

2010-05-20 Thread Julian Burgos
me4=strptime("2009 06 01 00 57 00",format="%Y %m %d %H %M") > > > > diff(c(time3,time4)) > Time difference of 10 mins > > > > I have version 2.10.1 > > On Thu, May 20, 2010 at 12:36 PM, Julian Burgos wrote: > > Dear list, > > > >

[R] Strange behaviour when using diff with POSIXt and POSIXlt objects

2010-05-20 Thread Julian Burgos
06 01 00 57 00",format="%Y %m %d %H %M") > > diff(c(time3,time4)) Time difference of NA secs ...which is not what I´m looking for. The difference should also be 10 minutes. I burned a few neurons (and searched the documentation) and I cannot figure why this happens. Any ide

[R] GAM for non-integer proportions

2010-02-14 Thread Julian Burgos
higher at intermediate values. It seems that the best way to go is to model this using the "quasi" family with a logit link and a mu(1-mu) variance. Or I am missing something obvious? I will be thankful for any input. All the best, Julian -- Julian Mariano Burgos Hafrannsóknastofnu

[R] linear models with variance in dependent and independent variables

2010-01-31 Thread Julian Burgos
(0.54) 7.4 (0.45) 1.4 (0.63)4.6 (0.37) What would be the best way to fit a multiple regression taking into account the variance of the individual values of the independent and dependent variables? Any references to methods or R packages will be greatly welcomed. Julian -- Julian

Re: [R] How to install old randomForest?

2010-01-13 Thread Julian Ramirez
microsoft. All that is free. I suggest you to read tutorials on how to build packages from sources in Windows. This website might be a good start point http://www.biostat.wisc.edu/~kbroman/Rintro/Rwinpack.html. Hope this helps, Julian Ramirez Research Assistant International Centre for Tropical

[R] package spam for R64-devel

2010-01-13 Thread Julian Ramirez
-tests flag Many thanks for any help you might be able to provide, Julian Ramirez Research Assistant International Centre for Tropical Agriculture, CIAT Colombia [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://

Re: [R] SOM library - where do I find it

2009-11-18 Thread Julian Burgos
You have to download it from CRAN and install it. From the GUI, do Packages->Install package(s). Pretty basic stuff...you should check the documentation before posting. Julian > > > R version 2.9.2 (2009-08-24) - for windows > >> library(SOM) > Error in library(SO

Re: [R] compute differences

2009-09-21 Thread Julian Burgos
rovide commented, minimal, self-contained, reproducible code. -- Julian Mariano Burgos Hafrannsóknastofnunin/Marine Research Institute Skúlagata 4, 121 Reykjavík, Iceland Sími/Telephone : +354-5752037 Bréfsími/Telefax: +354-5752001 Netfang/Email: jul...@hafro.is, jmbur...@u.washington.edu ___

[R] ggplot2 legend text....a basic question

2009-09-14 Thread Julian Burgos
end (for example, "number" instead of "n"). b) How to avoid having a legend for the polygon? Many thanks, Julian -- Julian Mariano Burgos Hafrannsóknastofnunin/Marine Research Institute Skúlagata 4, 121 Reykjavík, Iceland Sími/Telephone : +354-5752037 Bréfsími/Telefax: +354-

Re: [R] Out of memory issue

2009-04-24 Thread Julian Burgos
"PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code." Otherwise...the only other possible suggestion is to use a smaller file. Julian Neotropical bat risk assessments wrote: Hi all, I am

Re: [R] Forcing the extrapolation of loess through the origin

2009-04-14 Thread Julian Burgos
Hi Torsten, If you are fitting a line, why are you using "loess"? Why not simply use "lm" to fit a regression line that goes through the origin? (i.e. with no intercept). Julian jimm-pa...@gmx.de wrote: Hi all, I'm fitting a line to my dataset. Later I want to pr

[R] is there any IV probit estimator in R?

2009-02-23 Thread Julian TszKin Chan
Hi Is there any package in R that allow me to do IV probit, just like the ivprobit command in stata. or do i have to carry out the newey's two-step estimator myself? Regards, Julian [[alternative HTML version deleted]] __ R-h

Re: [R] coalesce columns within a data frame

2008-10-22 Thread Julian Burgos
You could do something like this: > Name.x=c('nx1','nx2',NA,NA) > Name.y=c('ny1','NA','ny3',NA) > Name=Name.x > Name[is.na(Name.x)]=Name.y[is.na(Name.x)] > Name [1] "nx1" "nx2" "ny3" NA Julian

Re: [R] 3 curves / 1 plot

2008-10-06 Thread Julian Burgos
You can use the points() and lines() functions to add points and lines to an existing plot. Julian Michel PETITJEAN wrote: I am a new user of R. Please does somebody knows how to plot 3 datasets (x1,a1),...,(xn,an), (x1,b1),...,(xn,bn), and (x1,c1),...,(xn,cn) on a single x,y plot, each of

Re: [R] Inverting data frame...row wise

2008-09-25 Thread Julian Burgos
How about something like my.data=my.data[,4:1] Julian milicic.marko wrote: Hi, I have the data.frame with 4 columns. I simply want to invert dataset so that last row becomes first... I tried with rev(my_data-frame) but I got my columns inverted... not my rows Thanks

Re: [R] Question about col.names in write.csv

2008-09-05 Thread Julian Burgos
primera fila, hice un vector con el nombre y nueve espacios en blanco (para completar las diez filas que tiene mi tabla). Al grabar la tabla usando write.csv, el comportamiento "default" es guardar los nombres de las columnas y filas. Saludos, Julian Luz Milena Zea Fernandez wrote:

Re: [R] Mclust - which cluster is each observation in?

2008-07-21 Thread Julian Burgos
by doing my.clusters$membership Hope this helps, Julian cnagy wrote: I'm trying to test a method of identifying individuals (birds) based on measured data (their calls). I have test data from known individual birds, and I am using the Mclust package to see if the program can corr

Re: [R] using which to identify a range of values

2008-07-18 Thread Julian Burgos
Do simply which(a<100 & a>=200) Julian sj wrote: Hello, I am trying to identify values that fall within a certain range. I thought that I might be able to use the which function to do this but I have been unable to figure out a way to do it. Perhaps a little code will illustrate

Re: [R] more columns that column names

2008-06-02 Thread Julian Burgos
Hi Paul, The easiest thing to to is to open the file using a text editor (Notepad will do) and examine the first few lines. You can add add a column name if needed. Julian Paul Adams wrote: Hello to everyone, I have gotten my file to print to screen but when I use read.table I am getting

Re: [R] Problemas usando jri0.4-1 y R 2.7.0

2008-06-02 Thread Julian Burgos
Hola Borja, Creo que vas a tener muy pocas respuestas a menos que escribas a la lista en inglés. Lamentablemente yo conozco poco de Java y no puedo responder tu pregunta. Saludos, Julian Borja Soto Varela wrote: Hola, es la primera vez que mando un correo a cualquiera de las listas de

Re: [R] Excluding/removing row and column names on text output files

2008-05-30 Thread Julian Burgos
Hi Steve, You can use write.table: write.table(x, file="/Users/Desktop/Data.txt", sep="",row.names=F,col.names=F) Cheers, Julian Stropharia wrote: Dear R users, I've had no joy finding a solution to this online or in any of my R books. Many thanks in advance

Re: [R] Plot colors

2008-05-30 Thread Julian Burgos
You could do something like this: mydata=c(1,2,1,1,6,7,-1,-1,5,-1) color= as.numeric(mydata== -1) +1 plot(mydata,col=color) This will give you a plot where the -1's are in red (color = 2) and the other numbers in black (color=1). Julian uv wrote: Hi. I am plotting graphs for v

Re: [R] Binning groups

2008-05-28 Thread Julian Burgos
see ?findIntervals francogrex wrote: Hi, this is probably quite simple but I can't seem to do it correctly. I have a data frame of counts of infections in different ages; something like: count=c(1, 1, 1, 2, 2, 2, 3, 3, 3, 3, 4, 4, 4, 4, 4, 7, 8, 8, 9, 9, 10, 11, 15, 17, 17, 17, 17, 19, 19, 19,

Re: [R] Help on nested FOR loops

2008-05-19 Thread Julian Burgos
rint (x) Var1 Var2 1 11 2 21 3 31 4 41 5 51 6 12 7 22 8 32 9 42 1052 1113 1223 1333 144 3 1553 Hope this helps, Julian Philip Twumasi-Ankrah wrote: I am new to more radical pro

Re: [R] Can R handle large dataset?

2008-05-14 Thread Julian Burgos
Depends on the RAM in your machine. And in your definition of 'handle'. You may be able to load a very large dataset into R, but won't be able to use some functions that require additional memory. A vague answer to a vague question... :) Julian Mingjun Huang wrote: Hello,

Re: [R] A very simple question

2008-05-14 Thread Julian Burgos
Try this: k=c(1,1,1,2,2,1,1,1) > k[(k!=1)] [1] 2 2 > k[(k!=2)] [1] 1 1 1 1 1 1 > k[(k!=3)] [1] 1 1 1 2 2 1 1 1 Julian Shubha Vishwanath Karanth wrote: Hi R, Suppose l=c(1,1,1,2,2,1,1,1) k[-which(k==1)] [1] 2 2 k[-which(k==2)] [1] 1 1 1 1 1 1 But, k[-w

Re: [R] wavlet analysis

2008-05-14 Thread Julian Burgos
Hi Stephen, Your link doesn't work. In any case, check out the wavCWT function in the wmtsa package. Julian stephen sefick wrote: http://ion.researchsystems.com/cgi-bin/ion-p I would like a continuous wavelet transform. I have downloaded wavethresh, Rwave, and waveslim. I would li

Re: [R] use list elements to subtract values from the dataframe

2008-05-07 Thread Julian Burgos
Try get(paste("wf$",fl[[1]],sep="")) See ?get Julian Dirkheld wrote: Hi, I have a dataframe wf existing of a header with different labels and beneath the values of those labels : wf: label1 label2 ... 0,450,21 0,100,45 I have a list fl &

Re: [R] How to count the overlapped in two vectors

2008-05-01 Thread Julian Burgos
See ?match ss wrote: Dear list, If I have two vector, t1 and t2, of different lengths. Is there an easy way to count the number of the overlapped in two vectors and show the result in the graph? Thanks much, Alex [[alternative HTML version deleted]]

Re: [R] How do you test for "consecutivity"?

2008-04-29 Thread Julian Burgos
rop=prop/length(x) return(prop)} This function first identifies which numbers in your original vector are part of a sequence of consecutive numbers. Julian Julian Burgos wrote: Hey Anthony, There must be many ways to do this. This is one of them: #First, define a function to calculate

Re: [R] How do you test for "consecutivity"?

2008-04-29 Thread Julian Burgos
ortions of consecutive numbers Julian Anthony28 wrote: I need to use R to model a large number of experiments (say, 1000). Each experiment involves the random selection of 5 numbers (without replacement) from a pool of numbers ranging between 1 and 30. What I need to know is what *proportion

Re: [R] means and variances of several groups in the matrix

2008-04-24 Thread Julian Burgos
x within each group. > means=tapply(x,group,mean) > vars=tapply(x,group,var) > means 1 2 3 -0.14711206 0.28314274 -0.07861427 > vars 1 2 3 1.4584971 0.3611996 0.6300624 Julian kathie wrote: > Dear R users, > > I have

Re: [R] Exporting of R dendrogram object in the format of E isen´s Cluster

2008-03-23 Thread Julian Lee
Hi You can try using the ctc package to export your dendrogram objects and use Cluster to open them. This package is part of biocondutor. >>source("http://bioconductor.org/biocLite.R";) >>biocLite('ctc') regards Julian Lee - Original Message - From: Ji

Re: [R] cumsum list..

2008-03-13 Thread Julian Burgos
In this case you can simply do cumsum(a[x,]+a[y,]) Julian yoo wrote: > Hi all, i have the following.. > > a <- data.frame(data = seq(1,10)) > > i have indices: > x <- c(1, 5, 3, 9) > y <- c(2, 7, 4, 10) > > I want the cumsum of a[1:2], a[5:7], a[3:4]

Re: [R] I need to buy a book in R

2008-03-03 Thread Julian Burgos
Hello Kayj, There are very good tutorials at the R website. See here: http://cran.r-project.org/other-docs.html Julian kayj wrote: > Hi All, > > I am a new user in R and I would like to buy a book that teaches me how to > use R. In addition, I may nees to do some advanced

Re: [R] using an element of an array as a new object

2008-01-16 Thread Julian Burgos
See ?assign. ppaarrkk wrote: > I have an array called filesBox. I want to take each element of the first > column and assign a dataset to it. > > For example : > > filesBox[4,1] returns > > [1] "fileR" > > Then I want to assign "fileR" which exists as a text file to the R object > "fileR" lik

[R] exact method in coxph

2008-01-16 Thread Julian Wucherpfennig
elp says that the exact method is computationally demanding, but even after days the computing it won't finish. Also, if I include a frailty-term, the exact method gives me results in no time. Is my setup incorrect? Many thanks in advance! - Julian Wucherpfennig PhD Student Politica

Re: [R] bootstrap sampling

2008-01-15 Thread Julian Burgos
See rnorm(). If you are sampling from a continuous normal distribution, it makes no sense to define a sample with replacement, because the probability of sampling twice the same number is zero. Julian sigalit mangut-leiba wrote: > Hello, > How do I sample observations with replacement

Re: [R] how to go to a line in R

2008-01-15 Thread Julian Burgos
Tinn-R also has this option. I suspect most editors will also do. Julian -Halcyon- wrote: > RWinEdt has line indication. You might want to try that. > > > Uwe Ligges-3 wrote: >> This depends on the editor you use for writing R code rather than on R. >> >> U

Re: [R] Great looking plot - but what does it mean?

2008-01-07 Thread Julian Burgos
Hi Mika03, It would be useful to know what function you used to create your plot. Assuming you used boxplot, do this: ?boxplot ?boxplot.stats Julian mika03 wrote: > > http://www.nabble.com/file/p14668788/paragraphs.png > > Hi, > > R is is world full of wonders... I cr

Re: [R] Calculate remainer

2007-12-20 Thread Julian Burgos
Hi Livia, There are several ways to do this. Try: a=50/12 floor(a) will give you the entire portion, and a-floor(a) will give you the remainder. Julian livia wrote: > Hello everyone, > > I have got a question about a simple calculation. If I would like to > calculate 50/12 an

Re: [R] detailed calculation of two way anova with unbalanced design

2007-12-14 Thread Julian Burgos
Hi Jack, Any intro stats book should have them. See for example chapter 11 in Sokal and Rohlf (2nd ed., 1981). Julian Jack Luo wrote: > Dear list, > > Could someone show me where can I find the detailed formula on how to > calculate the two way anova with unbalanced design? Say, i

Re: [R] [OT] vernacular names for circular diagrams

2007-12-13 Thread Julian Burgos
I should say that the name of this chart varies even among Spanish-speaking countries. In Argentina is "diagrama de torta" which is something like "cake-chart". Julian ahimsa campos-arceiz wrote: > Two non-eatable examples from Spain and Japan: > > in Span

Re: [R] Importing Large Dataset into Excel

2007-12-12 Thread Julian Burgos
large) and review line 528, instead of reviewing the spreadsheet in Excel. You should be able to spot the problem right away. Julian Wayne Aldo Gavioli wrote: > Hello all, > > > I seem to be having a problem importing a data set from Excel into R. I'm > using > the &q

Re: [R] combine variables to matrix

2007-12-12 Thread Julian Burgos
imply data.matrix=as.matrix(data) Julian Andre Jung wrote: > I just got stuck with a quite simple question. I've just read in an > ASCII table from a plain text file with read.table(). It's a 1200x1200 > table. R has assigned variables for each column: V1,V2,V3,V4,... > Fo

Re: [R] matrix graph

2007-12-12 Thread Julian Burgos
The basic functions you need are image() contour() although I like better the plot.surface() function in the 'fields' package. Julian threshold wrote: > Hi All, simple question: > do you know how to graph the following object/matrix in a 'surface manner': > >

Re: [R] Help rewriting looping structure?

2007-12-10 Thread Julian Burgos
it of clarification I may help you some more. Julian TLowe wrote: > Hey Folks, > > Could somebody help me rewrite the following code? > > I am looping through all records across 5 fields to calculate the cumulative > percentage of each record (relative to each individual field).

Re: [R] if/else for plot/lines?

2007-12-07 Thread Julian Burgos
previous plot)} Julian Roger Levy wrote: > I'm interested in writing a function that constructs a new plot on the > current graphics device if no plot exists there yet, but adds lines to > the existing plot if a plot is already there. How can I do this? It > seems to me that the

Re: [R] image() plot with z not in matrix format

2007-11-29 Thread Julian Burgos
nction. I hope that helps, Julian [EMAIL PROTECTED] wrote: > Hello R cracks > > The image() function requires strictly increasing x and y values and z as a > matrix. > > Actually, I don't have equally spaced variables, but anyway want to plot an > colored image() (wit

Re: [R] Filling in a Zero Matrix

2007-11-26 Thread Julian Burgos
clearly what you are trying to achieve, we can help you to simplify your code. Julian Julian Burgos wrote: > Hi Amy, > Many (perhaps most) of the people on the list do not receive emails with > html...so we can't see colored text. Also it would be helpfully to have > a bit of you

Re: [R] Filling in a Zero Matrix

2007-11-26 Thread Julian Burgos
lease explain what are the 'xlabel' and the 'ylabel' that you mentioned. I don't see those variables in your code. Julian Amit Patel wrote: > Hi > I am very new to R and statistical programming in general. I am trying to > reorder data from a .csv file. I have

Re: [R] Reconstruct array dataset

2007-11-21 Thread Julian Burgos
(as.matrix(x)) x=matrix(x,ncol=11,byrow=T) Julian marcg wrote: > This was my intention, but I'm not able to read it in as a vector (because i > don't know the function, neither I can convert the read in table to a vector > and then to matrix or directly. > > What did I miss o

Re: [R] better curve

2007-11-21 Thread Julian Burgos
analytical work, it is probably a good idea to learn about the different methods, its assumptions and limitations. Julian mysimbaa wrote: > http://www.nabble.com/file/p13880048/Fluctuation.jpeg > > Hi R users, > > I have collected data which I plot(x,y).The problem it has oscillat

Re: [R] Logarithmic axis

2007-11-20 Thread Julian Burgos
Like always, there is much to be learned from the R-help list! Another message had a much simpler approach. plot(xy.coords(reizstaerke, kennlinie1, log="x"), log="x") Julian Julian Burgos wrote: > Hey John, > > You can do simply > > plot(reizstaerke, k

Re: [R] Logarithmic axis

2007-11-20 Thread Julian Burgos
=F,xlim=log10(c(1,1000))) box() axis(1,at=log10(c(1,10,100,1000)),label=c(1,10,100,1000)) axis(2) Julian John Wiedenhoeft wrote: > Hi there, > > I guess this must be a standard issue, but I'm starting to go crazy with it. > I > simply want a plot with the x axis being log

Re: [R] Trying to get around R

2007-11-20 Thread Julian Burgos
e at the bottom of every message send through this list. The posting guide is not hard to find, and (in my opinion) it isn't long or difficult to understand (in particular for anyone taking college level statistics). Julian Loren Engrav wrote: > I am a newbie to R and Bio emails and > &g

Re: [R] delete a row in a matrix

2007-11-19 Thread Julian Burgos
columns, and you'll get everything except the rows/columns referenced. Julian Barb, Jennifer (NIH/CIT) [E] wrote: > Can anyone tell me how to delete a row in a matrix? I have searched > around and couldn't find a straightforward way to do this. > > Thanks, any help will be

Re: [R] Trying to get around R

2007-11-19 Thread Julian Burgos
should not expect for people on the list to do your homework for you. That is a big no-no. Cheers, Julian Epselon wrote: > I have three problems I am trying to simulate, that I am having difficulty > getting around with. > > Problem 1. > I want to determine the 85 percentile (

Re: [R] About print a label in plot

2007-11-16 Thread Julian Burgos
If you don't know the class of the object, you can do this: class(Disease.FL) Julian affy snp wrote: > Dear list, > > Hello! I have a question about how to print a label in the plot. > I am using the following code: > > width=5);plot(Disease.FL, index=i, type="Single&qu

Re: [R] About print a label in plot

2007-11-16 Thread Julian Burgos
posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > -- Julian M. Burgos Fisheries Acoustics Research Lab School of Aquatic and Fishery Science University of Washington 1122 NE Boat Street Seattle, WA 98105 Phone: 206-221-6864 __

Re: [R] Plot problem

2007-11-15 Thread Julian Burgos
self-contained, reproducible code' so we can find the problem. As a general comment, you can do any number of plots on a device (a window or a pdf file). The limit is only given by the number and size of the plots and the size of the device. Julian affy snp wrote: > Dear list, >

Re: [R] a repetition of simulation

2007-11-15 Thread Julian Burgos
summary(log_v) Julian sigalit mangut-leiba wrote: > Hello, > In addition to my question a few days ago, > Now I have a matrix of the coefficients, > how can I see all the P.Values (Pr(>|z|)) of the covariates from the 1000 > iterations? > I tried names(log_v) and couldn&#

Re: [R] Log random number

2007-11-14 Thread Julian Burgos
ogarithm is not defined for negative numbers. Julian Wensui Liu wrote: > dont think so, unless i miss something here. > please do check the range of normal random number and the domain of > log function. > > On 11/14/07, Tobias Schlottmann <[EMAIL PROTECTED]> wrote: >

Re: [R] How to get row numbers of a subset of rows

2007-11-14 Thread Julian Burgos
One way to do this is range(which(B[,2]==1)) Julian affy snp wrote: > Hello list, > > I read in a txt file using > > > by specifying the row.names=NULL so that the rows are numbered. > Below is an example after how the table looks like using > > >

  1   2   >