[R] CRAN Packages By Date

Hi, why the list of "Available CRAN Packages By Date of Publication" is stopped at 2022-12-25 ? Regards Marco Scazzocchio [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more,

[R] Error with stepwise models selection using NADA library

Hi, I'm working with environmental data which are left-censored and I found the R package NADA which seems to do the job. After fitting a complete model, using the cenreg function,I'd like to do a stepwise procedure to minimize the AIC using the stepAIC function. Unfortunately, after fitting the m

[R] clusplot() and fviz_cluster() functions

t=50)#library(cluster)windows()clusplot(z, fit$cluster, labels=3, lines=0, cex=0.6, col.txt="black", col.p="black")#library(factoextra)windows()fviz_cluster(fit, data=z)#* Thanks if someone has a suggestion! ----

Re: [R] Siegel nonparametric regression / mblm package

that we haven’t overlooked > something more basic, is that such > high breakdown estimates sacrifice some efficiency, that is to say, they > are more variable than other methods > when the data is well behaved, and of course, the Galton data is famously > “almost Gaussian”. > > On Feb

[R] Siegel nonparametric regression / mblm package

standard regression and Theil-Sen regression are practically identical. But the Siegel regression seems to have a bias that I cannot understand. May I ask for a possible explanation? The bias may be related to the number of ties in the set of data? Here's the code and the image. Best regards. M

[R] help

I am trying to write a function to make a matrix of precipitation with this secuencie; PRIMER PERIODO cordex1 <- nc_open("pr_CAM-44i_ICHEC-EC-EARTH_rcp45_r12i1p1_SMHI-RCA4_v1_mon_200601-201012.nc") fullmon1<-ncvar_get(cordex1,"pr") lat<-ncvar_get(cordex1,"lat", start=c(35.5),count=c(20))

Re: [R] I have been unable to run R from within QGIS on my Mac

din your system, or it is not correctly configured to be > used from QGIS > Click here to know more about how to install and configure R to be > used with QGIS" > > [[alternative HTML version deleted]] > Maybe you should ask in a QGIS forum, your QGIS question. But, j

Re: [R] bnlearn and cpquery

distributions from the data or 2) use posterior estimates for the parameters. Cheers, Marco On 13 July 2017 at 06:29, Ross Chapman wrote: > Hi all > > > > I have built a Bayesian network using discrete data using the bnlearn > package. > > > > When I try to run the cp

[R] [R-pkgs] lmvar

arry out cross-validations. The vignette 'Introduction to the package' provides an introduction, the vignette 'Math details' provides mathematical background. Version 1.2.1 is now available. Please send bugs, requests and comments to the package maintainer. Regards,

Re: [R] bnlearn impute problem

ve to upgrade to use it. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of Statistics University of Oxford, United Kingdom __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman

Re: [R] using a regular expression

s for deletion of the role pattern space, so if you are trying to substitute use: s/pattern//g # effectively removing pattern from the text Best Luck, Marco -- Marco Arthur @ (M)arco Creatives __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and

Re: [R] cpquery problem

an average of ~= 13 without considering FFB. If I assume FFB is positive, then I can easily see E(y) ~= 15 and E(y) - 1.96 * 0.96 s.d. ~= 13.5. So ABW < 11 has zero or almost zero probability mass. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of Statistic

[R] type of objects in last session

t;"character""character" op rain t t2 ts_ "character""character""character""character""character" And I noticed everything is type "cha

Re: [R] cpquery problem

7;y', TR = c(9, max(data$TR)), BU = c(15819, max(data$BU)), RF = c(2989, max(data$RF)), n=10^6, method = "lw") 3) look at the parameters in your fitted network and diagnose why this is happening. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department

Re: [R] cpquery problem

that these classes should both typically > return a value for ABW that is very much higher than the threshold value. That may be, but much depends on the specific sample the model was fitted from. How does the fitted network look like? Cheers, Marco -- Marco Scutari, Ph.D. L

[R] bnlearn cpquery

