[R] Calling Fortran from C++

2009-06-01 Thread Giuseppe
Hi, can anybody point me to a package with C++ code that call Fortran subroutines? I am trying to do the same thing but we scarce success. Error in dyn.load("utils.so") : unable to load shared library 'utils.so': dlopen(utils.so, 6): Symbol not found: _robcovf Referenced from: utils.so E

Re: [R] Calling Fortran from C++

2009-06-01 Thread Giuseppe
> From: r-help-boun...@r-project.org *> [mailto:r-help-bounces_at_r-project.org] On Behalf Of Thomas Lumley * *> Sent: Monday, June 01, 2009 9:33 AM * *> To: cls59 * *> Cc: r-help_at_r-project.org * *> Subject: Re: [R] Calling

[R] Lyapunov Discrete Time Equation

2010-01-29 Thread Giuseppe
Dear all, I need to solve the following Lyapunov Matrix equation: C=ACA' + B, with A and B given square symmetric matrices. Does anyone knows of a package that can solve the lyapunov matrix equation in R? Or even a C/Fortran implementation? I did not find one on netlib. Thank you. __

[R] liblapack

2010-01-29 Thread Giuseppe
Dear all, I am trying to put together in a package some functions that maybe of interest to other. A particular function of the package uses "zgges" from Lapack with a .Fortran call. Since R does not load all symbols from liblapack, I load all the liblapack symbols (lazyload, it is) with a .First.

[R] ERROR need finite 'ylim' values

2010-06-13 Thread Giuseppe
;,xlab="Surgical Pathology") My question is two fold: Why does not the subset command work within the stripchart command? Why the subsetted table cannot be used in the stripchart command, when the plotting variable is the same previously used in the subsetting process? Thank you G

Re: [R] ERROR need finite 'ylim' values

2010-06-13 Thread Giuseppe
Peter Ehlers wrote: > > Giuseppe, > > See comments below. > > On 2010-06-13 10:24, David Winsemius wrote: >> >> >> Giuseppe wrote: >>> >>> Hello: >>> I use R with MAC >>> I have a simple data table, numeric and text colu

[R] Help with multidimensional array

2011-03-28 Thread Giuseppe
I have the following situation. H12 is an array of dimension (n,k,m) and hessian_lambda is a numeric of length m. I need to multiply each matrix H12[,,1], H12[,,2], ..., H12[,m] by hessian_lambda[1], hessian_lambda[2], ..., hessian_lambda[m], respectively, and then add the resulting (n,k) matrice

Re: [R] Help with multidimensional array

2011-03-28 Thread Giuseppe
[2,]   11   14   17 > [3,]   12   15   18 > > , , 3 > > [,1] [,2] [,3] > [1,]   19   22   25 > [2,]   20   23   26 > [3,]   21   24   27 > >> xx <- 1:3 >> Reduce('+', lapply(xx, function(k) xx[k] * m[, , k])) >      [,1] [,2] [,3] > [1,] 

[R] Sweave and accented letters

2011-11-14 Thread Giuseppe
I often use Lyx/Sweave and I typically write in english. Today I had to write a document in Italian and, as many of you know, many italian popular words use è, ù, é. ò, etc. I discovered that if I type in Italian (that is there is at least one letter with accent) with the Sweave module selected

Re: [R] Sweave and accented letters

2011-11-15 Thread Giuseppe
> address this issue better. > > Regards, > Yihui > -- > Yihui Xie > Phone: 515-294-2465 Web: http://yihui.name > Department of Statistics, Iowa State University > 2215 Snedecor Hall, Ames, IA > > > > On Mon, Nov 14, 2011 at 7:21 PM, Giuseppe wrote: >> I often us

[R] Fwd: Using Geographic Information Systems and Remote Sensing to study disease vector habitat

2015-04-16 Thread Giuseppe Amatulli
ll-for-applications Best Rigards -- Giuseppe Amatulli, Ph.D. Department of Ecology and Evolutionary Biology, Yale University. Jetz Lab, OML Room 405 P.O. Box 208106 New Haven, CT 06520-8106 Teaching: spatial-ecology.net Work: http://sbsc.yale.edu/giuseppe-amatulli <http://www.spatial-e

[R] Problem with Kruskal–Wallis test

