[R] recovering PLSDA equations for prediction (outside R)

Hi all, I'm trying to use the result of a PLSDA model outside R, but I'm having a really hard time finding documentation on how to write the model from the results in the fit object. In order to provide a good starting point I believe that this code creates a good model that is stored in the fit va

Re: [R] Soft Question: Where to find this reference.

See book: Statistical Models in S (ISBN-10: 041283040X; ISBN-13: 978-0412830402) Chapter 5, p. 145. • Analysis of Variance; Designed Experiments, J. M. Chambers, A. E. Freeny and R. M. Heiberger, Statistical Models in S, J. M. Chambers and T. J. Hastie (editors), Wadsworth & Brooks/Cole, P

Re: [R] ncdf size error

name="GLW") varPBLH = ncvar_def("PBLH","m", list(dim1,dim2), -9, longname="PBLH") varTSLB = ncvar_def("TSLB","Celsius", list(dim1,dim2), -9, longname="TSLB") nc.ex = nc_create("/media/D/output.nc",list(varT2,

Re: [R] ncdf size error

), -9, longname="TSLB") # associate the netcdf variable with a netcdf file # put the variable into the file, and # close nc.ex = create.ncdf("/media/D/output.nc",list(varT2, varQ2, varQVAPOR, varQCLOUD, varPSFC, varW10, varW, varVEGFRA, varMAXVEGFRA, varTPREC, varFRAC_F

[R] ncdf size error

Hi I am using both ncdf and ncdf4 libraries and with both I keep getting the same error: Error in R_nc_enddef: NetCDF: One or more variable sizes violate format constraints. Error in R_nc_sync: NetCDF: Operation not allowed in define mode. This happens when I try to create.ncdf() a file with more

Re: [R] writing spdiags function for R

I am posting here the brilliant solutions, gently provided by Prof JC Nash , to me; so that people struggling in the future with the same issue can find a way through. FYI, compared to the original Matlab implementation: 1) "it does not handle the case with more than one input, and 2) (m > n)

Re: [R] writing spdiags function for R

al efficiency would still be there, even fixing the dimensionality issue. I really hope somebody out there can help me to figure this out. thanks a lot for your precious time and attention, Moreno -- View this message in context: http://r.789695.n4.nabble.com/writing-spdiags-function-for-R-tp4

[R] choose many files

Hello everyone on the list I want to run a FindU function on as many files as in a directory. Due of my limitations on data management I am calling the function, wait until finish the process, call for the next file and so on. I have almost 120 CSV files on the same directory. This the example

[R] problem on reading many files

Hi, I have to run almost 120 stations files of temperatura (mx and min), separated, and rain. I am following the instructions on RHtests tutorial as on http://www.cmc.org.ve/mediawiki/index.php?title=Preparando_los_datos link. Bit I have no success on running on multiple files. I have the .ls fi

Re: [R] changing cex pointwise in lattice

quot;group", "z"), row.names = c(NA, 50L ), class = "data.frame") -- --- José M. Blanco-Moreno Dept. de Biologia Vegetal (Botànica) Facultat de Biologia Universitat de Barcelona Av. Diagonal 643 08028 Barcelona SPAIN

[R] changing cex pointwise in lattice

11171L, 11219L, 11267L, 11315L, 11363L, 11411L, 11459L, 11507L, 11555L, 11603L, 11651L, 11699L, 11747L, 11795L, 11843L, 11891L, 11939L, 11987L, 12035L, 12083L, 12131L, 12179L, 12227L, 12275L)) -- --- José M. Blanco-Moreno Dept. de Biologia Vegetal (Botànica

Re: [R] writing spdiags function for R

ced by n <- nrow(A) p <- 2*n-1 d <- seq(-(n-1),(n-1)) yes, all good points, I have acted the changes you suggested, and it works smooth. thanks, Moreno -- The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336. ___

Re: [R] writing spdiags function for R

Hi Ben, thank you soo much:) your code worked, and we managed to halven computing time. I am copying below the final version of the function, if anybody is willing to use/improve it. thanks again, Moreno ### spdiags = function(A

[R] writing spdiags function for R

you spot other places where the function could be optimized I would be very glad to hear it! thank you very much in advance for your kind help, Best, Moreno ### ## it works only for square matrices ## it could work with spar

[R] XML package: modify info of internal node

result of the error... what am i missing here? is there any simpler way to do it? thanks a lot for your kind help, Moreno P.S. I also try: xmlSApply(xpathApply( current, "/annotation/filename")[[1]], xmlValue) = newname obtaining the same result. -- The University of Edinburgh i

Re: [R] modulo operation

or my hasty question when facing this behaviour. -- --- José M. Blanco-Moreno Dept. de Biologia Vegetal (Botànica) Facultat de Biologia Universitat de Barcelona Av. Diagonal 645 08028 Barcelona SPAIN --- phone: (+34) 934 039 86

[R] modulo operation