rvals in the evidence. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of Statistics University of Oxford, United Kingdom [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIB

Re: [R] web scraping tables generated in multiple server pages

> 4700 Keele StreetWeb:http://www.datavis.ca > Toronto, ONT M3J 1P3 CANADA > > # what is missing to you ?gsub # aliasing df <- sp.tab[, wanted] # convert to double as.double( # convert to double gsub(',', '.',

[R] Expand 80 columns terminal output

Hi R'ers, I noted that R breaks lines for any output greater than 80 chars, Is there a way to expand this limit ? Thks -- Marco Arthur @ (M)arco Creatives __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.et

Re: [R] bnlearn and TAN network

aining? TAN does not do any sort of feature selection, so all the nodes should be there. As for "showing the network after training", what kind of plot are you looking for? Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of S

Re: [R] bnlearn and TAN network

produce a graph + barplots plot from Rgraphviz. I have never managed to make it work, though. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of Statistics University of Oxford, United Kingdom __ R-help@r-project.org mai

Re: [R] Query in R.

als first. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of Statistics University of Oxford, United Kingdom __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEAS

Re: [R] problem in R

gt; > to programming in R. I see no future in this enterprise. > > > > *Gotta* be a fortune!!! > > cheers, > > Rolf > Jim++ Rolf++ must be. -- Marco Arthur @ (M)arco Creatives __ R-help@r-project.org mailing list -- To

Re: [R] replace NA's with row means for specific columns

tc..) and the row id itself (so consider pass the left hand of assignment entirely). I would use: fmean <- function(row, col_selection) { # homework for you here } Best Regards, -- Marco Arthur @ (M)arco Creatives __ R-help@r-project.org mailing

Re: [R] find all unit vectors which are orthogonal to a given vector

I think that's just what I needed. Thanks! On 2015-10-29 09:55 PM, peter dalgaard wrote: On 30 Oct 2015, at 00:03 , Duncan Murdoch wrote: On 29/10/2015 6:38 PM, Marco Inacio wrote: Is there a function in R to get the set of all unit vectors which are orthogonal to a given vector

Re: [R] find all unit vectors which are orthogonal to a given vector

Ok, thanks! On 2015-10-29 09:03 PM, Duncan Murdoch wrote: On 29/10/2015 6:38 PM, Marco Inacio wrote: Is there a function in R to get the set of all unit vectors which are orthogonal to a given vector? No. Duncan Murdoch __ R-help@r-project.org

[R] find all unit vectors which are orthogonal to a given vector

Is there a function in R to get the set of all unit vectors which are orthogonal to a given vector? __ 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:

[R] help with problem

*When i tried to read a table i**n .csv or .txt format R i get the following message and give some problem in orden to run the data o make test, etc* *Warning message* *In read.table(file = file, header = header, sep =";") *>* incomplete final line found by readTableHeader on 'test.csv* *I trie

[R] Running R on a hosting o server

(8) - daemon to execute scheduled commands (Vixie Cron) So you schedule your task. -- Marco Arthur @ (M)arco Creatives __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEA

Re: [R] Metafor and forest(); not showing 'ilab' and text

and correct that with struct="HAR" ? Sorry for the long mail, I wanted to explain as better as possible, and thank you again for your help, It is incredibly appreciated. Marco On 31 August 2015 at 19:28, Viechtbauer Wolfgang (STAT) wrote: > Have you read help(rma.mv)? It describes i

Re: [R] Metafor and forest(); not showing 'ilab' and text

The solution that you proposed works perfectly, thank you very much. I'll wait for Wolfgang answer as I'm having few doubts about the models. Thanks On 31 August 2015 at 18:34, Michael Dewey wrote: > Comments in line > > On 31/08/2015 16:08, Marco Colagrossi wrote: >>

Re: [R] Metafor and forest(); not showing 'ilab' and text

e there the formulas for it somewhere? So far I only found the formulas for the estimators included in the metafor package. 3. metareg1 <- rma.mv(pc, var, random = ~ 1 | author, mods = ~ pub + SIMiv, data=codebook) Again, if I'm correct this should be a multilevel meta regression (correct me

Re: [R] Metafor and forest(); not showing 'ilab' and text

ssage- >> From: Michael Dewey [mailto:li...@dewey.myzen.co.uk] >> Sent: Tuesday, August 25, 2015 13:23 >> To: Marco Colagrossi; Viechtbauer Wolfgang (STAT) >> Cc: r-help@r-project.org >> Subject: Re: [R] Metafor and forest(); not showing 'ilab' and text &

Re: [R] Metafor and forest(); not showing 'ilab' and text

ce and your time, and many compliments for the package, is guiding me through the use of R for the first time - as you might have guessed. Marco On 24 August 2015 at 16:50, Viechtbauer Wolfgang (STAT) wrote: > I cannot reproduce the issue with 'ilab' not being shown when using '

[R] Metafor and forest(); not showing 'ilab' and text

Hello folks, I have a couple of issues with the metafor package, specifically with the forest graphs. I am currently conducting a Meta-Analysis in economics throughout the metafor package. My meta-analysis has the specific of having different cases from single studies, and this proven to be chall

Re: [R] Mixed-effects model for pre-post randomization design

Thank you very much, Ben. Your answer has been very useful. __ 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

Re: [R] help in anova

I suspect that your "treatment" levels are perfectly nested inside the "soilgroup" levels. If this were the case you can either use one of the factor in your analysis or the other, depending on what you want to analyze. Il giorno Wed, 11 Feb 2015 14:47:18 + Mir Salam ha scritto: > Dear all,

[R] Mixed-effects model for pre-post randomization design

ite where it is explicitly reported to be "confidential"), so I have to review a bit of likelihood theory along with some implementation details. In the meantime, I wonder if anyone here could point out any related documentation to me. Thank you. Marco. __

