Displaying 20 results from an estimated 110 matches similar to: "Fitting a model"
2011 Jun 23
3
help- subtitles for multiple charts
Hello, I have a problem with my script. I don'y know how to apply subtitles.
I have 9 charts per page (for combination of 3 blocks and 3 treatments). I
want to have subtitles for this interaction (e.g. Block A Trt 1, Block A Trt
2, ...)
MBE$bt<- interaction(MBE$Block,MBE$trt)
par(mfrow=c(3,3))
for(i in unique(MBE$bt)){
ss <- MBE$bt==i
plot(MBE$Year[ss], MBE$DBH[ss])
2001 Mar 13
3
Invalid ICMP error
To anyone who is well versed with samba,
I am currently setup with ip masquerading with a cable modem connection to
my ISP. I have an error message that continues to show up on my linux box.
It says:
10.58.0.1 sent an invalid ICMP error to a broadcast.
This means it is my linux box sending an invalid packet out through my
modem. This started to happen when I installed samba. So I hav come
2010 May 19
3
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 Vall?s- SPAIN
Tlf. +34 93 581 1910
Fax: +34 93 581 1573
2008 Jul 09
3
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 Saub?
Grup de Biologia Molecular de Llevats
Facultat de Veterin?ria de la Universitat Aut?noma de
2008 Aug 27
3
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 Universitat Aut?noma de Barcelona
Edifici V, Campus UAB
08193 Cerdanyola del Vall?s- SPAIN
Centro de Investigaci?n
2009 Feb 17
3
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 Veterin?ria de la Universitat Aut?noma de Barcelona
Edifici V, Campus UAB
08193 Cerdanyola del Vall?s-
2008 Aug 28
2
Read a file
Hello,
I have a text file with this structure:
# File created = Thursday, August 28, 2008 3:33:02 PM GMT
# Data set = 373 2 1 C:\Bruker\TOPSPIN GABRMN
# Spectral Region:
# LEFT = 4.5 ppm. RIGHT = 0.5 ppm.
#
# SIZE = 13111 ( = number of points)
#
# In the following ordering is from the 'left' to the 'right' limits!
# Lines beginning with '#' must be considered as
2008 Feb 20
3
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 Barcelona
Edifici V, Campus UAB
08193 Cerdanyola del Vall?s- SPAIN
Centro de Investigaci?n Biom?dica en
2007 Nov 22
2
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 de Biologia Molecular de Llevats
Facultat de Veterin?ria de la Universitat Aut?noma de Barcelona
2009 Aug 19
3
Fitting a logistic regression
Hello,
I have this data:
Time AMP
0 0.2000000
10 0.1958350
20 0.2914560
40 0.6763628
60 0.8494534
90 0.9874526
120 1.0477692
where AMP is the concentration of this metabolite with time. If you plot
the data, you can see that it could be fitted using a logistic
regression. For this purpose, I used this code:
AMP.nls <- nls(AMP~SSlogis(Time,Asym, xmid, scal), data
2008 Feb 05
0
[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 <daniel.valverde at uab.cat>
A: R Help <r-help at r-project.org>
Hello,
I am tryining to build a package called NMRTools and when I check it
using R CMD
2008 Feb 19
1
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 Veterin?ria de la Universitat Aut?noma de Barcelona
Edifici V, Campus UAB
08193 Cerdanyola del Vall?s-
2008 Jun 06
1
lsmeans
Hello,
I have the next function call:
lme(fixed=Error ~ Temperature * Tumour ,random = ~1|ID, data=error_DB)
which returns an lme object. I am interested on carrying out some kind
of lsmeans on the data returned, but I cannot find any function to do
this in R. I'have seen the effect() function, but it does not work with
lme objects. Any idea?
Best,
Dani
--
Daniel Valverde Saub?
Grup
2008 Sep 01
2
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 Veterin?ria de la Universitat Aut?noma de Barcelona
Edifici V, Campus UAB
08193 Cerdanyola del Vall?s- SPAIN
2008 Feb 20
1
insert() function
Hello,
I am trying to insert a certain number of points into a certain position
of a vector with this code:
x <- seq(1:10909)
x1 <- c(13112-10909)
spect1 <- rnorm(13112)
interpol <- approx(x,spect1,xout=c(seq(from=1, by=((10909 - 1)/(x1 -
1)), length.out=x1)))
pos <- round(interpol$x,0)
intensities <- interpol$y
spect2 <- insert(spect1,ats=pos,values=intensities)
2008 Aug 27
1
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 Cerdanyola del Vall?s- SPAIN
Centro de Investigaci?n Biom?dica en Red
en Bioingenier?a, Biomateriales y
2008 Jun 02
1
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
--
Daniel Valverde Saub?
Grup de Biologia Molecular de Llevats
Facultat de Veterin?ria de la
2008 Sep 03
1
Problems with playwith()
Hello,
I am trying to execute an example provided in the help files for the
playwith() function. I try to execute the first example:
library(playwith)
if (interactive()) {
options(device.ask.default = FALSE)
## Scatterplot (Lattice graphics).
## Labels are taken from rownames of data.
## Just click on the plot to identify points.
playwith(xyplot(Income ~ log(Population / Area),
data =
2007 Aug 20
1
LDA decission boundaries
Hello,
I would like to plot the results of a LDA analysis plotting the
discriminant scores with the decission boundaries on it with rggobi. I
have GGobi already installed on my computer. I have three classes, so
the plot would be LD1xLD2 plus the decission boundaries. Here there is
the code I use make the plot:
library(MASS)
data <- zgcppr273K.pca$x[,1:7]
Tumor <-
2008 Feb 26
1
GLM
Hello,
I am trying to perform a glm analysis on a 68x13113 matrix (named
data.spect). The first column corresponds to the predictor
(data.spect[,1]) and the rest to the response variables
(data.spect[,2:13113]). When I try this code
glmObject <- glm(data.spect[,2:13113]~data.spect[,1])
I get the following error:
Error: (subscript) logical subscript too long
Could anyone help me on