2018-12-21 Thread Giuseppe Cillis
Dear all, I am a beginner with R (and also with the statistics) for which I hope to be clear. I should do this non-parametric test on data I extracted from maps. In practice I have a column that represents the landscape Dynamics of a certain time period (there are 3 dynamics, each of them marked by

[R] Fwd: Problem with Kruskal–Wallis test

2018-12-25 Thread Giuseppe Cillis
Dear all, I am a beginner with R (and also with the statistics) for which I hope to be clear. I should do this non-parametric test on data I extracted from maps. In practice I have a column that represents the landscape Dynamics of a certain time period (there are 3 dynamics, each of them marked by

Re: [R] Problem with Kruskal–Wallis test

2018-12-25 Thread Giuseppe Cillis
Dear Michael, Thanks for your answer. So, I'm not an expert in R and statistics, how can I create this interval of confidence of groups? Thanks Gc Il giorno sab 22 dic 2018, 13:34 Michael Dewey ha scritto: > Dear Giuseppe > > If I understand you correctly you have a very large sam

[R] Sankey Diagram

2019-03-07 Thread Giuseppe Cillis
Hi, I have a .csv hat is organized in this way: code88;code13;sup_88_13 1;1;77.67576251 1;2;10.1 1;3;2.4 1;4;0.0 1;5;0.2 1;6;0.1 .. I would create a simple Sankey diagram, how can I do it? I have little experience with R unfortunately ! thanks Giuseppe [[alternative HTML

[R] Union of columns of two matrices

2008-08-06 Thread Giuseppe Paleologo
re x and y are vectors. I could not come with anything nice. Any ideas? Giuseppe -- Giuseppe A. Paleologo :: Email: [EMAIL PROTECTED] :: AOL: gappy3000 :: Skype :: gappy3000 :: Gtalk: paleologo :: Mobile: 917.331.3497 fact: 2^32,582,657-1 is a prime [[alternative HTML v

Re: [R] Union of columns of two matrices

2008-08-07 Thread Giuseppe Paleologo
Aug 2008, Dan Davison wrote: > > On Wed, Aug 06, 2008 at 06:32:43PM -0400, Giuseppe Paleologo wrote: >> >>> I was posed the following problem/teaser: >>> >>> given two matrices, come up with an "elegant" (=fast & short) function >>> th

[R] A manipulation problem for a large data set in R

2008-08-27 Thread Giuseppe Paleologo
I have two questions for the group. One is very concrete, and is dangerously close to a "please do my homework" posting. The second follows from the first one but is more general. I would welcome the advice of experienced R users. As for the first one: I have a data frame with two variables X Y

[R] Strange zoo behaviour, possible bug?

2010-04-28 Thread Giuseppe Milicia
06 00:00:00 -100 2010-01-07 00:00:00 -100 2010-01-08 00:00:00 -100 2010-01-09 00:00:00 -100 2010-01-10 00:00:00 -100 Cheers, //Giuseppe MAKO This email and any files transmitted with it are confide...{{dropped:14}} __ R-help@r-project.org mail

Re: [R] Strange zoo behaviour, possible bug?

2010-04-28 Thread Giuseppe Milicia
't bump on this earlier... I have a function which can take as input zoo object whose coredata is either a matrix or a vector, to do that I'll have to handle the two cases slightly different. Not a bit deal, but kind of inelegant. Cheers, //Giuseppe ___

[R] Rcpp & static variables

2010-04-30 Thread Giuseppe Milicia
... The question is, can we use static variable to reliably preserve state between Rcpp calls? Cheers, //Giuseppe MAKO This email and any files transmitted with it are confide...{{dropped:17}} __ R-help@r

[R] corrgram (package corrgram): how to plot multiple correlograms in the same page?

2013-08-22 Thread Giuseppe Pagnoni
know a workaround for this? Many thanks in advance for any advice best giuseppe -- Giuseppe Pagnoni, PhD Dip. Scienze Biomediche, Metaboliche e Neuroscienze Sezione Fisiologia e Neuroscienze Univ. di Modena e Reggio Emilia Via Campi 287 I-41125 Modena, Italy Tel: +39-059-205-5742 Fax: +39-059-205

Re: [R] confidence / prediction ellipse

2013-01-27 Thread Giuseppe Amatulli
prediction ellipse should be more appropriate. But i'm not able to track back the function radius <- sqrt(dfn * qf(level, dfn, dfd)) in order to change it and draw a prediction ellipses. Regards Giuseppe On 26 January 2013 17:19, Bert Gunter wrote: > Well, I'd guess you have to f

Re: [R] confidence / prediction ellipse

2013-01-27 Thread Giuseppe Amatulli
% probability that my prediction fall inside the ellipse. Thanks again Giuseppe On 27 January 2013 17:26, Bert Gunter wrote: > All: > > Aha! -- The light dawneth, methinks (maybe...) > > Giuseppe: I re-read the SAS link you sent and if I have parsed it > correctly, what S

Re: [R] confidence / prediction ellipse

2013-02-07 Thread Giuseppe Amatulli
redictions derived by not-linear model (such as MARS and RF). Is it correct? or should i use an alternative solution ? Moreover, I was expecting that the abline (lm(b,a)) would be correspond to the main axis of the prediction ellipse, but is not this the case. why? Thanks in advance Giuseppe On

[R] fit an nls model

2013-02-12 Thread Giuseppe Amatulli
of-squares decrease. Should i just guess and change the main function (B + A * log(x06Veg)) with another one (e.g. B + A * sqrt(x06Veg)) and re-run or there is a better way? Thanks in Advance Giuseppe -- Giuseppe Amatulli Web: www.spatial-ecology.net [[alternative HTML version delet

[R] Function nlme::lme in Ubuntu (but not Win or OS X): "Non-positive definite approximate variance-covariance"

2013-06-07 Thread Giuseppe Pagnoni
s in advance for any suggestions, best giuseppe -- Giuseppe Pagnoni Dip. Scienze Biomediche, Metaboliche e Neuroscienze Sezione Fisiologia e Neuroscienze Univ. di Modena e Reggio Emilia Via Campi 287 I-41125 Modena, Italy Tel: +39-059-205-5742 Fax: +39-059-205-5363

Re: [R] curve fitting an exponential function

2012-10-08 Thread giuseppe calamita
Hi Waltenegus, you should consider to show us your code and the data you used to fit the curve. If you can't or if you prefer you can simply show the code and some sample data on which the code can be run. You'll find nice and useful tool in nlstools package. Hope this help

[R] Repeated-measures anova with a within-subject covariate (or varying slopes random-effects?)

2012-12-02 Thread Giuseppe Pagnoni
or(subj /(run + stim * cond)), data=rt.df) Perhaps should lme() be used instead (and with which formula?)? Many thanks in advance to anybody who'd be so kind to offer their advice or tip on this. I have scoured the web and some textbooks for a few days now, but to little avail. very best

[R] boxplot with x-axis time

2014-03-13 Thread giuseppe calamita
Dear R-users, I want to plot boxplots of a single variable collected a few times during almost one year and I would like the x-axis to recognize the date-class of the variable. I found some topics in the archive but: - some questions were poorly posed (http://r.789695.n4.nabble.com/Boxplot-with

[R] SPATIO-TEMPORAL ANALYSIS AND BIG DATA PROCESSING USING FREE AND OPEN SOURCE SW

2013-09-20 Thread Giuseppe Amatulli
Staff* Dr. Raul Zurita-Milla (University of Twente, NL) Dr. Giuseppe Amatulli (Yale University, USA) Dr. Stefano Casalegno (University of Exeter, UK) Dr. Pieter Kempeneers (VITO, BE) -- Giuseppe Amatulli Web: www.spatial-ecology.net [[alternat

Re: [R] Hausman test for endogeneity

2010-10-09 Thread Giuseppe Marinelli
e looking for. Cheers Giuseppe [1] http://cran.r-project.org/web/packages/systemfit/systemfit.pdf __ 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] Multivariate GARCH

2010-10-20 Thread Giuseppe Grillo
HI! I'm student of university and i need an example (an application) of Multivariate GARCH model for replicate it in my thesis. It's better with model CCC e DCC. Best Regards Giuseppe for answer: bepperozz...@hotmail.com [[alternative HTML versi

[R] Change intervals in apply.weekly

2010-12-06 Thread Giuseppe Vittucci
006-07-17 36210.00 ... p <- apply.weekly(Close,mean) [,1] 2006-07-03 3.00 2006-07-10 36500.20 2006-07-17 35802.60 ... Is there any command to set the week starting from Monday? Thanks Giuseppe __ R-help@r-project.org mailing list

[R] From ts to xts

2010-12-06 Thread Giuseppe Vittucci
I have the following problem. After having created an xts importing data from a csv with: data <- read.csv(mib.csv, header = TRUE, dec = ".", sep="\t") dates <- as.POSIXct(strptime(data[,1],format="%m/%d/%Y")) mib <- xts(data[,c(2:6)],order.by=dates)) I work out weekly log returns: p <- apply.we

[R] linear model with similar response predictor

2010-08-05 Thread Giuseppe Amatulli
pts is wrong it should be equal to 0. I compile the R 2.11.1 in the Fedora 13. Should i download a library? Some configure options are missing during the installation? or is a bug? Thanks in advance Giuseppe Amatulli __ R-help@r-project.org mailing list

[R] Regression Error Characteristic Curves

2010-08-11 Thread Giuseppe Amatulli
simple example all the tuning parameter settings. Thanks in advance Giuseppe Amatulli __ 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] Please help

2012-07-04 Thread giuseppe calamita
Dear Jun I think you should consider to submit your question to R-SIG (special interest group) about spatial data http://r-sig-geo.2731867.n2.nabble.com/ This will improve the probability you get some help. Cheers Giuseppe Calamita - Giuseppe Calamita PhD at CNR-IMAA Italian National

Re: [R] problem in lodging package

2012-07-04 Thread giuseppe calamita
you need to specify both the package containing the function and then the function E.G. pckg:function(parameters) hope this helps - Giuseppe Calamita PhD at CNR-IMAA Italian National Council of Research - Institute of Methodologies for Environmental Analysis, Tito Scalo -Potenza ITALY

Re: [R] svm with GRASS GIS

2012-03-05 Thread giuseppe calamita
Hi Ethienne, we (me and Gab) would like to thank you ! Finally we got what we were looking for but we did it without use the "raster" packagebut we are going to try also with it to see if it allow to have faster computation or data manipulation, If you're interested we can show you the code we

Re: [R] svm with GRASS GIS

2012-02-18 Thread giuseppe calamita
Dear Etienne, I'm a colleauge of Gabriele and I'm more into R (but he is more into GRASS). I'll try to explain you what we didi so far 1) Our ASTER images, (B1, B2 and B3) have 8363134 pixels; we made a subset in order to have training data sets: that is, for each band (B1,B2 and B3) 916 pixels w

Re: [R] svm with GRASS GIS

2012-02-26 Thread giuseppe calamita
s84)) #now the issue is that the "pred" object is str(pred) Factor w/ 4 levels "2","3","4","5": 3 3 3 3 3 3 3 3 3 3 ... - attr(*, "names")= chr [1:920591] "24389" "24390" "24391" "25729" ... that is, it con

Re: [R] How to remove $ (Dollar sign) from string

2012-04-22 Thread Giuseppe Marinelli
In data martedì 10 aprile 2012 13:34:13, Nevil Amos ha scritto: > How do I remove a "$" character from a string sub() and gsub() with "$" or > "\$" as pattern do not work. > > > sub("$","","ABC$DEF") > > [1] "ABC$DEF" > > > sub("\$","","ABC$DEF") > > Error: '\$' is an unrecognized escape in cha

[R] speeding up a special product of three arrays

2008-05-08 Thread Giuseppe Paleologo
ing the columns won't do. Is there any way to implement this efficiently in R or should I resign to do this in C? Thanks, Giuseppe [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinf

[R] mice package

2016-03-07 Thread Giuseppe Biondi Zoccai
nable to impute with mice the stress features among also those who did not undergo any stress test? Or should I best create a factor variable such as stress_status (0- no stress, 1-stress with low tolerance, 2-stress with high tolerance, and so forth)? Thanks for the help Giuseppe [[altern

[R] R: daisy(): space allocation issue

2010-08-26 Thread Giuseppe De Vitis
advance, shouldn't I? Thanks Giuseppe -Messaggio originale- Da: Gavin Simpson [mailto:gavin.simp...@ucl.ac.uk] Inviato: 26 August 2010 17:47 A: abanero Cc: r-help@r-project.org Oggetto: Re: [R] daisy(): space allocation issue On Thu, 2010-08-26 at 07:35 -0700, abanero wrote: >