Re: [R] Strange error while passing string as an argument to the function in bnlearn package

ot;method")= chr "lw" To compute P(A = a | whatever you conditioned on), just sum the corresponding weights over the total weight mass. Since no language trickery is involved, this works reliably. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of Sta

Re: [R] Strange error while passing string as an argument to the function in bnlearn package

e(text="(M=='s')")), evidence=list(lag1.M1='s'), method = "lw") passing str2 as a list. Cheers, Marco -- Marco Scutari, Ph.D. Lecturer in Statistics, Department of S

[R] coplot: extract regression output

cex=1, number=3, columns=3, col="black", pch=16, overlap = 0.1, panel=function(x,y,...) { points(x,y,pch=16) panel.smooth(x,y,span=.8,iter=5,...) abline(lm(y ~ x), col="blue")} ) Thank you. Marco __ R-help@r-project.org mailing li

Re: [R] Median expected survival

Hi, Lars. I don't understand well your question. why don't you simply type pred_leuk Call: survfit(formula = leuk.cox, newdata = leuk_new) records n.max n.start events median 0.95LCL 0.95UCL 1 4242 42 30 1 1 23 2 4242 42 30 7

Re: [R] survplot invert number at risk labels

my version), and this is why I put time.inc=1 in the previous call (This other issue is still present in 3.0.2) I'm going to upgrade my R and check that this little error is actually gone. Thank you Marco 2014-07-03 0:17 GMT+02:00 Frank Harrell : > Please try hard to cr

[R] survplot invert number at risk labels

