× Aidez la recherche contre le COVID-19 avec votre ordi ! Rejoignez l'équipe PC Astuces Folding@home
 > Tous les forums > Forum Bureautique
 Créer un fichier xls suivant des données de scripts
Ajouter un message à la discussion
Page : [1] 
Page 1 sur 1
williamsss
  Posté le 24/01/2015 @ 15:32 
Aller en bas de la page 
Petit astucien

Bonjour,

Pour réaliser ce graphique http://meteo.besse83.free.fr/imfix/trendtest-Moy%205%20s%E9ries-1998-30%202034.png une personne utilise un srcipt à voir ci-dessous suivant son message.

Le script qui permet de calculer l'IP à 95% de ce graphe est :

t <- read.table("C:/r/moy5.csv", sep=";", dec=",", header=T)

# note: this does not yet handle NAs in the trend well

# be careful with the bounderies - 1979 for sats, pre-1900 for ground

# need boundary tests

y1 <- 1998 # end of trend

t0 <- 30 # years of trend

y0 <- y1 - t0 # start of trend

t1 <- 3 * t0 # number of years to display

y2 <- y0 + t1 # max of x-axis

ymin <- -1.0 # min of y-axis

ymax <- 1.5 # max of y-axis

# column idx

gis <- 3

cru <- 4

noa <- 5

uah <- 6

rss <- 7

moy5 <- 8

idx <- gis

src = c("GISTEMP", "HADCRU", "NOAA", "UAH", "RSS", "Moy 5 séries")

# annualize monthly data

t$year <- as.integer(t[,1])

# get the record of interest, monthly to annual

Tann <- aggregate(t[,idx],list(t$year), FUN=mean, na.rm=T)

target <- Tann[Tann[,1]> y0 & Tann[,1]<= y1,]

# a series for x-axis

tyear <- y0:y2

# some set setup

x <- target[,1]

y <- target[,2]

new <- data.frame(x = seq(y0+1,y1))

# get trend; hi and lo for 95% ci

trend_lm <- lm(y ~ x)

trend_lo = cbind(x,predict(lm(y ~ x), interval="prediction")[,2])

trend_hi = cbind(x,predict(lm(y ~ x), interval="prediction")[,3])

lm_lo <- lm(trend_lo[,2] ~ trend_lo[,1])

lm_hi <- lm(trend_hi[,2] ~ trend_hi[,1])

# extend trend +/- ci into future years

trend_lm_lo <- tyear*lm_lo$coef[2] + lm_lo$coef[1]

trend_lm_hi <- tyear*lm_hi$coef[2] + lm_hi$coef[1]

trend_lm_mid <- tyear*trend_lm$coef[2] + trend_lm$coef[1]

# set y axis

ymax = max(trend_lm_hi) + 0.5

ymin = min(trend_lm_lo) - 0.5

# file name to write plot

png(paste("trendtest-",src[idx-2],"-",y1,"-",t0,".png",sep=""),height=480, width=540, bg="white")

# plot

dstamp=Sys.Date()

plot(tyear,rep(0,length(tyear)), ylim=c(ymin,ymax), col="white", xlab="Années", ylab="Anomalie (°C)") # hidden - it just sets the graph axis

title(main=paste("Le réchauffement s'est-il arrêté en",y1,"?\n","Test tendance/",t0,"ans",src[idx-2]),

sub=paste("http://rhinohide.wordpress.com+ modif ChristianP ", dstamp), cex.lab=0.75)

# warm, cold, no trend boxes

lines(c(y1,y1),c(ymin,ymax), type="l", col="gray") # beginning year of trend test

rect(y1,trend_lm_hi[t0+1],y2,ymax, col="lightpink")

rect(y1,trend_lm_lo[t0+1],y2,ymin, col="lightskyblue")

vtx <- min(min(tyear[trend_lm_lo > trend_lm_hi[t0+1]]), y2)

polygon(rbind(c(y1,trend_lm_lo[t0+1]),c(vtx,trend_lm_lo[vtx-y0+1]),c(y2,trend_lm_hi[t0+1]),c(y2,trend_lm_lo[t0+1])), col="light gray", lty="solid")

# label trend boxes

text(y2,(trend_lm_lo[t0+1]+trend_lm_hi[t0+1])/2, "Stabilité", pos=2, cex=0.8)

text(y1+1,(ymin+trend_lm_lo[t0+1])/2, "Refroidissement", pos=4, cex=0.8)

text(y1+1,(ymax+trend_lm_hi[t0+1])/2, "Réchauffement", pos=4, cex=0.8)

text(y0+1,(ymax+trend_lm_hi[t0+1])/3, paste("Tendance/",t0," ans \n95% intervalle de prédiction"), pos=4, cex=0.8)

# trend slope data

a1 <- format(trend_lm$coeff[2]*10,digits=2)

a2 <- format(trend_lm_hi[t0+1] - trend_lm_mid[t0+1],digits=2)

a3 <- format(-trend_lm_lo[t0+1] + trend_lm_mid[t0+1],digits=2)

text(y0+1,ymin+0.1, paste(a1,"°C par décennie\n IP +",a2,"/ -",a3), pos=4, cex=0.8)

# points in the trend

points(target[,1],target[,2], pch=16) # data points for annual global T

# trend lines and +/- prédiction interval

