Re: [R] standardized regression coefficients in GAM

Dear Dr Fox, Thank you so much for your response! I will report the results of the regression models only based on unstandardized coefficients then, it seems to be best. Thanks for taking the time to respond and clarify this issue for me! Best regards, Dani Sent from Outlook<h

[R] standardized regression coefficients in GAM

? I do not see how I can standardize dummy variables. Thank you! Best, Dani <http://aka.ms/weboutlook> [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/m

[R] How to calculate Confidence Interval for a prediction using Partial Regression?

y_0 values. Right now, the only way I can see to find that is to calculate a Confidence Interval for my partial regression-based predictions (the prediction interval), and I don't know how to. Does anyone know the formula to calculate it? Any help will be welcome! Thanks in advance!

[R] Pattern matching and replacement in matrix

Hello, Is there any function like gsub(), that can match and replace patterns in a matrix or a data frame? Cheers! Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària de la Universitat Autònoma de Barcelona Edifici V, Campus UAB 08193 Cerdanyola del

[R] Fitting a function

Hello, It is a very naive question, but here it is. I have this values: x: 0.5 4 6 8 12 y: 0.021 0.021 0.020 0.018 0.012 I need to fit a function to them. How can I do it with R? Thank you so much! Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària

Re: [R] Fitting a logistic regression