o solve the problem or to understand it. Thank you Marco. [[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/posti

Re: [R] NA/NaN values in bnlearn package R

owever, you can estimate the CPTs from incomplete data using table() and prop.table() and assemble them in a fitted BN with custom.fit(). On the other hand, maybe it would be better to write an EM wrapper around bn.fit() to make the best of the dependence structure of the data? Cheers,

[R] recall variable in a loop

what I have to write in place of "y" in the last loop in order to refer to xx1, xx2 or xx3?? I want to write something like: for (i in 1:3) { xxi<-1} Thank you for any advice. Marco __ R-help@r-project.org mailing list https

[R] [OT] Advice for medium size data management

me to bother at all. I'm not asking for prepackaged solutions, rather for help in documentation seeking and links to useful documentation or other threads (for example: is it worthwhile using parallel computing?) Thank you to anyone for reading this email. Marco Barbàra. P.S.: I work on a Debia

[R] Exposures in bernoulli glm()

Hi, is it possible to add exposures to a glm with family=binomial()? It's easy to do it for a Poisson/negative binomial: just multiply the mean by the exposure, that is, offset(log(exposure)): but this obviously wrong for a binomial/Bernoulli since the mean must be no bigger than 1. My goal w

Re: [R] proportional weights

No, you are perfectly fine using WLS. The constant of proportionality is the estimated error variance, i.e., the square of the residual standard error (as I think I said earlier). John You're right. That was a little hard for me to grasp. Thanks for the patience.

Re: [R] proportional weights

I think we can blame Tim Hesterberg for the confusion: He writes " I'll add: * inverse-variance weights, where var(y for observation) = 1/weight (as opposed to just being inversely proportional to the weight) * " And, although I'm not a native English speaker, I think there's a spurious c

Re: [R] proportional weights

Thanks for the answers. Dear Marco and Goran, Perhaps the documentation could be clearer, but it is after all a brief help page. Using weights of 2 to lm() is *not* equivalent to entering the observation twice. The weights are variance weights, not case weights. According to your post here

[R] proportional weights

Hello all, can help clarify something? According to R's lm() doc: Non-NULL weights can be used to indicate that different observations have different variances (with the values in weights being inversely *proportional* to the variances); or equivalently, when the elements of weights are positiv

Re: [R] Using "cpquery" function from bnlearn package inside loop

lihood weighting method will work. > Don't > hesitate to correct me if I'm wrong. >From your description, likelihood weighting should be fine. Cheers, Marco -- Marco Scutari, Ph.D. Research Associate, Genetics Institute (UGI) University College London (UCL), United Ki

Re: [R] Using "cpquery" function from bnlearn package inside loop

is easier to use in a programmatic way, but it is much more limited in the events you can use for conditioning. Logic sampling (method = "ls") accepts an arbitrary R expression as a conditioning event; likelihood weighting at the moment accepts only a single instantiation of the conditionin

Re: [R] Using "cpquery" function from bnlearn package inside loop

do so I suggest you should install the latest bugfix snapshot from bnlearn.com to avoid a few other bugs in cpquery(..., method = "lw"). Cheers, Marco -- Marco Scutari, Ph.D. Research Associate, Genetics Institute (UGI) University College London (UCL), United Kingdom __

[R] Besag endive data

ataset, but I can no longer find it. Does anybody know where I can find the dataset? Thanks, Marco -- ***** Marco Bee University of Trento - Department of Economics and Management via Inama, 5 - 38122 Trento, Italy Phone: +39-0461-282296, Fa

[R] IRanges::unlist in package

tried using only Imports: as suggested by Chambers but it breaks. Using Depends does not help. Am I having clashing namespace? Here is my Depends: (or Imports:) line: Depends: Rsamtools, GenomicFeatures, parallel, rtracklayer, edgeR Am I simply missing something? Thanks -- Marco Blanchette,

[R] Help with unzip error -1

I'm running Rx64 3.01 on Windows 7 and until the other day I could easily unzip zip files with unzip (package utils 3.01). Since yesterday however I get an error. Let's say my zipfile file path is called zfile (I'm omitting the output directory for this example). When I try to unzip it I get the fo

[R] colour coded dotchart

ealt with a similar task and found a solution? Thank you very much in advance for any help. Marco [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] Goodness-of-fit tests for Complex Survey Logistic Regression

~dnum+snum, weights=~pw, data=apiclus2) a<-svyglm(sch.wide~1, design=dclus2,family=quasibinomial()) #Null model c<-svyglm(sch.wide~ell*meals, design=dclus2,family=quasibinomial()) anova(a,c) Many thanks for your help Marco [[alternative HTML version deleted]]

Re: [R] Cross validation for Naive Bayes and Bayes Networks

to reassemble observed and predicted class labels and compute your metrics. > I also tried the *e1071* package, but I could not find a way to do > cross-validation. You might be able to trick the tune() function to do it, but I am not sure. Marco -- Marco Scutari, Ph.D. Research Assoc

Re: [R] bnlearn: how to compute boot strength with mmhc and a blacklist

different network every time it is run on > the same data set? This is not surprising, because mmhc() does not have a "start" argument, so it's starting from the same network over and over. There is no way to provide a random seed to mmhc(), so the only way

[R] creating a loop if to create a column of 0 and 1.

Hi, I am new with rstudio and i have a trouble with this program. I store 17 listsCPOD<-list()CPOD[[1]]<- GB1ACPOD[[2]]<- GB1CCPOD[[3]]<- GB1DCPOD[[4]]<- GB2ACPOD[[5]]<- GB2BCPOD[[6]]<- GB2CCPOD[[7]]<- GB2DCPOD[[8]]<- GB3ACPOD[[9]]<- GB3CCPOD[[10]]<- GB3DCPOD[[11]]<- GB4ACPOD[[12]]<- GB4CCPOD[[

[R] best ordination method for binary variables

are > 0, wascores) && any(comm < 0)) { : missing value where TRUE/FALSE needed* *In addition: Warning message:* *In Ops.factor(left, right) : < not meaningful for factor"* Do you think Principal coordinate analysis on a binary distance matrix is a decent strategy? Thanks

Re: [R] If cycle takes to much time...

mmetric matrix, where raws and colums are artists and value I get 1 (or true) if they had an exhibition in the same and in the same place. My unelegant code is working, but for 96000 observation is requiring months and months. Gerrit is very elegant, but i run out of memory... the problem is t

[R] clustering of binary data

1 1 0 1 0 0 0 case12 0 0 0 1 1 0 1 1 0 1 0 0 0 . So, my questions are the following: Is the Jaccard index a good strategy for my kind of data? Is binary distance used in pvclust is theoretically more correct? Is there any alternative to pvclust for testing the significance of my clusters? Tha

[R] fitting a gamma frailty model (coxph)

or in if (history[2, 3] < (history[1, 3] + 1)) theta <- mean(history[1:2, : missing value where TRUE/FALSE needed Does anyone have an idea on how to debug? Yours sincerely, Marco [[alternative HTML version deleted]] __ R-help@r-project.

[R] lapack routines cannot be loaded [Help request]

.4 [4] affy_1.34.0 Biobase_2.16.0BiocGenerics_0.2.0 loaded via a namespace (and not attached): [1] affyio_1.24.0 BiocInstaller_1.4.9 DBI_0.2-5 [4] IRanges_1.14.4preprocessCore_1.18.0 RSQLite_0.11.2 [7] stats4_2.15.2 tools_2.15.2

[R] plotting principal components on geographic map

that the dataset includes relatively few points. Therefore, I would also like to interpolate the values for the uncovered geographical areas. Thanks in advance for any feedback, best marco [[alternative HTML version deleted]] __ R-help@r-projec

[R] problem with the installation of r commander on a mac

Tcl/Tk and XQuark as recommended. I’ve also tried uninstalling R, reinstalling the latest version (2.15.1 signed) and all packages, and following the same official steps again. It didn’t work neither. Could you please help me? R Commander always worked fine on my iMac, until I updated to Mountain Lion.

Re: [R] Specifying a response variable in a Bayesian network

ne of the BN classifiers in bnlearn, naive.bayes()/tree.bayes(), which handle the concept of a response variable more naturally than general-purpose BNs. Hope it helps, Marco -- Marco Scutari, Ph.D. Research Associate, Genetics Institute (UGI) University C

[R] [R-pkgs] lqmm package 1.02

available on https://r-forge.r-project.org/R/?group_id=1396 prior to submission to CRAN. -- Marco Geraci PhD CStat Lecturer in Biostatistics Institute of Child Health, University College London London WC1N 1EH, UK Web: http://marcogeraci.wordpress.com [[alternative HTML version deleted

[R] Partial correlation analysis for censored data in R?

Hi does anyone know if there is any R function that preforms partial correlations analysis for censored data? I found a fortran program on the PSU website at http://www2.astro.psu.edu/statcodes/cens_tau.f but I was wondering if there is anything like that in any R package. Thanks m. _

Re: [R] bayesian gene network construction

kage from Korbinian Strimmer targets exactly that kind of appilication: http://cran.r-project.org/web/packages/GeneNet/ Regards, Marco -- Marco Scutari, Ph.D. Research Associate, Genetics Institute (UGI) University College London (UCL), United Kingdom ___

[R] Type II and III sum of squares (R and SPSS)

and III) fB42. fB * fattA 20. Error(fattB) 10.6667 between factor SourceSum of squares (type II and III) intercept1080.041667 fattA26.0417 Error 9.58333 The most stran

Re: [R] reshaping a dataset for a network

to pajek. I don't any other way to have them like this and save them as text file. Probably is my lack of knowledge... software? Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com -Original Message- From: Marco Guerzoni [mailto:marco.guerz...@unito.it] Sent: Wednesday, Mar

Re: [R] reshaping a dataset for a network

thank you very much I had managed with df <- data.frame(a,b) m <- lapply(split(df, df$a), function(x) x$b) n <- max(sapply(m, length)) a <- t(sapply(m, function (x) c(x, rep(NA, n - length(x) but your solution is much more elegant. best regards Marco Il 3/14/2012 1:54 PM

[R] How to extend a slot of a class?

ax" ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slots .. .. ..@ projargs: chr NA data is a data.frame including the information of long and lat location. I would like to have there a third information: time (but as POSIXct). Is there an elegant way to do this?

Re: [R] reshaping a dataset for a network

I would like: [,1] [,2] [,3] 1 11 27 34 4983 5 124 A solution could be rbind.fill , which does not seem to work with list. thanks Marco Guerzoni, Department of Economics University of Turin Note that your df<-cbind(a,b) produces a matrix, not the data.frame t

[R] reshaping a dataset for a network

8 3 5 12 4 a are vertex of a network, b the edges. In the data the lenght of a is about 5 I read several posts about reshape, reshape2, split, ldply but I couldn't manage to do it. The problem seems to be that the is not a real panel. Any help would be really appreciated, my

[R] SSA in R

Hi everyone I'm here with a brief question: Does anyone have some idea on how to perform Smallest Space Analysis in R? I'm aware of several packages/methods for non metric multidimensional scaling, but no clue about this specific technique. Thanks for any feedback Marco -- View th

[R] suggestions for ANOVA which includes the "year" as a factor

you think a mixed system could be suitable? In case which texts, paper, references, manuals, etc. could you suggest me? How can I compare means in a mixed model (for example LSD test) ? Thank you very much and sorry for bothering you. Sincerely. Marco ___

[R] Strange behaviour with Homals

d in knowing if someone has ever compared the homals solutions with the ones from the categorical quantizations performed with the Aspect package. And, yes, I wrote to package maintainer but I didn't get an answer. Any help is welcome. Thanks in advance. Marco -- "Anyone who lives within t

[R] Strange behaviour with the Homals Package

omeone has ever compared the homals solutions with the ones from the categorical quantizations performed with the Aspect package. And, yes, I wrote to package maintainer but I didn't get an answer. Thanks in advance. Marco -- "Anyone who lives within their means suffers f

[R] Strange behaviour with the Homals Package

omeone has ever compared the homals solutions with the ones from the categorical quantizations performed with the Aspect package. And, yes, I wrote to package maintainer but I didn't get an answer. Thanks in advance. Marco -- "Anyone who lives within their means suffers f

[R] converting matrix in array

Hi everyone, have a small problem trying to converting a dataset in matrix form to an array. Specifically: data include 3D measurement -x,y,z of 59 points in 36 objects. They are stored as a matrix (x) of 2124 rows and 3 columns. What I want to do is to extract each subject's dataset using an arra

Re: [R] Dichotomous variables

oes not solve the problem: same message. So: this is the blind alley. thanks for any feedback marco -- View this message in context: http://r.789695.n4.nabble.com/Dichotomous-variables-tp3721550p3724639.html Sent from the R help mailing list archive at Nabble.com. __

[R] Dichotomous variables

y for this long post marco -- View this message in context: http://r.789695.n4.nabble.com/Dichotomous-variables-tp3721550p3721550.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/ma

Re: [R] smooth scatterplot and geo map

thanks a lot!! -- View this message in context: http://r.789695.n4.nabble.com/smooth-scatterplot-and-geo-map-tp3702374p3702424.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list https://stat.ethz.ch/mail

[R] smooth scatterplot and geo map

n the map, I would like to have a "smoothed" scatterplot of the same superimposed on the map. Tried with the smoothScatter function, but really have no idea how to combine the map with the resulting graph. tx in adv marco -- View this message in context: http://r.789695.n4.nabble.com/smo

[R] smallest space analysis

Hi everyone just a brief question: I'm trying to perform a smallest space analysis on nonmetric variables with R. It seems there are at lear a couple of packages to perform non metric MDS but still not found any info concerning SSA. Any clue? thanks in adv. Marco -- View this message in co

[R] How to merge a matrix and a dataframe with different types of columns

7-15" [14,] "2011-07-16" [15,] "2011-07-17" [16,] "2011-07-18" [17,] "2011-07-19" [18,] "2011-07-20" [19,] "2011-07-21" [20,] "2011-07-22" [21,] NA [22,] NA [23,] NA

[R] How to merge a matrix and a dataframe with different types of columns

7-15" [14,] "2011-07-16" [15,] "2011-07-17" [16,] "2011-07-18" [17,] "2011-07-19" [18,] "2011-07-20" [19,] "2011-07-21" [20,] "2011-07-22" [21,] NA [22,] NA [23,] NA

Re: [R] summing array elements

Thanks a lot for your help marco 2011/5/27 Thomas Lumley > On Fri, May 27, 2011 at 10:06 AM, marco milella wrote: > > Hi to everybody > > I have an array with dimensions 2,4,3,3. > > Wanting to sum the matrices in the first two dimensions, I'm trying to > use

Re: [R] summing array elements

Thanks a lot for your hint marco 2011/5/27 Erich Neuwirth > > arr<-1:72 > > dim(arr)<-c(2,4,3,3) > > apply(arr,1:2,sum) > [,1] [,2] [,3] [,4] > [1,] 297 315 333 351 > [2,] 306 324 342 360 > > apply(arr,3:4,sum) > [,1] [,2] [,3] > [

[R] Why do not prevent users from creating a local copy of c() ?

doing such dangerous things? Thanks, Marco. __ 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, re

[R] summing array elements

Hi to everybody I have an array with dimensions 2,4,3,3. Wanting to sum the matrices in the first two dimensions, I'm trying to use the "apply" function, but with no results. Have to say I'm quite new with R syntax. tx in advance marco [[alternative

Re: [R] help with a survplot

Il giorno Mon, 2 May 2011 07:30:10 -0700 (PDT) Frank Harrell ha scritto: > Please elaborate. It is simply that, generally speaking, i don't like adding numbers to a plot. I eventually realized that riskset cardinality may be a useful indication. However, do not discriminate between events and ce

Re: [R] help with a survplot

Thank you very much. Despite prof. Harrell's support (for whom I feel great esteem) I still remain doubtful about this feature. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://w

[R] R: [BioC] samr - extract genes from siggenes.table

e more than enough... Anyway, most likely other people can give you better hints. All the best, Marco -- Marco Manca, MD University of Maastricht Faculty of Health, Medicine and Life Sciences (FHML) Cardiovascular Research Institute (CARIM) Mailing address: PO Box 616, 6200 MD Maastricht (The Net

Re: [R] Bayesian Belief Networks

On Fri, Dec 24, 2010 at 12:01 AM, Michael Bedward wrote: > Hello Walt, > > Have a look at the bnlearn and deal packages. > > Michael Dear Walt, take a look also to the catnet/mugnet and pcalg packages, and if you have questions about bnlearn feel free to ask. Regards, Ma

  1   2   >