> 25/0.2 [1] 125 > 25 - floor(25/0.25)*0.25 [1] 0 (At least this second one is what I would expect from the code in arithmetic.c, lines 168 to 178) -- --- José M. Blanco-Moreno Dept. de Biologia Vegetal (Botànica) Facultat de Biologia Universitat de

[R] Neyman-Scott in 1-dimension?

ckage can deal with linear point processes? If yes, how? Thank you very much for your help, Moreno -- The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336. __ R-help@r-project.org mailing list

Re: [R] updating arguments of formulae

am putting in copy the group you suggested me... Thanks again, Moreno Quoting "Meyners,Michael,LAUSANNE,AppliedMathematics" : Moreno, I leave the discussion on the mixed models to others (you might consider the SIG group on mixed models as well for this), but try a few hints to

[R] updating arguments of formulae

red only once. If the second version of the formula is the correct one, I don't understand how I can update arguments inside the formula rather than adding things on his right side... thus, in the ideal case, how do I go from something like this: depM ~ OS + (1 + OS | Sb2) to something like t

[R] updating specific arguments of formula

-- The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org

[R] updating specific arguments of formula

-- The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org

[R] updating specific arguments of formula

-- The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336. __ 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] Combine vectors in order to form matrixes with combn

> Moreno, > > I don't understand exactly what it is you are trying to do. Can you explain what you want your matrices to look like? Perhaps give an example. > > Bryan Dear Bryan, let's try with an example: i have 8 column vectors v1<-matrix(c(2,5,6,8)) v2<-

[R] Combine vectors in order to form matrixes with combn

invalid 'ncol' value (too large or NA) Inoltre: Warning message: In combn(c(coeff1, coeff2, coeff3, coeff4, coeff5, coeff6, coeff7, : si è prodotto un NA per coercizione Some ideas? Thanks in advance, Moreno Tutto Incluso Light: chiami e navighi senza limiti a soli 9,95 euro al mese.

Re: [R] saving large matrices with decimal numbers

be the problem? Thanks a lot, Moreno write.csv (image, "/afs/inf.ed.ac.uk/user/s07/s0788806/R/r_data/ProdEyeTrack/DataFrame/Clustered/Clut/Cluster1/Cl1Verb-tr.csv.1") Quoting Uwe Ligges : Moreno Ignazio Coco wrote: Hi everybody, I am struggling to save txt. or csv. large s

[R] saving large matrices with decimal numbers

. It seems that R manages to save up around 800 X 600 decimal matrices but for bigger sizes the output is an empty file. Does anyone know whether there is a limit in saving such matrices and how can I work around this problem? Thanks a lot, Regards, Moreno -- The University of Edinburgh

Re: [R] Solutions of equation systems

Ravi Varadhan wrote: [cut] It seems to work! Thanks a lot Ravi and Berend. Moreno. __ 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

Re: [R] Solutions of equation systems

n variables, solve works perfectly. Thanks in advance! Moreno __ 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.

Re: [R] Solutions of equation systems

? Thanks in advance, Moreno __ 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] Solutions of equation systems

Any idea? Thanks in advance, Moreno __ 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.

Re: [R] [R-sig-Geo] error in trmesh (alphahull package)

ist r-sig-...@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo -- --- José M. Blanco Moreno Dept. de Biologia Vegetal (Botànica) Universitat de Barcelona Av. Diagonal 645 08028 Barcelona SPAIN -

Re: [R] heatmap-changing column or row names

it is more where in the plot, zero is centered. If I say abline(v=0) it produces a vertical line much more right on the x axis than the actual 0. It is a problem of where the axis of the plots are set and they don't seem to match the layout of the heatmap. Any suggestion on this? Moren

[R] heatmap-changing column or row names

the heatmap some crucial timepoints. I wanted to do it by adding a straight abline on the heatmap plot abline(v=80,col=green) but also this doesn't seem to have any result It seems that it is quite difficult to add any information to the heatmap, but perhaps some of you might know how t

[R] 3dscatterplot -different colors for different factors of the same variable

points3d(pinkpts[,3],pinkpts[,2],pinkpts[,4],col="pink",type="o") R has plotted only the first set and wasn't able to overlay other points. 2) I have generated a "for" loop that should have taken points information ,arranged into a list object, plotting them on

[R] kalman filter random walk

Hi, I'm trying to use the kalman filter to estimate the variable drift of a random walk, given that I have a vector of time series data. Anyone have any thoughts on how to do this in R? Thanks, Alex [[alternative HTML version deleted]] __ R-h

[R] install problem 2.6 mac os x 10.5

I am having the very SAME problem. I removed Version 2.5.1 before trying to install 2.6.0. Please help. TIA, Oscar __ 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.

[R] standard errors of regression coefficients

Hi, If I have two vectors x and y and I do lm(y~x) and now I want to define variables that are the standard errors of the slope and intercept, how do I do that? Thanks, Alex [[alternative HTML version deleted]] __ R-help@r-project.org mailing