abline(trend_lm,col="black", lwd=2) # trendline for annual global T

lines(trend_lm_hi ~ tyear,col="red", lty=2, lwd=2) # 2 SD upper limit of monthly global T

lines(trend_lm_lo ~ tyear,col="blue", lty=2, lwd=2) # 2 SD lower limit of monthly global T

# add the trend test points

t_test <- Tann[Tann[,1]> y1,]

points(t_test[,1],t_test[,2], col="green", pch=16) # data points beyond the trend setters

dev.off()

Le fichier csv pour le script dessus est : http://meteo.besse83.free.fr/imfix/moy5.csv

Suivant ce script d'après-vous dans un fichier Excel avec les données et années de l'anomalie de la température comme ici http://cjoint.com/?0AypKsDGKy0 quels formules doit on y mettre pour pouvoir réaliser le graphique du 1er lien ci-dessus qui est réalisé avec ce script ???

Merci

Williams



Modifié par williamsss le 24/01/2015 15:34
Publicité
w36xb2w
 Posté le 24/01/2015 à 16:09 
Aller en bas de la page Revenir au message précédent Revenir en haut de la page
  Astucien

Bonjour

Déjà ici il y a une adresse qui ne fonctionnera peut être pas

sub=paste("http://rhinohide.wordpress.com+ modif ChristianP ", dstamp), cex.lab=0.75)

http://meteo.besse83.free.fr/

williamsss
 Posté le 24/01/2015 à 16:48 
Aller en bas de la page Revenir au message précédent Revenir en haut de la page
Petit astucien
w36xb2w a écrit :

Bonjour

Déjà ici il y a une adresse qui ne fonctionnera peut être pas

sub=paste("http://rhinohide.wordpress.com+ modif ChristianP ", dstamp), cex.lab=0.75)

http://meteo.besse83.free.fr/

Bonjour W36xb2w,

L'adresse http://rhinohide.wordpress.com fonctionne bien mais ci-dessus c'est le fait qu'il y a le "+" qui s'est incorporer dans l'adresse du lien comme tout est rattaché dans ce que la personne m'a donné concernant le code du script.

merci

Williams

Page : [1] 
Page 1 sur 1

Vous devez être connecté pour poster des messages. Cliquez ici pour vous identifier.

Vous n'avez pas de compte ? Créez-en un gratuitement !


Les bons plans du moment PC Astuces

Tous les Bons Plans
179,99 €Joystick Logitech G Saitek X52 Pro Flight Control System + Microsoft Flight Simulator à 179,99 €
Valable jusqu'au 25 Septembre

Cdiscount fait une promotion sur un lot comprenant le tout nouveau Microsoft Flight Simulator sur PC accompagné du joystick Logitech G Saitek X52 Pro Flight Control System qui passe à 179,99 € au lieu de 250 €. Ce système H.O.T.A.S. (3M) X52 Pro Pro Flight de Saitek est un contrôleur de vol main sur manche et manette entièrement intégré qui satisfait, voire dépasse, les demandes des pilotes sur simulateur de vol les plus aguerris.


> Voir l'offre
99,99 €Ecouteurs intra-auriculaires sans fil Bluetooth Jabra Elite 65t à 99,99 €
Valable jusqu'au 25 Septembre

Amazon fait une superbe promotion sur les écouteurs intra-auriculaires sans fil Bluetooth Jabra Elite 65t qui passent à 99,99 € livrés gratuitement. On les trouve ailleurs à partir de 150 €. Pensés pour vous permettre de vivre votre musique comme vous l'entendez et pour assurer des conversations d'une grande clarté, les écouteurs Jabra Elite 65t ont tout pour vous séduire ! Côté son, l'égaliseur personnalisable vous permet d'ajuster le rendu sonore tandis que la technologie avancée à quatre microphones supprime les bruits de fond lors des appels pour des conversations limpides. Au quotidien, ces écouteurs Jabra restent parfaitement en place et proposent même une fonction de commande vocale simplifiée. Enfin, vous apprécierez les 5 heures d'autonomie avec une seule charge, qui se transforment en 15 heures avec l'étui de recharge compact.  Le Jabra Elite 65t est certifié IP55 et offre une garantie de 2 ans contre l’humidité et la poussière. Où que vous soyez et quelle que soit votre activité, vous profitez d’une liberté absolue pour les appels et la musique. 


> Voir l'offre
21,91 €Souris sans fil Logitech M280 à 21,91 €
Valable jusqu'au 25 Septembre

Amazon fait une petite promotion sur la souris sans fil Logitech M280 qui passe à 21,91 €. On la trouve ailleurs à partir de 30 €. Sa forme profilée est adaptée aux droitiers. Son autonomie est de 18 mois (pile AA que vous pourrez remplacer).


> Voir l'offre

Sujets relatifs
Comment créer fichier client avec carte de fidélité
Créer un livret à paritr d'un fichier word ou pdf
Données d'un Fichier txt avec EXCEL
Créer texte dans une fichier JPG
créer fichier compressé avec pages word
CREATION BASE DE DONNEES AVEC 1 AUTRE FICHIER EXCE
creer un fichier pour envoyer un mailling
Comment créer un fichier .reg et l'installer
Fichier de base de données
créer un fichier pdf
Plus de sujets relatifs à Créer un fichier xls suivant des données de scripts
 > Tous les forums > Forum Bureautique