Thank you very much, Jun. This is what I was looking for. Best! Dani On Wed, 2009-08-19 at 09:52 -0500, Jun Shen wrote: > I would suggest a model with a baseline level, something like > > nls(AMP~E0+(Emax-E0)*Time**gamma/(EC50**gamma+Time**gamma),data=your > data, > start=list(

[R] Fitting a logistic regression

eve a better fit? Forgive me if it is a stupid question, but I am just starting with non linear regression. Thank you, Dani -- [?] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide h

Re: [R] Package HDF5

Hi Thomas, Thanks for responding. You are right about bioconductor, however their solution was to implement RHDF5 package for linux only, which allows selecting fields, but not ranges of the data. Any idea about ranges of data? My datafiles are as big as 20GB. thx, Dani On Fri, Mar 20, 2009 at

[R] Package HDF5

The package works fine, but it seems not to provide access to fields, but only to the whole data in a certain file (function hdf5load(file, load = TRUE, verbosity = 0, tidy = FALSE)) Since hdf5 organizes the data and metadata in a hierarchical structure, we must explore in our problem (file>7GB),

[R] R parser for If-else

Hi list, I don't know if somebody has spent a lot of time debugging strange problems with if else positioning - the parser seems to recognize only the syntax bellow - this is the only way of making these pieces of code to work. As far as i'm concerned, no examples were available (it would be so a

[R] Combination

Hello, I have a sequence of numbers: seq(1:50) and I would like to have all the possible combinations with this numbers without repeating any combination: 11, 12, 13, ... ,22,23,24,... How can I do it? Best, Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de

Re: [R] Fitting a model

N +34 93 5814126 En/na Eik Vettorazzi ha escrit: Actually "drm" as posted before fits a sigmoid curve (a generalized logistic function with 4 parameters, see ?LL.4), so I didn't get the point of your new question. Dani Valverde schrieb: Thank you all for your answers. If

Re: [R] Fitting a model

Thank you all for your answers. If you look at the plot resulting from my data, it seems that it is some kind of sigmoid function, not only polynomial. How could I fit it? Best, Dani Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària de la Universitat

[R] Fitting a model

ph <- c(4.4,4.6,4.8,5,5.2,5.4,5.6,5.8,6,6.2,6.4,6.6,6.8,7,7.2,7.4) plot(ph,delta,ylab=c(expression(Delta*delta)),xlab="pH") Which kind of model can I fit on these, so that can I predict for a given delta the pH of my sample? Once the model is fitted, how can I plot it on the graph? Best rega

[R] movie3d function in the rgl package

Dear list, I'm using the function "movie3d" in the package "rgl" to create a .gif animation of a 3d graphic. The program "ImageMagik" is working properly, R packages are working, basic examples available in the manual also working fine. Problem Solved: when I tried to create more complex movies,

[R] Problems with playwith()

GLib, and finally it says that the runtime environment asks to close the application in an appropiate way. Any idea on how to solve thi? I have all packages installed. I use R 2.7.2 and GTK+ 2.12.9. I have all the required packages installed. Best, Dani -- Daniel Valverde Saubí Grup de

Re: [R] Text

algaard ha escrit: Dani Valverde wrote: Hello, I would like to place some text that should appear in the same position of the graphic device, preferably outside the plotting area, regardless the x and y axes limits. How can I do this? I think you can use title() with an adj= argument. Othe

[R] Text

Hello, I would like to place some text that should appear in the same position of the graphic device, preferably outside the plotting area, regardless the x and y axes limits. How can I do this? Best, Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de

[R] Read a file

;left' to the 'right' limits! # Lines beginning with '#' must be considered as comment lines. # 1628.40625 1628.40625 1964.40625 2242.0625 2533.5 2937.90625 3448.0 3923.96875 Is it possible to read it with R, something like scan() but keeping the structure? Best, Dani

[R] Combine x-axis labels size

Hello, I would like to combine different label sizes in the same x-axis, for example 1 2 3 4 5 6 7 8 9 How can I do it? Best, Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària de la Universitat Autònoma de Barcelona Edifici V, Campus UAB 08193

[R] Closest value

Hello, I have this command: x.axis <- seq(from=0.5, to=4.5, length.out=13112) How can I which of the x.axis components is the closest to a given value, for example 3.2? Best, Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària de la Universi

[R] Anova

t;numeric" with a factor response will be ignored/ How can I solve this? Best, Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària de la Universitat Autònoma de Barcelona Edifici V, Campus UAB 08193 Cerdanyola del Vallès- SPAIN Centro de Inve

[R] Expression in axis

Hello, I am creating a plot and I would like to know how to put this expression to the y axis µmol/10^6 cells I've tried some combinations using the expression() function, but none of them worked. Any idea? Best, Dani -- Daniel Valverde

[R] ANOVA error

ike to perform an ANOVA on different variables. Is the call aov(Cr+mIb+...~Treatment, data=myMatrix) o correct way of doing it? What I would like to have is a p value for the variance for each response variable. Best, Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Ve

[R] fit.contrast error

or_DB) but I got this error Error in `contrasts<-`(`*tmp*`, value = c(-0.5, 0.5)) : contrasts apply only to factors My database is a dataframe, very similar to that of the Orthodont. Could anyone give me some advise on how to solve the problem? Best, Dani -- Daniel Valverde Saubí

[R] lsmeans

ffect() function, but it does not work with lme objects. Any idea? Best, Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària de la Universitat Autònoma de Barcelona Edifici V, Campus UAB 08193 Cerdanyola del Vallès- SPAIN Centro de Investigación Biomédica en R

[R] LDA and centroids

Hello, I have carried out an lda analysis using the lda function of MASS package. I have plotted the LD1xLD2 to represent the data. Now I would like to get the centroids for each group of data and plot it on the LD1xLD2 graph. How can I get the centroid value from the lda object? Best, Dani

[R] Problems with lme

me? Best, Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària de la Universitat Autònoma de Barcelona Edifici V, Campus UAB 08193 Cerdanyola del Vallès- SPAIN Centro de Investigación Biomédica en Red en Bioingeniería, Biomateriales y Nanomedicina (CIBER

[R] Clustering large data matrix

not just two principal components? Many thanks. Best, Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària de la Universitat Autònoma de Barcelona Edifici V, Campus UAB 08193 Cerdanyola del Vallès- SPAIN Centro de Investigación Biomédica en Red en

[R] GLM

get the following error: Error: (subscript) logical subscript too long Could anyone help me on solving this problem? Best, Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària de la Universitat Autònoma de Barcelona Edifici V, Campus UAB 08193 Cerdanyola

Re: [R] insert() function

Great! Now it works. Many thanks. Dani Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària de la Universitat Autònoma de Barcelona Edifici V, Campus UAB 08193 Cerdanyola del Vallès- SPAIN Centro de Investigación Biomédica en Red en Bioingeniería, Biomateriales y

Re: [R] insert() function

L_1.58 svMisc_0.9-5 svIDE_0.9-5 loaded via a namespace (and not attached): [1] cluster_1.11.9 grid_2.6.2 lattice_0.17-6 tools_2.6.2 Can you help? Best, Dani Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària de la Universitat Autòno

Re: [R] Rd to latex

Hello, I've tried with R CMD Rdeconv and it works, but, when I want to paste it into my main latex file, there are lots of undefined environments. I loaded the package Rd.sty in the latex document preface, but it does not seem to work. Any idea? Best, Dani Daniel Valverde Saubí Gr

[R] Rd to latex

Hello, I would like to convert an Rd object to a latex file, so that I can put it in my thesis. How can I do it? I tryed latex(), but it only works for code... Best, Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària de la Universitat Autònoma de

[R] insert() function

os, values = intensities) at insert(spect1, ats = pos, values = intensities) Can anyone help me on how can I solve it? I know that "pos" and "intensities" have the same length, so I don't know why I get this error. Any ideas will be welcome. Best, Dani -- Daniel Val

Re: [R] Interpolation between 2 vectors

code: spect1 <- c(1:10909) for(i in 1:length(interpol$x)){ append(spect1,interpol$y[i],interpol$x[i]) } but it didn't work. Any idea? Best, Dani Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària de la Universitat Autònoma de Barcelona Edifici V, Campus UAB

[R] Interpolation between 2 vectors

Hello, I have two vectors, one with 13112 points and the other one with 10909. I wonder if there is a way to interpolate the data so the shorter vectors has the same number of points as the longer one. Best, Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de

[R] Compile Rd file

Hello, I am trying to compile a file .Rd into .html using R CMD Rdconv -t=html /file.rd/>/file.html./ It seems that the process works ok but instead of having for example loo.cv(NMRTools) at the top of the html file I have loo.cv{unknown}. Any ideas on how to solve this? Best, D

[R] [Fwd: Problems when checking a package]

Forget my last message, it was a really stupid problem. Sorry about the mess. Best, Dani Missatge original Assumpte: Problems when checking a package Data: Tue, 05 Feb 2008 11:20:28 +0100 De: Dani Valverde <[EMAIL PROTECTED]> A: R Help Hello

[R] Problems when checking a package

otstrap.rd but not in lda.bootstrap.r. The problem is that everything appears in both files. Maybe I misundestand the error message. Any idea? Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària de la Universitat Autònoma de Barcelona Edifici V, Campus

[R] Package compilation

Hello, Is there an easy way to compile a package so that it can be installed as a binary file under Windows? Now I have the source code. I am not used in compiling, and the documentation seems too hard for me, so an easy way would be great. Best regards, Dani -- Daniel Valverde Saubí Grup

[R] Installing packages

suggestion to solve the problem? Best regards, Dani -- Daniel Valverde Saubí Grup de Biologia Molecular de Llevats Facultat de Veterinària de la Universitat Autònoma de Barcelona Edifici V, Campus UAB 08193 Cerdanyola del Vallès- SPAIN Centro de Investigación Biomédica en Red en Bioingeniería

Re: [R] Odp: Cutting a vector

I do not want to write breaks=c(4,8,12,...,13108, +Inf), but some expresion to avoid writing such a long argument. Best regards, Dani En/na Petr PIKAL ha escrit: > Hi > > you shall be more specific. Do you want to split your vector according > some pionts > > split(x, fi

[R] Cutting a vector

Hello, I have a vector, lets say x <- 1:50 I would like it to be cut at certain points, being for example 1:5, 6:11, 12:17, ... How can I do it? I have tried the cut() function, but I don not know how to place the cutting points properly. Best regards, Dani -- Daniel Valverde Saubí Grup